Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spherical harmonic analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Effects of initial conditions and transport on ram pressure, Mach number, and uniformity for plasma liner formation and implosion

Plasma jet driven magneto-inertial fusion involves the merging of high Mach number plasma jets in a spherically converging geometry to form an imploding plasma liner. A three-dimensional smooth particle hydrodynamic simulation tool was used to study the number of plasma jets, jet density, and implosion velocity on the effects of formation and subsequent implosion of a plasma liner. The effects of different transport physics are presented. Thermal radiation and two-temperature modeling were found to have the most significant impact on the ram pressure. Solid angle averaged ram pressure and Mach number degradation were of particular interest in this study. Higher jet numbers were found to result in higher peak ram pressure. A spherical harmonic analysis was conducted for each case, providing a quantitative measure of the non-uniformity of the liner during the implosion process and how non-uniformities are affected by the different parameters explored in this study.

Schillo, Kevin↗
DOE OSTI2022

Constraining low redshift [C II ] emission by cross-correlating FIRAS and BOSS data

ABSTRACT We perform a tomographic cross-correlation analysis of archival FIRAS data and the BOSS galaxy redshift survey to constrain the amplitude of [C II] 2P3/2 → 2P1/2 fine structure emission. Our analysis employs spherical harmonic tomography (SHT), which is based on the angular cross-power spectrum between FIRAS maps and BOSS galaxy over-densities at each pair of redshift bins, over a redshift range of 0.24 < z < 0.69. We develop the SHT approach for intensity mapping, where it has several advantages over existing power spectral estimators. Our analysis constrains the product of the [C II] bias and [C II] specific intensity, $b_{\rm [C \small{\rm II}]}I_{\rm [C \small{\rm II}]}$, to be <0.31 MJy/sr at z ≈ 0.35 and <0.28 MJy/sr at z ≈ 0.57 at $95{{\ \rm per\ cent}}$ confidence. These limits are consistent with most current models of the [C II] signal, as well as with higher-redshift [C II] cross-power spectrum measurements from the Planck satellite and BOSS quasars. We also show that our analysis, if applied to data from a more sensitive instrument such as the proposed PIXIE satellite, can detect pessimistic [C II] models at high significance.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Discrete spherical harmonic functions for texture representation and analysis

A basis of discrete harmonic functions for efficient representation and analysis of crystallographic texture is presented. Discrete harmonics are a numerical representation of the harmonics on the sphere. A finite element formulation is utilized to calculate these orthonormal basis functions, which provides several advantageous features for quantitative texture analysis. These include high-precision numerical integration, a simple implementation of the non-negativity constraint and computational efficiency. Simple examples of pole figure and texture interpolation and of Fourier filtering using these basis sets are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Towards a self-consistent analysis of the anisotropic galaxy two- and three-point correlation functions on large scales: application to mock galaxy catalogues

We establish a practical method for the joint analysis of anisotropic galaxy two- and three-point correlation functions (2PCF and 3PCF, respectively) on the basis of the decomposition formalism of the 3PCF using tripolar spherical harmonics. We perform such an analysis with MultiDark-Patchy mock catalogues to demonstrate and understand the benefit of the anisotropic 3PCF. We focus on scales above $80\, h^{-1}\, {\rm Mpc}$, and use information from the shape and the baryon acoustic oscillation (BAO) signals of the 2PCF and 3PCF. We also apply density field reconstruction to increase the signal-to-noise ratio of BAO in the 2PCF measurement, but not in the 3PCF measurement. In particular, we study in detail the constraints on the angular diameter distance and the Hubble parameter. We build a model of the bispectrum or 3PCF that includes the non-linear damping of the BAO signal in redshift space. We carefully account for various uncertainties in our analysis including theoretical models of the 3PCF, window function corrections, biases in estimated parameters from the fiducial values, the number of mock realizations to estimate the covariance matrix, and bin size. The joint analysis of the 2PCF and 3PCF monopole and quadrupole components shows a $30{{\ \rm per\ cent}}$ and $20{{\ \rm per\ cent}}$ improvement in Hubble parameter constraints before and after reconstruction of the 2PCF measurements, respectively, compared to the 2PCF analysis alone. Furthermore, this study clearly shows that the anisotropic 3PCF increases cosmological information from galaxy surveys and encourages further development of the modelling of the 3PCF on smaller scales than we consider.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Functional Principal Component Analysis of Radio–Optical Reference Frame Tie

The Gaia optical reference frame is intrinsically undefined with respect to global orientation and spin, so it needs to be anchored in the radio-based International Celestial Reference Frame (ICRF) to provide a referenced and quasi-inertial celestial coordinate system. The link between the two fundamental frames is realized through two samples of distant extragalactic sources, mostly active galactic nuclei and quasars, but only the smaller sample of radio-loud ICRF sources with optical counterparts is available to determine the mutual orientation. The robustness of this link can be mathematically formulated in the framework of functional principal component analysis using a set of vector spherical harmonics to represent the differences in celestial positions of the common objects. The weakest eigenvectors are computed, which describe the greatest deficiency of the link. The deficient or poorly determined terms are specific vector fields on the sphere that carry the largest errors of absolute astrometry using Gaia in reference to the ICRF. This analysis provides guidelines for the future development of the ICRF maximizing the accuracy of the link over the entire celestial sphere. A measure of the robustness of a least-squares solution, which can be applied to any linear model fitting problem, is introduced to help discriminate between reference frame tie models of different degrees.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Quantifying molecular deformation in polymer melts by a generalized Zimm plot approach

The Zimm plot has been widely used to characterize the molecular dimensions of polymers from small-angle scattering experiments, where the reciprocal intensity is analyzed as a function of the square of the magnitude of the scattering wavevector Q. Here, this work explores the benefits of analyzing the reciprocal scattering intensity from deformed polymers, extending the original Zimm plot to anisotropic materials. In the small-angle limit, a tensorial extension of the Guinier law is found for the gyration tensor and the reciprocal single-chain structure factor. In the high-Q limit, application of the spherical harmonic expansion technique to the reciprocal structure factor permits direct model-independent analysis of spatially dependent molecular deformation of polymers. Additionally, the contributions from high-order spherical harmonics become insignificant in the reciprocal-intensity representation. The proposed generalized Zimm plot approach is demonstrated computationally with the affine deformation model and the Rouse model, and experimentally with small-angle neutron scattering measurements of deformed polystyrene melts.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Quantification of Deformation-Induced Concentration Fluctuations in Polymeric Liquids by Small-Angle Neutron Scattering

This study demonstrates how the large concentration fluctuations of polymeric liquids under deformation and flow can be quantitatively studied by applying the spherical harmonic expansion technique to small-angle neutron scattering. Using a series of binary polymer blends, it is shown that the emergence of the so-called butterfly patterns is caused by the change of sign in the leading anisotropic component of the small-angle spectrum, when the scattering is dominated by intermolecular correlation associated with flow-induced demixing. The increasing spatial fluctuations of concentration are evidenced by the enhancement of the isotropic component of the scattering spectrum in the zero-angle limit and peak shift of the leading anisotropic coefficients toward low Q. Additionally, the spherical harmonic expansion framework permits real-space analysis of anisotropic scattering length density correlation in a convenient form. The methodology described in this work provides a concrete venue for quantitative studies of phase behavior of fluids under deformation and flow via small-angle scattering techniques, where mainly qualitative approaches were previously employed.

Chemical structure↗
DOE OSTI2022

Energy Exascale Earth System Model v2.0.1

First patch release of v2.0.0 Changes since v2.0.0 [Important change] Fix ocean threading bug seen in debug cases on Chrysalis with Intel 20.0.4. Was introduced around time of v2.0.0 tag. Does not change v2.0.0 answers on Chrysalis because those didn't use threading or debugging. [EAM] Add semi-lagrangian tracer transport for theta-l (F90 and C++), add new algorithm for finding tropopause, add DSCREAM to allow v2 and SCREAM settings in same code such as adjust_ps [EAM-MMF] 60L default, allow C++ back end of RRTMGP (EAM too). [EAMxx] add nu-top functionality, fix forcing functor, add ttype9 and dcmip2012 tests 2.1, 2.2, and 3 HOMME: remove obsolete remap algs, option to specify dynamics alg indep of tracer, new sponge layer, add imex tests [ELM] Add topography-based subgrid (topounits), add FATES-ELM Nitro., Phos. and CH4 coupling, add land-use ts for NARRM, add lulc for SSP3 RCP7, Fix nutrient fertilization exp test and carbon isotope flux, Fix xactive lnd dry deposition, add lake water storage option, fix plant hydraulics 2d params, fix carbon budget calc, fix soil nutrient conc. bug, fix mosart dam bug, add test for new ELM, MOSART features, fix bug in O3 dry dep stomatal resistances, fix plant hydraulics restart BFB error, update mkmapdata. [MOSART] fix bug for reading the latitude from an unstructured input file, fix oversat in bubble test. [MPAS-ocean] Add CFC11, CFC12 tracers, add 2D spherical transport tests, fix del4 tracer mixing, add MARBL ocean tracer mixing, modify harmonic analysis options, add GPU port of vmix routines, fix calc of ML-averaged BV freq. [MPAS-seaice] Change extents of initial polar disks for oRRS18to6v3 grid, fix ice BGC with MARBL, update spherical test cases, fix DON coupling, Remove Cf from sea ice constants. [MPAS-landice] add CRYO1850-4xCO2 compset [CIME] add GCP, ANL GCE, Spock, Perlmutter, deprecate config_compilers.xml, fix and clean-up cmake macros, fix slurm bindings, refactor CIME internal testing, cleanup SCORPIO perf data, allow position independent compset naming, [also] update v2 benchmarking suite, extend e3sm_prod with throughput and memory checks

E3SM Project, DOE↗
DOE OSTI2023

Coarse-grained methods for heterogeneous vesicles with phase-separated domains: Elastic mechanics of shape fluctuations, plate compression, and channel insertion

We develop coarse-grained particle approaches for studying the elastic mechanics of vesicles with heterogeneous membranes having phase-separated domains. We perform simulations both of passive shape fluctuations and of active systems where vesicles are subjected to compression between two plates or subjected to insertion into narrow channels. Analysis methods are developed for mapping particle configurations to continuum fields with spherical harmonics representations. Heterogeneous vesicles are found to exhibit rich behaviors where the heterogeneity can amplify surface two-point correlations, reduce resistance during compression, and augment vesicle transport times in channels. The developed methods provide general approaches for characterizing the mechanics of coarse-grained heterogeneous systems taking into account the roles of thermal fluctuations, geometry, and phase separation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Scattering insights into shear-induced scission of rod-like micelles

Understanding the scission of rod-like micelles under mechanical forces is crucial for optimizing their stability and behavior in industrial applications. This study investigates how micelle length, flexibility, and external forces interact, offering insights into the design of micellar systems in processes influenced by mechanical stress. Although significant, direct experimental observations of flow-induced micellar scission using scattering techniques remain scarce. Small angle neutron scattering (SANS) is used to explore the shear response of aqueous cetyltrimethylammonium bromide (CTAB) solutions with sodium nitrate. Rheological tests show shear thinning with no shear banding, ensuring a uniform flow field for reliable interpretation of scattering data. As shear rate increases, the scattering spectra show angular distortion, which is analyzed using spherical harmonic decomposition to characterize flow-induced scission and micelle orientation under shear. Two analysis steps are used: a model-independent spectral eigendecomposition reveals a decrease in micellar length, while regression analysis quantifies the evolution of the length distribution and mean length with shear rate. Additionally, micelle alignment increases with shear, quantified by the orientational distribution function. In conclusion, these findings provide experimental evidence for flow-induced alignment and scission, offering a new framework for understanding shear-induced phenomena in micellar systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Griffin: A MOOSE-based reactor physics application for multiphysics simulation of advanced nuclear reactors

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor physics application for multiphysics simulations of advanced reactor designs jointly developed by Idaho National Laboratory and Argonne National Laboratory. This paper summarizes the motivation, significance, architecture, design, and features of Griffin. Griffin offers flexible and extensible features to address the challenges associated with advanced reactor designs. These features range from fundamental particle transport to specific reactor physics tasks. The features cover a wide range including on-the-fly and traditional two-step cross-section generation methods, steady-state and transient transport solvers suitable for both heterogeneous and homogeneous models, high-fidelity depletion where thousands of isotopes can be tracked and low-fidelity depletion characterized by burnup, etc. The most fundamental aspect that sets Griffin apart from other reactor analysis codes is that it is developed based on the MOOSE framework. A modular development approach is strongly enforced, with multiphysics being an essential element considered since the beginning of Griffin’s development. Griffin links various MOOSE physics modules and couples to other MOOSE-based applications and non-MOOSE-based applications for multiphyiscs simulations. Griffin includes three modules: ISOXML for preparing and managing multigroup cross sections, radiation transport for solving the neutron transport equation, and reactor analysis for user-oriented reactor physics analysis functionalities. Griffin uses various finite element methods for spatial discretization, multigroup approximation for energy discretization and discrete ordinates method, spherical harmonics expansion method, and diffusion approximation for streaming direction discretization to solve the neutron transport equation. Griffin’s flexibility is evidenced through Griffin’s various applications to fast reactor, high-temperature reactor, pebble bed reactor, molten salt reactor, and microreactor designs. Griffin development follows the software quality assurance procedure for MOOSE-based applications and with software requirements consistent with the ASME NQA-1 standard. Griffin has been adopted into the reactor analysis system for the U.S. NRC and is in use at U.S. companies, universities and national laboratories.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Short and medium range structure in elastic deformation of metallic and covalent glasses

Here, we present a concise methodology to analyze structural response to the applied stress in amorphous solids, including metallic glasses (MG), glassy selenium, silica and polycarbonate, using high energy x-ray diffraction and atomic pair distribution function (PDF) analysis. To assess the structural anisotropy induced by applied axial stress, diffraction data were expanded into spherical harmonics. Using Bessel transformation, components of the structure function were converted into isotropic and anisotropic PDFs. The PDFs were compared to the expected model behavior for ideal elastic deformation to separate homogeneous affine strain from local non-affine strains. In metallic glass the range of non-affine deformation is limited to the nearest neighbor shell, suggesting local strain relaxation under stress that occurs even in the elastic regime. Beyond the second atomic shell strain is uniform. However, in glassy silica, polycarbonate and selenium strong local bonding inhibits local displacements and strain in short range order is accommodated by rotation of local units. Interestingly, beyond a molecular unit, deformation in covalent systems is similar to MG, and response of the medium range order scales with the macroscopic stress.

glassy structure↗
DOE OSTI2021

Real-time nuclear activation detectors for measuring neutron angular distributions at the National Ignition Facility (invited)

The Real Time Nuclear Activation Detector (RTNAD) array at NIF measures the distribution of 14 MeV neutrons emitted by deuterium-tritium (DT) fueled inertial confinement fusion implosions. The uniformity of the neutron distribution is an important indication of implosion symmetry and DT shell integrity. The array consists of 48 LaBr 3 (Ce) crystal gamma-ray spectrometers mounted outside the NIF target chamber, which continuously monitor the slow decay of the 909 keV gamma-ray line from activated 89 Zr located in Zr cups surrounding each crystal. The measured decay rate dramatically increases during a DT implosion in proportion to the number of 14 MeV neutrons striking each Zr cup. The neutrons produce activated 89 Zr through an (n, 2n) reaction on 90 Zr, which is insensitive to low energy neutrons. The neutron flux along the detector line-of-sight at shot time is determined by extrapolating the fitted 909 keV decay curve back to shot time. Automatic analysis algorithms were developed to handle the non-stop data stream. The large number of detectors and the high statistical accuracy of the array enable the spherical harmonic modes of the neutron angular distribution to be measured up to L ≤ 4 to provide a better understanding of implosion dynamics. In addition, these data combined with measurements of the down-scattered neutrons can be used to derive fuel areal density distributions. This paper will describe the RTNAD hardware and analysis procedures.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Implementation of hybrid finite element method based transport solver in GRIFFIN

A new transport solver option based on the hybrid FEM (HFEM) was implemented in GRIFFIN, the MOOSE-based reactor analysis code, as an effort to support routine core design calculations for advanced reactor applications. The HFEM formulation with P{sub N} (spherical harmonics expansion), akin to the variational nodal method, is effective for solving a spatially homogenized problem with strong transport effect. The residual and Jacobian evaluations of the HFEM weak form were derived and successfully implemented in GRIFFIN, having the diffusion and the PN options available in the new HFEM based transport solver. The performance was tested with the simplified ABTR benchmark problems. The results indicate that the HFEM-based transport solver is a feasible option for solving problems with spatially homogenized and strong streaming by providing superior accuracy with a proper p-refinement. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Precise and accurate cosmology with CMB×LSS power spectra and bispectra

With the advent of a new generation of cosmological experiments that will provide high-precision measurements of the cosmic microwave background (CMB) and galaxies in the large-scale structure, it is pertinent to examine the potential of performing a joint analysis of multiple cosmological probes. In this paper, we study the cosmological information content contained in the one-loop power spectra and tree bispectra of galaxies cross-correlated with CMB lensing. We use the FFTLog method to compute angular correlations in spherical harmonic space, applicable for wide angles that can be accessed by forthcoming galaxy surveys. Here, we find that adding the bispectra and cross-correlations with CMB lensing offers a significant improvement in parameter constraints, including those on the total neutrino mass, M ν , and local non-Gaussianity amplitude, ƒ NL . In particular, our results suggest that the combination of the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) and CMB-S4 will be able to achieve σ(M ν )=42 meV from galaxy and CMB lensing correlations, and σ(M ν )=12 meV when further combined with the CMB temperature and polarization data, without any prior on the optical depth.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗
DOE OSTI2022

Dark Energy Survey Year 3 results: A 2.7% measurement of baryon acoustic oscillation distance scale at redshift 0.835

We present angular diameter measurements obtained by measuring the position of baryon acoustic oscillations (BAO) in an optimized sample of galaxies from the first three years of Dark Energy Survey data (DES Y3). The sample consists of 7 million galaxies distributed over a footprint of 4100 deg 2 with 0.6 < z photo < 1.1 and a typical redshift uncertainty of 0.03 ( 1 + z ) . The sample selection is the same as in the BAO measurement with the first year of DES data, but the analysis presented here uses three times the area, extends to higher redshift, and makes a number of improvements, including a fully analytical BAO template, the use of covariances from both theory and simulations, and an extensive preunblinding protocol. We used two different statistics; angular correlation function and power spectrum, and validate our pipeline with an ensemble of over 1500 realistic simulations. Both statistics yield compatible results. We combine the likelihoods derived from angular correlations and spherical harmonics to constrain the ratio of comoving angular diameter distance D M at the effective redshift of our sample to the sound horizon scale at the drag epoch. We obtain D M ( z eff = 0.835 ) / r d = 18.92 ± 0.51 , which is consistent with, but smaller than, the Planck prediction assuming flat Λ CDM , at the level of 2.3 σ . The analysis was performed blind and is robust to changes in a number of analysis choices. It represents the most precise BAO distance measurement from imaging data to date, and is competitive with the latest transverse ones from spectroscopic samples at z > 0.75 . When combined with DES 3 x 2 pt + SNIa , they lead to improvements in H 0 and Ω m constraints by ∼ 20 % .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

BEYONDPLANCK VIII. Efficient sidelobe convolution and corrections through spin harmonics

We introduce a new formulation of the Conviqt convolution algorithm in terms of spin harmonics, and apply this to the problem of sidelobe correction for BEYONDPLANCK, the first end-to-end Bayesian Gibbs sampling framework for CMB analysis. We compare our implementation to the previous Planck LevelS implementation, and find good agreement between the two codes in terms of accuracy, but with a speed-up reaching a factor of 3–10, depending on the frequency bandlimits, l max and m max . The new algorithm is significantly simpler to implement and maintain, since all low-level calculations are handled through an external spherical harmonic transform library. We find that our mean sidelobe estimates for Planck LFI are in good agreement with previous efforts. Additionally, we present novel sidelobe rms maps that quantify the uncertainty in the sidelobe corrections due to variations in the sky model.

79 ASTRONOMY AND ASTROPHYSICS↗