Engineering Papers⌕ Search

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.

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Detection of the Pairwise Kinematic Sunyaev-Zel'dovich Effect with BOSS DR11 and the Atacama Cosmology Telescope

We present a new measurement of the kinematic Sunyaev-Zel'dovich effect using data from the Atacama Cosmology Telescope (ACT) and the Baryon Oscillation Spectroscopic Survey (BOSS). Using 600 square degrees of overlapping sky area, we evaluate the mean pairwise baryon momentum associated with the positions of 50,000 bright galaxies in the BOSS DR11 Large Scale Structure catalog. A non-zero signal arises from the large-scale motions of halos containing the sample galaxies. The data fits an analytical signal model well, with the optical depth to microwave photon scattering as a free parameter determining the overall signal amplitude. We estimate the covariance matrix of the mean pairwise momentum as a function of galaxy separation, using microwave sky simulations, jackknife evaluation, and bootstrap estimates. The most conservative simulation-based errors give signal-to-noise estimates between 3.6 and 4.1 for varying galaxy luminosity cuts. We discuss how the other error determinations can lead to higher signal-to-noise values, and consider the impact of several possible systematic errors. Estimates of the optical depth from the average thermal Sunyaev-Zel'dovich signal at the sample galaxy positions are broadly consistent with those obtained from the mean pairwise momentum signal.

De Bernardis, F.↗
NASA NTRS2019

Multi-Arc Filtering During the Navigation Campaign of the OSIRIS-REx Mission

The Navigation Campaign of the OSIRIS-REx mission consisted of three phases: Approach, Preliminary Survey and Orbital-A. These phases were designed to optimize the initial characterization of Bennu's mass, shape and spin-state to support a safe orbit insertion and a quick transition to landmark-based optical navigation tracking. The standard orbit determination filtering techniques used to navigate the spacecraft were unable to fit data from these three phases simultaneously due to numerical issues associated with the nonlinear dynamics and the long arc length. Consequently, a multi-arc filtering algorithm was implemented in order to combine the information from each of these arcs. Multi-arc solutions for Bennu's spin state and gravity field are presented here.

French, Andrew S.↗
NASA NTRS2019

An ALMA+ACA measurement of the shock in the Bullet Cluster

Context: The thermal Sunyaev-Zeldovich (SZ) effect presents a relatively new tool for characterizing galaxy cluster merger shocks, traditionally studied through X-ray observations. Widely regarded as the “textbook example” of a cluster merger bow shock, the western, most-prominent shock front in the Bullet Cluster (1E0657-56) represents the ideal test case for such an SZ study. Aims: We aim to characterize the shock properties using deep, high-resolution interferometric SZ effect observations in combination with priors from an independent X-ray analysis. Methods: Our analysis technique relies on the reconstruction of a parametric model for the SZ signal by directly and jointly fitting data from the Atacama Large Millimeter/submillimeter Array (ALMA) and Atacama Compact Array (ACA) in Fourier space. Results: The ALMA+ACA data are primarily sensitive to the electron pressure difference across the shock front. To estimate the shock Mach number M, this difference can be combined with the value for the upstream electron pressure derived from an independent Chandra X-ray analysis. In the case of instantaneous electron-ion temperature equilibration, we find M = 2.08 (+0.12,−0.12), in ≈ 2.4σ tension with the independent constraint from Chandra, M(X)= 2.74 ± 0.25. The assumption of purely adiabatic electron temperature change across the shock leads to M = 2.53 (+0.33,−0.25), in better agreement with the X-ray estimate M(X) = 2.57 ± 0.23 derived for the same heating scenario. Conclusion: We have demonstrated that interferometric observations of the thermal SZ effect provide constraints on the properties of the shock in the Bullet Cluster that are highly complementary to X-ray observations. The combination of X-ray and SZ data yields a powerful probe of the shock properties, capable of measuring M and addressing the question of electron-ion equilibration in cluster shocks. Our analysis is however limited by systematics related to the overall cluster geometry and the complexity of the post-shock gas distribution. To overcome these limitations, a simultaneous, joint-likelihood analysis of SZ and X-ray data is needed.

Luca Di Mascolo↗
DOE OSTI2020

Multiplexed binary assembly and quantitative tracking of bacterial populations

Compositions and methods for tracking combinations of mutations in populations by coupling DNA assembly and paired-end sequencing technology with high-throughput sequencing. This allows for an unlimited number of mutation sites in a population of microorganisms to be identified, and allows sites to be sequenced across several replicates. Unique identifiers (DNA barcodes) can be used to measure fitness data by sequencing of the barcodes multiplexed with the mutation sites with high-throughput short sequencing read technology, allowing selections to be rapidly performed on populations of known combinatorial genotypes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

High Temperature, High Pressure Stability of Aqueous Foams for Potential Application in Enhanced Geothermal System (EGS)

Enhanced Geothermal Systems (EGS) are typically constructed by injecting high-pressure water into deep hot dry rocks (HDR) under carefully controlled conditions to create new or re-open existing fractures, which usually uses an immense quantity of water. In lieu, a more sustainable technique is to utilize aqueous foams as fracturing fluid to reduce water consumption and waste-water treatments of conventional hydraulic fracturing. Although foam-based fracturing has shown promising results in oil and gas industries, its feasibility is not demonstrated in EGS conditions that usually involve high temperature and high pressures. One potential barrier of utilizing foams in EGS applications is that foams are thermodynamically unstable and will become more unstable with increasing temperature due to liquid drainage, bubble coarsening, and coalescence. This work focuses on evaluation of the stability of selected aqueous foams under high-temperature and high-pressure conditions. Specifically, foams generated with surfactant alfa olefin sulfonate (AOS) were studied at temperature up to 200oC, while the pressure ranged between 100 psi and 1000 psi. The effect of additional stabilizing agents was also examined, including guar gum, bentonite clay, borate salt crosslinker, silicon dioxide nanoparticles (SiO2), and graphene oxide (GO) dispersion. Results showed that the stabilizing agents can enhance the foam thermal stability. Foams made with AOS and the borate salt crosslinker exhibited the longest half-life of 20 min at 200°C when the pressure was at 1000 psi. Data fitting showed that foam stability decreased exponentially as temperature increased. On the other hand, pressure increased the foam half-life which followed a power model. This study indicates that it may be possible to obtain highly stable foams at high-temperature and high-pressure conditions with appropriate stabilizing agents.

Thakore, Viren↗
DOE OSTI2021

Analysis of Lobe Power Calculator and Indication System with Physics and Cycle Based Models

The Advanced Test Reactor (ATR) at INL measures reactor power through two methods, thermal and Nitrogen-16 (N-16) activity. Water power calculator (WPC) is a thermal power system that measures flow and temperature to determine the thermal quadrant powers. The N-16 system utilizes a beta chamber detector that outputs reactivity levels to calculate lobe power through an algorithm called lobe power calculation and indication system (LPCIS). The LPCIS utilizes the N-16 system and the WPC system to determine reactor core power levels. The WPC provides accurate calculations of quadrant and total reactor thermal power. With the use of WPC measurements, thermal-to-N-16 (T2N) power ratios are produced to determine if the two indication systems agree on core power. Relative magnitude equations are used to utilize N-16 coefficients and multipliers to improve the indications of the LPCIS. This is crucial for maintaining safety limits. Currently, the LPCIS system uses linear equations and matrices to calculate lobe power through multipliers and coefficients. Advancements in technology and system upgrades have assisted system engineers at INL with the objective to reach a more dynamic system. In return the system demonstrates an increase in the accuracy of power reading while maintaining safety margins. The new proposed coefficient and multiple method implements a cycle specific, and a physics-based model intended for changing coefficients during operation. This calibration experiment focuses on power splitting, outer shim and neck shim, as well as fuel burning into the RDAS weighting factors. Results show that the Physics learning method yields a smaller error margin inside of the desired power range for the data set 166-A. This proved true for both constrained and unconstrained testing. This is most likely due to the physics data fitting approach that resulted in favorable coefficients and multipliers for 166-A. Continuing to improve the physics-based model will help improve the power accuracy of the LPCIS system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

NOvA Dual-Baseline Search for Active-to-Sterile Neutrino Oscillations using Neutrino- and Antineutrino-Enriched Samples

We report a search for neutrino oscillations to sterile neutrinos in the NOvA detectors under a model with three active and one sterile neutrinos. This search simultaneously fits data in the two NOvA detectors and is the first from NOvA to use both neutrino- and antineutrino-mode beams, with exposures of $26.61\times10^{20}$ and $12.50\times10^{20}$ protons on target, respectively. There is no evidence for sterile neutrinos in the data and we are able to exclude regions of parameter space that were allowed by previous experiments, including most of the allowed region reported by IceCube.

Abubakar, S. [Erciyes U.]↗
NASA NTRS1995

Multiparameter linear least-squares fitting to Poisson data one count at a time

A standard problem in gamma-ray astronomy data analysis is the decomposition of a set of observed counts, described by Poisson statistics, according to a given multicomponent linear model, with underlying physical count rates or fluxes which are to be estimated from the data. Despite its conceptual simplicity, the linear least-squares (LLSQ) method for solving this problem has generally been limited to situations in which the number n(sub i) of counts in each bin i is not too small, conventionally more than 5-30. It seems to be widely believed that the failure of the LLSQ method for small counts is due to the failure of the Poisson distribution to be even approximately normal for small numbers. The cause is more accurately the strong anticorrelation between the data and the wieghts w(sub i) in the weighted LLSQ method when square root of n(sub i) instead of square root of bar-n(sub i) is used to approximate the uncertainties, sigma(sub i), in the data, where bar-n(sub i) = E(n(sub i)), the expected value of N(sub i). We show in an appendix that, avoiding this approximation, the correct equations for the Poisson LLSQ (PLLSQ) problems are actually identical to those for the maximum likelihood estimate using the exact Poisson distribution. We apply the method to solve a problem in high-resolution gamma-ray spectroscopy for the JPL High-Resolution Gamma-Ray Spectrometer flown on HEAO 3. Systematic error in subtracting the strong, highly variable background encountered in the low-energy gamma-ray region can be significantly reduced by closely pairing source and background data in short segments. Significant results can be built up by weighted averaging of the net fluxes obtained from the subtraction of many individual source/background pairs. Extension of the approach to complex situations, with multiple cosmic sources and realistic background parameterizations, requires a means of efficiently fitting to data from single scans in the narrow (approximately = 1.2 keV, HEAO 3) energy channels of a Ge spectrometer, where the expected number of counts obtained per scan may be very low. Such an analysis system is discussed and compared to the method previously used.

Wheaton, Wm. A.↗
NASA NTRS2021

Dealing with Ion LET Uncertainties: An Application of Generalized Linear Models

Although most SEE rate estimation methods presume a fit to SEE cross section vs. LET, fitting SEE data is challenging because the data are not compatible with the assumptions of many common fitting techniques (e.g. linear regression. The difficulty of fitting such data is compounded when the LET of the ion responsible for an SEE is uncertain. We modify a Generalized Linear Model SEE data fitting methodology to accommodate uncertain LET and apply the method to the problem of backside heavy-ion SEE testing to demonstrate the utility of the method, explore the dependence of systematic errors that arise from improper treatment of LET uncertainty and develop guidelines for minimizing such systematic errors when proper treatment is not possible. Additional applications are suggested and assessed for suitability of treatment by the model.

Single-event effects↗
NASA NTRS1995

Fitting Prony Series To Data On Viscoelastic Materials

Improved method of fitting Prony series to data on viscoelastic materials involves use of least-squares optimization techniques. Based on optimization techniques yields closer correlation with data than traditional method. Involves no assumptions regarding the gamma'(sub i)s and higher-order terms, and provides for as many Prony terms as needed to represent higher-order subtleties in data. Curve-fitting problem treated as design-optimization problem and solved by use of partially-constrained-optimization techniques.

Hill, S. A.↗
NASA NTRS1980

Cosmic ray studies with a gas Cherenkov counter in association with an ionization spectrometer

The results from a balloon-borne gas Cherenkov counter (threshold 16.5 GeV/nucleon) and an ionization spectrometer are presented. The gas Cherenkov counter provides an absolute energy distribution for the response of the calorimeter for 5 or = Z 26 nuclei of cosmic rays. The contribution of scintillation to the gas Cherenkov pulse height was obtained by independently selecting particles below the gas Cherenkov threshold using the ionization spectrometer. Energy spectra were derived by minimizing the chi squared between Monte Carlo simulted data and flight data. Best fit power laws, dN/dE = AE-gamma, were determined for C, N, O, Ne, Mg, and Si. The power laws, all consistent with E (-2.7) are not good fits to the data. A better fit is obtained using the spectrum derived from the spectrometer. The data from the ionization calorimeter and the gas Cherenkov are thus completely self-consistent.

Balasubrahmanyan, V. K.↗
DOE OSTI2023

Rate Equations for Reversible Disproportionation Reactions and Fitting to Time-Course Data

Integrated rate equations are straightforward to fit to experimental data to verify a proposed mechanism and to extract kinetic parameters. Such equations are derived for reversible disproportionation/comproportionation reactions with any set of initial concentrations. Extraction of forward and reverse rate constants from experimental data by fitting the rate law to the data is demonstrated for the disproportionation of 2,2,6,6-tetramethyl-1-piperidinyl-N-oxyl (TEMPO) under acidic conditions where the approach to equilibrium is observed.

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

Beyond traditional diagnostics: Identifying active galactic nuclei using spectral energy distribution fitting in DESI data

Active galactic nuclei (AGN) are typically identified through their distinctive X-ray or radio emissions, mid-infrared (MIR) colors, or emission lines. However, each method captures different subsets of AGN due to signal-to-noise (S/N) limitations, redshift coverage, and extinction effects, underscoring the necessity for a multiwavelength approach for comprehensive AGN samples. This study explores the effectiveness of spectral energy distribution (SED) fitting as a robust method for AGN identification. Using CIGALE optical-MIR SED fits on DESI Early Data Release galaxies, we compare SED-based AGN selection (AGNFRAC ≥ 0.1) with traditional methods including BPT diagrams, WISE colors, X-ray, and radio diagnostics. The SED fitting identifies ∼70% of narrow- and broad-line AGN and 87% of WISE-selected AGN. Incorporating high S/N WISE photometry reduces star-forming galaxy contamination from 62% to 15%. Initially, ∼50% of SED-AGN candidates are undetected by standard methods, but additional diagnostics classify ∼85% of these sources, revealing low-ionization nuclear emission-line regions and retired galaxies potentially representing evolved systems with weak AGN activity. Further spectroscopic and multiwavelength analysis will be essential to determine the true AGN nature of these sources. SED fitting provides complementary AGN identification, unifying multiwavelength AGN selections. This approach enables more complete – albeit somewhat contaminated – AGN samples, which are essential for upcoming large-scale surveys where spectroscopic diagnostics may be limited.

Seyfert↗
NASA NTRS1982

Relativistic precessing jets in quasars and radio galaxies - Models to fit high resolution data

The formulation of generalized models tracing the geometry and intensity of the synchrotron emission from precessing, twin, relativistic jets as projected on the plane of the sky is presented. It is shown that neither the shape of the image nor its relative intensities are altered by including the effects of a cosmological redshift and a relative velocity between the source and observer. The models are fitted to the available data for several quasars and radio galaxies and demonstrate the plausibility of the phenomenon. Probable selection effects are considered and diagnostics given for recognizing objects showing this behavior. In the radio galaxies considered, velocities up to about 0.2c and precession periods of 1,000,000 yr are deduced. In the QSOs investigated, velocities of 0.7c and greater are found and periods of order 10,000 yr. In some cases precession cone angles increase with time. Consequences in terms of lifetimes of QSO behavior and binary supermassive objects are discussed.

Gower, A. C.↗
NASA NTRS1980

Cosmic ray studies with a gas Cerenkov counter in association with an ionization spectrometer

The results from a balloon-borne gas Cerenkov counter (threshold 16.5 GeV/nuc) and an ionization spectrometer are presented. The gas Cerenkov counter provides an absolute energy calibration for the response of the calorimeter for the Z range of 5-26 nuclei of cosmic rays. The contribution of scintillation to the gas Cerenkov pulse height has been obtained by independently selecting particles below the gas Cerenkov threshold using the ionization spectrometer. Energy spectra were derived by minimizing the chi-squared between a Monte Carlo simulated data and flight data. Best fit power laws were determined for C, N, O, Ne, Mg, and Si. The power laws, all consistent with E exp-2.7, are not good fits to the data. A better fit is obtained using the spectrum derived from the spectrometer.

Balasubrahmanyan, V. K.↗
NASA NTRS1997

A Numerical Method for Calculating Stellar Occultation Light Curves from an Arbitrary Atmospheric Model

We present a method for speeding up numerical calculations of a light curve for a stellar occultation by a planetary atmosphere with an arbitrary atmospheric model that has spherical symmetry. This improved speed makes least-squares fitting for model parameters practical. Our method takes as input several sets of values for the first two radial derivatives of the refractivity at different values of model parameters, and interpolates to obtain the light curve at intermediate values of one or more model parameters. It was developed for small occulting bodies such as Pluto and Triton, but is applicable to planets of all sizes. We also present the results of a series of tests showing that our method calculates light curves that are correct to an accuracy of 10(exp -4) of the unocculted stellar flux. The test benchmarks are (i) an atmosphere with a l/r dependence of temperature, which yields an analytic solution for the light curve, (ii) an atmosphere that produces an exponential refraction angle, and (iii) a small-planet isothermal model. With our method, least-squares fits to noiseless data also converge to values of parameters with fractional errors of no more than 10(exp -4), with the largest errors occurring in small planets. These errors are well below the precision of the best stellar occultation data available. Fits to noisy data had formal errors consistent with the level of synthetic noise added to the light curve. We conclude: (i) one should interpolate refractivity derivatives and then form light curves from the interpolated values, rather than interpolating the light curves themselves; (ii) for the most accuracy, one must specify the atmospheric model for radii many scale heights above half light; and (iii) for atmospheres with smoothly varying refractivity with altitude, light curves can be sampled as coarsely as two points per scale height.

Chamberlain, D. M.↗
NASA NTRS1974

A Monte Carlo investigation of experimental data requirements for fitting polynomial functions

This report examines the extent to which sample size affects the accuracy of a low order polynomial approximation of an experimentally observed quantity and establishes a trend toward improvement in the accuracy of the approximation as a function of sample size. The task is made possible through a simulated analysis carried out by the Monte Carlo method, in which data are generated by using several transcendental or algebraic functions as models. Contaminated data of varying amounts are fitted to linear quadratic or cubic polynomials, and the behavior of the mean-squared error of the residual variance is determined as a function of sample size. Results indicate that the effect of the size of the sample is significant only for relatively small sample sizes and diminishes drastically for moderate and large amounts of experimental data.

Canavos, G. C.↗
NASA NTRS1995

Metafitting: Weight optimization for least-squares fitting of PTTI data

For precise time intercomparisons between a master frequency standard and a slave time scale, we have found it useful to quantitatively compare different fitting strategies by examining the standard uncertainty in time or average frequency. It is particularly useful when designing procedures which use intermittent intercomparisons, with some parameterized fit used to interpolate or extrapolate from the calibrating intercomparisons. We use the term 'metafitting' for the choices that are made before a fitting procedure is operationally adopted. We present methods for calculating the standard uncertainty for general, weighted least-squares fits and a method for optimizing these weights for a general noise model suitable for many PTTI applications. We present the results of the metafitting of procedures for the use of a regular schedule of (hypothetical) high-accuracy frequency calibration of a maser time scale. We have identified a cumulative series of improvements that give a significant reduction of the expected standard uncertainty, compared to the simplest procedure of resetting the maser synthesizer after each calibration. The metafitting improvements presented include the optimum choice of weights for the calibration runs, optimized over a period of a week or 10 days.

Douglas, Rob J.↗