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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Clustering Algorithm for AM Parts using GSH and EDT with Autoencoder

SAND2025-10103O The Clustering Algorithm for AM Parts Using GSH and (EDT With Autoencoder is a software tool. It uses a clustering algorithm for additive manufacturing (AM) parts using generalized spherical harmonics (GSH) and Euclidean distance transform (EDT) with an autoencoder to quantify material microstructure. The tool offers improved sensitivity to microstructural changes compared to traditional approaches. The tool integrates multiple microstructural properties, such as grain morphology, crystallographic orientation, and material phase information, to provide a comprehensive analysis of material microstructures. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Rodgers, Theron [Sandia National Lab. (SNL-CA), Li↗
NASA NTRS1985

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

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

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

encore : an O ( N g2) estimator for galaxy N -point correlation functions

ABSTRACT We present a new algorithm for efficiently computing the N-point correlation functions (NPCFs) of a 3D density field for arbitrary N. This can be applied both to a discrete spectroscopic galaxy survey and a continuous field. By expanding the statistics in a separable basis of isotropic functions built from spherical harmonics, the NPCFs can be estimated by counting pairs of particles in space, leading to an algorithm with complexity $\mathcal {O}(N_\mathrm{g}^2)$ for Ng particles, or $\mathcal {O}(N_\mathrm{FFT}\log N_\mathrm{FFT})$ when using a Fast Fourier Transform with NFFT grid-points. In practice, the rate-limiting step for N > 3 will often be the summation of the histogrammed spherical harmonic coefficients, particularly if the number of radial and angular bins is large. In this case, the algorithm scales linearly with Ng. The approach is implemented in the encore code, which can compute the 3PCF, 4PCF, 5PCF, and 6PCF of a BOSS-like galaxy survey in ${\sim}100$ CPU-hours, including the corrections necessary for non-uniform survey geometries. We discuss the implementation in depth, along with its GPU acceleration, and provide practical demonstration on realistic galaxy catalogues. Our approach can be straightforwardly applied to current and future data sets to unlock the potential of constraining cosmology from the higher point functions.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2001

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

Quantifying Microstructure Variability in Laser Powder Bed Fusion 316 L Stainless Steel Microstructures with Spatial Statistics

Here, we have explored data-driven methods for material microstructure quantification that improve sensitivity to microstructural changes compared to traditional approaches. The methods integrate multiple microstructural properties, including grain morphology, crystallographic orientation, and material phase information. The simpler method employs maps of the Euclidean distance transformation metric to evaluate the morphology of grain boundary networks. The more intensive approach employs generalized spherical harmonic mapping for crystallographic orientations, per-pixel phase information, and a variational auto-encoder for dimensionality reduction and results in a multidimensional clustering of by microstructure similarity. Applied to an experimental dataset of additively manufactured steel, both methods detected slight variations in samples produced under nominally identical processing conditions. Both methods were able to distinguish between samples from multiple (nominally identical) builds, while the generalized spherical harmonics-based method could additionally cluster data samples rotated at two orientations on the build plate. The improved sensitivity of the methods, demonstrated through comparison with traditional microstructure characterization techniques, offers advantages for microstructure quantification and comparisons in advanced manufacturing applications.

SS316L↗
NASA NTRS1990

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

Radial deformation of the earth by oceanic tidal loading

A high-degree spherical harmonic series is used to compute the radial deformation of the Earth by oceanic tidal loading. By exploiting fast numerical transforms, this approach is found to be much more efficient, but no less accurate, than the traditional Green's function approach. The method is used to derive an atlas of load tide maps for 10 constitutents of the NSWC ocean tide model.

Ray, R. D.↗
NASA NTRS1974

Finite-amplitude thermal convection in a spherical shell

The properties of finite-amplitude thermal convection for a Boussinesq fluid contained in a spherical shell are investigated. All nonlinear terms are retained in the equations, and both axisymmetric and nonaxisymmetric solutions are studied. The velocity is expanded in terms of poloidal and toroidal vectors. Spherical surface harmonics resolve the horizontal structure of the flow, but finite differences are used in the vertical. With a few modifications, the transform method developed by Orszag (1970) is used to calculate the nonlinear terms, while Green's function techniques are applied to the poloidal equation and diffusion terms.

Young, R. E.↗
NASA NTRS2019

A New Model of the Crustal Magnetic Field of Mars Using MGS and MAVEN

While devoid of an active magnetic dynamo field today, Mars possesses a remanent magnetic field that may reach several thousand nanoteslas locally. The exact origin and the events that have shaped the crustal magnetization remain largely enigmatic. Three magnetic field data sets from two spacecraft collected over 13 cumulative years have sampled the Martian magnetic field over a range of altitudes from 90 up to 6,000 km: (a) Mars Global Surveyor (MGS) magnetometer (1997–2006), (b) MGS Electron Reflectometer (1999–2006), and (c) Mars Atmosphere and Volatile EvolutioN (MAVEN) magnetometer (2014 to today). In this paper we combine these complementary data sets for the first time to build a new model of the Martian internal magnetic field. This new model improves upon previous ones in several aspects: comprehensive data coverage, refined data selection scheme, modified modeling scheme, discrete-to-continuous transformation of the model, and increased model resolution. The new model has a spatial resolution of ∼160 km at the surface, corresponding to spherical harmonic degree 134. It shows small scales and well-defined features, which can now be associated with geological signatures.

Magnetic↗
NASA NTRS1976

Structure of the atmosphere of Jupiter from Pioneer 10 and 11 radio occultation measurements

A technique for incorporating effects of the oblateness of Jupiter's atmosphere into the data analysis procedure for radio-occultation measurements is described which makes use of a spherical harmonic representation of the gravity field to compute the shape of the planet. With this technique, the center of refraction is located by the radius of curvature and the normal direction at the closest approach point of the ray. The present technique, a subsequent approximation involving the use of a constant center of refraction, and the Abel integral transform inversion method for spherical planets are applied to Pioneer 10 and 11 data. All the intermediate results obtained show a temperature inversion between 10 and 100 mb, with temperatures from 130 to 170 K at 10 mb and from 80 to 120 K at 100 mb. Comparison of the radio-occultation profiles with radiative-convective equilibrium models and the temperature profile based on Pioneer 10 IR radiometer data indicates good agreement between the models and the occultation results.

Kliore, A. J.↗
NASA NTRS2013

A T Matrix Method Based upon Scalar Basis Functions

A surface integral formulation is developed for the T matrix of a homogenous and isotropic particle of arbitrary shape, which employs scalar basis functions represented by the translation matrix elements of the vector spherical wave functions. The formulation begins with the volume integral equation for scattering by the particle, which is transformed so that the vector and dyadic components in the equation are replaced with associated dipole and multipole level scalar harmonic wave functions. The approach leads to a volume integral formulation for the T matrix, which can be extended, by use of Green's identities, to the surface integral formulation. The result is shown to be equivalent to the traditional surface integral formulas based on the VSWF basis.

wave functions↗
DOE OSTI2025

Directional finite difference method for directly solving 3D gyrokinetic field equations with enhanced accuracy

The gyrokinetic (GK) field equation is a three-dimensional (3D) elliptic equation, but it is often simplified to a set of two-dimensional (2D) equations by assuming that the field does not vary along a specific direction. However, this simplification can introduce inevitable 0th-order numerical errors, as nonlinear mode coupling in toroidal geometry can produce undesirable harmonic modes that violate the assumption. In this work, we propose a novel directional finite difference method (FDM) with a local coordinate transformation to better resolve the target field of interest. The directional FDM can accurately solve 3D GK field equations without simplifications, which can overcome the limitations of conventional methods. The accuracy and efficiency of different FDMs are analyzed in great detail for a variety of geometries, from simple 2D Cartesian coordinates to realistic 3D curvilinear coordinates. The 0th-order numerical errors of simplified 2D GK equations were found to be more problematic for low-harmonic modes and low aspect ratio geometries such as spherical tokamaks. On the other hand, the directional 3D FDM can accurately resolve a much wider range of harmonic modes aligned to the direction of interest, including the low-harmonic modes. In conclusion, we demonstrate that the directional 3D FDM is a highly effective algorithm for solving the 3D GK field equations, achieving accuracy improvements of 10 to 100 times or more, particularly for low-harmonic modes in spherical tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Sequence modeling of higher-order wave modes of quasi-circular, spinning, non-precessing binary black hole mergers

Higher-order gravitational wave modes from quasi-circular, spinning, non-precessing binary-black-hole (BBH) mergers encode rich information about the nonlinear dynamics of strong-field gravity. We present a transformer-based sequence-completion surrogate that, given an early-inspiral segment, forecasts the subsequent late inspiral, merger, and ringdown. The intended applications are (i) patching or completing expensive or interrupted numerical-relativity (NR) simulations and (ii) providing late-time cross-checks and rapid hybridization studies. The training set is built from the NRHybSur3dq8 surrogate, which provides spherical-harmonic modes up to $\ell$ ≤ 4 (excluding (4, 0) and (4,±1), and including (5, 5)) for mass ratios q ≤ 8, dimensionless spin components s$^{z}_{1,2}$ ϵ[–0.8, 0.8], and inclination angles θ ϵ [0, π]. Waveforms are supplied on the interval t ϵ [–5000M, –100 M) and the model autoregressively generates the plus and cross polarizations (h + , h x ) on t ϵ [–100 M, 130M]. Training on the Delta supercomputer with 16 NVIDIA A100 GPUs required ~15 h on more than 14 million hybrid waveforms. Evaluation on a held-out test set of 840,000 samples yields mean and median overlaps of 0.996 and 0.997, respectively, with respect to the surrogate ground truth.

black-hole merger↗