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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
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 NTRS2017

B-plane Evolution Under Highly Non-Keplerian Dynamics

Juno is the second of a series of New Frontiers missions and was launched in 2011. The spacecraft was set en route to Jupiter, and its Jupiter Orbit Insertion occurred on 2016-07-04 after a five-year cruise in deep space. The mission phase just prior to this large maneuver is called approach. As navigation team, we were responsible for precise orbit determination of the Juno spacecraft to ensure a successful orbit insertion. To evaluate the navigation performance we employed the B-plane mapping, where the incoming hyperbolic velocity defines the plane perpendicular to it, and the time to go to hit the plane will directly map into the uncertainties of the target point. The B-plane is a convenient tool when assessing the navigation performance as it is not affected by non-linearities due to the gravitational pull of the targeted body. However, as the Jupiter Orbit Insertion approached we started to realize that the B-plane mapping uncertainties were significantly larger than expected. We found that the B-plane mapping time had tremendous influence on the covariance inflation in the B-plane. Because of the close approach distance during the Jupiter Orbit Insertion and the strong interaction with J2 spherical harmonic coefficient, the Juno dynamics were too far from a Keplerian one assumed for the B-plane mapping. We first discuss the analytical approach of the B-plane mapping uncertainty and perform a numerical analysis to isolate the components of the Juno dynamics that cause the B-plane covariance inflation.

Bordi, John↗
NASA NTRS2019

Application of GRACE to the Estimation of Groundwater Storage Change in a Data-poor Region: A Case Study of Ngadda Catchment in the Lake Chad Basin

The present study is to explore the feasibility of GRACE-based estimation of a groundwater storage change in a data-poor region using a case study of the Ngadda catchment in the Lake Chad Basin. Although the Ngadda catchment has only one set of in situ time series data of groundwater from 2006 to 2009 and a limited number of groundwater measurements in 2005 and 2009, GRACE-based groundwater storage change can be evaluated against thein situ groundwater measurements combined with specific yield data. The cross-correlation analysis in the Ngadda catchment shows that maximum rainfall reached in July and August, whereas both the maximum total water storage anomaly and the maximum groundwater storage anomaly occurred 2months later. Whereas the mean annual amplitude of total water storage anomaly is about 17cm from both the average total water storage anomaly from three mascon products and the one from three spherical harmonic products, the mean annual amplitude of soil moisture storage anomaly is substantially varied from5.58cm for CLM to about 14cm for NOAH and Mosaic. The goodness-of-fit tests show that CLM soil moisture produces the closest estimation of groundwater storage anomaly to thein situ groundwater measurements. The present study shows that GRACE-based estimation for groundwater storage anomaly can be a cost-effective and alternative tool to observe how groundwater changes in a basin scale under the limitation of modelling and in situ data availability.

GRACE↗
NASA NTRS2020

Julia Language 1.1 Ephemeris Reader and Gravitational Modeling Program for Solar System Bodies

This paper analyzes the advancements to the Julia Language 1.1 Ephemeris and Physical constants Reader including the addition of gravitational modeling. Originally written in MATLAB, this Julia Language program is intended to be used in for trajectory design. Written in an open-source coding language, this ephemeris reader can output the state of planetary bodies including asteroids as well as other constants such as gravitational parameters. Two primary methods were chosen to calculate the gravitational potentials which include polyhedral modeling and spherical harmonics.

Gray, Brennan↗
NASA NTRS2021

Search for anisotropic gravitational-wave backgrounds using data from Advanced LIGO and Advanced Virgo's first three observing runs

We report results from searches for anisotropic stochastic gravitational-wave backgrounds using data from the first three observing runs of the Advanced LIGO and Advanced Virgo detectors. For the first time, we include Virgo data in our analysis and run our search with a new efficient pipeline called PyStochon data folded over one sidereal day. We use gravitational-wave radiometry(broadband and narrow band) to produce sky maps of stochastic gravitational-wave backgrounds and to search for gravitational waves from point sources. A spherical harmonic decomposition method is employed to look for gravitational-wave emission from spatially-extended sources. Neither technique found evidence of gravitational-wave signals. Hence we derive 95% confidence-level upper limit sky maps on the gravitational-wave energy flux from broadband point sources, ranging from F(α,Θ) < (0.013−7.6)×10^(−8)erg/sq. cm s Hz,and on the (normalized) gravitational-wave energy density spectrum from extended sources, ranging from Ω(α,Θ) < (0.57−9.3)×10^(−9) per sr, depending on direction (Θ) and spectral index (α). These limits improve upon previous limits by factors of 2.9−3.5. We also set 95% confidence level upper limits on the frequency-dependent strain amplitudes of quasimonochromatic gravitational waves coming from three interesting targets, Scorpius X-1, SN1987A and the Galactic Center, with best upper limits range fromh(0) < (1.7−2.1)×10^(−25), a factor of ≥ 2.0 improvement compared to previous stochastic radiometer searches.

R. Abbott↗
NASA NTRS2021

High-Spatial-Resolution Mass Rates From GRACE and GRACE-FO: Global and Ice Sheet Analyses

GRACE and GRACE-FO gravimetry data are unique in its ability to observe global water mass variability and has frequently been applied for the determination of global mass rates to quantify the impacts of climate change and human activity on global freshwater availability and land ice evolution. However, its relatively coarse spatial resolution of ∼300-500 km at monthly intervals precludes its direct use to measure mass changes in smaller drainage basins and individual glacier systems. We present a new global high resolution mascon product that improves spatial resolution and signal recovery by estimating trends directly from more than 16 years of Level 1B GRACE and GRACE-FO data. We compare this product to independent gravimetry spherical harmonic trend estimates and to ice sheet mass rates derived from ICESat and ICESat-2 laser altimetry. We demonstrate enhanced global signal recovery and improved agreement with ICESat/ ICESat-2 and discuss likely causes for the remaining discrepancies between the gravimetry and laser altimetry solutions.

Bryant D. Loomis↗
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↗
NASA NTRS2024

Sensitivity studies of nighttime top-of-atmosphere radiances from artificial light sources using a 3-D radiative transfer model for nighttime aerosol retrievals

By accounting for surface-based light source emissions and top-of-atmosphere (TOA) downward lunar fluxes, we adapted the spherical harmonics discrete ordinate method (SHDOM) 3-dimensional (3-D) radiative transfer model (RTM) to simulate nighttime 3-D TOA radiances as observed from the Visible Infrared Imaging Radiometer Suite (VIIRS) Day/Night Band (DNB) on board the Suomi-NPP satellite platform. Used previously for daytime 3-D applications, these new SHDOM enhancements allow for the study of the impacts of various observing conditions and aerosol properties on simulated VIIRS-DNB TOA radiances. Observations over Dakar, Senegal, selected for its bright city lights and a large range of aerosol optical depth (AOD), were investigated for potential applications and opportunities for using observed radiances containing VIIRS-DNB “bright pixels” from artificial light sources to conduct aerosol retrievals. We found that using the standard deviation (SD) of such bright pixels provided a more stable quantity for nighttime AOD retrievals than direct retrievals from TOA radiances. Further, both the mean TOA radiance and SD of TOA radiances over artificial sources are significantly impacted by satellite viewing angles. Light domes, the enhanced radiances adjacent to artificial light sources, are strong functions of aerosol properties and especially aerosol vertical distribution, which may be further utilized for retrieving aerosol layer height in future studies. Through inter-comparison with both day- and nighttime Aerosol Robotic Network (AERONET) data, the feasibility of retrieving nighttime AODs using 3-D RTM SHDOM over artificial light sources was demonstrated. Our study shows strong potential for using artificial light sources for nighttime AOD retrievals, while also highlighting larger uncertainties in quantifying surface light source emissions. This study underscores the need for surface light emission source characterizations as a key boundary condition, which is a complex task that requires enhanced input data and further research. We demonstrate how quality-controlled nighttime light data from NASA’s Black Marble product suite could serve as a primary input into estimations of surface light source emissions for nighttime aerosol retrievals.

Jianglong Zhang↗
NASA NTRS2023

The Impact of New Accelerometer Transplant Data (ACH) on GRACE Follow-On Along-Orbit Inter-Satellite Laser Ranging Observations and Monthly Time-Variable Gravity and Mascon Solutions

GRACE-D accelerometer data show significant bias jumps since one month after the launch of the GRACE Follow-On (GRACE-FO) satellites in May 2018, making them inapplicable for correcting GRACE-FO products for non-gravitational accelerations. The GRACE-FO Science Data System (SDS) compensated this issue by transplanting the GRACE-C accelerometer data toward that of GRACE-D. Recently, GRACE-FO SDS implemented an updated transplant method, used in the latest release of GRACE-FO data. Here, we examine the impact of updated accelerometer transplant data (ACH) on GRACE-FO measurements at all levels: (a) Level-1B inter-satellite laser ranging residuals measured along satellite orbit, (b) Level-2 time-variable gravity solutions from all SDS centers (JPL, CSR, and GFZ), and (c) Level-3 mascon solutions. We show that inter-satellite laser ranging residuals are modified at low frequencies below 1 mHz, affecting the along-orbit analysis of large-scale time-variable gravity signals. When mapped into monthly Level-2 spherical harmonic coefficients of geopotential, the low-frequency change in inter-satellite ranging residuals leads to substantial improvement of coefficients associated with resonant orders (i.e., 15, 30, 45, etc.) and C 30 . We also present an improved SLR-derived C 30 which significantly improves the agreement with updated GRACE-FO C 30 at seasonal and interannual timescales. Moreover, we demonstrate the noise reduction in mass change estimates from new GRACE-FO Level-2 data over oceans, Greenland, and Antarctica for all SDS solutions. GRACE-FO mascon solutions show a moderate change in the updated release. Our comprehensive analyses demonstrate high-quality estimates of non-gravitational accelerations by the updated transplant method, resulting in more accurate GRACE-FO time-variable gravity and mass change observations.

Laser ranging interferometer↗
NASA NTRS2025

Assessment of Radiation Scattering by Dust in Mars Entry

The impact of gas phase radiation scattering by dust particle presence in Mars entry scenarios is assessed using a suite of numerical tools. A hypersonic gas and dilute particle multiphase flow is simulated using a mixed continuum model and an efficient size distribution quadrature scheme. Radiative transport is simulated with a hybrid ray-tracing and P1 spherical harmonics model. The NASA HARA radiation code and a Lorentz-Mie scattering model are used to predict the gas-particle mixture radiative properties. Application of the numerical suite to a representative case based on the Mars 2020 entry shows that, even at conditions in excess of the strongest recorded dust storm, the augmentation in surface radiative heat flux due to scattering alone is lower than 1%, and therefore negligible. Predictions with even higher dust loading conditions confirm that even highly uncertain dust mass profiles are unlikely to be a concern for the scattering augmentation mechanism.

Andrew Hinkle↗
DOE OSTI2022

Improvements to the Griffin Transport Solvers

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of steady-state solvers for fixed-source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. This document summarizes the transport solver development efforts pursued during Fiscal Year 2022. We added the multiphysics transient capability for the coarse-mesh finite difference accelerated Richardson iteration for discontinuous finite element method discrete ordinates (DFEM-SN) scheme to support high-order heterogeneous transport simulations. HFEM (hybrid finite element method) - PN (spherical harmonics expansion) was completed and red-black iteration was added for solving the HFEM-PN system with both preconditioned Jacobian-free Newton Krylov and Richardson iteration solvers. The HFEM-PN scheme, as one of the low-order transport schemes, is expected for supporting routine design simulations. Pin power reconstruction capability was also designed and implemented with the Griffin ISOXML module to enhance all the low-order transport solvers for more accurate multiphysics simulations. Numerical results are presented for demonstrating the capabilities and verifying their performance, and future works are discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Formulation of the density eigenvalue problem in neutron transport for relevant engineering applications

A new formulation of the density eigenvalue problem for the neutron transport equation is presented. This new eigenvalue, named ζ eigenvalue can be introduced freely in the transport model, acting on a selected portion of the phase space. Despite its broader applications and its connection with the nature of the transport operator, the ζ eigenvalue has been presented here mainly as a design-oriented technique for the efficient evaluation of the critical concentration for a specific nuclide (or mixture of nuclides). This new eigenvalue is particularly adequate to study the definition of the material composition in the criticality design process of a multiplying system. The method is then applied for the study of classical problems such as the critical moderation ratio and the poison concentration to control the reactor. Some results are presented in one dimensional configuration using the multi- group spherical harmonics approach. This eigenvalue formulation proves to be a convenient and useful way to attain criticality, also for complex, realistic systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Frequency-domain vs time-domain SP{sub N} equations to simulate neutron noise

Inside the reactor core mechanical vibrations of fuel assemblies can produce high fluctuations around a steady-state configuration, known as neutron noise. This effect can cause the triggering of power reduction measures. Classically, diffusion theory has been used to simulate this behavior. However, this equation has some limitations if the materials of the reactor have strong variations. In this work, we use the diffusive time-dependent simplified spherical harmonics equations that improve the previous results without the necessity of using high computational requirements. In particular, two types of analyses with these equations (SP{sub 3}) are made: a frequency-domain and a time-domain. A numerical neutron noise benchmark tests the methodology and compare both formulations. First, numerical results show a good agreement between the amplitudes and phases of the SP3 equations computed with the frequency-domain and time-domain. Therefore, as the frequency-domain computation only requires to solve a linear system, it is a recommendable option for neutron noise computations. Second, one can conclude that for this type of nuclear systems, where the assemblies are not homogenized, the SP{sub 3} approximation results improve considerably the accuracy of the diffusion theory.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

FolpsD: combining EFT and phenomenological approaches for joint power spectrum and bispectrum analyses

We present a theoretical model for the power spectrum and bispectrum of galaxy clustering that exploits the complementarity between small-scale power spectrum information and large-scale bispectrum measurements. We extend the FOLPS code by combining its one-loop EFT galaxy power spectrum with a tree-level galaxy bispectrum projected onto the tripolar spherical harmonics (Sugiyama) basis. To access additional small-scale information, we also consider a line-of-sight damping factor in both statistics, mirroring approaches commonly used in studies of redshift-space distortions. We test the model using DESI DR2 galaxy mocks. Even without damping, the joint analysis of the EFT power spectrum and bispectrum significantly improves constraints and reduces parameter degeneracies relative to power spectrum analyses alone. For LRG-like samples, including the damping further extends the range beyond $k\sim 0.3 \,h \text{Mpc}^{-1}$ in the power spectrum and $k \sim 0.24 \,h \text{Mpc}^{-1}$ in the bispectrum without introducing statistically significant parameter biases. This leads to up to $\sim 30\%$ tighter constraints on $A_s$ and $ω_{cdm}$. For low signal-to-noise tracers such as QSOs, however, the damping parameters are weakly constrained and can absorb noise fluctuations, leading to shifts in inferred parameters. Similar limitations may arise in models where cosmological information is encoded in power-spectrum shape features degenerate with the damping, such as scenarios with massive neutrinos. In contrast, for $w_0w_a$CDM we obtain $15\%$ and $21\%$ tighter constraints on $w_0$ and $w_a$, respectively, yielding a deviation from constant dark energy at slightly more than the $1σ$ level using full-shape information alone. The code is publicly available at https://github.com/cosmodesi/FolpsD

Bansal, P. [Michigan U., MCTP; Michigan U.] (ORCID↗
DOE OSTI2021

GPU-accelerated multitiered iterative phasing algorithm for fluctuation X-ray scattering

The multitiered iterative phasing (MTIP) algorithm is used to determine the biological structures of macromolecules from fluctuation scattering data. It is an iterative algorithm that reconstructs the electron density of the sample by matching the computed fluctuation X-ray scattering data to the external observations, and by simultaneously enforcing constraints in real and Fourier space. This paper presents the first ever MTIP algorithm acceleration efforts on contemporary graphics processing units (GPUs). The Compute Unified Device Architecture (CUDA) programming model is used to accelerate the MTIP algorithm on NVIDIA GPUs. The computational performance of the CUDA-based MTIP algorithm implementation outperforms the CPU-based version by an order of magnitude. Furthermore, the Heterogeneous-Compute Interface for Portability (HIP) runtime APIs are used to demonstrate portability by accelerating the MTIP algorithm across NVIDIA and AMD GPUs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Efficient computation of N -point correlation functions in D dimensions

We present efficient algorithms for computing the N-point correlation functions (NPCFs) of random fields in arbitrary D-dimensional homogeneous and isotropic spaces. Such statistics appear throughout the physical sciences and provide a natural tool to describe stochastic processes. Typically, algorithms for computing the NPCF components have $\mathscr O$(n N ) complexity (for a dataset containing n particles); their application is thus computationally infeasible unless N is small. By projecting the statistic onto a suitably defined angular basis, we show that the estimators can be written in a separable form, with complexity $\mathscr O$(n 2 ) or $\mathscr O$(n g log n g ) if evaluated using a Fast Fourier Transform on a grid of size n g . Our decomposition is built upon the D-dimensional hyperspherical harmonics; these form a complete basis on the (D – 1) sphere and are intrinsically related to angular momentum operators. Concatenation of (N – 1) such harmonics gives states of definite combined angular momentum, forming a natural separable basis for the NPCF. As N and D grow, the number of basis components quickly becomes large, providing a practical limitation to this (and all other) approaches: However, the dimensionality is greatly reduced in the presence of symmetries; for example, isotropic correlation functions require only states of zero combined angular momentum. We provide a Julia package implementing our estimators and show how they can be applied to a variety of scenarios within cosmology and fluid dynamics. The efficiency of such estimators will allow higher-order correlators to become a standard tool in the analysis of random fields.

97 MATHEMATICS AND COMPUTING↗