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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Measurement and Interpretation of Temporal Variations of the Earths Gravity Field Using GPS and SLR Data

This investigation has had four main thrusts: 1) The analysis of seasonal variations of the Earth's gravitational field using Lageos 1 and 2 SLR data and comparisons to geophysical models. We have estimated the annual variation of the gravity field via a spherical harmonic expansion complete to degree and order 4. We have also constructed a similar model using models of the annual variation in the gravity field due to atmospheric, hydrologic, and ocean mass redistribution. These three models, when combined together, are in excellent agreement with the variations observed by satellite laser ranging. An article on these results was published in the journal Geophysical Research Letters. 2) The second thrust of our investigation has been to analyze the output of a Global Climate Model (GCM) to determine if the GRACE gravity mission can be expected to detect climate change signals. Working with Gary Russell at the Goddard Institute for Space Studies (GISS), we have determined that there are several large secular signals that GRACE might be able to detect, including secular changes in snow cover, sea ice, polar ice, ocean mass, and other variables. It is possible that some of these signals could be detected with 5 years of GRACE measurements - its hard to judge this because the interannual variability in the GCM, which could mask the climate signals, is unreliable. Certainly a follow-on GRACE mission could detect these signals when compared to the data from the initial GRACE mission.). An article on these results will be published in the journal Journal of Geophysical Research. 3) In the last year of the investigation, we developed a new technique for analyzing temporal gravity variations using "geophysical fingerprints", which was successfully demonstrated on 20 years of satellite laser ranging data [Nerem et al., 20031. 4]. We also participated in a workshop on future satellite gravity measurements, which resulted in paper on measuring ocean mass variations using GRACE [Nerem et al., 20031 and on using laser interferometry for future gravity missions [Bender et al., 20031].

Nerem, R. Steven↗
NASA NTRS2017

B-plane Evolution Under Highly Non-Keplerian Dynamics

Juno is the second of a series of New Frontiers missions and was launched in 2011. The spacecraft was set en route to Jupiter, and its Jupiter Orbit Insertion occurred on 2016-07-04 after a five-year cruise in deep space. The mission phase just prior to this large maneuver is called approach. As navigation team, we were responsible for precise orbit determination of the Juno spacecraft to ensure a successful orbit insertion. To evaluate the navigation performance we employed the B-plane mapping, where the incoming hyperbolic velocity defines the plane perpendicular to it, and the time to go to hit the plane will directly map into the uncertainties of the target point. The B-plane is a convenient tool when assessing the navigation performance as it is not affected by non-linearities due to the gravitational pull of the targeted body. However, as the Jupiter Orbit Insertion approached we started to realize that the B-plane mapping uncertainties were significantly larger than expected. We found that the B-plane mapping time had tremendous influence on the covariance inflation in the B-plane. Because of the close approach distance during the Jupiter Orbit Insertion and the strong interaction with J2 spherical harmonic coefficient, the Juno dynamics were too far from a Keplerian one assumed for the B-plane mapping. We first discuss the analytical approach of the B-plane mapping uncertainty and perform a numerical analysis to isolate the components of the Juno dynamics that cause the B-plane covariance inflation.

Bordi, John↗
DOE OSTI2023

BEYONDPLANCK VIII. Efficient sidelobe convolution and corrections through spin harmonics

We introduce a new formulation of the Conviqt convolution algorithm in terms of spin harmonics, and apply this to the problem of sidelobe correction for BEYONDPLANCK, the first end-to-end Bayesian Gibbs sampling framework for CMB analysis. We compare our implementation to the previous Planck LevelS implementation, and find good agreement between the two codes in terms of accuracy, but with a speed-up reaching a factor of 3–10, depending on the frequency bandlimits, l max and m max . The new algorithm is significantly simpler to implement and maintain, since all low-level calculations are handled through an external spherical harmonic transform library. We find that our mean sidelobe estimates for Planck LFI are in good agreement with previous efforts. Additionally, we present novel sidelobe rms maps that quantify the uncertainty in the sidelobe corrections due to variations in the sky model.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Free-air and Bouguer gravity anomalies and the Martian crustal dichotomy

Free-air and Bouguer gravity anomalies from a 50x50 field, derived from re-analysis of Viking Orbiter and Mariner 9 tracking data and using a 50x50 expansion of the current Mars topography and the GSFC degree 50 geoid as the equipotential reference surface, with the Martian crustal dichotomy are compared. The spherical harmonic topography used has zero mean elevation, and differs from the USGS maps by about 2 km. In this field the dichotomy boundary in eastern Mars lies mostly at -1 to -2 km elevation. Bouguer gravity anomalies are shown on a map of Noachian, Hesperian, and Amazonian age terrains, simplified from current geologic maps. The map is centered at 300 deg W to show the continuity of the dichotomy boundary. Contour interval is 100 mgals. Gravity and topography were compared along approximately 40 profiles oriented parallel to the dichotomy boundary topographic gradient, to determine how the geophysical character of the boundary changes along its length and what this implies for its origin and development.

Frey, Herbert↗
NASA NTRS1980

Deconvolution and analysis of wide-angle longwave radiation data from Nimbus 6 Earth radiation budget experiment for the first year

One year of longwave radiation data from July 1975 through June 1976 from the Nimbus 6 satellite Earth radiation budget experiment is analyzed by representing the radiation field by a spherical harmonic expansion. The data are from the wide field of view instrument. Contour maps of the longwave radiation field and spherical harmonic coefficients to degree 12 and order 12 are presented for a 12 month data period.

Bess, T. D.↗
NASA NTRS1992

Analysis of crustal structure of Venus utilizing residual Line-of-Sight (LOS) gravity acceleration and surface topography data. A trial of global modeling of Venus gravity field using harmonic spline method

To construct Venus' gravity disturbance field (or gravity anomaly) with the spacecraft-observer line of site (LOS) acceleration perturbation data, both a global and a local approach can be used. The global approach, e.g., spherical harmonic coefficients, and the local approach, e.g., the integral operator method, based on geodetic techniques are generally not the same, so that they must be used separately for mapping long wavelength features and short wavelength features. Harmonic spline, as an interpolation and extrapolation technique, is intrinsically flexible to both global and local mapping of a potential field. Theoretically, it preserves the information of the potential field up to the bound by sampling theorem regardless of whether it is global or local mapping, and is never bothered with truncation errors. The improvement of harmonic spline methodology for global mapping is reported. New basis functions, a singular value decomposition (SVD) based modification to Parker & Shure's numerical procedure, and preliminary results are presented.

Fang, Ming↗
NASA NTRS1976

Spherical harmonic representation of the gravity field of Mars using a short-arc technique

Planetary gravity fields are primarily determined from Doppler tracking of orbiting spacecraft which are typically in high eccentricity orbits. Thus, the non-central gravitational components perturb the orbit primarily in the vicinity of periapsis. To find a computationally efficient technique and to reduce the influence of unmodeled forces on the spacecraft, an analysis has been performed of the effectiveness of processing only data in the vicinity of periapsis using multiple short arcs. Mariner 9 data are processed for the study, using 1, 2, 3 and 4 hour arcs and estimating gravity fields from 4th through 7th degree. The resulting gravity fields are inter-compared and also compared with previously published results.

Daniels, E. F.↗
NASA NTRS1995

Upper Thermosphere Winds and Temperatures in the Geomagnetic Polar Cap: Solar Cycle, Geomagnetic Activity, and Interplanetary Magnetic Field Dependencies

Ground-based Fabry-Perot interferometers located at Thule, Greenland (76.5 deg. N, 69.0 deg. W, lambda = 86 deg.) and at Sondre Stromfjord, Greenland (67.0 deg. N, 50.9 deg. W, lambda = 74 deg.) have monitored the upper thermospheric (approx. 240-km altitude) neutral wind and temperature over the northern hemisphere geomagnetic polar cap since 1983 and 1985, respectively. The thermospheric observations are obtained by determining the Doppler characteristics of the (OI) 15,867-K (630.0-nm) emission of atomic oxygen. The instruments operate on a routine, automatic, (mostly) untended basis during the winter observing seasons, with data coverage limited only by cloud cover and (occasional) instrument failures. This unique database of geomagnetic polar cap measurements now extends over the complete range of solar activity. We present an analysis of the measurements made between 1985 (near solar minimum) and 1991 (near solar maximum), as part of a long-term study of geomagnetic polar cap thermospheric climatology. The measurements from a total of 902 nights of observations are compared with the predictions of two semiempirical models: the Vector Spherical Harmonic (VSH) model of Killeen et al. (1987) and the Horizontal Wind Model (HWM) of Hedin et al. (1991). The results are also analyzed using calculations of thermospheric momentum forcing terms from the Thermosphere-ionosphere General Circulation Model TGCM) of the National Center for Atmospheric Research (NCAR). The experimental results show that upper thermospheric winds in the geomagnetic polar cap have a fundamental diurnal character, with typical wind speeds of about 200 m/s at solar minimum, rising to up to about 800 m/s at solar maximum, depending on geomagnetic activity level. These winds generally blow in the antisunward direction, but are interrupted by episodes of modified wind velocity and altered direction often associated with changes in the orientation of the Interplanetary Magnetic Field (IMF). The central polar cap (greater than approx. 80 magnetic latitude) antisunward wind speed is found to be a strong function of both solar and geomagnetic activity. The polar cap temperatures show variations in both solar and geomagnetic activity, with temperatures near 800 K for low K(sub p) and F(sub 10.7) and greater than about 2000 K for high K(sub p) and F(sub 10.7). The observed temperatures are significantly greater than those predicted by the mass spectrometer/incoherent scatter model for high activity conditions. Theoretical analysis based on the NCAR TIGCM indicates that the antisunward upper thermospheric winds, driven by upstream ion drag, basically 'coast' across the polar cap. The relatively small changes in wind velocity and direction within the polar cap are induced by a combination of forcing terms of commensurate magnitude, including the nonlinear advection term, the Coriolis term, and the pressure gradient force term. The polar cap thennospheric thermal balance is dominated by horizontal advection, and adiabatic and thermal conduction terms.

Killeen, T. L.↗
DOE OSTI2021

A selection function toolbox for subsets of astronomical catalogues

ABSTRACT Large catalogues are ubiquitous throughout astronomy, but most scientific analyses are carried out on smaller samples selected from these catalogues by chosen cuts on catalogued quantities. The selection function of that scientific sample – the probability that a star in the catalogue will satisfy these cuts and so make it into the sample – is thus unique to each scientific analysis. We have created a general framework that can flexibly estimate the selection function of a sample drawn from a catalogue in terms of position, magnitude and colour. Our method is unique in using the binomial likelihood and accounting for correlations in the selection function across position, magnitude, and colour using Gaussian processes and spherical harmonics. We have created a new open-source python package selectionfunctiontoolbox that implements this framework and used it to make three different estimates of the APOGEE DR16 red giant sample selection function, as a subset of 2MASS, with each estimate using an increasing amount of technical knowledge of the APOGEE targeting. In a companion paper, we applied our methodology to derive estimates of the astrometric and spectroscopic selection functions of Gaia EDR3. Our framework will make it trivial for astrophysicists to estimate the selection function that they should be using with the custom sample of stars that they have chosen to answer their scientific question.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1990

Gravity field improvement using GPS data from Topex/Poseidon - A covariance analysis

A covariance analysis is performed using a realistic scenario for processing 10 days of GPS data, to obtain the expected improvement to the GEM-T2 gravity field. The gravity bin technique has been refined to compute the covariance matrix associated with the spherical harmonic gravity field. It is shown that the GPS data from one ten-day arc of Topex/Poseidon with no a priori can improve medium degree and order (3-26) sigmas for the parameters in the GEM-T2 gravity field by more than an order of magnitude.

Bertiger, Willy I.↗
NASA NTRS2011

Rapid Calculation of Spacecraft Trajectories Using Efficient Taylor Series Integration

A variable-order, variable-step Taylor series integration algorithm was implemented in NASA Glenn's SNAP (Spacecraft N-body Analysis Program) code. SNAP is a high-fidelity trajectory propagation program that can propagate the trajectory of a spacecraft about virtually any body in the solar system. The Taylor series algorithm's very high order accuracy and excellent stability properties lead to large reductions in computer time relative to the code's existing 8th order Runge-Kutta scheme. Head-to-head comparison on near-Earth, lunar, Mars, and Europa missions showed that Taylor series integration is 15.8 times faster than Runge- Kutta on average, and is more accurate. These speedups were obtained for calculations involving central body, other body, thrust, and drag forces. Similar speedups have been obtained for calculations that include J2 spherical harmonic for central body gravitation. The algorithm includes a step size selection method that directly calculates the step size and never requires a repeat step. High-order Taylor series integration algorithms have been shown to provide major reductions in computer time over conventional integration methods in numerous scientific applications. The objective here was to directly implement Taylor series integration in an existing trajectory analysis code and demonstrate that large reductions in computer time (order of magnitude) could be achieved while simultaneously maintaining high accuracy. This software greatly accelerates the calculation of spacecraft trajectories. At each time level, the spacecraft position, velocity, and mass are expanded in a high-order Taylor series whose coefficients are obtained through efficient differentiation arithmetic. This makes it possible to take very large time steps at minimal cost, resulting in large savings in computer time. The Taylor series algorithm is implemented primarily through three subroutines: (1) a driver routine that automatically introduces auxiliary variables and sets up initial conditions and integrates; (2) a routine that calculates system reduced derivatives using recurrence relations for quotients and products; and (3) a routine that determines the step size and sums the series. The order of accuracy used in a trajectory calculation is arbitrary and can be set by the user. The algorithm directly calculates the motion of other planetary bodies and does not require ephemeris files (except to start the calculation). The code also runs with Taylor series and Runge-Kutta used interchangeably for different phases of a mission.

Scott, James R.↗
DOE OSTI2020

Learning physics at future $e^–e^+$ colliders with machine

Information deformation and loss in jet clustering are one of the major limitations for precisely measuring hadronic events at future $e^–e^+$ colliders. Because of their dominance in data, the measurements of such events are crucial for advancing the precision frontier of Higgs and electroweak physics in the next decades. We show that this difficulty can be well-addressed by synergizing the event-level information into the data analysis, with the techniques of deep neutral network. In relation to this, we introduce a CMB-like observable scheme, where the event-level kinematics is encoded as Fox-Wolfram (FW) moments at leading order and multi-spectra of spherical harmonics at higher orders. Then we develop a series of jet-level (w/ and w/o the FW moments) and event-level classifiers, and analyze their sensitivity performance comparatively with two-jet and four-jet events. As an application, we analyze measuring Higgs decay width at $e^–e^+$ colliders with the data of 5ab –1 @240GeV. The precision obtained is significantly better than the baseline ones presented in documents. We expect this strategy to be applied to many other hadronic- event measurements at future $e^–e^+$ colliders, and to open a new angle for evaluating their physics capability.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2017

Drag-Free Control and Drag Force Recovery of Small Satellites

Drag-free satellites provide autonomous precision orbit determination, accurately map the static and time varying components of Earth's mass distribution, aid in our understanding of the fundamental force of gravity, and will ultimately open up a new window to our universe through the detection and observation of gravitational waves. At the heart of this technology is a gravitational reference sensor, which (a) contains and shields a free-floating proof mass from all non-gravitational forces, and (b) precisely measures the position of the test mass inside the sensor. Thus, both test mass and spacecraft follow a pure geodesic in spacetime. By tracking the position of a low Earth orbiting drag-free satellite we can directly determine the detailed shape of geodesics and through analysis, the higher order harmonics of the Earths geopotential. This paper explores two different drag-free control systems on small satellites. The first drag-free control system is a continuously compensated single thruster 3-unit CubeSat with a suspension-free spherical proof-mass. A feedback control system commands the thruster and Attitude and Determination Control System to fly the tender spacecraft with respect to the test mass. The spheres position is sensed with a LED-based differential optical shadow sensor, its electric charge controlled by photoemission using UV LEDs, and the spacecraft position is maintained with respect to the sphere using an ion electrospray propulsion system. This configuration is the most fuel-efficient drag-free system possible today. The second drag-free control system is an electro-statically suspended cubical proof-mass that is operated with a low duty cycle, limiting suspension force noise over brief, known time intervals on a small GRACE-II -like satellite. The readout is performed using a laser interferometer, which is immune to the dynamic range limitations of voltage references. This system eliminates the need for a thruster, enabling drag-free control systems for passive satellites. In both cases, the test mass position, GPS tracking data, and commanded actuation, either thrust or suspension system, can be analyzed to estimate the 3-axis drag forces acting on the satellite. The data produces the most precise maps of upper atmospheric drag forces and with additional information, detailed models that describe the dynamics of the upper atmosphere and its impact on all satellites that orbit the Earth. This paper highlights the history, applications, design, laboratory technology development and highly detailed simulation results of each control system.

Small Satellites↗
NASA NTRS1975

The complex magnetic field of Jupiter

An analysis of the characteristics of the magnetic field of the planet Jupiter is presented. The data were obtained during the flight of Pioneer 11 space probe, using a high field triaxial fluxgate magnetometer. The data are analyzed in terms of traditional Schmitt normalized spherical harmonic expansion fitted to the observations in a least squares sense. Tables of data and graphs are provided to summarize the findings.

Acuna, M. H.↗
NASA NTRS1982

An analysis of the structure of Saturn's magnetic field using charged particle absorption signatures

A new technique is derived for determining the structure of Saturn's magnetic field. This technique uses the observed positions of charged particle signatures due to the satellites and rings of Saturn to determine the parameters of an axially symmetric, spherical harmonic model of the magnetic field using the method of least squares. Absorption signatures observed along the Pioneer 11, Voyager 1, and Voyager 2 spacecraft trajectories are used to derive values for the orientation of the magnetic symmetry axis relative to Saturn's axis of rotation, the axial displacement of the center of the magnetic dipole from the center of Saturn, and the magnitude of the external field component. Comparing these results with the magnetic field model parameters deduced from analyses of magnetometer data leads to a preference for models that incorporate a northward offset of the dipole center by about 0.05 Saturn radii.

Chenette, D. L.↗
NASA NTRS1993

Theory and operation of a three-gate time-of-flight velocity analyzer

Theoretical considerations and test results are presented for a new-type velocity analyzer for incident fast neutral particles, positive ions, and negative ions. Velocity analysis is carried out by means of a pulsed, three-gate time-of-flight (TOF) technique capable of eliminating alias velocities (harmonics) to sixth order. In addition the design and operation are presented of a four-element ion lens system, with small spherical and chromatic aberrations, suitable for interfacing a large-diameter ion beam from the TOF section with a subsequent mass analyzer.

Martus, K. E.↗
NASA NTRS1981

Deconvolution of wide field-of-view radiometer measurements of earth-emitted radiation. II - Analysis of first year of Nimbus 6 ERB data

The theory of deconvolution considered by Smith and Green (1981) is applied to Nimbus 6 data in order to interpret the data with the maximum obtainable resolution. Attention is given to the data sample, sun-contaminated measurements, the measurement model, the deconvolution, the degree variance, the spherical harmonic coefficients, the geographical distribution of longwave radiation, time histories of zonal coefficients, and the effect of a grid system. Degree variance plots for 12 months of longwave radiation data show that the limit for a spherical harmonic representation of the Nimbus 6 wide field-of-view longwave data is 12th degree. The degree variance plots also reveal that most of the power is in the lower degree terms. Contour maps of the radiation field show the geographical distribution of earth emitted radiant exitance for each month.

Bess, T. D.↗