Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spherical harmonic method”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

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↗
NASA NTRS2012

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↗
DOE OSTI2020

Filtered Backprojection in Compton Imaging using a Spherical Harmonic Wiener Filter with Pixelated CdZnTe

Filtered backprojection is an image reconstruction technique for Compton imaging that provides reasonably high resolution at much lower computational costs when compared to iterative methods. Here, this work applies a Wiener filter that has been derived for spherical harmonics on Compton imaging using the OrionUM pixelated CdZnTe imaging-spectrometer. To regularize the filter, an investigation is made into the power spectral density of the signal and noise to develop an appropriate spectral signal-to-noise ratio model for the restoration process. Experimental measurements were conducted with two 228Th sources placed 30° apart. The resulting filtered image of the two sources have an average full-width-at-half-maximum (FWHM) of 9.8° or 7.5° when using a mean squared error and structural similarity optimization approach respectively; an improvement from the 29.0° FWHM image when using simple backprojection.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1979

Gravitational spectra from direct measurements

A simple rapid method is described for determining the spectrum of a surface field (in spherical harmonics) from harmonic analysis of direct (in situ) measurements along great circle arcs. The method is shown to give excellent overall trends (smoothed spectra) to very high degree from even a few short arcs of satellite data. Three examples are taken with perfect measurements of satellite tracking over a planet made up of hundreds of point masses using (1) altimetric heights from a low-orbiting spacecraft, (2) velocity (range rate) residuals between a low and a high satellite in circular orbits, and (3) range rate data between a station at infinity and a satellite in a highly eccentric orbit. In particular, the smoothed spectrum of the earth's gravitational field is determined to about degree 400(50-km half wavelength) from 1 x 1 deg gravimetry and the equivalent of 11 revolutions of GEOS 3 and Skylab altimetry. This measurement shows that there is about 46 cm of geoid height (rms worldwide) remaining in the field beyond degree 180.

Wagner, C. A.↗
NASA NTRS1990

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

Equilibrium shape and location of a liquid drop acoustically positioned in a resonant rectangular chamber

The effect of a standing wave field in a rectangular chamber on the shape and location of an acoustically positioned drop or bubble is calculated. The sample deformation and equilibrium position are obtained from an analysis of the spherical harmonic projections of the total surface stress tensor. The method of calculation relies on the assumed condition that the sample is only slightly distorted from a spherical form. The equilibrium location of a levitated drop is combined with a formula introduced by Hasegawa (1979) to calcualte the ka dependence of the radiation force function. The present theory is valid for large as well as small ka values. Calculations in the small ka limit agree with previous theories and experimental results. Examples are presented for nonplane-wave modes as well as plane-wave rectangular modes.

Jackson, H. W.↗
DOE OSTI2024

Solidification and crystallographic texture modeling of laser powder bed fusion Ti-6Al-4V using finite difference-monte carlo method

Laser powder bed fusion (LPBF) additive manufacturing makes near-net-shaped parts with reduced material cost and time, rising as a promising technology to fabricate Ti-6Al-4V, a widely used titanium alloy in aerospace and medical industries. However, LPBF Ti-6Al-4V parts produced with 67° rotation between layers, a scan strategy commonly used to reduce microstructure and property inhomogeneity, have varying grain morphologies and weak crystallographic textures that change depending on processing parameters. Here, this study predicts LPBF Ti-6Al-4V solidification at three energy levels using a finite difference-Monte Carlo method and validates the simulations with large-area electron backscatter diffraction (EBSD) scans. The developed model accurately shows that a <001> texture forms at low energy and a <111> texture occurs at higher energies parallel to the build direction but with a lower strength than the textures observed from EBSD. A validated and well-established method of combining spatial correlation and general spherical harmonics representation of texture is developed to calculate a difference score between simulations and experiments. The quantitative comparison enables effective fine-tuning of nucleation density (N 0 ) input, which shows a nonlinear relationship with increasing energy level. Future improvements in texture prediction code and a more comprehensive study of N 0 with different energy levels will further advance the optimization of LPBF Ti-6Al-4V components. These developments contribute a novel understanding of crystallographic texture formation in LPBF Ti-6Al-4V, the development of robust model validation and calibration pipeline methodologies, and provide a platform for mechanical property prediction and process parameter optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

Limits on coronal reflection using high-frequency solar oscillations

Acoustic waves in the Sun with frequencies above about 5.3 mHz can propagate in the chromosphere. We examine imaged solar intensity data for evidence of reflection of these waves in the upper chromosphere, where the temperature increases by a large factor over a short distance. Our method is to compare the observed and theoretically derived frequency spacings between peaks in the power spectrum. We find that our theoretical frequencies provide the best fit to the data when the reflection in the upper atmosphere is eliminated. In particular, the model of Kumar (1993b), which includes the source depth, and radiative damping, in the calculation of power spectra but ignores chromospheric reflection, gives peak frequencies that are in good agreement with the observations. For acoustic waves of frequency greater than 6 mHz we put an upper limit to the reflectivity of chromosphere and corona, using our method, of about 10%. At a given spherical harmonic degree, the frequency spacing between peaks in the data generally decreases with increasing frequency, because the lower turning point of the waves is moving inward. However, between 5 and 5.5 mHz the frequency spacing increases slightly. This feature is probably associated with the acoustic cutoff frequency in the solar atmosphere, i.e., it indicates a transition from trapped waves to propagating waves. We are able to reproduce the observed behavior by a crude modeling of the solar atmosphere. Further study of these peaks should provide an independent way of exploring the mean structure of the solar atmosphere, particularly around the temperature minimum region.

Kumar, P.↗
NASA NTRS2015

Swarm Flyby Gravimetry

This study describes a new technology for discerning the gravity fields and mass distribution of a solar system small body, without requiring dedicated orbiters or landers. Instead of a lander, a spacecraft releases a collection of small, simple probes during a flyby past an asteroid or comet. By tracking those probes from the host spacecraft, one can estimate the asteroid's gravity field and infer its underlying composition and structure. This approach offers a diverse measurement set,equivalent to planning and executing many independent and unique flyby encounters of a single spacecraft. This report assesses a feasible hardware implementation, derives the underlying models,and analyzes the performance of this concept via simulation.In terms of hardware, a small, low mass, low cost implementation is presented, which consists of a dispenser and probes. The dispenser contains roughly 12 probes in a tube and has a total size commensurate with a 6U P-Pod. The probes are housed in disc shaped sabots. When commanded,the dispenser ejects the top-most probe using a linear motor. The ejected probe separates from its sabots and unfolds using internal springs. There are two types of probes, each designed for a particular tracking modality. The reflective probe type, tracked by a telescope, unfolds to forma diffusely reflective sphere. The retroreflector probe type, tracked by a lidar, unfolds to form a corner-cube retroreflector assembly. Both types are designed to spherical so that their attitude doesn't affect the spacecraft's tracking performance.This analysis indicates that the point-mass term of small bodies larger than roughly 500 m in diameter can be observed from a host spacecraft that tracks locally deployed probes throughout a flyby to an uncertainty of better than 5%. The conditions by which this measurement is possible depends on the characteristics of the asteroid (size, type), the flyby velocity, and the type of tracking available (angles-only or angles+ranging). For most encounters, a few (1-3) well placed probes can be very effective, with marginal improvement for additional probes. Given realistic deployment errors, an encounter may require roughly 10-12 probes to ensure that 1-3 achieve their target. Long duration tracking of probes flying by large asteroids (>5 km diameter) can sometimes provide observability of the gravity field's first spherical harmonic, J( sub 2). In summary, this method offers a feasible, affordable approach to enabling or augmenting flyby science.

Atchison, Justin↗
NASA NTRS1989

Bidirectional canopy reflectance and its relationship to vegetation characteristics

A method for determining the relationship between bidirectional canopy reflectance (CR) of a vegetated area and its characteristics (type and growth stage) is presented. In the method, the three-dimensional CR surface is expressed in terms of spherical harmonics which are weighted by appropriate coefficients, and the coefficients are analyzed by a cluster technique. The method is applied to four vegetation types (soybeans, corn, shinnery oak, and orchard grass) at various stages of growth.

Goel, Narendra S.↗
NASA NTRS2023

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↗
DOE OSTI2021

Crystal mechanics-based thermo-elastic constitutive modeling of orthorhombic uranium using generalized spherical harmonics and first-order bounding theories

In earlier works, a mathematical procedure for invertible microstructure-property linkages was developed using computationally efficient spectral methods for polycrystalline cubic and hexagonal metals. This paper formulates such invertible microstructure–property linkages for orthorhombic polycrystalline metals relying on the generalized spherical harmonics (GSH) spectral basis. The procedure is used to compute property closures of orthorhombic polycrystals. The closures represent the complete set of theoretically possible combinations of effective properties for a selected material. The procedure relies on the first-order bounding theories and considers orientation distribution functions (ODFs) as the main microstructural descriptor influencing homogenized properties. Numerous examples of these closures involving second-rank thermal expansion and fourth-rank elastic stiffness tensorial properties over a broad range of temperatures are presented for α-uranium (α-U). In doing so, certain key properties of these closures are exploited to facilitate their computation with drastically reduced computational effort. Along with the recently developed GSH-based interpolation procedure for ODFs from coarsely spaced experimental measurement grids to finely spaced finite element mesh resolution grids presented in Barrett et al., the developed computationally efficient ODF-effective property linkages are used to establish a crystal mechanics-based simulation framework coupled with the finite element method (FEM). The ODF dependent thermal expansion and elastic stiffness tensors are efficiently calculated at every integration point and used by the FEM to predict the overall distortion of a hemispherical part made of α-U during heating. In conclusion, it is shown that the developed framework can be used to simulate microstructurally heterogeneous components under thermo-mechanical loadings in a computationally efficient manner.

36 MATERIALS SCIENCE↗
NASA NTRS2019

Efficient computation of optimal low thrust gravity perturbed orbit transfers

We have developed a new method for solving low-thrust fuel-optimal orbit transfer problems in the vicinity of a large body (planet or asteroid), considering a high-fidelity spherical harmonic gravity model. The algorithm is formulated via the indirect optimization method, leading to a two-point boundary value problem (TPBVP). We make use of a hyperbolic tangent smoothing law for performing continuation on the thrust magnitude to reduce the sharpness of the control switches in early iterations and thus promote convergence. The TPBVP is solved using the method of particular solutions (MPS) shooting method and Picard-Chebyshev numerical integration. Application of Picard-Chebyshev integration affords an avenue for increased efficiency that is not available with step-by-step integrators. We demonstrate that computing the particular solutions with only a low-fidelity force model greatly increases the efficiency of the algorithm while ultimately achieving near machine precision accuracy. A salient feature of the MPS is that it is parallelizable, and thus further speedups are available. It is also shown that, for near Earth orbits and over a small number of en-route revolutions around the Earth, only the zonal perturbation terms are required in the costate equations to obtain a solution that is accurate to machine precision and optimal to engineering precision. The proposed framework can be used for trajectory design around small asteroids and also for orbit debris rendezvous and removal tasks.

Junkins, John L.↗
NASA NTRS1979

Combination solution for gravity field including altimetry

A previous global solution for the earth's gravitational field and station positions employed a combination of 21 weeks of Doppler tracking data from seven satellites, a complete set of 5-deg mean gravity anomalies, and geometric constraints from the worldwide satellite triangulation network on 17 of the stations. To this is now added geoid heights over the major part of the oceans obtained from GEOS 3 altimeter observations. The solution has a form based on the simple density-layer method and is equivalent to a solution in spherical harmonics to the tenth degree and order. The rms difference of 5-deg mean geoid heights of the new solution with the GEM 10 solution is 5.6 m.

Chovitz, B. 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

Venus' center of mass - center of figure displacement and implications

Earth, Moon, Mars, and Venus all have centers of mass (C.M.) that are displaced from their centers of figure (C.F.) by amounts which range from 340 meters (Venus) to 2.5 km (Mars). These offsets have all been calculated from the first degree terms in spherical harmonic expansions of topography. We describe an alternate method for calculating C.M. - C.F. offsets directly from a global topographic data set and apply it to Venus. Using Magellan altimetry, we find that Venus' C.F. is displaced approximately 280 meters from its C.M. in the direction of Western Aphrodite Terra (4.4 deg S, 135.8 deg E). We investigate several simple models for this offset and find that it is most consistent with thickened crust in Ovda and Thetis Regiones (which constitute most of W. Aphrodite). The location of the C.F. offset also places constraints on the degree of crustal thickening in Western Ishtar Terra and/or this highland's mode of origin. We favor a model in which offset due to thick crust in Western Ishtar Terra is balanced by an opposing offset due to cold, downwelling mantle material beneath the highland.

Bindschadler, D. L.↗
DOE OSTI2023

A finite element method for angular discretization of the radiation transport equation on spherical geodesic grids

Discrete ordinate (S N ) and filtered spherical harmonics (FP N ) based schemes have been proven to be robust and accurate in solving the Boltzmann transport equation but they have their own strengths and weaknesses in different physical scenarios. We present a new method based on a finite element approach in angle that combines the strengths of both methods and mitigates their disadvantages. The angular variables are specified on a spherical geodesic grid with functions on the sphere being represented using a finite element basis. A positivity-preserving limiting strategy is employed to prevent non-physical values from appearing in the solutions. Here, the resulting method is then compared with both S N and FP N schemes using four test problems and is found to perform well when one of the other methods fail.

79 ASTRONOMY AND ASTROPHYSICS↗