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

Estimation of time- and state-dependent delays and other parameters in functional differential equations

A parameter estimation algorithm is developed which can be used to estimate unknown time- or state-dependent delays and other parameters (e.g., initial condition) appearing within a nonlinear nonautonomous functional differential equation. The original infinite dimensional differential equation is approximated using linear splines, which are allowed to move with the variable delay. The variable delays are approximated using linear splines as well. The approximation scheme produces a system of ordinary differential equations with nice computational properties. The unknown parameters are estimated within the approximating systems by minimizing a least-squares fit-to-data criterion. Convergence theorems are proved for time-dependent delays and state-dependent delays within two classes, which say essentially that fitting the data by using approximations will, in the limit, provide a fit to the data using the original system. Numerical test examples are presented which illustrate the method for all types of delay.

Murphy, K. A.↗
NASA NTRS1992

Magellan ephemeris improvement using synthetic aperture radar landmark measurements

A technique is described for measuring the positions of landmarks in multiple SAR images of the surface of Venus taken aboard the Magellan spacecraft. These measurements are then used to improve the spacecraft orbit estimate. The Venus-fixed coordinates of the landmarks are also estimated, as are the low-order coefficients of the gravitational field. Sample results are shown for five-orbit and 13-orbit data arcs using hundreds of landmark measurements. Reasonably good fits to the data are obtained for the short-arc solutions, while the data fits over long arcs are poorer, possibly due to higher-order uncertainties in the gravitational field. A comparison of post-fit orbit uncertainties shows that the SAR data significantly improves the orbit estimate.

Chodas, Paul W.↗
NASA NTRS1995

Determination of plasma temperatures and luminosities using multiple extreme-ultraviolet and X-ray filters

We carefully examine the techniques used to infer temperatures of stellar coronal plasmas from the count rates of several broadband instruments in the X-ray and extreme-ultraviolet spectral ranges. In particular, we determine to what extent temperatures can be constrained and the corresponding uncertainties in the luminosities and emission measures lowered by fitting simultaneously count rates from the Einstein imaging proportional counter (IPC), the ROSAT Position Sensitive Proportional Counter (PSPC), the ROSAT Wide Field Camera (WFC) (both filters), and the EXOSAT Low Energy Telescope (LET) with the 3-Lex filter. We use published plasma emissivities with solar photospheric abundances. Since it has been found that single-temperature plasmas do not fit IPC data well, we assume a two-temperature plasma model. We find that, even with count rates from all of the above filters and overly optimistic error estimates, it is still not possible to determine a unique two-temperature solution. However, since the use of count rates from many filters can reduce substantially the number of possible solutions, temperature solutions determined by other means can be tested. We carry out such an analysis on a set of 18 nearby late-type stars to determine possible two-temperature solutions using multifilter photometry, and we compare these results with the temperature solutions derived by Schmitt et al. (1990) using IPC spectral data. In general, the two-temperature fits derived from the IPC spectral data are inconsistent with our results, with our data implying that, for many stars, the two temperatures derived by the IPC may be too low by about a factor of 2. The EXOSAT transmission grating Spectrometer (TGS) spectra of capella and sigma(exp 2) CrB support this conclusion. For Procyon and 70 Oph, though, the presence of a temperature component cooler than a million degress (not detected by the IPC) is deduced. While our analysis suggests the existence of more than one temperature in the coronae of late-type stars, in many instances our WFC data appear to be inconsistent with the presence of significant emission measure over a broad temperature distribution. This, together with the success of two-temperature plasmas in fitting IPC and TGS data, implies that for many stars, the coronal emission measure distribution may in fact be dominated by two distinct temperature regimes.

Wood, Brian E.↗
NASA NTRS2025

An Efficient GPU-Accelerated Multi-Source Global Fit Pipeline for LISA Data Analysis

The large-scale analysis task of deciphering gravitational wave signals in the LISA data stream will be difficult, requiring a large amount of computational resources and extensive development of computational methods. Its high dimensionality, multiple model types, and complicated noise profile require a global fit to all parameters and input models simultaneously. In this work, we detail our global fit algorithm, called “Erebor,” designed to accomplish this challenging task. It is capable of analysing current state-of-the-art datasets and then growing into the future as more pieces of the pipeline are completed and added. We describe our pipeline strategy, the algorithmic setup, and the results from our analysis of the LDC2A Sangria dataset, which contains Massive Black Hole Binaries, compact Galactic Binaries, and a parameterized noise spectrum whose parameters are unknown to the user. The Erebor algorithm includes three unique and very useful contributions: GPU acceleration for enhanced computational efficiency; ensemble MCMC sampling with multiple MCMC walkers per temperature for better mixing and parallelized sample creation; and special online updates to reversible-jump (or trans-dimensional) sampling distributions to ensure sampler mixing and accurate initial estimates for detectable sources in the data. We recover posterior distributions for all 15 (6) of the injected MBHBs in the LDC2A training (hidden) dataset. We catalog ∼12000 Galactic Binaries (∼8000 as high confidence detections) for both the training and hidden datasets. All of the sources and their posterior distributions are provided in publicly available catalogs.

LISA global fit↗
DOE OSTI2022

Apodization Specific Fitting for Improved Resolution, Charge Measurement, and Data Analysis Speed in Charge Detection Mass Spectrometry

Short-time Fourier transforms with short segment lengths are typically used to analyze single ion charge detection mass spectrometry (CDMS) data either to overcome effects of frequency shifts that may occur during the trapping period or to more precisely determine the time at which an ion changes mass or charge, or enters an unstable orbit. The short segment lengths can lead to scalloping loss unless a large number of zero-fills are used, making computational time a significant factor in real-time analysis of data. Apodization specific fitting leads to a 9-fold reduction in computation time compared to zero-filling to a similar extent of accuracy. This makes possible real-time data analysis using a standard desktop computer. Rectangular apodization leads to higher resolution than the more commonly used Gaussian or Hann apodization and makes it possible to separate ions with similar frequencies, a significant advantage for experiments in which the masses of many individual ions are measured simultaneously. Equally important is a >20% increase in S/N obtained with rectangular apodization compared to Gaussian or Hann, which directly translates to a corresponding improvement in accuracy of both charge measurements and ion energy measurements that rely on the amplitudes of the fundamental and harmonic frequencies. Finally, combined with computing the fast Fourier transform in a lower-level language, this fitting procedure eliminates computational barriers and should enable real-time processing of CDMS data on a laptop computer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

Spectroscopic analysis of global tide gauge sea level data

Yearly and monthly global tide-gage sea-level data are fitted to numerically generated tidal data in order to search for the 18.6-yr lunar nodal tide and 14-month pole tide. Both of these tides are clearly evident, with amplitudes and phases that are consistent with a global equilibrium response. The ocean's response to atmospheric pressure is studied with the least-squares fit technique. Consideration is given to the global rise in sea level, the effects of postglacial rebound, and the possible causes of the enhanced pole tides in the North Sea, the Baltic Sea, and the Gulf of Bothnia. The results support O'Connor's (1986) suggestion that the enhanced pole tide in these regions is due to meteorological forcing rather than a basin-scale resonance. Also, the global average of the tide-gage data show an increase in sea level over tha last 80 yr of between 1.1 and 1.9 mm/yr.

Trupin, A.↗
NASA NTRS2017

Raising the Albedo of 2010 GY6: Fitting ATPM to Wise Data

Near-Earth Asteroid 462775 (2010 GY6) is in the Apollo orbit-family with a 1.46 year orbital period. 2010 GY6 was measured by WISE and fitted with NEATM, yielding NEATM model parameters of D=1.1 km, pv=0.028 and eta=2.3.The NEATM-derived geometric albedo of 2010 GY6 is lower than the surface of comet 67P/C-G. The eta value is considerably higher than typical for its phase angle of 33 deg, indicating a cooler surface due to non-zero thermal inertia and/or surface roughness are important. If the thermal inertia and surface roughness are constrained by fitting the Advanced Thermophysical Model (ATPM) to the WISE data, what would the resulting geometric albedo? We find pv=0.06-0.08, in the same range as B- or C-type NEAs like Bennu or JU3.

Near-Earth Asteroid 462775↗
DOE OSTI2023

Thermal expansion of 4H and 6H SiC from 5 K to 340 K

The first thermal expansion measurements of the 4H and 6H polytypes of SiC below room temperature are reported. The measurements were carried out on single-crystal specimens using high-resolution capacitive-based dilatometry. For both polytypes, the thermal expansion coefficient is below 2.4 × 10 -6 1/K near room temperature. No phase transitions are observed over the 5 K to 340 K temperature range of the measurements. The thermal expansion coefficient α of 4H SiC is slightly anisotropic for measurements parallel and perpendicular to the crystallographic c axis with α ∥ about 2.2 × 10 -7 1/K larger than α ⟂ near room temperature. For 6H SiC no discernible anisotropy is observed. The differences in anisotropy can be understood by considering the ratio of hexagonal to cubic bonds of each polytype. Narrow regions with negative thermal expansion that are within the limits of our resolution (~ 1 x 10 -8 ) are observed in the vicinity of 30 K for both specimens. In conclusion, tabulated data, polynomial fits, fit parameters, and comparison to data based on lattice-parameter measurements are provided.

36 MATERIALS SCIENCE↗
NASA NTRS2006

Using Quasiparticle Poisoning To Detect Photons

According to a proposal, a phenomenon associated with excitation of quasiparticles in certain superconducting quantum devices would be exploited as a means of detecting photons with exquisite sensitivity. The phenomenon could also be exploited to perform medium-resolution spectroscopy. The proposal was inspired by the observation that Coulomb blockade devices upon which some quantum logic gates are based are extremely sensitive to quasiparticles excited above the superconducting gaps in their leads. The presence of quasiparticles in the leads can be easily detected via the charge states. If quasiparticles could be generated in the leads by absorption of photons, then the devices could be used as very sensitive detectors of electromagnetic radiation over the spectral range from x-rays to submillimeter waves. The devices in question are single-Cooper-pair boxes (SCBs), which are mesoscopic superconducting devices developed for quantum computing. An SCB consists of a small superconducting island connected to a reservoir via a small tunnel junction and connected to a voltage source through a gate capacitor. An SCB is an artificial two-level quantum system, the Hamiltonian of which can be controlled by the gate voltage. One measures the expected value of the charge of the eigenvectors of this quantum system by use of a radio-frequency single-electron transistor. A plot of this expected value of charge as a function of gate voltage resembles a staircase that, in the ideal case, consists of steps of height 2 e (where e is the charge of one electron). Experiments have shown that depending on the parameters of the device, quasiparticles in the form of "broken" Cooper pairs present in the reservoir can tunnel to the island, giving rise to steps of 1 e. This effect is sometimes called "poisoning." Simulations have shown that an extremely small average number of quasiparticles can generate a 1-e periodic signal. In a device according to the proposal, this poisoning would be turned to advantage. Depending on the wavelength, an antenna or other component would be used to couple radiation into the reservoir, wherein the absorption of photons would break Cooper pairs, thereby creating quasiparticles that, in turn, would tunnel to the island, creating a 1-e signal. On the basis of conservative estimates of device parameters derived from experimental data and computational simulations that fit the data, it has been estimated that the noise equivalent power of a device according to the proposal could be as low as 6 10(exp -22) W/Hz(exp 1/2). It has also been estimated that the spectroscopic resolution (photon energy divided by increment of photon energy) of such a device in visible light would exceed 100.

Echternach, Pierre↗
NASA NTRS2005

Spectral Photometric Properties of the Moon

We modeled the solar phase curves of the moon at a series of wavelengths using the full disk telescopic observations [1]. We endeavored to keep the database self-contained, that is, to use the values derived for the solar magnitude and phase curves of the disk-integrated [1]. These observations were made in a suite of 10 narrowband filters between 0.315 microns and 1.06 microns, and in the broad band Johnson UBV filters, as part of a larger program to obtain photoelectric photometry of the larger planets. Two aspects of the lunar observations are unique. First, the observations cover phase angles from 6deg through 120deg. More importantly, the observers used a special 20-mm diameter f/15 fused quartz lens constructed solely for this purpose. The lens reduced the whole lunar image in the focal plane to a size comparable to the planets observed as part of the same program. This image was fed directly into the photometer. Thus, these observations constitute the only existing set of phase curves of the entire lunar disk over a range of wavelengths. Table 1 lists the values of the Hapke model parameters which fit the data. Figure 1 is an example of the model fits to the data.

Dominque, D.↗
NASA NTRS2025

Efficient GPU-Accelerated MultiSource Global Fit Pipeline for LISA Data Analysis

The large-scale analysis task of deciphering gravitational-wave signals in the LISA data stream will be difficult, requiring a large amount of computational resources and extensive development of computational methods. Its high dimensionality, multiple model types, and complicated noise profile require a global fit to all parameters and input models simultaneously. In this work, we detail our global fit algorithm, called “Erebor,” designed to accomplish this challenging task. It is capable of analyzing current state-of-the-art datasets and then growing into the future as more pieces of the pipeline are completed and added. We describe our pipeline strategy, the algorithmic setup, and the results from our analysis of the LDC2A Sangria dataset, which contains massive black hole binaries, compact galactic binaries, and a parametrized noise spectrum whose parameters are unknown to the user. The Erebor algorithm includes three unique and very useful contributions: GPU acceleration for enhanced computational efficiency; ensemble Markov Chain Monte Carlo (MCMC) sampling with multiple MCMC walkers per temperature for better mixing and parallelized sample creation; and special online updates to reversible-jump (or transdimensional) sampling distributions to ensure sampler mixing and accurate initial estimates for detectable sources in the data.We recover posterior distributions for all 15 (6) of the injected massive black hole binaries (MBHB) in the LDC2A training (hidden) dataset. We catalog ∼12000 galactic binaries (∼8000 as high confidence detections) for both the training and hidden datasets. All of the sources and their posterior distributions are provided in publicly available catalogs.

LISA↗
NASA NTRS1994

Detection of the SO2 atmosphere on Io with the Hubble Space Telescope

Observations of the trailing hemisphere of Io made with the Faint Object Spectrograph of the Hubble Space Telescope (HST) in March 1992 have resulted in the first detection of atmospheric SO2 absorption bands in the ultraviolet. These observations represent only the third positive means of detection of what is widely believed to be Io's primary atmospheric constituent. Below approximately 2130 A the geometric albedo of the satellite is dominated by SO2 gas absorption band signatures, which have been analyzed using models that include the effects of optical thickness, temperature, and spatial distribution. The disk-intergrated HST data cannot resolve the spatial distribution, but it is possible to define basic properties and set constraints on the atmosphere at the time of the observations. Hemispheric atmospheres with average column density N = 6 - 10 x 10(exp 15)/sq sm and T(gas) = 110 - 500 K fit the data, with preference for temperatures of approximately 200 - 250 K. Better fits are found as the atmosphere is spatially confined, with a limit of approximately 8% hemispheric areal coverage and N approximately equal to 3 x 10(exp 17)/sq cm with colder 110 - 250 K temepratures. A dense (N greater than or equal to 10(exp 16)/sq cm), localized component of SO2 gas, such as that possibly associated with active volcanoes, can generate the observed spectral constrast only when the atmosphere is cold (110 K) and an extended component such as Pele is included. The combination of a dense, localized atmosphere with a tenuous component (N less than 10(exp 16)/sq cm, either patchy or extended) also fits the data. In all cases the best fit models imply a disk-averaged column density larger than exospheric but approximately 10 - 30 times less than the previous upper limit from near-UV observations.

Ballester, G. E.↗
DOE OSTI2026

Measurements of three-flavor neutrino oscillations from a PISCES two-detector fit to the NOvA Experiment data

NOvA is a long-baseline neutrino oscillation experiment with two functionally identical detectors: a Near Detector (ND) at Fermilab, placed 1 km from the neutrino source, and a Far Detector (FD) located 810 km away from the ND in Minnesota. NOvA s primary physics goals are to measure the neutrino oscillation parameters $\theta_{23}$ and $\Delta m^2_{32}$ with high precision, determine the neutrino mass hierarchy, and constrain the value of $\delta_{CP}$, primarily via the study of muon neutrino to electron neutrino oscillation. Extracting values for oscillation parameters from fits to data usually relies on treating systematic uncertainties as nuisance parameters, a strategy that suffers from poor scalability as the number of uncertainties becomes larger. This work introduces PISCES (Parameter Inference with Systematic Covariance and Exact Statistics), a novel method that circumvents this scalability problem by encoding systematic uncertainties into a covariance matrix. PISCES utilizes a nested minimization in which optimal systematic pulls are first computed using the covariance matrix in an inner minimization step, then the oscillation parameters are profiled over in the outer minimization. PISCES also uses a Poisson Likelihood term, making it ideal for the inclusion of low-statistic samples in the fits. PISCES is a flexible framework that also supports complex fits, such as a joint Near and Far detector fit. In the standard NOvA analysis, oscillation parameters are extracted using an extrapolation technique in which the ND data indirectly constrain the FD prediction via a ratio method. PISCES, on the other hand, enables a simultaneous ND+FD fit, allowing the high-statistics ND data to directly constrain systematic uncertainties across all samples. This thesis presents the full PISCES joint ND+FD fit for the NOvA three-flavor analysis, details its implementation, and evaluates its performance through extensive robustness tests and fake data studies. It also provides a comparison between the PISCES joint ND+FD results and the standard NOvA extrapolation method using the full NOvA 10-year data set. The results demonstrate that PISCES can successfully fit NOvA data while incorporating the constraints from the ND detectors consistently, using physically motivated systematic uncertainties to account for data/MC discrepancies.

Rajaoalisoa, Miriama [Cincinnati U.]↗
DOE OSTI2022

Consistent lensing and clustering in a low- S 8 Universe with BOSS, DES Year 3, HSC Year 1, and KiDS-1000

ABSTRACT We evaluate the consistency between lensing and clustering based on measurements from Baryon Oscillation Spectroscopic Survey combined with galaxy–galaxy lensing from Dark Energy Survey (DES) Year 3, Hyper Suprime-Cam Subaru Strategic Program (HSC) Year 1, and Kilo-Degree Survey (KiDS)-1000. We find good agreement between these lensing data sets. We model the observations using the Dark Emulator and fit the data at two fixed cosmologies: Planck (S8 = 0.83), and a Lensing cosmology (S8 = 0.76). For a joint analysis limited to large scales, we find that both cosmologies provide an acceptable fit to the data. Full utilization of the higher signal-to-noise small-scale measurements is hindered by uncertainty in the impact of baryon feedback and assembly bias, which we account for with a reasoned theoretical error budget. We incorporate a systematic inconsistency parameter for each redshift bin, A, that decouples the lensing and clustering. With a wide range of scales, we find different results for the consistency between the two cosmologies. Limiting the analysis to the bins for which the impact of the lens sample selection is expected to be minimal, for the Lensing cosmology, the measurements are consistent with A = 1; A = 0.91 ± 0.04 (A = 0.97 ± 0.06) using DES+KiDS (HSC). For the Planck case, we find a discrepancy: A = 0.79 ± 0.03 (A = 0.84 ± 0.05) using DES+KiDS (HSC). We demonstrate that a kinematic Sunyaev–Zeldovich-based estimate for baryonic effects alleviates some of the discrepancy in the Planck cosmology. This analysis demonstrates the statistical power of small-scale measurements; however, caution is still warranted given modelling uncertainties and foreground sample selection effects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Mirror twin Higgs cosmology: constraints and a possible resolution to the H$_{0}$ and S$_{8}$ tensions

The mirror twin Higgs model (MTH) is a solution to the Higgs hierarchy problem that provides well-predicted cosmological signatures with only three extra parameters: the temperature of the twin sector, the abundance of twin baryons, and the vacuum expectation value (VEV) of twin electroweak symmetry breaking. These parameters specify the behavior of twin radiation and the acoustic oscillations of twin baryons, which lead to testable effects on the cosmic microwave background (CMB) and large-scale structure (LSS). While collider searches can only probe the twin VEV, through a fit to cosmological data we show that the existing CMB (Planck18 TTTEEE+lowE+lowT+lensing) and LSS (KV450) data already provide useful constraints on the remaining MTH parameters. Additionally, we show that the presence of twin radiation in this model can raise the Hubble constant H$_{0}$ while the scattering twin baryons can reduce the matter fluctuations S$_{8}$, which helps to relax the observed H$_{0}$ and S$_{8}$ tensions simultaneously. This scenario is different from the typical ΛCDM + ΔN$_{eff}$ model, in which extra radiation helps with the Hubble tension but worsens the S$_{8}$ tension. For instance, when including the SH0ES and 2013 Planck SZ data in the fit, we find that a universe with ≳ 20% of the dark matter comprised of twin baryons is preferred over ΛCDM by ~ 4σ. If the twin sector is indeed responsible for resolving the H$_{0}$ and S$_{8}$ tensions, future measurements from the Euclid satellite and CMB Stage 4 experiment will further measure the twin parameters to O(1 - 10%)-level precision. Our study demonstrates how models with hidden naturalness can potentially be probed using precision cosmological data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1985

Improved mapping of radio sources from VLBI data by least-square fit

A method is described for producing improved mapping of radio sources from Very Long Base Interferometry (VLBI) data. The method described is more direct than existing Fourier methods, is often more accurate, and runs at least as fast. The visibility data is modeled here, as in existing methods, as a function of the unknown brightness distribution and the unknown antenna gains and phases. These unknowns are chosen so that the resulting function values are as near as possible to the observed values. If researchers use the radio mapping source deviation to measure the closeness of this fit to the observed values, they are led to the problem of minimizing a certain function of all the unknown parameters. This minimization problem cannot be solved directly, but it can be attacked by iterative methods which we show converge automatically to the minimum with no user intervention. The resulting brightness distribution will furnish the best fit to the data among all brightness distributions of given resolution.

Rodemich, E. R.↗
DOE OSTI2024

Generating synthetic signaling networks for in silico modeling studies

Predictive models of signaling pathways have proven to be difficult to develop. Reasons include the uncertainty in the number of species, the complexity in species’ interactions, and the sparseness and uncertainty in experimental data. Traditional approaches to developing mechanistic models rely on collecting experimental data and fitting a single model to that data. This approach works for simple systems but has proven unreliable for complex systems such as biological signaling networks. For example, uncertainty and sparseness of the data often result in overfitted models that have little predictive value beyond recapitulating the experimental data itself. Thus, there is a need to develop new approaches to create predictive mechanistic models of complex systems. However, to determine the effectiveness of any new algorithm, a baseline model is needed to test its performance. To meet this need, we developed a method for generating artificial synthetic networks that are reasonably realistic and thus can be treated as ground truth models. These synthetic models can then be used to generate synthetic data for developing and testing algorithms designed to recover the underlying network topology and associated parameters. Here, we describe a simple approach for generating synthetic signaling networks that can be used for this purpose.

42 ENGINEERING↗