Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Tensor factorizations”

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 307 records · Page 17

Post-Newtonian cosmological dynamics of plane-parallel perturbations and back-reaction

We study the general relativistic non-linear dynamics of self-gravitating irrotational dust in a cosmological setting, adopting the comoving and synchronous gauge, where all the equations can be written in terms of the metric tensor of spatial hyper-surfaces orthogonal to the fluid flow. Performing an expansion in inverse powers of the speed of light, we obtain the post-Newtonian equations, which yield the lowest-order relativistic effects arising during the non-linear evolution. We then specialize our analysis to globally plane-parallel configurations, i.e. to the case where the initial perturbation field depends on a single coordinate. The leading order of our expansion, corresponding to the ''Newtonian background'', is the Zel'dovich approximation, which, for plane-parallel perturbations in the Newtonian limit, represents an exact solution. This allows us to find the exact analytical form for the post-Newtonian metric, thereby providing the post-Newtonian extension of the Zel'dovich solution: this accounts for some relativistic effects, such as the non-Gaussianity of primordial perturbations. An application of our solution in the context of the back-reaction proposal is eventually given, providing a post-Newtonian estimation of kinematical back-reaction, mean spatial curvature, average scale-factor and expansion rate.

79 ASTRONOMY AND ASTROPHYSICS↗

Quark and gluon contributions to the QCD trace anomaly

We show that, in dimensional reguralization in the minimal subtraction scheme, the QCD trace anomaly can be unambiguously decomposed into two parts coming from the renormalized quark and gluon energy momentum tensors. Here, we carry out this decomposition at the two-loop level. The result can be used to constrain the renormalization group properties of the nucleon’s twist-four gravitational form factor $\bar{C}$ q,g .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Constraints on Inflationary Gravitational Waves with Two Years of SPT-3G Data

We present a measurement of the $B$-mode polarization power spectrum of the cosmic microwave background anisotropies at 32 $\le$$\ell$$<$ 502 for three bands centered at 95, 150, and 220 GHz using data from the SPT-3G receiver on the South Pole Telescope. This work uses SPT-3G observations from the 2019 and 2020 winter observing seasons of a $\sim$1500 deg$^2$ patch of sky that directly overlaps with fields observed with the BICEP/Keck family of telescopes, and covers part of the proposed Simons Observatory and CMB-S4 deep fields. Employing new techniques for mitigating polarized atmospheric noise, the SPT-3G data demonstrates a white noise level of 9.3 (6.7) $\mu$K-arcmin at $\ell \sim 500$ for the 95 GHz (150 GHz) data, with a $1/\ell$ noise knee at $\ell$=128 (182). We fit the observed six auto- and cross-frequency $B$-mode power spectra to a model including lensed $\Lambda$CDM $B$-modes and a combination of Galactic and extragalactic foregrounds. This work characterizes foregrounds in the vicinity of the BICEP/Keck survey area, finding foreground power consistent with that reported by the BICEP/Keck collaboration within the same region, and a factor of $\sim$ 3 higher power over the full SPT-3G survey area. Using SPT-3G data over the BICEP/Keck survey area, we place a 95% upper limit on the tensor-to-scalar ratio of $r < 0.25$ and find the statistical uncertainty on $r$ to be $\sigma(r) = 0.067$.

Zebrowski, J. A. [Chicago U., KICP; Chicago U., As↗

Prediction of absolute infrared intensities for the fundamental vibrations of H2O2

Absolute infrared intensities are predicted for the vibrational bands of gas-phase H2O2 by the use of a hydrogen atomic polar tensor transferred from the hydroxyl hydrogen atom of CH3OH. These predicted intensities are compared with intensities predicted by the use of a hydrogen atomic polar tensor transferred from H2O. The predicted relative intensities agree well with published spectra of gas-phase H2O2, and the predicted absolute intensities are expected to be accurate to within at least a factor of two. Among the vibrational degrees of freedom, the antisymmetric O-H bending mode nu(6) is found to be the strongest with a calculated intensity of 60.5 km/mole. The torsional band, a consequence of hindered rotation, is found to be the most intense fundamental with a predicted intensity of 120 km/mole. These results are compared with the recent absolute intensity determinations for the nu(6) band.

Rogers, J. D.↗

Ab initio calculation of infrared intensities for the linear isoelectronic series HCN, HNC, CO, HCO/+/, and HOC/+/

Ab initio infrared intensities and dipole moment derivatives expressed in atomic polar tensor form are calculated using the 4-31 and 6-31G(double asterisk) basis sets for the isoelectronic HCN, HNC, CO, HCO(+), and HOC(+) series of molecules. The calculated atomic polar tensors are analyzed in terms of the charge-charge flux-overlap model, which is found to be useful in explaining some of the trends observed in the dipole moment derivatives for this series of molecules. A detailed examination of the dipole moment derivatives for the structural isomers indicates some of the ways in which experimental atomic polar tensors for one isomer should be modified to predict infrared intensities for the other isomer. The absolute intensities calculated for the HCO(+) and HOC(+) ions are believed to be accurate to within a factor of 2 and thus should be useful in astrophysical applications.

Rogers, J. D.↗

Cosmic microwave background probes models of inflation

Inflation creates both scalar (density) and tensor (gravity wave) metric perturbations. We find that the tensor-mode contribution to the cosmic microwave background anisotropy on large-angular scales can only exceed that of the scalar mode in models where the spectrum of perturbations deviates significantly from scale invariance. If the tensor mode dominates at large-angular scales, then the value of DeltaT/T predicted on 1 deg is less than if the scalar mode dominates, and, for cold-dark-matter models, bias factors greater than 1 can be made consistent with Cosmic Background Explorer (COBE) DMR results.

Davis, Richard L.↗

Optimizing Tensor Contraction Expressions for Hybrid CPU-GPU Execution

Tensor contractions are generalized multidimensional matrix multiplication operations that widely occur in quantum chemistry. Efficient execution of tensor contractions on Graphics Processing Units (GPUs) requires several challenges to be addressed, including index permutation and small dimension-sizes reducing thread block utilization. Moreover, to apply the same optimizations to various expressions, we need a code generation tool. In this paper, we present our approach to automatically generate CUDA code to execute tensor contractions on GPUs, including management of data movement between CPU and GPU. To evaluate our tool, GPU-enabled code is generated for the most expensive contractions in CCSD(T), a key coupled cluster method, and incorporated into NWChem, a popular computational chemistry suite. For this method, we demonstrate speedup over a factor of 8.4 using one GPU (instead of one core per node) and over 2.6 when utilizing the entire system using hybrid CPU+GPU solution with 2 GPUs and 5 cores (instead of 7 cores per node). Finally, we analyze the implementation behavior on future GPU systems.

CUDA↗

Explosive Yield Estimation Using Regional Seismic Moment Tensors

Here, we use the Pasyanos and Chiang (2022) data set to calculate the seismic moment M 0 for each explosion and use the measured explosive yield W to validate the W~M 0 relationship in Denny and Johnson (1991; hereafter, DJ91). The M 0 is corrected by transforming to a potency tensor and applying more appropriate near-source geophysical parameter values in the moment estimate. The mean residual between observed and predicted yield is near zero; however, the standard deviation of the residuals results in an F-value (a 95% confidence factor) of about 5. We re-estimate the coefficients in the DJ91 model and find similar values and only a slight improvement in the F-value. Next, we embark on a similar model selection process as DJ91, allowing for non-cube-root yield scaling and other plausible near-source elastic moduli. As was found by DJ91, the yield dependence is not significantly different from unity, and a cube root assumption is valid. Therefore, we yield scale the seismic moment and test the significance of all plausible explanatory variables. Isotropic moment performs better in the response variable than total moment. The preference for isotropic moment could be due to its relationship to volume change, which would be more directly affected by explosive yield. Surprisingly, we find that the overburden pressure, which is a function of depth, is not a significant parameter in the model. We hypothesize that this is due to the competing depth effects on source asymmetry and the incorporation of depth in the Green’s functions used to calculate the seismic moment tensors. Importantly, this emphasizes that only seismic moment tensor-derived moments should be used in these models. After removing insignificant model parameters, we are left with a simple model to predict explosive yield $\widehat{W}$ in kt from isotropic moment M I in N·m, $\widehat{W}$=κ –1.4132 10 0.035626GP M I , in which κ and GP are the near-source bulk modulus and gas porosity in Pa and %, respectively. The F-value for this model is approximately 3.

58 GEOSCIENCES↗

Multimodality in the Search for New Physics in Pulsar Timing Data and the Case of Kination-amplified Gravitational-wave Background from Inflation

We investigate the kination-amplified inflationary gravitational-wave background (GWB) interpretation of the signal recently reported by various pulsar timing array (PTA) experiments. Kination is a post-inflationary phase in the expansion history dominated by the kinetic energy of some scalar field, characterized by a stiff equation of state w = 1. Within the inflationary GWB model, we identify two modes that can fit the current data sets (NANOGrav and EPTA) with equal likelihood: the kination-amplification (KA) mode and the ordinary, no-kination-amplification (no-KA) mode. The multimodality of the likelihood motivates a Bayesian analysis with nested sampling. We analyze the free spectra of current PTA data and mock free spectra constructed with higher signal-to-noise ratios using nested sampling. The analysis of the mock spectrum designed to be consistent with the best fit to the NANOGrav 15 yr (NG15) data successfully reveals the expected bimodal posterior for the first time while excluding the reheating mode that appears in the fit to the current NG15 data, making a case for our correct and comprehensive treatment of potential multimodal posteriors arising from future PTA data sets. The resultant Bayes factor is $\mathcal{B}$ $\equiv$ Z no–KA /Z KA = 2.9 ± 1.9, indicating comparable statistical significance between the two modes. Given the theoretical model-building challenges of producing highly blue-tilted primordial tensor spectra, the KA mode has the advantage of requiring less blue primordial spectra, compared with the no-KA mode. The synergy between future cosmic microwave background polarization, pulsar timing, and laser interferometer measurements of gravitational waves will help resolve the ambiguity implied by the multimodal posterior in PTA-only searches.

Cosmology↗

Measurements of B-mode polarization of the cosmic microwave background from 500 square degrees of SPTpol data

We report a B-mode power spectrum measurement from the cosmic microwave background (CMB) polarization anisotropy observations made using the SPTpol instrument on the South Pole Telescope. This work uses 500 deg 2 of SPTpol data, a five-fold increase over the last SPTpol B-mode release. As a result, the bandpower uncertainties have been reduced by more than a factor of two, and the measurement extends to lower multipoles: 52 < ` < 2301. Data from both 95 and 150 GHz are used, allowing for three cross-spectra: 95 GHz x 95 GHz, 95 GHz x 150 GHz, and 150 GHz x 150 GHz. B-mode power is detected at very high significance; we find P(BB < 0) = 5.8 × 10 -71 , corresponding to a 18.1 σ detection of power. With a prior on the galactic dust from Planck, WMAP and BICEP2/Keck observations, the SPTpol B-mode data can be used to set an upper limit on the tensor-to-scalar ratio, r < 0.44 at 95% confidence (the expected 1 σ constraint on r given the measurement uncertainties is 0.22). We find the measured B-mode power is consistent with the Planck best-fit ΛCDM model predictions. Scaling the predicted lensing B-mode power in this model by a factor A lens , the data prefer A lens = 1.17 ± 0.13. These data are currently the most precise measurements of B-mode power at ℓ > 320.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamical Theory of Segmental Relaxation and Emergent Elasticity in Supercooled Polymer Melts

We generalize the force-level Elastically Collective Nonlinear Langevin Equation theory of supercooled molecular liquid dynamics to polymer melts based on mapping chains to disconnected and noninterpenetrating Kuhnsized spheres. This allows first-principles, no adjustable parameter calculations to be performed for the temperature-dependent mean segmental relaxation time of chemically diverse van der Waals polymers over a wide range of molecular weights. Despite the simplicity of the mapping, the theory does a good job of a priori predicting the glass transition temperature (T g ), the dynamic fragility, and full temperature dependence of the -relaxation time for some high molecular weight polymers and the chain length dependence of T g as the consequence of the molecular weight dependence of backbone stiffness. The minimalist model does not capture the unusually low and high fragilities of certain long chain polymers which are not typical of van der Waals molecular liquids. This seems likely due to the simple coarse graining adopted which ignores longer range chain connectivity and nonuniversal factors on the sub-Kuhn length scale. Elasticity, not of an entropic single chain origin, emerges in deeply supercooled polymer liquids due to transient segmental localization and is studied at the microscopic stress-tensor level. Calculations of the frequency-dependent dynamic storage modulus close to T g appear to be qualitatively consistent with recent measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Physiological Observations and Omics to Develop Personalized Sensormotor Adaptability Countermeasures Using Bed Rest and Space Flight Data

Astronauts experience sensorimotor disturbances during the initial exposure to microgravity and during the re-adapation phase following a return to an earth-gravitational environment. These alterations may disrupt the ability to perform mission critical functional tasks requiring ambulation, manual control and gaze stability. Interestingly, astronauts who return from space flight show substantial differences in their abilities to readapt to a gravitational environment. The ability to predict the manner and degree to which individual astronauts would be affected would improve the effectiveness of countermeasure training programs designed to enhance sensorimotor adaptability. For such an approach to succeed, we must develop predictive measures of sensorimotor adaptability that will allow us to foresee, before actual space flight, which crewmembers are likely to experience the greatest challenges to their adaptive capacities. The goals of this project are to identify and characterize this set of predictive measures that include: 1) behavioral tests to assess sensory bias and adaptability quantified using both strategic and plastic-adaptive responses; 2) imaging to determine individual brain morphological and functional features using structural magnetic resonance imaging (MRI), diffusion tensor imaging, resting state functional connectivity MRI, and sensorimotor adaptation task-related functional brain activation; 3) genotype markers for genetic polymorphisms in Catechol-O-Methyl Transferase, Dopamine Receptor D2, Brain-derived neurotrophic factor and genetic polymorphism of alpha2-adrenergic receptor that play a role in the neural pathways underlying sensorimotor adaptation. We anticipate these predictive measures will be significantly correlated with individual differences in sensorimotor adaptability after long-duration space flight and an analog bed rest environment. We will be conducting a retrospective study leveraging data already collected from relevant ongoing/completed bed rest and space flight studies. These data will be combined with predictor metrics that will be collected prospectively - behavioral, brain imaging and genomic measures; from these returning subjects to build models for predicting post-mission (bed rest - non-astronauts or space flight - astronauts) adaptive capability as manifested in their outcome measures. Comparisons of model performance will allow us to better design and implement sensorimotor adaptability training countermeasures that are customized for each crewmember's sensory biases, adaptive capacity, brain structure and functional capacities, and genetic predispositions against decrements in post-mission adaptive capability. This ability will allow more efficient use of crew time during training and will optimize training prescriptions for astronauts to ensure expected outcomes.

Mulavara, A. P.↗

Nucleon form factors in dispersively improved chiral effective field theory: Scalar form factor

We propose a method for calculating the nucleon form factors (FFs) of $G$-parity-even operators by combining Chiral Effective Field Theory ($\chi$EFT) and dispersion analysis. The FFs are expressed as dispersive integrals over the two-pion cut at $t > 4 M_\pi^2$. The spectral functions are obtained from the elastic unitarity condition and expressed as products of the complex $\pi\pi \rightarrow N\bar N$ partial-wave amplitudes and the timelike pion FF. $\chi$EFT is used to calculate the ratio of the partial-wave amplitudes and the pion FF, which is real and free of $\pi\pi$ rescattering in the $t$-channel ($N/D$ method). The rescattering effects are then incorporated by multiplying with the squared modulus of the empirical pion FF. The procedure results in a marked improvement compared to conventional $\chi$EFT calculations of the spectral functions. We apply the method to the nucleon scalar FF and compute the scalar spectral function, the scalar radius, the $t$-dependent FF, and the Cheng-Dashen discrepancy. Higher-order chiral corrections are estimated through the $\pi N$ low-energy constants. Results are in excellent agreement with dispersion-theoretical calculations. We elaborate several other interesting aspects of our method. The results show proper scaling behavior in the large-$N_c$ limit of QCD because the $\chi$EFT includes $N$ and $\Delta$ intermediate states. The squared modulus of the timelike pion FF required by our method can be extracted from Lattice QCD calculations of vacuum correlation functions of the operator at large Euclidean distances. Our method can be applied to the nucleon FFs of other operators of interest, such as the isovector-vector current, the energy-momentum tensor, and twist-2 QCD operators (moments of generalized parton distributions).

1=Nc expansion↗

Field validation of data-driven BSDF and peak extraction models for light-scattering fabric shades

Shading and daylighting systems affect cooling, heating, and lighting energy use by modulating solar radiation through the building façade. Characterizing shading systems holistically and accurately helps designers and engineers evaluate shading systems to achieve energy and non-energy performance goals. These complex fenestration systems can be modeled using Bidirectional Scattering Distribution Functions (BSDF), which map incident radiation to hemispherical distributions of outgoing radiation. Data-driven, tabulated BSDFs are derived from interpolated goniophotometer measured data, then sampled during the raytracing calculation. A peak extraction (PE) algorithm was developed to circumvent limits in BSDF angular resolution, where the specular peak is extracted during simulation by evaluating the BSDF in the through direction and surrounding region. The objective of this study was to validate this measurement and modeling workflow using field monitored data from a full scale testbed with eleven installed fabrics of different weaves, openness factors, and colors and assess the accuracy of the workflow under different adaptation and contrast conditions. Test conditions were limited to clear sky conditions with the sun in the field of view. Results showed that, for tensor tree datasets, vertical illuminance, solar luminance (2.5° apex), and daylight glare probability (DGP) were predicted to within a mean bias error (MBE) error of -456 lx (-12.3%), -3.46e5 (-38.4%), and -0.042 (-7.8%) when full PE occurred. With a binary classification of glare/ no glare, DGP was predicted accurately with a true positive rate of 0.98 and true negative rate of 1.0 using tensor tree data and less accurately with Klems BSDF data, particularly for cases of no glare. The workflow may be of insufficient accuracy to distinguish borderline performance between fabrics using the four-point glare scale, particularly under low adaptation, high contrast daylit conditions. Errors were due to reductions in peak shape and intensity across the BSDF interpolation and data reduction workflow. Future work is needed to better preserve measurement fidelity during interpolation and sampling, which in turn will improve PE performance.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Impact of anisotropy on TRISO fuel performance

Manufacturing of tristructural isotropic (TRISO) particles involves the deposition of pyrolytic carbon (PyC) and silicon carbide (SiC) layers using the fluidized bed chemical vapor deposition (CVD) process. The CVD process is known to generate polycrystalline layers with crystallographic textures, which imparts anisotropic thermophysical properties to the layers. Past studies have shown the risk for particle failure increases with an increase in anisotropy. The limit beyond which the anisotropy of PyC layers becomes unacceptable due to failure risk has been identified as a high-priority knowledge gap. This work presents a first systematic study on the effects of anisotropic thermal and mechanical properties on TRISO fuel performance. This computational study, performed using the fuel performance code BISON, investigates how the anisotropy in elasticity and thermal properties affect the stresses, temperature, and failure of a TRISO particle. The influence of other factors, such as operating temperature and particle geometry on the anisotropy effects, also has been analyzed. The studies utilize the recently published anisotropic elasticity and thermal behavior models for TRISO PyC and SiC layers implemented using tensors with full anisotropic capability. The spherical TRISO particles with anisotropic properties were found to have greater maximum tensile stress and significantly higher failure probability than the spherical particles with isotropic properties. In conclusion, the fuel performance predicted using these recently developed models was found to be comparable with the performance obtained using the historical models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Adiabatic regularization of power spectra in nonminimally coupled chaotic inflation

We investigate the effect of adiabatic regularization on both the tensor- and scalar-perturbation power spectra in nonminimally coupled chaotic inflation. Similar to that of the minimally coupled general single-field inflation, we find that the subtraction term is suppressed by an exponentially decaying factor involving the number of e -folds. By following the subtraction term long enough beyond horizon crossing, the regularized power spectrum tends to the ''bare'' power spectrum. This study justifies the use of the unregularized (''bare'') power spectrum in standard calculations.

79 ASTRONOMY AND ASTROPHYSICS↗

Logamediate Inflation in f ( T ) Teleparallel Gravity

We study logamediate inflation in the context of f ( T ) teleparallel gravity. f ( T )-gravity is a generalization of the teleparallel gravity which is formulated on the Weitzenbock spacetime, characterized by the vanishing curvature tensor (absolute parallelism) and the non-vanishing torsion tensor. We consider an f ( T )-gravity model which is sourced by a canonical scalar field. Assuming a power-law f ( T ) function in the action, we investigate an inflationary universe with a logamediate scale factor. Our results show that, although logamediate inflation is completely ruled out by observational data in the standard inflationary scenario based on Einstein gravity, it can be compatible with the 68% confidence limit joint region of Planck 2015 TT,TE,EE+lowP data in the framework of f ( T )-gravity.

79 ASTRONOMY AND ASTROPHYSICS↗