SEARCH · Search NASA

Results for “generalized spherical harmonics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Statistics of the geomagnetic secular variation for the past 5 m.y

A new statistical model is proposed for the geomagnetic secular variation over the past 5Ma. Unlike previous models, the model makes use of statistical characteristics of the present day geomagnetic field. The spatial power spectrum of the non-dipole field is consistent with a white source near the core-mantle boundary with Gaussian distribution. After a suitable scaling, the spherical harmonic coefficients may be regarded as statistical samples from a single giant Gaussian process; this is the model of the non-dipole field. The model can be combined with an arbitrary statistical description of the dipole and probability density functions and cumulative distribution functions can be computed for declination and inclination that would be observed at any site on Earth's surface. Global paleomagnetic data spanning the past 5Ma are used to constrain the statistics of the dipole part of the field. A simple model is found to be consistent with the available data. An advantage of specifying the model in terms of the spherical harmonic coefficients is that it is a complete statistical description of the geomagnetic field, enabling us to test specific properties for a general description. Both intensity and directional data distributions may be tested to see if they satisfy the expected model distributions.

Constable, C. G.↗

NASA NTRS2004

Magnetic Probing of Core Geodynamics

To better understand geomagnetic theory and observation, we can use spatial magnetic spectra for the main field and secular variation to test core dynamical hypotheses against seismology. The hypotheses lead to theoretical spectra which are fitted to observational spectra. Each fit yields an estimate of the radius of Earth's core and uncertainty. If this agrees with the seismologic value, then the hypothesis passes the test. A new way to obtain theoretical spectra extends the hydromagnetic scale analysis of Benton to scale-variant field and flow. For narrow scale flow and a dynamically weak field by the top of Earth's core, this yields a generalized Stevenson-McLeod spectrum for the core-source field, and a secular variation spectrum modulated by a cubic polynomial in spherical harmonic degree n. The former passes the tests. The latter passes many tests, but does not describe rapid dipole decline and quadrupole rebound; some tests suggest it is a bit hard, or rich in narrow scale range. In a core geodynamo, motion of the fluid conductor does work against the Lorentz force. This converts kinetic into magnetic energy which, in turn, is lost to heat via Ohmic dissipation. In the analysis at length-scale 1/k, if one presumes kinetic energy is converted in either eddy-overturning or magnetic free-decay time-scales, then Kolmogorov or other spectra in conflict with observational spectra can result. Instead, the rate work is done roughly balances the dissipation rate, which is consistent with small-scale flow. The conversion time-scale depends on dynamical constraints. These are summarized by the magnetogeostrophic vertical vorticity balance by the top of the core, which includes anisotropic effects of rotation, the magnetic field, and the core-mantle boundary. The resulting theoretical spectra for the core-source field and its SV are far more compatible with observation. The conversion time-scale of order 120 years is pseudo-scale-invariant. Magnetic spectra of other planets may differ; however, if a transition to non-conducting fluid hydrogen in Jupiter acts as a barrier to vertical flow, as well as current, then the shape of the jovi-magnetic spectrum could be remarkably Earth-like.

Voorhies, Coerte V.↗

NASA NTRS2004

Magnetic Probing of Core Geodynamics

To better understand geomagnetic theory and observation, we can use spatial magnetic spectra for the main field and secular variation to test core dynmcal hypotheses against seismology. The hypotheses lead to theoretical spectra which are fitted to observational spectra. Each fit yields an estimate of the radius of Earth's core and uncertainty. If this agrees with the seismologic value, then the hypothes pass the test. A new way to obtain theoretical spectra extends the hydromagnetic scale analysis of Benton to scale-variant field and flow. For narrow scale flow and a dynamically weak field by the top of Earth's core, this yields a generalized Stevenson-McLeod spectrum for the core-source field, and a secular variation spectrum modulated by a cubic polynomial in spherical harmonic degree n. The former passes the tests. The latter passes many tests, but does not describe rapid dipole decline and quadrupole rebound; some tests suggest it is a bit hard, or rich in narrow scale change. In a core geodynamo, motion of the fluid conductor does work against the Lorentz force. This converts kinetic into magnetic energy which, in turn, is lost to heat via Ohmic dissipation. In the analysis at lentgh-scale l/k, if one presumes kinetic energy is converted in either eddy- overturning or magnetic free-decay time-scales, then Kolmogorov or other spectra in conflict with observational spectra can result. Instead, the rate work is done roughly balances the dissipation rate, which is consistent with small scale flow. The conversion time-scale depends on dynamical constraints. These are summarized by the magneto-geostrophic vertical vorticity balance by the top of the core, which includes anisotropic effects of rotation, the magnetic field, and the core- mantle boundary. The resulting theoretical spectra for the core-source field and its SV are far more compatible with observation. The conversion time-scale of order l20 years is pseudo-scale-invarient. Magnetic spectra of other planets may differ; however, if a transition to non-conductmg fluid hydrogen in Jupiter acts as barrier to vertical flow, as well as current, then the shape of the jovi-magnetic spectrum could be remarkably Earth-like.

Voorhies, Coerte V.↗

NASA NTRS1976

On the question of the energy of the precessional dynamo

The various estimates for the coupling mechanism by which precession transfers rotational, kinetic energy of earth into the energy of its magnetic field are generally considering hydromagnetic stresses that originate between mantle and core. Estimates of the energy of the geomagnetic field calculated from the data of spherical harmonic analysis derive precession energy values in reasonable agreement with the observed external energy of the geomagnetic field and with the rate of ohmic dissipation of energy in the core.

Dolginov, S. S.↗

DOE OSTI2025

Analytical expression of a finite, long, conical canted-cosine-theta coil for particle collider interaction regions

Magnets in the accelerator interaction region (IR) present significant challenges because of high field requirements and limited available space. Conical-shaped magnets offer advantages in these environments by allowing closer placement to the interaction point while maintaining clearance from synchrotron radiation. Interestingly, numerical studies have shown that conical canted-cosine-theta (CCT) designs produce a constant field distribution along the axial direction in the IR quadrupoles for the Electron-Ion Collider (EIC) at Brookhaven National Laboratory. However, the field harmonics generated by conical CCT windings are not yet fully understood. This paper presents an analytical approach to describe the magnetic field produced by a conical surface current and proposes a method for designing conical CCT magnets for accelerator applications. First, we begin with a surface current sheet having a general cosine-theta distribution in spherical coordinates and solve the vector potential using the Green’s function. The magnetic fields generated by the conical current sheet are expressed using associated Legendre polynomials. These results are then related to circular field harmonics and integral field harmonics for designing a coil that produces a pure multipole field. Next, a single layer of the conical CCT winding path is produced based on the cosine-theta current distribution. Finally, the magnetic field quality of dipole and quadrupole conical CCT coils with multiple layers is verified using the Biot-Savart law.

Yang, Ye↗

NASA NTRS1994

Radiation from advanced solid rocket motor plumes

The overall objective of this study was to develop an understanding of solid rocket motor (SRM) plumes in sufficient detail to accurately explain the majority of plume radiation test data. Improved flowfield and radiation analysis codes were developed to accurately and efficiently account for all the factors which effect radiation heating from rocket plumes. These codes were verified by comparing predicted plume behavior with measured NASA/MSFC ASRM test data. Upon conducting a thorough review of the current state-of-the-art of SRM plume flowfield and radiation prediction methodology and the pertinent data base, the following analyses were developed for future design use. The NOZZRAD code was developed for preliminary base heating design and Al2O3 particle optical property data evaluation using a generalized two-flux solution to the radiative transfer equation. The IDARAD code was developed for rapid evaluation of plume radiation effects using the spherical harmonics method of differential approximation to the radiative transfer equation. The FDNS CFD code with fully coupled Euler-Lagrange particle tracking was validated by comparison to predictions made with the industry standard RAMP code for SRM nozzle flowfield analysis. The FDNS code provides the ability to analyze not only rocket nozzle flow, but also axisymmetric and three-dimensional plume flowfields with state-of-the-art CFD methodology. Procedures for conducting meaningful thermo-vision camera studies were developed.

Farmer, Richard C.↗

NASA NTRS2004

Estimating Antarctic Near-Surface Magnetic Anomalies from Oersted and CHAMP Satellite Magnetometer Observations

Significant improvement in predicting near-surface magnetic anomalies can result from the highly accurate magnetic observations of the CHAMP satellite that is orbiting at about 400 km altitude. In general, regional magnetic signals of the crust are strongly masked by the core field and its secular variations due to wavelength coupling in the spherical harmonic representation and thus are difficult to isolate in the satellite measurements. However, efforts to isolate the regional lithospheric from core field components can exploit the correlations between the CHAMP magnetic anomalies and the pseudo magnetic effects inferred from gravity-derived crustal thickness variations. In addition, we can use spectral correlation theory to filter the static lithospheric field components from the dynamic external field effects. Employing these procedures, we processed the CHAMP magnetic conservations for an improved magnetic anomaly map of the Antarctic crust. Relative to the much higher altitude Oersted and noisier Magsat observations, CHAMP magnetic anomalies at 400 km altitude reveal new details on the effects of intra-crustal magnetic features and crustal thickness variations of the Antarctic. Moreover, these results greatly facilitate predicting magnetic anomalies in the regional coverage gaps of the ADMAP compilation of Antarctic magnetic anomalies from shipborne, airborne and ground surveys. Our analysis suggests that considerable new insights on the magnetic properties of the lithosphere may be revealed by a further order-of-magnitude improvement in the accuracy of the magnetometer.

vonFrese, Ralph R. B.↗

DOE OSTI2021

An exact inversion method for extracting orientation ordering by small-angle scattering

Here, we outline a nonparametric inversion strategy for determining the orientation distribution function (ODF) of sheared interacting rods using small-angle scattering techniques. With the presence of direct inter-rod interaction and fluid mechanical forces, the scattering spectra are no longer characterized by the azimuthal symmetry in the coordinates defined by the principal directions of simple shear conditions, which severely compounds the reconstruction of ODFs based on currently available methods developed for dilute systems. Using a real spherical harmonic expansion scheme, the real-space ODFs are uniquely determined from the anisotropic scattering spectra and their numerical accuracy is verified computationally. Our method can be generalized to extract ODFs of uniaxially anisotropic objects under different flow conditions in a properly transformed reference frame with suitable basis vectors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1985

Hopf bifurcation in the presence of symmetry

Group theory is applied to obtain generalized differential equations from the Hopf bifurcation theory on branching to periodic solutions. The conditions under which the symmetry group will admit imaginary eigenvalues are delimited. The action of the symmetry group on the circle group are explored and the Liapunov-Schmidt reduction is used to prove the Hopf theorem in the symmetric case. The emphasis is on simplifying calculations of the stability of bifurcating branches. The resulting general theory is demonstrated in terms of O(2) acting on a plane, O(n) in n-space, and O(3) and an irreducible model for spherical harmonics.

Golubitsky, M.↗

NASA NTRS1991

Use of non-spherical gravity harmonics for relative motion GN/C (category 2)

Detailed analysis of the Automatic Rendezvous and Capture problem indicate a need for three different regions of mathematical description for the GN&C algorithms: (1) multi-vehicle orbital mechanics to the rendezvous interface point, i.e., within 100 nm; (2) relative motion solutions (such as Clohessy-Wiltshire type) from the far-field to the near-field interface, i.e., within 1 nm and; (3) close proximity motion - the near-field motion where the relative differences in the gravitational and orbit inertial accelerations can be neglected from the equations of motion. Limit boundaries to these regions can be precisely defined by further analysis and will be functions of the tracking measurement accuracies and the computer resources available for the solution of the algorithms. This paper analyzes the relative motion in Regions 2 and 3 above and present the derivation and discussion of the general case of non-spherical gravitational perturbed relative motion. Mathematical deviations from the numerically integrated spherical gravity case and solutions from the Clohessy-Wiltshire equations are presented in the analysis. Based upon this preliminary analysis, it is recommended that further efforts be used to assess the relative position and velocity differences in Region 2 due to non-spherical gravity harmonics and that viable GN&C algorithms be developed to include these gravity perturbations (especially the effects of the first gravity harmonic, J2).

Henderson, David M.↗

DOE OSTI2022

Full-potential KKR within the removed-sphere method: A practical and accurate solution to the Poisson equation

Here an efficient and accurate generalization of the removed-sphere method (RSM) to solve the Poisson equation for total charge density in a solid with space-filling convex Voronoi polyhedra (VPs) and any symmetry is presented. The generalized RSM avoids the use of multipoles and VP shape functions for cellular integrals, which have associated ill-convergent large, double-internal L sums in spherical-harmonic expansions, so that fast convergence in single-L sums is reached. Our RSM adopts full Ewald formulation to work for all configurations or when symmetry breaking occurs, such as for atomic displacements or elastic constant calculations. The structure-dependent coefficients A L that define RSM can be calculated once for a fixed structure and speed up the whole self-consistent-field procedure. The accuracy and rapid convergence properties are confirmed using two analytic models, including the Coulomb potential and energy. We then implement the full-potential RSM using the Green's function Korringa-Kohn-Rostoker (KKR) method for real applications and compare the results with other first-principle methods and experimental data, showing that they are equally as accurate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2007

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

kinematics↗

NASA NTRS2013

Global Propagation of Gravity Waves Generated with the Whole Atmosphere Transfer Function Model

A brief review is presented of the Transfer Function Model (TFM) [e.g., Mayr et al., Space Science Reviews, 1990], which describes acoustic gravity waves (AGW) that propagate across the globe in a dissipative and static (no winds) background atmosphere with globally uniform temperature and density variations extending from the ground to 700 km. Unique among existing models, the TFM can be placed between the analytical approach on one end, and the rigorous numerical approach of general circulation models (GCM). The time consuming numerical integration of the conservation equations is restricted to compute the transfer function (TF) for a broad range of frequencies and spherical harmonics. Given TF, the atmospheric response for a chosen source configuration is then obtained in short order. Computationally efficient, the model is well suited to serve as experimental and educational tool for simulating propagating wave patterns across the globe. By design, the TFM is also semi-analytical and therefore well suited to explore the different wave modes that can be generated under different dynamical conditions.

Transfer Function Model↗

NASA NTRS1994

Euler potentials of current-free fields expressed in spherical harmonics

Given a magnetic field B = -del(vector differential operator)(sub gamma) with gamma expanded in spherical harmonics, it is shown that analytic Euler potentials may be derived for B if gamma is asymmetrical but contains only the contribution of a single index n. This work generalizes a result for sectorial harmonics with n = m, derived by Willis and Gardiner (1988).

Stern, David P.↗

NASA NTRS1986

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↗

DOE OSTI2020

SCEPTRE 2.1 Quick Start Guide

This report provides a summary of notes for building and running the Sandia Computational Engine for Particle Transport for Radiation Effects (SCEPTRE) code. SCEPTRE is a general purpose C++ code for solving the Boltzmann transport equation in serial or parallel using unstructured spatial finite elements, multigroup energy treatment, and a variety of angular treatments including discrete ordinates and spherical harmonics. Either the first-order form of the Boltzmann equation or one of the second-order forms may be solved. SCEPTRE requires a small number of open-source Third Party Libraries (TPL) to be available, and example scripts for building these TPL's are provided. The TPL's needed by SCEPTRE are Trilinos, boost, and netcdf. SCEPTRE uses an autoconf build system, and a sample configure script is provided. Running the SCEPTRE code requires that the user provide a spatial finite-elements mesh in Exodus format and a cross section library in a format that will be described. SCEPTRE uses an xml-based input, and several examples will be provided.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1974

A general circulation model of stratospheric ozone

The paper reports on a general circulation model being developed with the ultimate objective of assessing the effect of SST operations on the atmospheric ozone distribution. The model variables are represented in the spectral domain using 79 spherical harmonics and 26 vertical levels between the ground and 70 km. Initial calculations have attempted to simulate the seasonal variations of the unperturbed global ozone distribution. The model shows the presence of a westerly stratospheric jet in the winter hemisphere with amplitudes similar to those observed in the polar night jet. Easterlies dominate in the summer hemisphere. The model also shows ozone to be transported polewards and downwards. Furthermore the annual cycle in ozone concentrations at mid-latitudes has been simulated. Preliminary comparisons are also made with other features of the dynamics and the ozone distribution in the stratosphere.

Cunnold, D. M.↗

NASA NTRS1994

Gravity model development for TOPEX/POSEIDON: Joint gravity models 1 and 2

The TOPEX/POSEIDON (T/P) prelaunch Joint Gravity Model-1 (JGM-1) and the postlaunch JGM-2 Earth gravitational models have been developed to support precision orbit determination for T/P. Each of these models is complete to degree 70 in spherical harmonics and was computed from a combination of satellite tracking data, satellite altimetry, and surface gravimetry. While improved orbit determination accuracies for T/P have driven the improvements in the models, the models are general in application and also provide an improved geoid for oceanographic computations. The postlaunch model, JGM-2, which includes T/P satellite laser ranging (SLR) and Doppler orbitography and radiopositioning integrated by satellite (DORIS) tracking data, introduces radial orbit errors for T/P that are only 2 cm RMS with the commission errors of the marine geoid for terms to degree 70 being +/- 25 cm. Errors in modeling the nonconservative forces acting on T/P increase the total radial errors to only 3-4 cm root mean square (RMS), a result much better than premission goals. While the orbit accuracy goal for T/P has been far surpassed geoid errors still prevent the absolute determination of the ocean dynamic topography for wavelengths shorter than about 2500 km. Only a dedicated gravitational field satellite mission will likely provide the necessary improvement in the geoid.

Nerem, R. S.↗