Search NASA⌕ Search

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.

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

122 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The r -modes of slowly rotating, stratified neutron stars

ABSTRACT The only r-modes that exist in a globally barotropic, rotating, Newtonian star are the fundamental l = |m| solutions, where l and m are the indices of the spherical harmonic $Y_l^m$ that describe the mode’s angular dependence. This is in stark contrast to a stellar model that is non-barotropic throughout its interior, which hosts all the l ≥ |m| perturbations including radial overtones. In reality, neutron stars are stratified with locally barotropic regions. Therefore, we explore how stratification alters a star’s ability to support r-modes. We consider the globally stratified case and examine the behaviour of the modes as the star gets close to barotropicity. In this limit, we find that all but the fundamental l = |m| perturbations change character and become generic inertial modes. Restricting the analysis to l = |m| perturbations, we develop the r-mode equations in order to consider stellar models that exhibit local barotropicity. Our results for such models show that the r-mode overtones diverge and join the inertial modes. In order to see which r-modes persist and retain their character in realistic neutron stars, these calculations will need to be brought into full general relativity.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Preliminary assessment of the accuracy and precision of TOPEX/POSEIDON altimeter data with respect to the large-scale ocean circulation

TOPEX/POSEIDON sea surface height measurements are examined for quantitative consistency with known elements of the oceanic general circulation and its variability. Project-provided corrections were accepted but are at tested as part of the overall results. The ocean was treated as static over each 10-day repeat cycle and maps constructed of the absolute sea surface topography from simple averages in 2 deg x 2 deg bins. A hybrid geoid model formed from a combination of the recent Joint Gravity Model-2 and the project-provided Ohio State University geoid was used to estimate the absolute topography in each 10-day period. Results are examined in terms of the annual average, seasonal average, seasonal variations, and variations near the repeat period. Conclusion are as follows: the orbit error is now difficult to observe, having been reduced to a level at or below the level of other error sources; the geoid dominates the error budget of the estimates of the absolute topography; the estimated seasonal cycle is consistent with prior estimates; shorter-period variability is dominated on the largest scales by an oscillation near 50 days in spherical harmonics Y(sup m)(sub 1)(theta, lambda) with an amplitude near 10 cm, close to the simplest alias of the M(sub 2) tide. This spectral peak and others visible in the periodograms support the hypothesis that the largest remaining time-dependent errors lie in the tidal models. Though discrepancies attribute to the geoid are within the formal uncertainties of the good estimates, removal of them is urgent for circulation studies. Current gross accuracy of the TOPEX/POSEIDON mission is in the range of 5-10 cm, distributed overbroad band of frequencies and wavenumbers. In finite bands, accuracies approach the 1-cm level, and expected improvements arising from extended mission duration should reduce these numbers by nearly an order of magnitude.

Wunsch, Carl↗
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 OSTI2022

Near-field radiative heat transfer between irregularly shaped dielectric particles modeled with the discrete system Green's function method

Near-field radiative heat transfer (NFRHT) between irregularly shaped dielectric particles made of SiO 2 and morphology characterized by Gaussian random spheres is studied. Particles are modeled using the discrete system Green's function (DSGF) approach, which is a volume integral numerical method based on fluctuational electrodynamics. This method is applicable to finite, three-dimensional objects, and all system interactions are defined independent of thermal excitation by a generalized system Green's function. The DSGF method is deemed suitable to model NFRHT between irregularly shaped particles after verification against the analytical solution for chains of two and three SiO 2 spheres. The NFRHT results reveal that geometric irregularity in particles leads to a reduction of the total conductance from that of comparable perfect spheres at vacuum separation distances smaller than the particle size, a regime in which NFRHT is a surface phenomenon. At vacuum separation distances larger than the particle size, NFRHT becomes a volumetric process, and the total conductance between irregularly shaped particles converges to that of comparable perfect spheres. Spectral analysis reveals, however, that particle irregularity leads to damping and broadening of resonances at all separation distances, thereby highlighting the importance of the DSGF method for spectral engineering in the near field. The reduced spectral coherence when particle size is larger than the vacuum separation distance is attributed to coupling of surface phonon-polaritons within the randomly generated, distorted particle features. For particle size smaller than the vacuum separation distance, resonance broadening and damping are linked with the multiple localized surface phonon modes supported by the composite spherical harmonic morphologies of the Gaussian random spheres. In conclusion, this paper has direct implications for thermal management of packed particle systems, with applications in radiative property control, electronics, energy conversion, and nanomanufacturing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2001

Satellite Remote Sensing of Tropical Precipitation and Ice Clouds for GCM Verification

This project, supported by the NASA New Investigator Program, has primarily been funding a graduate student, Darren McKague. Since August 1999 Darren has been working part time at Raytheon, while continuing his PhD research. Darren is planning to finish his thesis work in May 2001, thus some of the work described here is ongoing. The proposed research was to use GOES visible and infrared imager data and SSM/I microwave data to obtain joint distributions of cirrus cloud ice mass and precipitation for a study region in the Eastern Tropical Pacific. These joint distributions of cirrus cloud and rainfall were to be compared to those from the CSU general circulation model to evaluate the cloud microphysical amd cumulus parameterizations in the GCM. Existing algorithms were to be used for the retrieval of cloud ice water path from GOES (Minnis) and rainfall from SSM/I (Wilheit). A theoretical study using radiative transfer models and realistic variations in cloud and precipitation profiles was to be used to estimate the retrieval errors. Due to the unavailability of the GOES satellite cloud retrieval algorithm from Dr. Minnis (a co-PI), there was a change in the approach and emphasis of the project. The new approach was to develop a completely new type of remote sensing algorithm - one to directly retrieve joint probability density functions (pdf's) of cloud properties from multi-dimensional histograms of satellite radiances. The usual approach is to retrieve individual pixels of variables (i.e. cloud optical depth), and then aggregate the information. Only statistical information is actually needed, however, and so a more direct method is desirable. We developed forward radiative transfer models for the SSM/I and GOES channels, originally for testing the retrieval algorithms. The visible and near infrared ice scattering information is obtained from geometric ray tracing of fractal ice crystals (Andreas Macke), while the mid-infrared and microwave scattering is computed with Mie scattering. The radiative transfer is performed with the Spherical Harmonic Discrete Ordinate Method (developed by the PI), and infrared molecular absorption is included with the correlated k-distribution method. The SHDOM radiances have been validated by comparison to version 2 of DISORT (the community "standard" discrete-ordinates radiative transfer model), however we use SHDOM since it is computationally more efficient.

Evans, K. Franklin↗
DOE OSTI2023

Unsupervised Learning Based Interaction Force Model for Nonspherical Particles in Incompressible Flows

This project provides a neural network-based interaction force model for gas-solid flows from low to intermediate Reynolds numbers and concentration, which can be linked to MFiX-DEM. We have constructed a database of the interaction force between the irregular-shaped particles using a spherical harmonic method and the fluid phase based on the particle-resolved direct numerical simulation (PR-DNS) with immersed boundary-based gas kinetic scheme. Unsupervised learning method, i.e., variational auto-encoder (VAE) has been applied to extract the primitive shape factors determining the drag force, lifting forces, and torque. The interaction force model has been trained and validated with a simple but effective multi-layer feed-forward neural network: multi-layer perceptron (MLP), which will be concatenated after the encoder of the previously trained VAE for geometry feature extraction for single, irregular particles. We have trained transpose convolutional neural networks with the PR-DNS data to predict the velocity and pressure gradient of the single particle systems and utilized them to calculate drag force of multi-particle systems. This model can provide high computational efficiency because it does not require collecting multiparticle system data from PR-DNS.

99 GENERAL AND MISCELLANEOUS↗
NASA NTRS2010

High-Speed Solution of Spacecraft Trajectory Problems Using Taylor Series Integration

It has been known for some time that Taylor series (TS) integration is among the most efficient and accurate numerical methods in solving differential equations. However, the full benefit of the method has yet to be realized in calculating spacecraft trajectories, for two main reasons. First, most applications of Taylor series to trajectory propagation have focused on relatively simple problems of orbital motion or on specific problems and have not provided general applicability. Second, applications that have been more general have required use of a preprocessor, which inevitably imposes constraints on computational efficiency. The latter approach includes the work of Berryman et al., who solved the planetary n-body problem with relativistic effects. Their work specifically noted the computational inefficiencies arising from use of a preprocessor and pointed out the potential benefit of manually coding derivative routines. In this Engineering Note, we report on a systematic effort to directly implement Taylor series integration in an operational trajectory propagation code: the Spacecraft N-Body Analysis Program (SNAP). The present Taylor series implementation is unique in that it applies to spacecraft virtually anywhere in the solar system and can be used interchangeably with another integration method. SNAP is a high-fidelity trajectory propagator that includes force models for central body gravitation with N X N harmonics, other body gravitation with N X N harmonics, solar radiation pressure, atmospheric drag (for Earth orbits), and spacecraft thrusting (including shadowing). The governing equations are solved using an eighth-order Runge-Kutta Fehlberg (RKF) single-step method with variable step size control. In the present effort, TS is implemented by way of highly integrated subroutines that can be used interchangeably with RKF. This makes it possible to turn TS on or off during various phases of a mission. Current TS force models include central body gravitation with the J2 spherical harmonic, other body gravitation, thrust, constant atmospheric drag from Earth's atmosphere, and solar radiation pressure for a sphere under constant illumination. The purpose of this Engineering Note is to demonstrate the performance of TS integration in an operational trajectory analysis code and to compare it with a standard method, eighth-order RKF. Results show that TS is 16.6 times faster on average and is more accurate in 87.5% of the cases presented.

Scott, James R.↗
DOE OSTI2023

Cosmic topology. Part I. Limits on orientable Euclidean manifolds from circle searches

The Einstein field equations of general relativity constrain the local curvature at every point in spacetime, but say nothing about the global topology of the Universe. Cosmic microwave background anisotropies have proven to be the most powerful probe of non-trivial topology since, within ΛCDM, these anisotropies have well-characterized statistical properties, the signal is principally from a thin spherical shell centered on the observer (the last scattering surface), and space-based observations nearly cover the full sky. The most generic signature of cosmic topology in the microwave background is pairs of circles with matching temperature and polarization patterns. No such circle pairs have been seen above noise in the WMAP or Planck temperature data, implying that the shortest non-contractible loop around the Universe through our location is longer than 98.5% of the comoving diameter of the last scattering surface. We translate this generic constraint into limits on the parameters that characterize manifolds with each of the nine possible non-trivial orientable Euclidean topologies, and provide a code which computes these constraints. In all but the simplest cases, the shortest non-contractible loop in the space can avoid us, and be shorter than the diameter of the last scattering surface by a factor ranging from 2 to at least 6. Further, this result implies that a broader range of manifolds is observationally allowed than widely appreciated. Probing these manifolds will require more subtle statistical signatures than matched circles, such as off-diagonal correlations of harmonic coefficients.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

A spectral formalism for computing three-dimensional deformations due to surface loads. 1: Theory

We outline a complete spectral formalism for computing high spatial resolution three-dimensional deformations arising from the surface mass loading of a spherically symmetric planet. The main advantages of the formalism are that all surface mass loads are always described using a consistent mathematical representation and that calculations of deformation fields for various spatial resolutions can be performed by simpley altering the spherical harmonic degree truncation level of the procedure. The latter may be important when incorporating improved observational constraints on a particular surface mass load, when considering potential errors in the computed field associated with mass loading having a spatial scale unresolved by the observational constraints, or when treating a number of global surface mass loads constrained with different spatial resolutions. The advantages do not extend to traditional 'Green's function' approaches which involve surface element discretizations of the global mass loads. Another advantage of the spectral formalism, over the Green's function approach, is that a posteriori analyses of the computed deformation fields are easily performed. In developing the spectral formalism, we consider specific cases where the Earth's mantle is assumed to respond as an elastic, slightly anelastic, or linear viscoelastic medium. In the case of an elastic or slightly anelastic mantle rheology the spectral response equations incorporate frequency dependent Love numbers. The formalism can therefore be used, for example, to compute the potentially resonant deformational response associated with the free core nutation and Chandler wobble eigenfunctions. For completeness, the spectral response equations include both body forces, as arise from the gravitational attraction of the Sun and the Moon, and surface mass loads. In either case, and for both elastic and anelastic mantle rheologies, we outline a pseudo-spectral technique for computing the ocean adjustment associated with the total gravitational perturbation induced by the external forcing. Three-dimensional deformations computed using the usual Love number approach are generally referenced to an origin at the center of mass of the undeformed planet. We derive a spectral technique for transforming the results to an origin located at the center of mass of the deformed planet.

Mitrovica, J. X.↗
NASA NTRS2019

Quantifying Water Storage Change and Land Subsidence Induced by Reservoir Impoundment Using GRACE, Landsat, and GPS Data

The construction of hydropower dams is a common strategy to support a country's increasing need for electricityand river water management for industry and agriculture. Although the hydrological and geophysical impacts ofwater relocation are usually assessed prior to impoundment, their accuracy is generally limited due to the lack ofin situ observations, especially in a remote area. This study presents a workflow to quantify the terrestrial waterstorage change (TWS) and land subsidence induced by a reservoir's water impoundment using multiple satelliteobservations (GRACE, Landsat), land surface models (CABLE, GLDAS, NCEP, ECMWF), and GPS data. The studysite is the Bakun Dam, located in Sarawak, Malaysia, which is the largest hydropower dam in Southeast Asia.Commencing operation in late 2010, the dam induced a change of water mass and lake surface area that wasclearly observed by GRACE and Landsat observations, respectively. During the 17-month impounding period(from August 2010 to December 2011), GRACE observed a dramatic increase of approximately 200mmequivalent water height, while Landsat detected an increased lake extent of around 600 km2. In this paper, aforward model is developed to determine the increased water surface level corresponding to GRACE observations,estimated to be about 120 m. In contrast to GRACE, the TWS derived from land surface models cannotcapture the increased TWS, due to the lack of reservoir routing algorithms in the models. In addition, the landsubsidence was calculated using the disk load model constructed based on the GRACE-derived lake level andLandsat-derived lake extent; the result is validated with the GPS data from BIN1 station, located at the westerncoast of Borneo. The commencement stage of the Bakun Dam induces the large-scale land subsidence, whichcauses the GPS-BIN1 station to subside by ~9 mm, and move toward the Bakun Lake by ~4 mm. Computation ofthe surface displacements directly from GRACE spherical harmonic coefficient data fails to capture the subsidencefeature, mainly due to the truncation error. Overall, this study demonstrates that evaluating GRACE inconjunction with Landsat, LSMs, and GPS data allows the exploitation of the gravity signal at a much smallerspatial scale than its intrinsic resolution. Benefiting from global coverage, the newly developed satellite-basedalgorithm is a valuable tool for assessing the impacts of reservoir operation on hydrological and geophysicalchanges from local to regional scales.

Tangdamrongsub, Natthachet↗
NASA NTRS2019

Quantifying Water Storage Change and Land Subsidence Induced by Reservoir Impoundment Using GRACE, Landsat, and GPS Data

The construction of hydropower dams is a common strategy to support a country's increasing need for electricity and river water management for industry and agriculture. Although the hydrological and geophysical impacts of water relocation are usually assessed prior to impoundment, their accuracy is generally limited due to the lack of in situ observations, especially in a remote area. This study presents a workflow to quantify the terrestrial water storage change (TWS) and land subsidence induced by a reservoir's water impoundment using multiple satellite observations (GRACE, Landsat), land surface models (CABLE, GLDAS, NCEP, ECMWF), and GPS data. The study site is the Bakun Dam, located in Sarawak, Malaysia, which is the largest hydropower dam in Southeast Asia. Commencing operation in late 2010, the dam induced a change of water mass and lake surface area that was clearly observed by GRACE and Landsat observations, respectively. During the 17-month impounding period (from August 2010 to December 2011), GRACE observed a dramatic increase of approximately 200 mm equivalent water height, while Landsat detected an increased lake extent of around 600 km2. In this paper, a forward model is developed to determine the increased water surface level corresponding to GRACE observations, estimated to be about 120 m. In contrast to GRACE, the TWS derived from land surface models cannot capture the increased TWS, due to the lack of reservoir routing algorithms in the models. In addition, the land subsidence was calculated using the disk load model constructed based on the GRACE-derived lake level and Landsat-derived lake extent; the result is validated with the GPS data from BIN1 station, located at the western coast of Borneo. The commencement stage of the Bakun Dam induces the large-scale land subsidence, which causes the GPS-BIN1 station to subside by ~9 mm, and move toward the Bakun Lake by ~4 mm. Computation of the surface displacements directly from GRACE spherical harmonic coefficient data fails to capture the subsidence feature, mainly due to the truncation error. Overall, this study demonstrates that evaluating GRACE in conjunction with Landsat, LSMs, and GPS data allows the exploitation of the gravity signal at a much smaller spatial scale than its intrinsic resolution. Benefiting from global coverage, the newly developed satellite-based algorithm is a valuable tool for assessing the impacts of reservoir operation on hydrological and geophysical changes from local to regional scales.

Natthachet Tangdamrongsub↗
DOE OSTI2020

Synthesis of TiO{sub 2} nanoparticles decorated with gold nanoclusters using pulsed laser ablation in liquid

Amorphous spherical titanium dioxide (TiO{sub 2}) nanoparticles, decorated gold nanoparticles, are shown to be formed under irradiation of commercial TiO{sub 2} nanopowders dispersed in aqueous solution of hydrogen tetrachloroaurate (HAuCl{sub 4}) by second-harmonic nanosecond pulses of an Nd : YAG laser. It is found that these hybrid nanostructures are formed during laser-induced remelting of the initial nanoparticles, stimulated by gold nanoclusters (providing more efficient absorption of visible laser radiation by titanium dioxide) recovered on their surface. The morphology and chemical composition of the newly formed hybrid nanomaterials are investigated in detail by electron microscopy, Raman spectroscopy, and energy-dispersive analysis. It is shown that the average size and number of gold nanoparticles recovered on the TiO{sub 2} surface can be controlled by varying the concentration ratio of the initial nanomaterial and HAuCl{sub 4}. A spectroscopic analysis of light scattering by single hybrid nanoparticles and the results of numerical calculation of the structure of electromagnetic fields near their surface indicate good prospects of these hybrid nanomaterials for various applications in modern optics, optoelectronics, and nanophotonics; for example, for designing chemo- and biosensor platforms and new-generation solar cells. (nanoparticles)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Gravity of two photon decay and its quantum coherence

A linear analytical solution is derived for the gravitational shock wave produced by a particle of mass M that decays into a pair of null particles. The resulting space-time is shown to be unperturbed and isotropic, except for a discontinuous perturbation on a spherical null shell. Formulae are derived for the perturbation as a function of polar angle, as measured by an observer at the origin observing clocks on a sphere at distance R. The effect of the shock is interpreted physically as an instantaneous displacement in time and velocity when the shock passes the clocks. The time displacement is shown to be anisotropic, dominated by a quadrupole harmonic aligned with the particle-decay axis, with a magnitude δτ ~ GM/c 3 , independent of R. The velocity displacement is isotropic. The solution is used to derive the gravitational effect of a quantum state with a superposition of a large number of randomly oriented, statistically isotropic particle decays. This approach is shown to provide a well-controlled approximation to estimate the magnitude of gravitational fluctuations in systems composed of null point particles up to the Planck energy in a causal diamond of duration τ = 2R/c, as well as quantum-gravitational fluctuations of black holes and cosmological horizons. Coherent large-angle quantum distortions of macroscopic geometry from fluctuations up to the Planck scale are shown to grow linearly with the duration, with a variance $\langle$δτ 2 $\rangle$ ~ τtP much larger than that produced in models without causal quantum coherence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Performance Improvements for the Griffin Transport Solvers

Griffin is a Multiphysics Object-Oriented Simulation Environment based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. Griffin includes a variety of deterministic radiation transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. A code assessment performed in FY-20 identified two significant issues with the transport solvers in Griffin: first, the primary heterogeneous SN (discrete ordinates) transport solver based on continuous finite element methods required significant mesh refinement and higher memory usage compared to solvers based on the method of characteristic for equivalent accuracy. Second, the homogeneous PN (spherical harmonics expansion) transport solver did not adequately support polynomial refinement, which is a feature usually required for problems with spatial homogenization and pronounced streaming, typical in fast or gas-cooled reactor systems. To address the first issue, the development effort focused on the more promising discontinuous finite element method (DFEM)-based SN transport solver in Griffin. The addition of an asynchronous parallel transport sweeper and coarse mesh finite difference (CMFD) acceleration have rendered a superior heterogeneous SN transport capability for multiphysics problems that requires far less computing resources in terms of both CPU time and memory usage. This is demonstrated with typical thermal- and fast-spectrum reactor benchmark problems, including 2D Transient Reactor Test, 3D Advanced Burner Test Reactor (ABTR), and 2D and 3D Empire microreactor. For the second issue, the development effort focused on a new transport solver based on the hybrid finite element PN method (HFEM-PN), equivalent to the variational nodal method, as well as a new diffusion solver based on HFEM-Diffusion. This solver is intended for homogenized domains with multiphysics coupling (i.e., supports mesh displacement, seamless temperature feedback, etc.). Initial calculations with the HFEM-Diffusion implementation show very good parallel efficiency for the residual evaluations with the 2D ABTR benchmark. A future development effort will be centered on further improvements to the CMFD, HFEM-PN, and DFEM diffusion solvers to ensure Griffin meets performance and software quality assurance requirements for advanced reactor design and analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗