Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spherical harmonic analysis”

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 gravity field of Mars

A complete sixth-degree spherical-harmonic model of the gravitational potential of Mars has been produced from an analysis of Doppler data obtained from the radio tracking of the Mariner-9 spacecraft for over 200 orbital revolutions. The equipotential surface defined by this model is estimated to be accurate to within 100 m below 60 deg north latitude and to within 300 m above. The dominant feature of this surface is a Tharsis-related bulge which rises nearly two kilometers above the surrounding areas. Gravitational 'lows' are found to correspond to the Hellas depression and the Valles Marineris.

Reasenberg, R. D.↗

Summary of initial results from the GSFC fluxgate magnetometer on Pioneer 11

The main magnetic field of Jupiter was measured by the Fluxgate Magnetometer on Pioneer 11 and analysis reveals it to be relatively more complex than expected. In a centered spherical harmonic representation with a maximum order of n = 3 (designated GSFC model 04), the dipole term (with opposite polarity to the Earth's) has a moment of 4.28 Gauss x (Jupiter radius cubed), tilted by 9.6 deg towards a system 111 longitude of 232. The quadrupole and octupole moments are significant, 24% and 21% of the dipole moment respectively, and this leads to deviations of the planetary magnetic field from a simple offset tilted dipole for distances smaller than three Jupiter radii. The GSFC model shows a north polar field strength of 14 Gauss and a south polar field strength of 10.4 Gauss. Enhanced absorption effects in the radiation belts may be predicted as a result of field distortion.

Acuna, M. H.↗

Structure of the atmosphere of Jupiter from Pioneer 10 and 11 radio occultation measurements

A technique for incorporating effects of the oblateness of Jupiter's atmosphere into the data analysis procedure for radio-occultation measurements is described which makes use of a spherical harmonic representation of the gravity field to compute the shape of the planet. With this technique, the center of refraction is located by the radius of curvature and the normal direction at the closest approach point of the ray. The present technique, a subsequent approximation involving the use of a constant center of refraction, and the Abel integral transform inversion method for spherical planets are applied to Pioneer 10 and 11 data. All the intermediate results obtained show a temperature inversion between 10 and 100 mb, with temperatures from 130 to 170 K at 10 mb and from 80 to 120 K at 100 mb. Comparison of the radio-occultation profiles with radiative-convective equilibrium models and the temperature profile based on Pioneer 10 IR radiometer data indicates good agreement between the models and the occultation results.

Kliore, A. J.↗

Gaussian statistics of the cosmic microwave background: Correlation of temperature extrema in the COBE DMR two-year sky maps

We use the two-point correlation function of the extrema points (peaks and valleys) in the Cosmic Background Explorer (COBE) Differential Microwave Radiometers (DMR) 2 year sky maps as a test for non-Gaussian temperature distribution in the cosmic microwave background anisotropy. A maximum-likelihood analysis compares the DMR data to n = 1 toy models whose random-phase spherical harmonic components a(sub lm) are drawn from either Gaussian, chi-square, or log-normal parent populations. The likelihood of the 53 GHz (A+B)/2 data is greatest for the exact Gaussian model. There is less than 10% chance that the non-Gaussian models tested describe the DMR data, limited primarily by type II errors in the statistical inference. The extrema correlation function is a stronger test for this class of non-Gaussian models than topological statistics such as the genus.

Kogut, A.↗

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

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

The Local Stellar Velocity Field via Vector Spherical Harmonics

We analyze the local field of stellar tangential velocities for a sample of 42,339 nonbinary Hipparcos stars with accurate parallaxes, using a vector spherical harmonic formalism. We derive simple relations between the parameters of the classical linear model (Ogorodnikov-Milne) of the local systemic field and low-degree terms of the general vector harmonic decomposition. Taking advantage of these relationships, we determine the solar velocity with respect to the local stars of (V(sub X), V(sub Y), V(sub Z)) (10.5, 18.5, 7.3) +/- 0.1 km s(exp -1) not corrected for the asymmetric drift with respect to the local standard of rest. If only stars more distant than 100 pc are considered, the peculiar solar motion is (V(sub X), V(sub Y), V(sub Z)) (9.9, 15.6, 6.9) +/- 0.2 km s(exp -1). The adverse effects of harmonic leakage, which occurs between the reflex solar motion represented by the three electric vector harmonics in the velocity space and higher degree harmonics in the proper-motion space, are eliminated in our analysis by direct subtraction of the reflex solar velocity in its tangential components for each star. The Oort parameters determined by a straightforward least-squares adjustment in vector spherical harmonics are A=14.0 +/- 1.4, B=13.1 +/- 1.2, K=1.1 +/- 1.8, and C=2.9 +/- 1.4 km s(exp -1) kpc(exp -1). The physical meaning and the implications of these parameters are discussed in the framework of a general linear model of the velocity field. We find a few statistically significant higher degree harmonic terms that do not correspond to any parameters in the classical linear model. One of them, a third-degree electric harmonic, is tentatively explained as the response to a negative linear gradient of rotation velocity with distance from the Galactic plane, which we estimate at approximately -20 km s(exp -1) kpc(exp -1). A similar vertical gradient of rotation velocity has been detected for more distant stars representing the thick disk (z greater than 1 kpc), but here we surmise its existence in the thin disk at z less than 200 pc. The most unexpected and unexplained term within the Ogorodnikov-Milne model is the first-degree magnetic harmonic, representing a rigid rotation of the stellar field about the axis -Y pointing opposite to the direction of rotation. This harmonic comes out with a statistically robust coefficient of 6.2 +/- 0.9 km s(exp -1) kpc(exp -1) and is also present in the velocity field of more distant stars. The ensuing upward vertical motion of stars in the general direction of the Galactic center and the downward motion in the anticenter direction are opposite to the vector field expected from the stationary Galactic warp model.

Galaxy↗

Spatial and temporal variations of monthly mean total columnar ozone derived from 7 years of BUV data

Ninth-degree spherical harmonic models are used to represent monthly mean total columnar ozone fields, with spatial spectra defined in terms of simple analytic expressions representing the rms value of each spherical harmonic coefficient. The spectra are combined with estimation theory techniques to develop a realistic error analysis accounting for errors due to nonglobal coverage. Amplitude and phase variations for longitudinal waves one, two and three are presented, and it is shown that these waves reinforce in space and time in the period of high poleward ozone transport in the northern hemisphere. In the southern hemisphere, the wave two and three components are much smaller, with wave one maximum amplitude occurring in September and October.

Tolson, R. H.↗

Large-scale, near-Earth, magnetic fields from external sources and the corresponding induced internal field

Data from MAGSAT analyzed as a function of the Dst index to determine the first degree/order spherical harmonic description of the near-Earth external field and its corresponding induced field. The analysis was done separately for data from dawn and dusk. The MAGSAT data was compared with POGO data. A local time variation of the external field persists even during very quiet magnetic conditions; both a diurnal and 8-hour period are present. A crude estimate of Sq current in the 45 deg geomagnetic latitude range is obtained for 1966 to 1970. The current strength, located in the ionosphere and induced in the Earth, is typical of earlier determinations from surface data, although its maximum is displaced in local time from previous results.

Langel, R. A.↗

Large-scale, near-field magnetic fields from external sources and the corresponding induced internal field

Data from Magsat analyzed as a function of the Dst index to determine the first degree/order spherical harmonic description of the near-earth external field and its corresponding induced field. The analysis was done separately for data from dawn and dusk. The Magsat data was compared with POGO data. A local time variation of the external field persists even during very quiet magnetic conditions; both a diurnal and 8-hour period are present. A crude estimate of Sq current in the 45 deg geomagnetic latitude range is obtained for 1966 to 1970. The current strength, located in the ionosphere and induced in the earth, is typical of earlier determinations from surface data, although its maximum is displaced in local time from previous results.

Langel, R. A.↗

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

Griffin: A MOOSE-based reactor physics application for multiphysics simulation of advanced nuclear reactors

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor physics application for multiphysics simulations of advanced reactor designs jointly developed by Idaho National Laboratory and Argonne National Laboratory. This paper summarizes the motivation, significance, architecture, design, and features of Griffin. Griffin offers flexible and extensible features to address the challenges associated with advanced reactor designs. These features range from fundamental particle transport to specific reactor physics tasks. The features cover a wide range including on-the-fly and traditional two-step cross-section generation methods, steady-state and transient transport solvers suitable for both heterogeneous and homogeneous models, high-fidelity depletion where thousands of isotopes can be tracked and low-fidelity depletion characterized by burnup, etc. The most fundamental aspect that sets Griffin apart from other reactor analysis codes is that it is developed based on the MOOSE framework. A modular development approach is strongly enforced, with multiphysics being an essential element considered since the beginning of Griffin’s development. Griffin links various MOOSE physics modules and couples to other MOOSE-based applications and non-MOOSE-based applications for multiphyiscs simulations. Griffin includes three modules: ISOXML for preparing and managing multigroup cross sections, radiation transport for solving the neutron transport equation, and reactor analysis for user-oriented reactor physics analysis functionalities. Griffin uses various finite element methods for spatial discretization, multigroup approximation for energy discretization and discrete ordinates method, spherical harmonics expansion method, and diffusion approximation for streaming direction discretization to solve the neutron transport equation. Griffin’s flexibility is evidenced through Griffin’s various applications to fast reactor, high-temperature reactor, pebble bed reactor, molten salt reactor, and microreactor designs. Griffin development follows the software quality assurance procedure for MOOSE-based applications and with software requirements consistent with the ASME NQA-1 standard. Griffin has been adopted into the reactor analysis system for the U.S. NRC and is in use at U.S. companies, universities and national laboratories.

97 MATHEMATICS AND COMPUTING↗

Short and medium range structure in elastic deformation of metallic and covalent glasses

Here, we present a concise methodology to analyze structural response to the applied stress in amorphous solids, including metallic glasses (MG), glassy selenium, silica and polycarbonate, using high energy x-ray diffraction and atomic pair distribution function (PDF) analysis. To assess the structural anisotropy induced by applied axial stress, diffraction data were expanded into spherical harmonics. Using Bessel transformation, components of the structure function were converted into isotropic and anisotropic PDFs. The PDFs were compared to the expected model behavior for ideal elastic deformation to separate homogeneous affine strain from local non-affine strains. In metallic glass the range of non-affine deformation is limited to the nearest neighbor shell, suggesting local strain relaxation under stress that occurs even in the elastic regime. Beyond the second atomic shell strain is uniform. However, in glassy silica, polycarbonate and selenium strong local bonding inhibits local displacements and strain in short range order is accommodated by rotation of local units. Interestingly, beyond a molecular unit, deformation in covalent systems is similar to MG, and response of the medium range order scales with the macroscopic stress.

glassy structure↗

Influence of spatial filtering on possible anisotropies in solar oscillations

Full disk Doppler observations of solar oscillations were used to compare the amplitudes of sectoral modes propagating along the equator with those of similar modes propagating along a great circle aligned with the poles. It is found that the amplitudes are generally not equal for the two classes of modes, but the results are sensitive to analysis procedures attempting to isolate the different modes of oscillation. Spatial filtering of the data using spherical harmonics suggests that greater amplitudes are associated with 'polar' sectoral modes than with 'equatorial' sectoral modes.

Hill, Frank↗

The Impact of New Accelerometer Transplant Data (ACH) on GRACE Follow-On Along-Orbit Inter-Satellite Laser Ranging Observations and Monthly Time-Variable Gravity and Mascon Solutions

GRACE-D accelerometer data show significant bias jumps since one month after the launch of the GRACE Follow-On (GRACE-FO) satellites in May 2018, making them inapplicable for correcting GRACE-FO products for non-gravitational accelerations. The GRACE-FO Science Data System (SDS) compensated this issue by transplanting the GRACE-C accelerometer data toward that of GRACE-D. Recently, GRACE-FO SDS implemented an updated transplant method, used in the latest release of GRACE-FO data. Here, we examine the impact of updated accelerometer transplant data (ACH) on GRACE-FO measurements at all levels: (a) Level-1B inter-satellite laser ranging residuals measured along satellite orbit, (b) Level-2 time-variable gravity solutions from all SDS centers (JPL, CSR, and GFZ), and (c) Level-3 mascon solutions. We show that inter-satellite laser ranging residuals are modified at low frequencies below 1 mHz, affecting the along-orbit analysis of large-scale time-variable gravity signals. When mapped into monthly Level-2 spherical harmonic coefficients of geopotential, the low-frequency change in inter-satellite ranging residuals leads to substantial improvement of coefficients associated with resonant orders (i.e., 15, 30, 45, etc.) and C 30 . We also present an improved SLR-derived C 30 which significantly improves the agreement with updated GRACE-FO C 30 at seasonal and interannual timescales. Moreover, we demonstrate the noise reduction in mass change estimates from new GRACE-FO Level-2 data over oceans, Greenland, and Antarctica for all SDS solutions. GRACE-FO mascon solutions show a moderate change in the updated release. Our comprehensive analyses demonstrate high-quality estimates of non-gravitational accelerations by the updated transplant method, resulting in more accurate GRACE-FO time-variable gravity and mass change observations.

Laser ranging interferometer↗

Spherical harmonic representation of the main geomagnetic field for world charting and investigations of some fundamental problems of physics and geophysics

Progress in the harmonic analysis of MAGSAT data is reported. Single-day data sets were subdivided into information on the sunrise side of the Earth and information on the sunset side of the Earth. Data for the main and external fields each demonstrate a clear and consistent systematic difference between the sets of data which was determined to be, due to ionospheric currents which differ from the sunset to the sunrise terminator. A toroidal field was analyzed for and determined to be an apparent toroidal field resulting from electric currents concentrated in the two terminators. Progressive elimination of auroral zone data demonstrates that the information presented does not arise from complications due to Birkeland currents.

Barraclough, D. R.↗

Stellar convection. II - A multimode numerical solution for convection in spheres

A stable, equilibrium solution for convection in a self-gravitating sphere of Boussinesq fluid was computed by using a modal analysis in which the theta, phi dependence of the fluid is expanded in a set of 168 spherical harmonics. The temperature, velocity, and convective flux of the fluid as well as the kinetic and thermal energy spectra as functions of wavelength are computed. The spectra are found to be in agreement with both experimental observations and analytic scaling laws. The dynamics of the energy cascade is examined by computing the ratio of the amount of energy dissipation at a particular wavelength to the amount of energy produced at that same wavelength. It is found that there is only a slight cascade of kinetic energy to smaller wavelengths but a large cascade of thermal energy.

Marcus, P. S.↗

Speculation on a Solar Chronometer for Climate

Solar activity has been correlated to climatic fluctuations and has been postulated as a major factor in quasi-periodic global climatic change. However, correlations are not explanations of physical mechanisms and do not couple cause with effect. A mechanism for a chronometer for solar output variability is proposed based on relations between properties of thermonuclear fusion, nuclear magnetic moment, and nuclear magnetic resonance. A fundamental oscillation of a nucleus with a net nuclear magnetic moment (NMM) is the precession of its axis of rotation when subjected to a magnetic field. Nuclear magnetic resonance (NMR) is the preferred frequency of precession for a nucleus of a particular isotope when placed in a magnetic field of specific intensity. The NMM for those isotopes involved in the proton-proton (p-p) chain pathway for solar fusion varies from strong positive to strong negative. Individual fusion events, for hydrogen and helium isotopes which release varying amounts of energy, may be controlled by NMR frequencies. The pulses of energy from fusion events occurring at NMR frequencies in the solar interior may be transformed into pressure or gravity waves that emerge as gravity or acoustic waves at the surface. Dictated by spherical harmonics, certain wavelengths may be reinforced and reenter the solar interior to modulate the fusion process. Qualitative analysis of solar and climatic data support the interaction of the three basic components of the chronometer, magnetic activity, oscillation frequency, and solar energy output.

Perry, Charles A.↗