Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vector fitting”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Comparative Study of Wind Energy Potential Estimation Methods for Wind Sites in Togo and Benin (West Sub-Saharan Africa)

The characterization of wind speed distribution and the optimal assessment of wind energy potential are critical factors in selecting a suitable site for wind power plants (WPP). The Weibull distribution law has been used extensively to analyze the wind characteristics of candidate WPP sites, and to estimate the available and deliverable energy. This paper presents a comparative study of five wind energy resource assessment methods as they applied to the context of wind sites in West Sub-Saharan Africa. We investigated three numerical approaches, namely, the adaptive neuro-fuzzy inference system (ANFIS), the multilayer perceptron method (MLP), and support vector regression (SVR), to derive the distribution law of wind speeds and to optimally quantify the corresponding wind energy potential. Next, we compared these three approaches to two well-known Weibull distribution law-based methods: the empirical method of Justus (EMJ) and the maximum likelihood method (MLM). Case study results indicated that the neural network-based methods, ANFIS and MLP, yielded the most accurate distribution fits and wind energy potential estimates, and consequently, are the most recommended methods for the wind sites in Togo and Benin. The orders of magnitude of the root mean squared error (RMSE) in estimating the recoverable energy using ANFIS were, respectively, 10-4 and 10-5 for Lomé and Cotonou, while MLP achieved an RMSE order of magnitude of 10-3 for both sites.

17 WIND ENERGY↗
DOE OSTI2024

Efficiently evaluating loop integrals in the EFTofLSS using QFT integrals with massive propagators

We develop a new way to analytically calculate loop integrals in the Effective Field Theory of Large Scale-Structure. Previous available methods show severe limitations beyond the one-loop power spectrum due to analytical challenges and computational and memory costs. Our new method is based on fitting the linear power spectrum with cosmology-independent functions that resemble integer powers of quantum field theory massive propagators with complex masses. A remarkably small number of them is sufficient to reach enough accuracy. Similarly to former approaches, the cosmology dependence is encoded in the coordinate vector of the expansion of the linear power spectrum in our basis. We first produce cosmology-independent tensors where each entry is the loop integral evaluated on a given combination of basis vectors. For each cosmology, the evaluation of a loop integral amounts to contracting this tensor with the coordinate vector of the linear power spectrum. The 3-dimensional loop integrals for our basis functions can be evaluated using techniques familiar to particle physics, such as recursion relations and Feynman parametrization. We apply our formalism to evaluate the one-loop bispectrum of galaxies in redshift space. The final analytical expressions are quite simple and can be evaluated with little computational and memory cost. We show that the same expressions resolve the integration of all one-loop N-point function in the EFTofLSS. This method, which is originally presented here, has already been applied in the first one-loop bispectrum analysis of the BOSS data to constraint ΛCDM parameters and primordial non-Gaussianities [1, 2].

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Speculations on the W-mass measurement at CDF*

Abstract The W mass determination at the Tevatron CDF experiment reported a deviation from the SM expectation at the 7 σ level. We discuss a few possible interpretations and their collider implications. We perform electroweak global fits under various frameworks and assumptions. We consider three types of electroweak global fits in the effective-field-theory framework: the S - T , S - T - , and eight-parameter flavor-universal one. We discuss the amounts of tensions between different measurements reflected in these fits and the corresponding shifts in central values of these parameters. With these electroweak fit pictures in hand, we present a few different classes of models and discuss their compatibility with these results. We find that while explaining the discrepancy, the single gauge boson extensions face strong LHC direct search constraints unless the is fermiophobic (leptophobic), which can be realized if extra vector fermions (leptons) mix with the SM fermions (leptons). Vector-like top partners can partially generate the needed shift to the electroweak observables. The compatibility with the top squark is also studied in detail. We find that the non-degenerate top squark soft masses enhance the needed operator coefficients, enabling an allowed explanation compatible with current LHC measurements. Overall, more theoretical and experimental developments are highly in demand to reveal the physics behind this discrepancy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Modeling Stochastic Variability in Multiband Time-series Data

In preparation for the era of time-domain astronomy with upcoming large-scale surveys, we propose a state-space representation of a multivariate damped random walk process as a tool to analyze irregularly-spaced multifilter light curves with heteroscedastic measurement errors. We adopt a computationally efficient and scalable Kalman filtering approach to evaluate the likelihood function, leading to maximum O(k 3 n) complexity, where k is the number of available bands and n is the number of unique observation times across the k bands. This is a significant computational advantage over a commonly used univariate Gaussian process that can stack up all multiband light curves in one vector with maximum O(k 3 n 3 ) complexity. Using such efficient likelihood computation, we provide both maximum likelihood estimates and Bayesian posterior samples of the model parameters. Three numerical illustrations are presented: (i) analyzing simulated five-band light curves for a comparison with independent single-band fits; (ii) analyzing five-band light curves of a quasar obtained from the Sloan Digital Sky Survey Stripe 82 to estimate short-term variability and timescale; (iii) analyzing gravitationally lensed g- and r-band light curves of Q0957+561 to infer the time delay. Two R packages, Rdrw and timedelay, are publicly available to fit the proposed models.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1985

Development of an explicit multigrid algorithm for quasi-three-dimensional viscous flows in turbo-machinery

A rapid quasi three-dimensional analysis was developed for blade-to-blade flows in turbomachinery. The analysis solves the unsteady Euler or thin layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stream-surface thickness. The Baldwin-Lomax eddy-viscosity model is used for turbulent flows. The equations which are solved by a two-stage Runge-Kutta scheme made efficient by use of vectorization, a variable time-step, and a flux-based multigrid scheme, are described. A stability analysis is presented for the two-stage scheme. Results for a flat-plate model problem show the applicability of the method to axial, radial, and rotating geometries. Results for a centrifugal impeller and a radial diffuser show that the quasi three-dimensional viscous analysis can be a practical design tool.

Chima, R. V.↗
NASA NTRS1986

Development of an explicit multigrid algorithm for quasi-three-dimensional viscous flows in turbomachinery

A rapid quasi three-dimensional analysis was developed for blade-to-blade flows in turbomachinery. The analysis solves the unsteady Euler or thin layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stream-surface thickness. The Baldwin-Lomax eddy-viscosity model is used for turbulent flows. The equations which are solved b a two-stage Runge-Kutta scheme made efficient by use of vectorization, a variable time-step, and a flux-based multigrid scheme, are described. A stability analysis is presented for the two-stage. Results for a flat-plate model problem show the applicability of the method to axial, radial, and rotating geometries. Results for a centrifugal impeller and a radial diffuser show that the quasi three-dimensional viscous analysis can be a practical design tool.

Chima, R. V.↗
NASA NTRS1996

Evaluation of IGS Orbits with Satellite Laser Ranging

The accuracy with which orbits for the Global Positioning System (GPS) spacecraft, can be computed directly affects the accuracy of the resulting site coordinates and polar motion. Several groups routinely analyze GPS ground tracking data to compute precise orbits and terrestrial reference frame solutions. In this paper, we infer the accuracy of the orbits of two of the GPS satellites by comparing to independent laser ranges of subcentimeter accuracy obtained by a small but reasonably well distributed network of tracking sites. We find that all seven International GPS Service for Geodynamics (IGS) analysis centers achieve range residual root mean square (rms) errors at or below the 100 mm level. The best orbit solutions, from JPL, CODE, and the IGS combined product, yield a residual rms of about 50 mm. These residuals are consistent with three dimensional orbit errors of less than 150 mm. Estimating yaw rates for the spacecraft during shadow events, and using these estimates to compute the laser residual, significantly improves the fit. A small mean residual value of -15 to -30 mm seems to exist for most centers and laser sites which is not fully explained at present, but may be due to uncertainties in the corrections to the laser data, such as the reflector to spacecraft center of mass vector or small reference frame differences between the SLR sites and the GPS orbits.

Watkins, M. M.↗
DOE OSTI2026

Two-pion contribution to the hadronic vacuum polarization with staggered quarks

We present results from the first lattice QCD calculation of the two-pion contributions to the light-quark connected vector-current correlation function obtained from staggered-quark operators. We employ the MILC Collaboration’s gauge-field ensemble with 2 + 1 + 1 flavors of highly improved staggered sea quarks at a lattice spacing of a ≈ 0.15 fm with a light sea-quark mass at its physical value. The two-pion contributions allow for a refined determination of the noisy long-distance tail of the vector-current correlation function, which we use to compute the light-quark connected contribution to hadronic vacuum polarization (HVP) with improved statistical precision. We compare our results with traditional noise-reduction techniques used in lattice QCD calculations of the light-quark connected HVP, namely, the so-called fit and bounding methods. We observe a factor of roughly 3 improvement in the statistical precision in the determination of the HVP contribution to the muon’s anomalous magnetic moment over these approaches. We also lay the group theoretical groundwork for extending this calculation to finer lattice spacings with increased numbers of staggered two-pion taste states.

Lahert, Shaun [Utah U.; Illinois U., Urbana] (ORCI↗
DOE OSTI2021

Band mixing in 96,98 Mo isotopes *

Abstract We use the two lowest weight states to fit E 2 strengths connecting the and transitions in Mo. Our results confirm that the and states are maximally mixed, and that the states are weakly mixed in both nuclei. An appropriate Hamiltonian to represent the band mixing is found to be exactly solvable, and its eigenstates can be expressed as the basis vectors in the configuration mixing scheme and interacting boson model. The interacting boson model and coexistence mixing configuration under the solvable methods are suitable models for analyzing the band mixing with high accuracy.

Physics↗
NASA NTRS1975

Time-controlled descent guidance in uncertain winds

A procedure has been developed for constructing a statistical model of the altitude-dependent mean wind profile from the historical record of wind measurements at particular locations. The model is constructed by fitting a Markov process, with altitude as the stage variable, to the historical wind data. The wind model, together with the aircraft dynamics and the error characteristics of the navigation system, are incorporated in the design of a state estimator, which gives the minimum variance estimate of the aircraft state and the wind vector. The state and wind estimates are used as inputs to a linear feedback law for guiding the aircraft along the nominal trajectory. An example design of a time-constrained (4D RNAV) descent guidance system is presented, showing tracking accuracy, control activity, and probability of arrival time with and without the wind estimator.

Menga, G.↗
DOE OSTI2022

In Silico Prediction of the Toxicity of Nitroaromatic Compounds: Application of Ensemble Learning QSAR Approach

In this work, a dataset of more than 200 nitroaromatic compounds is used to develop Quantitative Structure–Activity Relationship (QSAR) models for the estimation of in vivo toxicity based on 50% lethal dose to rats (LD 50 ). An initial set of 4885 molecular descriptors was generated and applied to build Support Vector Regression (SVR) models. The best two SVR models, SVR_A and SVR_B, were selected to build an Ensemble Model by means of Multiple Linear Regression (MLR). The obtained Ensemble Model showed improved performance over the base SVR models in the training set (R 2 = 0.88), validation set (R 2 = 0.95), and true external test set (R 2 = 0.92). The models were also internally validated by 5-fold cross-validation and Y-scrambling experiments, showing that the models have high levels of goodness-of-fit, robustness and predictivity. The contribution of descriptors to the toxicity in the models was assessed using the Accumulated Local Effect (ALE) technique. The proposed approach provides an important tool to assess toxicity of nitroaromatic compounds, based on the ensemble QSAR model and the structural relationship to toxicity by analyzed contribution of the involved descriptors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Impact of W and Z Production Data and Compatibility of Neutrino DIS Data in Nuclear Parton Distribution Functions

Vector boson production and neutrino deep-inelastic scattering (DIS) data are crucial for constraining the strange quark parton distribution function (PDF) and more generally for flavor decomposition in PDF extractions. We extend the nCTEQ15 nuclear PDFs (nPDFs) by adding the recent W and Z production data from the LHC in a global nPDF fit. The new nPDF set, referred to as nCTEQ15WZ, is used as a starting point for a follow-up study in which we assess the compatibility of neutrino DIS data with charged lepton DIS data. Specifically, we re-analyze neutrino DIS data from NuTeV, Chorus, and CDHSW, as well as dimuon data from CCFR and NuTeV. To scrutinize the level of compatibility, different kinematic regions of the neutrino data are investigated. Fits to the neutrino data alone and a preliminary global fit are performed and compared to nCTEQ15WZ.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Dynamic logistic regression and variable selection: Forecasting and contextualizing civil unrest

Civil unrest can range from peaceful protest to violent furor, and researchers are working to monitor, forecast, and assess such events to allocate resources better. Twitter has become a real-time data source for forecasting civil unrest because millions of people use the platform as a social outlet. Additionally, daily word counts are used as model features, and predictive terms contextualize the reasons for the protest. To forecast civil unrest and infer the reasons for the protest, we consider the problem of Bayesian variable selection for the dynamic logistic regression model and propose using penalized credible regions to select parameters of the updated state vector. This method avoids the need for shrinkage priors, is scalable to high-dimensional dynamic data, and allows the importance of variables to vary in time as new information becomes available. A substantial improvement in both precision and F1-score using this approach is demonstrated through simulation. Finally, we apply the proposed model fitting and variable selection methodology to the problem of forecasting civil unrest in Latin America. Our dynamic logistic regression approach shows improved accuracy compared to the static approach currently used in event prediction and feature selection.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Structure Sensitive Reaction Kinetics of Chiral Molecules on Intrinsically Chiral Surfaces

Enantiospecific heterogeneous catalysis utilizes chiral surfaces to resolve enantiomers via structure sensitive surface chemistry. The catalyst design challenge is the identification of chiral surface structures that maximize enantiospecificity. Herein, we develop data driven models for the enantiospecificity of tartaric acid reactions on chiral Cu(hkl) R&S surfaces. Measurements of enantiospecific rate constants were obtained by using curved Cu(hkl) R&S surfaces that enable kinetic measurements on hundreds of chiral surface orientations. One model uses feature vectors derived from generalized coordination numbers to capture the local structure around Cu atoms exposed by the Cu(hkl) R&S surfaces. The second model introduces the use of chiral cubic harmonic functions to capture the symmetry constraints of the face-centered cubic Cu structure. The model using 58 generalized coordination numbers has a fitting error similar to that of the model using only 5 cubic harmonic functions. The two models predict maxima in the enantiospecificity on surfaces with very similar surface orientations. The models developed in this work are applicable for any enantiospecific reaction happening on any chiral material with a cubic lattice structure, opening the way to understanding the surface structure sensitivity of the enantiospecific reaction kinetics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Study of High-Transverse-Momentum Higgs Boson Production in Association with a Vector Boson in the qqbb Final State with the ATLAS Detector

This Letter presents the first study of Higgs boson production in association with a vector boson (V=W or Z) in the fully hadronic qqbb final state using data recorded by the ATLAS detector at the LHC in proton-proton collisions at $\sqrt{s}$=13 TeV and corresponding to an integrated luminosity of 137 fb –1 . The vector bosons and Higgs bosons are each reconstructed as large-radius jets and tagged using jet substructure techniques. Dedicated tagging algorithms exploiting b-tagging properties are used to identify jets consistent with Higgs bosons decaying into $b\bar{b}$. Dominant backgrounds from multijet production are determined directly from the data, and a likelihood fit to the jet mass distribution of Higgs boson candidates is used to extract the number of signal events. The VH production cross section is measured inclusively and differentially in several ranges of Higgs boson transverse momentum: 250–450, 450–650, and greater than 650 GeV. The inclusive signal yield relative to the standard model expectation is observed to be μ=1.4$^{+1.0}_{–0.9}$ and the corresponding cross section is 3.1±1.3⁢(stat)$^{+1.8}_{–1.4}$(syst.) pb.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Computing rank‐revealing factorizations of matrices stored out‐of‐core

This paper describes efficient algorithms for computing rank-revealing factorizations of matrices that are too large to fit in main memory (RAM), and must instead be stored on slow external memory devices such as disks (out-of-core or out-of-memory). Traditional algorithms for computing rank-revealing factorizations (such as the column pivoted QR factorization and the singular value decomposition) are very communication intensive as they require many vector-vector and matrix-vector operations, which become prohibitively expensive when data is not in RAM. Randomization allows to reformulate new methods so that large contiguous blocks of the matrix are processed in bulk. The paper describes two distinct methods. The first is a blocked version of column pivoted Householder QR, organized as a “left-looking” method to minimize the number of the expensive write operations. The second method results employs a UTV factorization. It is organized as an algorithm-by-blocks to overlap computations and I/O operations. As it incorporates power iterations, it is much better at revealing the numerical rank. Numerical experiments on several computers demonstrate that the new algorithms are almost as fast when processing data stored on slow memory devices as traditional algorithms are for data stored in RAM.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Improving resbos for the precision needs of the LHC

The resummation calculation (esos) is a widely used tool for the simulation of single vector boson production at colliders. In this work, we develop a significant improvement over the esos code by increasing the accuracy from NNLL + NLO to N 3 LL + NNLO and release the esos v2.0 code. Furthermore, we propose a new nonperturbative function that includes information about the rapidity of the system (IFY). The IFY functional form was fitted to data from fixed target experiments, the Tevatron, and the LHC. We find that the nonperturbative function has mild rapidity dependence based on the results of the fit. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1988

North America-Pacific plate boundary, an elastic-plastic megashear - Evidence from very long baseline interferometry

Data obtained by Mark III VLBI measurements of radio signals from permanent and mobile VLBI sites for 5.5 years of observations, starting in October 1982, were used to derive a picture of the earth crust deformation near the North America-Pacific plate boundary. The data, which included the vector positions of the VLBI sites and their rate of change, were used for comparison with a number of lithospheric deformation models based upon the concept that the motions of points near the North America-Pacific plate boundary are a linear combination of North America and Pacific velocities. The best of these models were found to fit 95 percent of the variance in 139 VLBI length and transverse velocity observations. Instantaneous shear deformation associated with plate tectonics is apparently developing in a zone 450 km wide paralleling the San Andreas Fault; some of this deformation will be recovered through elastic rebound, while the rest will be permanently set through plastic processes. Because the VLBI data have not been collected for a significant fraction of the earthquake cycle, they cannot discriminate between elastic and plastic behaviors.

Ward, Steven N.↗