Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spherical harmonic transform”

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 19 records

Tests of a simple data merging algorithm for the GONG project

The GONG (Global Oscillation Network Group) project proposes to reduce the impact of diurnal variations on helioseismic measurements by making long-term observations of solar images from six sites placed around the globe. The sun will be observed nearly constantly for three years, resulting in the acquisition of l+ terabyte of image data. To use the solar network to maximum advantage, the images from the sites must be combined into a single time series to determine mode frequencies, amplitudes, and line widths. Initial versions of combined, i.e., merged, time series were made using a simple weighted average of data from different sites taken simultaneously. In order to accurately assess the impact of the data merge on the helioseismic measurements, a set of artificial solar disk images was made using a standard solar model and containing a well known set of oscillation modes and frequencies. This undegraded data set and data products computed from it were used to judge the relative merits of various data merging schemes. The artificial solar disk images were subjected to various instrumental and atmospheric degradations, dependent on site and time, in order to create a set of images simulating those likely to be taken at the site. The degraded artificial solar disk images for the six observing sites were combined in various ways to form merged time series of images and mode coefficients. Various forms of a weighted average were used, including an equally-weighted average, an average with weights dependent upon air mass and averages with weights dependent on various quality assurance parameters. Both the undegraded solar disk image time series and several time series made up of various combinations of the degraded solar disk images from the six sites were subjected to standard helioseismic measurement processing. This processing consisted of coordinate remapping, detrending, spherical harmonic transformation, computation of power series for the oscillation mode coefficients, and mode frequency identification. Visual and statistical evaluation of the merged data sets themselves and differences between the merged and undegraded data set shows good agreement between the two data sets. Some slight differences in image scale and registration appear between the undegraded data set and the various merged data sets. In the set of power series made from the mode coefficients of the merged data sets, some power leakage is observed into the background and into slightly lower l-value modes, especially at higher l-value mode frequencies. The results of the comparison of time series and mode oscillation frequencies of the undegraded data with those of the data merged using weighted averages indicate that, at least for p-mode solar oscillations, a weighted average of either the detrended remapped images or the mode coefficients gives good determination of the mode frequencies and adequate-to-good determination of the amplitudes and widths of the mode frequency lines. This conclusion is most advantageous to the analysis of the massive amounts of data to be received by the GONG network.

Williams, W. E.↗

An analytic solution for the orbital perturbations of the Venus Radar Mapper due to gravitational harmonics

Hill's variational equations are solved analytically for the orbital perturbations of a spacecraft nominally in an elliptic orbit around a non-spherical body. The rotation of the central planet about its spin-axis is not considered in the analysis. The perturbations are restricted to the planetary gravitational harmonics only. An extremely simple algorithm is derived to transform the spherical harmonic potentials to the orbital coordinate system, and the resulting accelerations are shown to be simply trigonometric functions of the true anomaly. With the principal matrix solution for the differential equations of the adjoint system given in closed form, the orthogonality of the trigonometric functions makes it possible to obtain an analytic solution for the non-homogeneous problem, at intervals of 2 pi in true anomaly. The solution for orbital perturbations can be extended over several revolutions by applying well-known results from Floquet's theory. The technique is demonstrated with results presented on the spacecraft periapsis altitude for the forthcoming Venus Radar Mapper Mission.

Vijayaraghavan, A.↗

Retrieving the Balanced Winds on the Globe as a Generalized Inverse Problem

A generalized inverse technique is applied to retrieve two types of balanced winds that characterize the large-scale dynamics of the atmosphere: rotational winds based upon the linear balance equation, and divergent winds based upon the vorticity budget equation. Both balance equations are singular at or near the equator. The balance equations are transformed in spherical harmonic function space to an under-determined system, for which the scale-weighed least-squares solution consists of a sum of principal and singular components. The principal components represent the response to the source function for the regular eigenmodes, while the singular components are determined by the projection of an independent measurement on the singular eigenmodes. The method was tested with the NCEP/NCAR reanalysis data in which a quasi-balance condition exists. A realistic balanced wind field is retrievable when the singular components are computed based upon the reanalyzed wind data.

Lu, Huei-Iin↗

Sampling functions for geophysics

A set of spherical sampling functions is defined such that they are related to spherical-harmonic functions in the same way that the sampling functions of information theory are related to sine and cosine functions. An orderly distribution of (N + 1) squared sampling points on a sphere is given, for which the (N + 1) squared spherical sampling functions span the same linear manifold as do the spherical-harmonic functions through degree N. The transformations between the spherical sampling functions and the spherical-harmonic functions are given by recurrence relations. The spherical sampling functions of two arguments are extended to three arguments and to nonspherical reference surfaces. Typical applications of this formalism to geophysical topics are sketched.

Giacaglia, G. E. O.↗

The Impact of Dimensionality Reduction of Ion Counts Distributions on Preserving Moments, With Applications to Data Compression

The field of space physics has a long history of utilizing dimensionality reduction methods to distill data, including but not limited to spherical harmonics, the Fourier Transform, and the wavelet transform. Here, we present a technique for performing dimensionality reduction on ion counts distributions from the Multiscale Mission/Fast Plasma Investigation (MMS/FPI) instrument using a data-adaptive method powered by neural networks. This has applications to both feeding low-dimensional parameterizations of the counts distributions into other machine learning algorithms, and the problem of data compression to reduce transmission volume for space missions. The algorithm presented here is lossy, and in this work, we present the technique of validating the reconstruction performance with calculated plasma moments under the argument that preserving the moments also preserves fluid-level physics, and in turn a degree of scientific validity. The method presented here is an improvement over other lossy compressions in loss-tolerant scenarios like the Multiscale Mission/Fast Plasma Investigation Fast Survey or in non-research space weather applications.

D. da Silva↗

Equipotential transformation of multipole systems to dipole systems

The vector sum of fields produced by a set of magnetic dipoles with different magnetic moments that are both tilted from the planetary spin axis and offset from the planetary center by different amounts can be used to completely represent a planetary magnetic field in the conventional form of spherical harmonic expansion. The scalar equipotential transformation analytically results in 24 equations that may subsequently be solved for the 24 adjustable parameters in dipole systems with the predetermined main dipole. Attention is given to this method's application to the Jovian magnetic field.

Wang, W. X.↗

The Hadley and Rossby regimes in a spherical atmosphere

The properties of the steady Hadley and Rossby regimes for a thermally forced rotating fluid on a sphere are studied. The two layer modified geostrophic model is employed which allows for thermal advection by the divergent wind and time dependent static stability. Heating processes are parameterized using the Newtonian approximation and Rayleigh friction is accounted for. The equations are transformed to spectral form using spherical harmonics and then truncated retaining a simple axisymmetric state and initial, one wave. A time independent Hadley circulation is obtained which is neutral to axisymmetric disturbances but unstable to wave like perturbations for intermediate values of the meridional temperature gradient, indicating the existence of both an upper and lower symmetric Hadley regime. An analytical solution for the steady Rossby circulation is determined for values of the meridional temperature gradient where the Hadley regime is unstable. Linear perturbation theory is used to show that within the steady Rossby regime two or more waves cannot exist simultaneously.

Feldstein, S. B.↗

Spectral Spherical Harmonics Discrete Ordinate Method

A new method for modeling the radiative transfer in inhomogeneous three-dimensional media illuminated by a Gaussian beam is described. This approach, called the Spectral Spherical Harmonics Discrete Ordinate Method (SSHDOM), uses the Fourier expansion method to transform the three-dimensional radiative transfer into an one-dimensional equation in the spectral domain, and the Spherical Harmonics Discrete Ordinate Method (SHDOM) for its solution. Specifically, (i) the source function is represented in the spectral domain through a spherical harmonic expansion, (ii) the spectral one-dimensional radiative transfer equation is integrated along discrete ordinates through a spatial grid, and (iii) the solution method is based on the Picard iteration. Both SSHDOM and SHDOM algorithms are implemented in a common computer code.

Gaussian beam↗

Geopotential Field Anomaly Continuation with Multi-Altitude Observations

Conventional gravity and magnetic anomaly continuation invokes the standard Poisson boundary condition of a zero anomaly at an infinite vertical distance from the observation surface. This simple continuation is limited, however, where multiple altitude slices of the anomaly field have been observed. Increasingly, areas are becoming available constrained by multiple boundary conditions from surface, airborne, and satellite surveys. This paper describes the implementation of continuation with multi-altitude boundary conditions in Cartesian and spherical coordinates and investigates the advantages and limitations of these applications. Continuations by EPS (Equivalent Point Source) inversion and the FT (Fourier Transform), as well as by SCHA (Spherical Cap Harmonic Analysis) are considered. These methods were selected because they are especially well suited for analyzing multi-altitude data over finite patches of the earth such as covered by the ADMAP database. In general, continuations constrained by multi-altitude data surfaces are invariably superior to those constrained by a single altitude data surface due to anomaly measurement errors and the non-uniqueness of continuation.

Kim, Jeong Woo↗

Geopotential Field Anomaly Continuation with Multi-Altitude Observations

Conventional gravity and magnetic anomaly continuation invokes the standard Poisson boundary condition of a zero anomaly at an infinite vertical distance from the observation surface. This simple continuation is limited, however, where multiple altitude slices of the anomaly field have been observed. Increasingly, areas are becoming available constrained by multiple boundary conditions from surface, airborne, and satellite surveys. This paper describes the implementation of continuation with multi-altitude boundary conditions in Cartesian and spherical coordinates and investigates the advantages and limitations of these applications. Continuations by EPS (Equivalent Point Source) inversion and the FT (Fourier Transform), as well as by SCHA (Spherical Cap Harmonic Analysis) are considered. These methods were selected because they are especially well suited for analyzing multi-altitude data over finite patches of the earth such as covered by the ADMAP database. In general, continuations constrained by multi-altitude data surfaces are invariably superior to those constrained by a single altitude data surface due to anomaly measurement errors and the non-uniqueness of continuation.

Kim, Jeong Woo↗

Converting gravity bins to spherical harmonic coefficients

The gravity bin technique as originally formulated recovers the local gravity field from the bin parameters by finite differencing. The spherical harmonic coefficients of the gravity field are then computed by an orthogonal transformation of the local gravity field. The result differs from that of the traditional method. This paper discusses the difference and proposes a new algorithm to convert the bin parameters to spherical harmonic coefficients. It is shown that the new method produces the same gravity field as the traditional method and maintains the high computational efficiency of the basic gravity bin technique.

Wu, J. T.↗

Second-order Compton-Getting effect on arbitrary intensity distribution

Theoretical studies of energetic particles in space are often referred to a special frame of reference. To compare theory with experiment, one has to transform the particle distribution from the special frame to the observer's frame, or vice versa. Various methods have been derived to obtain the directional distribution in the comoving frame from the directional fluxes measured on a spacecraft. These methods have become progressively complicated as increasingly detailed directional particle data become available. A set of 2nd order correct formulae for the transformation of an arbitrary differential intensity distribution, expressed as a series of spherical harmonics, between any two frames in constant relative motion is presented. These formulae greatly simplify the complicated procedures currently in use for the determination of the differential intensity distribution in a comoving frame.

Ng, C. K.↗

A Spectral Algorithm for Solving the Relativistic Vlasov-Maxwell Equations

A spectral method algorithm is developed for the numerical solution of the full six-dimensional Vlasov-Maxwell system of equations. Here, the focus is on the electron distribution function, with positive ions providing a constant background. The algorithm consists of a Jacobi polynomial-spherical harmonic formulation in velocity space and a trigonometric formulation in position space. A transform procedure is used to evaluate nonlinear terms. The algorithm is suitable for performing moderate resolution simulations on currently available supercomputers for both scientific and engineering applications.

Shebalin, John V.↗

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