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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Coupled Groups of g-Modes in a Sun with Mixed Core

Groups of linear g-modes can sum to create long-lived nonlinear oscillations in small "hot volumes" very deep in the Sun that help drive the modes. In these volumes (dimensions -10 Mm), the time average rate of He-3 burning doubles as temperature fluctuations exceed 10% and rises by an order of magnitude for fluctuations of 25%. To be consistent with locally large motions, we impose a mixed shell on an otherwise standard solar model before computing g-mode solutions. Mixing in the assumed shell r = (0.10+/-0.03) R(sub sun) is rapid (<<10(exp 6) yr) with slower mixing somewhat beyond. If l is the principal spherical harmonic index, a set of g-modes for any single l less than or equal to 15 with five consecutive radial harmonics can be excited with nearly linear thermal amplitudes, A(sub T) less than or equal to 0.053, throughout the star and a fractional temperature fluctuation in its hot volume of (Delta)T/T less than or equal to 0.18. These thresholds for excitation will become smaller when sets for several values of l are computed simultaneously. There is some evidence for the rotation of g-mode sets in the long solar activity record and g-mode upward wave flux has been suggested to explain the 1.3 yr reversing flows tentatively detected below the Sun's convective envelope (CE). The large local amplitudes needed for excitation implies that g-modes may transport a non-negligible fraction of the solar luminosity, yet their near linear amplitudes outside the hot volume suggests amplitudes over most of the solar surface that would be barely detectable for l > 3. A formalism is presented for summing the g-modes and estimating growth rates under the approximation that modes are strictly linear except in a hot volume which holds only a few percent of mode kinetic energy. Finally over the range 2 less than or equal to l less than or equal to 30, we summed all zonal harmonics, m, for a given l and computed the relative angular orientations that would release the most nuclear energy. This should be close to the physically preferred angular state of such a family and a few examples were displayed.

Wolff, Charles L.↗
NASA NTRS2005

Re-Entry Simulation and Landing Area for YES2

The REST simulator includes many parameters: a) Inertial <-> Fix to Earth reference system; b) Geodetic <-> Geocentric coordinates; c) Rotational velocity of the Atmosphere; d) Effect of the rotation of the Earth; e) Bulge effect of the Earth; f) Spherical harmonic expansion for the Earth s gravitational potential, J2 (zonal); g) Heat flux, temperature in the wall; h) Drag coefficient for different regimes; i) Flow regime status; j) Density model NRLMSISE-00; k) Wind model HWM-93; l) G2S atmospheric model with the latest meteorological conditions and m) Landing area (Monte Carlo Simulations)

Calzada, Silvia↗
NASA NTRS2000

Computer-implemented method and apparatus for autonomous position determination using magnetic field data

A computer-implemented method and apparatus for determining position of a vehicle within 100 km autonomously from magnetic field measurements and attitude data without a priori knowledge of position. An inverted dipole solution of two possible position solutions for each measurement of magnetic field data are deterministically calculated by a program controlled processor solving the inverted first order spherical harmonic representation of the geomagnetic field for two unit position vectors 180 degrees apart and a vehicle distance from the center of the earth. Correction schemes such as a successive substitutions and a Newton-Raphson method are applied to each dipole. The two position solutions for each measurement are saved separately. Velocity vectors for the position solutions are calculated so that a total energy difference for each of the two resultant position paths is computed. The position path with the smaller absolute total energy difference is chosen as the true position path of the vehicle.

Ketchum, Eleanor A.↗
NASA NTRS1998

Multipath calibration in GPS pseudorange measurements

Novel techniques are disclosed for eliminating multipath errors, including mean bias errors, in pseudorange measurements made by conventional global positioning system receivers. By correlating the multipath signals of different satellites at their cross-over points in the sky, multipath mean bias errors are effectively eliminated. By then taking advantage of the geometrical dependence of multipath, a linear combination of spherical harmonics are fit to the satellite multipath data to create a hemispherical model of the multipath. This calibration model can then be used to compensate for multipath in subsequent measurements and thereby obtain GPS positioning to centimeter accuracy.

Kee, Changdon↗
NASA NTRS2008

Fast Plasma Instrument for MMS: Simulation Results

Magnetospheric Multiscale (MMS) mission will study small-scale reconnection structures and their rapid motions from closely spaced platforms using instruments capable of high angular, energy, and time resolution measurements. The Dual Electron Spectrometer (DES) of the Fast Plasma Instrument (FPI) for MMS meets these demanding requirements by acquiring the electron velocity distribution functions (VDFs) for the full sky with high-resolution angular measurements every 30 ms. This will provide unprecedented access to electron scale dynamics within the reconnection diffusion region. The DES consists of eight half-top-hat energy analyzers. Each analyzer has a 6 deg. x 11.25 deg. Full-sky coverage is achieved by electrostatically stepping the FOV of each of the eight sensors through four discrete deflection look directions. Data compression and burst memory management will provide approximately 30 minutes of high time resolution data during each orbit of the four MMS spacecraft. Each spacecraft will intelligently downlink the data sequences that contain the greatest amount of temporal structure. Here we present the results of a simulation of the DES analyzer measurements, data compression and decompression, as well as ground-based analysis using as a seed re-processed Cluster/PEACE electron measurements. The Cluster/PEACE electron measurements have been reprocessed through virtual DES analyzers with their proper geometrical, energy, and timing scale factors and re-mapped via interpolation to the DES angular and energy phase-space sampling measurements. The results of the simulated DES measurements are analyzed and the full moments of the simulated VDFs are compared with those obtained from the Cluster/PEACE spectrometer using a standard quadrature moment, a newly implemented spectral spherical harmonic method, and a singular value decomposition method. Our preliminary moment calculations show a remarkable agreement within the uncertainties of the measurements, with the results obtained by the Cluster/PEACE electron spectrometers. The data analyzed was selected because it represented a potential reconnection event as currently published.

Figueroa-Vinas, Adolfo↗
NASA NTRS2008

Cluster/Peace Electrons Velocity Distribution Function: Modeling the Strahl in the Solar Wind

We present a study of kinetic properties of the strahl electron velocity distribution functions (VDF's) in the solar wind. These are used to investigate the pitch-angle scattering and stability of the population to interactions with electromagnetic (whistler) fluctuations. The study is based on high time resolution data from the Cluster/PEACE electron spectrometer. Our study focuses on the mechanisms that control and regulate the pitch-angle and stability of strahl electrons in the solar wind; mechanisms that are not yet well understood. Various parameters are investigated such as the electron heat-flux and temperature anisotropy. The goal is to check whether the strahl electrons are constrained by some instability (e.g., the whistler instability), or are maintained by other types of processes. The electron heat-flux and temperature anisotropy are determined by fitting the VDF's to a spectral spherical harmonic model from which the moments are derived directly from the model coefficients.

Figueroa-Vinas, Adolfo↗
NASA NTRS2009

Modeling the Strahl in the Solar Wind Using CLUSTER/PEACE Electron Distribution Functions

We present a study of kinetic properties of the strahl electron velocity distribution functions (VDF's) in the solar wind. These are used to investigate the pitch-angle scattering and stability of the population to interactions with electromagnetic (whistler) fluctuations. The study is based on high time resolution data from the Cluster/PEACE electron spectrometer. Our study focuses on the mechanisms that control and regulate the pitch-angle and stability of strahl electrons in the solar wind; mechanisms that are not yet well understood. Various parameters are investigated such as the strahl-electron density, temperature anisotropy, and electron heat-flux. The goal is to check whether the strahl electrons are constrained by some instability (e.g., the whistler instability), or are maintained by other types of processes. The electron heat-flux and temperature anisotropy are determined by modeling the 3D-VDFs to a spectral spherical harmonic model from which the moments are obtained directly from the spectral coefficients.

Vinas, Adolfo F.↗
NASA NTRS2010

A New Representative for Small-Body Gravity Estimation

One of the options proposed in the 2009 Review of the U. S. Human Spaceflight Plans Committee calls for the consideration of robotic and human missions to asteroids. Beginning in 2000, the NEAR Shoemaker spacecraft inserted into orbit about the asteroid 433 Eros, and studied it for a year before landing. The Hayabusa satellite journeyed to 25143 Itokawa, and operated in a proximity of the asteroid for several months before landing to collect a return sample. Additionally, the DAWN satellite, which was launched in 2007, will explore the dwarf-planet Ceres and the asteroid 4 Vesta in the next three years. These are just several examples of current interests in exploring asteroids. In this paper, we explore an alternative for estimating the gravity field of a small body using observations of an orbiting satellite. Our approach uses two new mathematical tools, the quadratures (on a sphere) invariant under the icosahedral group, and a multiresolution representation of the gravity potential. With the near optimal quadratures on the sphere, we minimize the number of parameters needed to recover a function. Also, instead of estimating the coefficients of the spherical harmonics, the new quadratures allow us to estimate directly the values of the gravity field since they offer an analogue of the Lagrange-type interpolation.

Jones, Brandon A.↗
NASA NTRS2010

Seismic Constraints on Interior Solar Convection

We constrain the velocity spectral distribution of global-scale solar convective cells at depth using techniques of local helioseismology. We calibrate the sensitivity of helioseismic waves to large-scale convective cells in the interior by analyzing simulations of waves propagating through a velocity snapshot of global solar convection via methods of time-distance helioseismology. Applying identical analysis techniques to observations of the Sun, we are able to bound from above the magnitudes of solar convective cells as a function of spatial convective scale. We find that convection at a depth of r/R(solar) = 0.95 with spatial extent l < 30, where l is the spherical harmonic degree, comprise weak flow systems, on the order of 15 m/s or less. Convective features deeper than r/R(solar) = 0.95 are more difficult to image due to the rapidly decreasing sensitivity of helioseismic waves.

Hanasoge, Shravan M.↗
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.↗
NASA NTRS2011

Time-Distance Analysis of Deep Solar Convection

Recently it was shown by Hanasoge, Duvall, and DeRosa (2010) that the upper limit to convective flows for spherical harmonic degrees l<? is considerably smaller than the flows predicted by the ASH simulations (Miesch et a7. ref) at the depth r/R=0.95 ' The deep-focusing Lime-distance technique used to develop the upper limit was applied to linear acoustic simulations of a solar interior perturbed by convective flows in order to calibrate the technique. This technique has been applied to other depths in the convection zone and the results will be presented. The deep-focusing technique has considerable sensitivity to the flow ' signals at the desired subsurface location ' However, as shown by Birch {ref}, there is remaining much sensitivity to near-surface signals. Modifications to the technique using multiple bounce signals have been examined in a search for a more refined sensitivity, or kernel function. Initial results are encouraging and results will be presented'

Duvall, T. L., Jr.↗
NASA NTRS2011

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

Effect of Cross-Correlation on Geomagnetic Forecast Accuracies

Surface geomagnetic observation can determine up to degree L = 14 time-varying spherical harmonic coefficients of the poloidal magnetic field. Assimilation of these coefficients to numerical dynamo simulation could help us understand better the dynamical processes in the Earth's outer core, and to provide more accurate forecast of geomagnetic secular variations (SV). In our previous assimilation studies, only the poloidal magnetic field in the core is corrected by the observations in the analysis. Unobservable core state variables (the toroidal magnetic field and the core velocity field) are corrected via the dynamical equations of the geodynamo. Our assimilation experiments show that the assimilated core state converges near the CMB, implying that the dynamo state is strongly constrained by surface geomagnetic observations, and is pulled closer to the truth by the data. We are now carrying out an ensemble of assimilation runs with 1000 years of geomagnetic and archeo/paleo magnetic record. In these runs the cross correlation between the toroidal and the poloidal magnetic fields is incorporated into the analysis. This correlation is derived from the physical boundary conditions of the toroidal field at the core-mantle boundary (CMB). The assimilation results are then compared with those of the ensemble runs without the cross-correlation, aiming at understanding two fundamental issues: the effect of the crosscorrelation on (1) the convergence of the core state, and (2) the SV prediction accuracies. The constrained dynamo solutions will provide valuable insights on interpreting the observed SV, e.g. the near-equator magnetic flux patches, the core-mantle interactions, and possibly other geodynamic observables.

Kuang, Weijia↗
NASA NTRS2011

Simulations of Recovery of Time-Varying Gravity from DECIGO Pathfinder

We simulated time-varying Earth's gravity field recovered from DPF to evaluate an impact of DPF and future satellite gradiometry mission on earth science. From hydrological water movement data and orbit information, gravity gradients to be measured at altitude about ~500km were generated. Errors caused by atmospheric and oceanic variations and instrumental noise were added. Monthly gravity fields were estimated solving normal equations between spherical harmonic coefficients and simulated gravity gradient data. Simulation results show that DPF likely provides monthly hydrological water storage change with spatial scale between 400 and 1000km. Sensitivities to large scale estimates depends on long-term stability of gravity gradient measurement, and errors in short scale estimates are caused by instrumental noise and imperfections in atmospheric and ocean model. With acceleration noise level is lower than ~5 x 10(exp -14) [m/s2/sqrtHz] at frequency higher than 3mHz, water storage changes at limited small basins will be provided by DPF. To monitor continental scale hydrological water movement, noise level must be lower than ~5 x 10(exp -14) [m/s2/sqrtHz] at frequency higher than 1mHz.

Hasegawa, Takashi↗
NASA NTRS2011

Toward Realistic Dynamics of Rotating Orbital Debris, and Implications for Lightcurve Interpretation

Optical observations of rotating space debris near GEO contain important information on size, shape, composition, and rotational states, but these aspects are difficult to extract due to data limitations and the high number of degrees of freedom in the modeling process. For tri-axial rigid debris objects created by satellite fragmentations, the most likely initial rotation state has a large component of initial angular velocity directed along the intermediate axis of inertia, leading to large angular reorientations of the body on the timescale of the rotation period. This lends some support to the simplest possible interpretation of light curves -- that they represent sets of random orientations of the objects of study. However, effects of internal friction and solar radiation are likely to cause significant modification of rotation states within a time as short as a few orbital periods. In order to examine the rotational dynamics of debris objects under the influences of these effects, a set of seven first-order coupled equations of motion were assembled in state form: three are Euler equations describing the rates of change of the components of angular velocity in the body frame, and four describe the rates of change of the components of the unit quaternion. Quaternions are a four-dimensional extension of complex numbers that form a seamless, singularity-free representation of body orientation on S3. The Euler equations contain explicit terms describing torque from solar radiation in terms of spherical harmonics, and terms representing effects of a prescribed rate of internal friction. Numerical integrations of these equations of motion are being performed, and results will be presented. Initial tests show that internal friction without solar radiation torque leads to rotation about the maximum principal axis of inertia, as required, and solar radiation torque is expected to lead to spin-up of objects. Because the axis of maximum rotational inertia tends to be roughly coincident with the normal to the largest projected cross-sectional area, internal friction is expected to lead to reduced variation of light curve amplitudes at a given phase angle, but a large dependence of the same on phase angle. At a given phase angle, databases are generated which contain reflected intensities for comprehensive sets of equally-likely orientations, represented as unit quaternions. When projected onto three dimensions (S2) and color-coded by intensity, the set is depicted as points within a solid, semi-transparent unit sphere, within which all possible reflected intensities for an object at a given phase angle may be inspected simultaneously. Rotational sequences are represented by trajectories through the sphere. Databases are generated for each of a set of phase angles separately, forming a comprehensive dataset of reflected intensities spanning all object orientations and solar phase angles. Symmetries in the problem suggest that preferred rotation states are likely, defined relative to the object-sun direction in inertial space and relative to the maximum principal axis of inertia in the body coordinate system. Such rotation states may greatly simplify the problem of light curve interpretation by reducing the number of degrees of freedom in the problem.

Ojakangas, Gregory W.↗
NASA NTRS2011

Investigation of the Crust of the Pannonian Basin, Hungary Using Low-Altitude CHAMP Horizontal Gradient Magnetic Anomalies

The Pannonian Basin is a deep intra-continental basin that formed as part of the Alpine orogeny. It is some 600 by 500 km in area and centered on Hungary. This area was chosen since it has one of the thinnest continental crusts in Europe and is the region of complex tectonic structures. In order to study the nature of the crustal basement we used the long-wavelength magnetic anomalies acquired by the CHAMP satellite. The SWARM constellation, scheduled to be launched next year, will have two lower altitude satellites flying abreast, with a separation of between ca. 150 to 200 km. to record the horizontal magnetic gradient. Since the CHAMP satellite has been in orbit for eight years and has obtained an extensive range of data, both vertically and horizontally there is a large enough data base to compute the horizontal magnetic gradients over the Pannonian Basin region using these many CHAMP orbits. We recomputed a satellite magnetic anomaly map, using the spherical-cap method of Haines (1985), the technique of Alsdorf et al. (1994) and from spherical harmonic coefficients of MF6 (Maus et aI., 2008) employing the latest and lowest altitude CHAMP data. We then computed the horizontal magnetic anomaly gradients (Kis and Puszta, 2006) in order to determine how these component data will improve our interpretation and to preview what the SW ARM mission will reveal with reference to the horizontal gradient anomalies. The gradient amplitude of an 1000 km northeast-southwest profile through our horizontal component anomaly map varied from 0 to 0.025 nT/km with twin positive anomalies (0.025 and 0.023 nT/km) separated by a sharp anomaly negative at o nT/km. Horizontal gradient indicate major magnetization boundaries in the crust (Dole and Jordan, 1978 and Cordell and Grauch, 1985). Our gradient anomaly was modeled with a twodimensional body and the anomaly, of some 200 km, correlates with a 200 km area of crustal thinning in the southwestern Pannonian Basin.

Taylor, Patrick T.↗
NASA NTRS2010

Lunar Gravity Field Determination Using SELENE Same-Beam Differential VLBI Tracking Data

A lunar gravity field model up to degree and order 100 in spherical harmonics, named SGM 100i, has been determined from SELENE and historical tracking data, with an emphasis on using same-beam S-band differential VLBI data obtained in the SELENE mission between January 2008 and February 2009. Orbit consistency throughout the entire mission period of SELENE as determined from orbit overlaps for the two sub-satellites of SELENE involved in the VLBI tracking improved consistently from several hundreds of metres to several tens of metres by including differential VLBI data. Through orbits that are better determined, the gravity field model is also improved by including these data. Orbit determination performance for the new model shows improvements over earlier 100th degree and order models, especially for edge-on orbits over the deep far side. Lunar Prospector orbit determination shows an improvement of orbit consistency from I-day predictions for 2-day arcs of 6 m in a total sense, with most improvement in the along and cross-track directions. Data fit for the types and satellites involved is also improved. Formal errors for the lower degrees are smaller, and the new model also shows increased correlations with topography over the far side. The estimated value for the lunar GM for this model equals 4902.80080 +/- 0.0009 cu km/sq s (10 sigma). The lunar degree 2 potential Love number k2 was also estimated, and has a value of 0.0255 +/- 0.0016 (10 sigma as well).

Goossens, S.↗
NASA NTRS2011

Improved Systematic Pointing Error Model for the DSN Antennas

New pointing models have been developed for large reflector antennas whose construction is founded on elevation over azimuth mount. At JPL, the new models were applied to the Deep Space Network (DSN) 34-meter antenna s subnet for corrections of their systematic pointing errors; it achieved significant improvement in performance at Ka-band (32-GHz) and X-band (8.4-GHz). The new models provide pointing improvements relative to the traditional models by a factor of two to three, which translate to approximately 3-dB performance improvement at Ka-band. For radio science experiments where blind pointing performance is critical, the new innovation provides a new enabling technology. The model extends the traditional physical models with higher-order mathematical terms, thereby increasing the resolution of the model for a better fit to the underlying systematic imperfections that are the cause of antenna pointing errors. The philosophy of the traditional model was that all mathematical terms in the model must be traced to a physical phenomenon causing antenna pointing errors. The traditional physical terms are: antenna axis tilts, gravitational flexure, azimuth collimation, azimuth encoder fixed offset, azimuth and elevation skew, elevation encoder fixed offset, residual refraction, azimuth encoder scale error, and antenna pointing de-rotation terms for beam waveguide (BWG) antennas. Besides the addition of spherical harmonics terms, the new models differ from the traditional ones in that the coefficients for the cross-elevation and elevation corrections are completely independent and may be different, while in the traditional model, some of the terms are identical. In addition, the new software allows for all-sky or mission-specific model development, and can utilize the previously used model as an a priori estimate for the development of the updated models.

Rochblatt, David J.↗