SEARCH · Engineering Papers
Results for “Data 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.
Removing flat directions in standard model EFT fits: How polarized electron-ion collider data can complement the LHC
We study the potential of future electron-ion collider (EIC) data to probe four-fermion operators in the Standard Model effective field theory (SMEFT). The ability to perform measurements with both polarized electron and proton beams at the EIC provides a powerful tool that can disentangle the effects from different SMEFT operators. We compare the potential constraints from an EIC with those obtained from Drell-Yan data at the Large Hadron Collider. We show that EIC data play an important complementary role since they probe combinations of Wilson coefficients not accessible through available Drell-Yan measurements.
Measurements of polarization and spin correlation and observation of entanglement in top quark pairs using lepton + jets events from proton-proton collisions at $\sqrt{𝑠}$ = 13 TeV
Measurements of the polarization and spin correlation in top quark pairs ($t\bar{t}$) are presented using events with a single electron or muon and jets in the final state. The measurements are based on proton-proton collision data from the LHC at $\sqrt{𝑠}$ =13 TeV collected by the CMS experiment, corresponding to an integrated luminosity of 138 fb −1 . All coefficients of the polarization vectors and the spin correlation matrix are extracted simultaneously by performing a binned likelihood fit to the data. The measurement is performed inclusively and in bins of additional observables, such as the mass of the $t\bar{t}$ system and the top quark scattering angle in the $t\bar{t}$ rest frame. The measured polarization and spin correlation are in agreement with the standard model. From the measured spin correlation, conclusions on the $t\bar{t}$ spin entanglement are drawn by applying the Peres-Horodecki criterion. The standard model predicts entangled spins for $t\bar{t}$ states at the production threshold and at high masses of the $t\bar{t}$ system. Entanglement is observed for the first time in events at high $t\bar{t}$ mass, where a large fraction of the $t\bar{t}$ decays are spacelike separated, with an expected and observed significance of above 5 standard deviations.
The relative line strength and intensity of the N II 2143 doublet
The doublet emission from N II at 2139.7 A and 2143.6 A was observed by a 1.4-m scanning spectrometer with 3.1 A resolution in the daytime, high-altitude thermosphere during moderate levels of solar activity. The spectrometer viewed the earth's limb 5 deg below the local horizontal to give a nominal tangent height of 152 km. Both sub band heads of the nitric oxide gamma band system were resolved in the data at the resolution used. The emission features from N II are clearly evident on the short wavelength shoulder of the (1, 0) band. Synthetic profiles of the (1, 0) gamma band and the (0, 3) delta band of nitric oxide were fitted to the data using a chi-square analysis. These contributions were removed from the data leaving a residual emission, considered to be the N II doublet. A chi-square minimization of the data relative to a synthetic intensity profile was done. The minimum was for a line strength ratio between the 2139 A and 2143 A lines of 0.58 + or - 0.08. The mean solar EUV flux deduced from the intensity of the N II emission in this experiment is lower than other reported observations, consistent with a lower solar activity level.
Identifying anomalous and non-sulfate aerosol layers through Angstrom exponent goodness-of-fit: proof of concept from in situ data
Explore the source record for details and available documents.
SparTen: Leveraging Kokkos for On-node Parallelism in a Second-Order Method for Fitting Canonical Polyadic Tensor Models to Poisson Data.
Abstract not provided.
SparTen: Leveraging Kokkos for On-node Parallelism in a Second-Order Method for Fitting Canonical Polyadic Tensor Models to Poisson Data.
Abstract not provided.
Uranium Retardation Capacity of Lithologies from the Negev Desert, Israel—Rock Characterization and Sorption Experiments
A series of batch experiments were performed to assess the uranium sorption capacity of four mineralogically distinct lithologies from the Negev Desert, Israel, to evaluate the suitability of a potential site for subsurface radioactive waste disposal. The rock specimens consisted of an organic-rich phosphorite, a bituminous marl, a chalk, and a sandstone. The sorption data for each lithology were fitted using a general composite surface complexation model (GC SCM) implemented in PHREEQC. Sorption data were also fitted by a non-mechanistic Langmuir sorption isotherm, which can be used as an alternative to the GC SCM to provide a more computationally efficient method for uranium sorption. This is because all the rocks tested have high pH/alkalinity/calcium buffering capacities that restrict groundwater chemistry variations, so that the use of a GC SCM is not advantageous. The mineralogy of the rocks points to several dominant sorption phases for uranyl (UO 2 2+ ), including apatite, organic carbon, clays, and iron-bearing phases. The surface complexation parameters based on literature values for the minerals identified overestimate the uranium sorption capacities, so that for our application, an empirical approach that makes direct use of the experimental data to estimate mineral-specific sorption parameters appears to be more practical for predicting uranium sorption.
Updating High-Order Aeroservoelastic Models from Low-Order System Identification Results
Estimating aircraft models from test data requires several simplifying assumptions that introduce biases into the parameter estimates. In this paper, these biases are defined and a method for estimating the biases is discussed. Having an estimate of the bias allows the parameters estimated from test to be integrated into a high-order model. The Integrated Adaptive Wing Technology Maturation (IAWTM) wind tunnel model is discussed and the bias is demonstrated for one of the testing configurations. The methodology was able to estimate these biases and apply corrections to high-order aeroelastic models to improve the fit to test data. The consideration of the biases allows more meaningful comparisons and avoids the erroneous differences between pretest predictions and the fitted model.
Updating High-Order Aeroservoelastic Models from Low-Order System Identification Results
Estimating aircraft models from test data requires several simplifying assumptions that introduce biases into the parameter estimates. In this paper, these biases are defined and a method for estimating the biases is discussed. Having an estimate of the bias allows the parameters estimated from test to be integrated into a high-order model. The Integrated Adaptive Wing Technology Maturation (IAWTM) wind tunnel model is discussed and the bias is demonstrated for one of the testing configurations. The methodology was able to estimate these biases and apply corrections to high-order aeroelastic models to improve the fit to test data. The consideration of the biases allows more meaningful comparisons and avoids the erroneous differences between pretest predictions and the fitted model.
PolyFit: A C++ code for Polynomial Curve Fit with Calculation of Error Bars
In radiobiology, many dose-response results are modeled using the so-called linear-quadratic (LQ)model, which means that results are modeled as a function of dose Das 𝑅(𝐷)=𝛽0+𝛽1𝐷+𝛽2𝐷2. The coefficients 𝛽0, 𝛽1and 𝛽2are obtained from fitting a series of data points (𝑥𝑖, 𝑦𝑖), which is usually done using a least-square method. The LQ and more generally the polynomial fit capability is implemented in many software that analyzes data. However, there are some instances where the fitting needs to be done programmatically. Furthermore, depending on the software used, some features may not be implemented. In this mini-review, I discuss the basis of polynomial fitting, including the calculation of errors on the coefficients and results, use of weighting and fixing the intercept value (the coefficient 𝛽0). A simple C++ code to perform the polynomial curve fitting is also provided. This code should be useful not only in radiobiology but in other fields of science as well.
A Markov chain Monte Carlo (MCMC) Bayesian inference approach to analyze apparent activation barriers and reaction orders from microreactor data
Statistical analysis of steady-state catalytic kinetic data is often limited by data sparsity due to the slow pace at which the data is collected. Data sparsity and limitations in statistical analysis make it difficult to differentiate between mechanistic models and catalytic sites. A Bayesian inference tool is reported for catalysis researchers to estimate error in the determination of reaction orders from steady state microreactor data. The benefits of a Bayesian inference approach are discussed, as an alternative to the more common frequentist approach. The approach incorporates prior knowledge of the system and the data collected to form an error estimate on reaction orders. We investigated the effects of three distinct data treatments—individual fitting of trials, pooled analysis, and constrained regression methods—on the precision and uncertainty of reaction order determinations. To assess the robustness of our findings, we conducted sensitivity analyses to evaluate the influence of Bayesian parameters on uncertainty estimation. Additionally, we utilized synthetic data to illustrate how data quality impacts the precision of uncertainty assessments. We show Bayesian analysis can obtain a more precise estimation of error with a sparse data set than a frequentist analysis. Finally, this work provides strong evidence that the adoption of Bayesian analysis of kinetic data may help researchers make more precise arguments as to the strength of their evidence for a particular mechanistic hypothesis, or in comparing across different catalysts.
Recent Advances in Small Angle X-ray Scattering for Superlattice Study
Small-angle x-ray scattering is used for the structure determination of superlattice for its superior resolution, nondestructive nature, and high penetration power of x rays. With the advent of high brilliance x-ray sources and innovative computing algorithms, there have been notable advances in small angle x-ray scattering analysis of superlattices. High brilliance x-ray beams have made data analyses less model-dependent. Additionally, novel data acquisition systems are faster and more competitive than ever before, enabling a more accurate mapping of the superlattices' reciprocal space. Fast and high-throughput computing systems and algorithms also make possible advanced analysis methods, including iterative phasing algorithms, non-parameterized fitting of scattering data with molecular dynamics simulations, and the use of machine learning algorithms. As a result, solving nanoscale structures with high resolutions has become an attainable task. In this review, we highlight new developments in the field and introduce their applications for the analysis of nanoscale ordered structures, including nanoparticle supercrystals, nanoscale lithography patterns, and supramolecular self-assemblies. Particularly, we highlight the reciprocal space mapping techniques and the use of iterative phase retrieval algorithms. We also cover coherent-beam-based small angle x-ray scattering techniques such as ptychography and ptycho-tomography in view of the traditional small angle x-ray scattering perspective.
A method for determining crack opening load from load-displacement data
The fatigue crack opening load is determined as the tangent point on the ascending load-displacement data between the curved portion and the upper linear region. A model for the 'unzipping' behavior of the crack indicates that the curved portion of the load-deflection curve is second order. The opening load is determined by a nonlinear, least squares fit of the data to the model, which optimally locates the tangent point of the two curves. The method provided consistent results for determining opening load P(op) for 7475-T731 aluminum using data from a crack tip opening gage.
Remote and Ground Truth Spectral Measurement Comparisons of FORMOSAT III
FORMOSAT III are a set of six research satellites from Taiwan that were launched in April 2006. The satellites are in 800 km, 71 degree inclination orbits and separated by 24 degrees in ascending node. Laboratory spectral measurements were taken of outer surface materials on FORMOSAT III. From those measurements, a computer model was built to predict the spectral reflectance accounting for both solar phase angle and orientation of the spacecraft relative to the observer. However, materials exposed to the space environment have exhibited spectral changes including a darkening and a "reddening" of the spectra. This "reddening" is characterized by an increase in slope of the reflectance as the wavelength increases. Therefore, the model of pre-flight materials was augmented to include the presumed causative agent: space weathering effects. Remote data were collected on two of the six FORMOSAT satellites using the 1.6 meter telescope at the AMOS (Air Force Maui Optical and Supercomputing) site with the Spica spectrometer. Due to the separation in ascending node, observations were acquired of whichever one of the six satellites was visible on that specific night. Three nights of data were collected using the red (6000 - 9500 angstroms) filter and five nights of data were collected using the blue (3200 - 6600 angstroms) filter. A comparison of the data showed a good match to the pre-flight models for the blue filter region. The absorption feature near 5500 angstroms due to the copper colored Kapton multi-layer insulation (MLI) was very apparent in the remote samples and a good fit to the data was seen in all satellites observed. The features in the red filter regime agreed with the pre-flight model up through 7000 angstroms where the reddening begins and the slope of the remote sample increases. A comparison of the satellites showed similar features in the red and blue filter regions, i.e. the satellite surfaces were aging at the same rate. A comparison of the pre-flight model to the first month of remote measurements showed the amount by which the satellite had reddened. The second month of data observed a satellite at a higher altitude and was therefore, not compared to the first month. A third month of data was collected but of satellites at the lower altitude regime and can only be compared to the first month. One cause of the reddening that was ruled out in early papers was a possible calibration issue.
Nonmonotonic-potential description of polarization effects, fusion, and nuclear rainbows in elastic scattering of 6 Li + 12 C at 4.5–600 MeV
The experimental differential cross-section (CS) and analyzing power (AP) data of the 6 Li + 12 C elastic scattering over a wide laboratory energy scale (4.5MeV≤𝐸 lab ≤600MeV) are analyzed within the framework of the optical model (OPM) using nonmonotonic (NM) nucleus–nucleus potentials. The real part of the NM potentials is derived from the Pauli-embodied energy density-functional (EDF) formalism with the sudden approximation. The real part of the noncentral spin-orbit and tensor terms, as well as the imaginary parts, are treated phenomenologically. The effect of the radius of sensitivity on the CS and AP data is found to be more important at lower energies. The diffractive and refractive scattering with Airy structures in the whole angular region of the elastic scattering across the studied energy range is successfully described within the OPM using the NM 6 Li + 12 C potential. The near- and far-side (N and F) decomposition of the total elastic-scattering amplitudes has also been studied using our NM potentials. The evolution of the Airy minima in the angular distributions, coupled with the fitting of the AP data, provides an accurate description of Airy minima of different orders. The OPM calculations with the NM potentials describe exceptionally well the CS, vector analyzing power (VAP), and tensor analyzing power data at 𝐸 lab =9.0,19.24,20,30, and 50MeV. In agreement with our past successful descriptions of CS and the opposite signs of the VAP data for the 6 Li and 7 Li elastic scattering using NM potentials in OPM, the present results appear to provide a better fit, so far, than those obtained from the coupled-channels method. The fusion cross sections of 6 Li + 12 C have been predicted in the energy range (4.5MeV≤𝐸 lab ≤20MeV), fitting the experimental data well in the range 𝐸 lab =2.97–11.77MeV. The EDF potential without any energy dependence and renormalization is also found to describe satisfactorily the experimental CS and AP data at energies up to several hundreds of MeV.
Modified 316H constitutive model updated with high temperature stress relaxation test data
This report describes the calibration of a new high temperature constitutive model for 316H stainless steel, suitable for use with the ASME Boiler & Pressure Vessel Section III, Division 5, Class A rules for design by inelastic analysis. The model retains the same mathematical form used by the reference model included in Nonmandatory Appendix Z of the Code, but refits the model to an expanded dataset including all the data used to fit the original model plus seven new stress relaxation tests. The addition of these high temperature stress relaxation tests improves the model's accuracy in predicting relaxation at temperatures greater than 700 ⁰C, without compromising the accuracy of the model versus the original calibration data.
Recovery of atmospheric refractivity profiles from simulated satellite-to-satellite tracking data
Techniques for recovering atmospheric refractivity profiles from simulated satellite-to-satellite tracking data are documented. Examples are given using the geometric configuration of the ATS-6/NIMBUS-6 Tracking Experiment. The underlying refractivity model for the lower atmosphere has the spherically symmetric form N = exp P(s) where P(s) is a polynomial in the normalized height s. For the simulation used, the Herglotz-Wiechert technique recovered values which were 0.4% and 40% different from the input values at the surface and at a height of 33 kilometers, respectively. Using the same input data, the model fitting technique recovered refractivity values 0.05% and 1% different from the input values at the surface and at a height of 50 kilometers, respectively. It is also shown that if ionospheric and water vapor effects can be properly modelled or effectively removed from the data, pressure and temperature distributions can be obtained.