Search NASA⌕ Search

SEARCH · Search NASA

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 109 records · Page 6

Paired autoencoders for likelihood-free estimation in inverse problems

Abstract We consider the solution of nonlinear inverse problems where the forward problem is a discretization of a partial differential equation. Such problems are notoriously difficult to solve in practice and require minimizing a combination of a data-fit term and a regularization term. The main computational bottleneck of typical algorithms is the direct estimation of the data misfit. Therefore, likelihood-free approaches have become appealing alternatives. Nonetheless, difficulties in generalization and limitations in accuracy have hindered their broader utility and applicability. In this work, we use a paired autoencoder framework as a likelihood-free estimator (LFE) for inverse problems. We show that the use of such an architecture allows us to construct a solution efficiently and to overcome some known open problems when using LFEs. In particular, our framework can assess the quality of the solution and improve on it if needed. We demonstrate the viability of our approach using examples from full waveform inversion and inverse electromagnetic imaging.

Chung, Matthias (ORCID:0000000178224539)↗

Acclimation of leaf respiration consistent with optimal photosynthetic capacity

Abstract Plant respiration is an important contributor to the proposed positive global carbon‐cycle feedback to climate change. However, as a major component, leaf mitochondrial (‘dark’) respiration ( R d ) differs among species adapted to contrasting environments and is known to acclimate to sustained changes in temperature. No accepted theory explains these phenomena or predicts its magnitude. Here we propose that the acclimation of R d follows an optimal behaviour related to the need to maintain long‐term average photosynthetic capacity ( V cmax ) so that available environmental resources can be most efficiently used for photosynthesis. To test this hypothesis, we extend photosynthetic co‐ordination theory to predict the acclimation of R d to growth temperature via a link to V cmax , and compare predictions to a global set of measurements from 112 sites spanning all terrestrial biomes. This extended co‐ordination theory predicts that field‐measured R d and V cmax accessed at growth temperature ( R d,tg and V cmax,tg ) should increase by 3.7% and 5.5% per degree increase in growth temperature. These acclimated responses to growth temperature are less steep than the corresponding instantaneous responses, which increase 8.1% and 9.9% per degree of measurement temperature for R d and V cmax respectively. Data‐fitted responses proof indistinguishable from the values predicted by our theory, and smaller than the instantaneous responses. Theory and data are also shown to agree that the basal rates of both R d and V cmax assessed at 25°C ( R d,25 and V cmax,25 ) decline by ~4.4% per degree increase in growth temperature. These results provide a parsimonious general theory for R d acclimation to temperature that is simpler—and potentially more reliable—than the plant functional type‐based leaf respiration schemes currently employed in most ecosystem and land‐surface models.

Wang, Han↗

Improved Parametric Models for Explosion Pressure Signals Derived From Large Datasets

Accurate recording and characterization of explosion-induced pressure signals are key components of the forensic analysis of explosion events in the atmosphere. Parametric overpressure models based on several key waveform features (peak overpressure, positive pulse duration, and impulse) are widely used to estimate explosion energy in terms of trinitrotoluene equivalent yield. However, those models are often developed by a limited dataset, including only a few events or recordings at relatively short propagation distances. Here, we develop empirical waveform-parameter models based on a regression analysis of a large set of data curated from four chemical explosion experiments including 16 detonations. We measured peak overpressure and impulse for positive and negative phases from 1000 pressure signals recorded at local ranges (<20 km) with scaled distance up to 8000 m/kg 1/3 . Additionally, the measured waveform parameters showed large variation with respect to observing distances indicating the effects of atmospheric propagation. In this study, a second-order polynomial model was used in a least-squares regression to account for those propagation effects and to improve data fitting. In addition to model parameters for waveform features, we also determined range-dependent model uncertainties based on data variance. The model uncertainty represents the prediction error of our models and can be critical to evaluating the uncertainty of yield estimate.

58 GEOSCIENCES↗

Real-time data reduction at 100 Tbps: Challenge and opportunity for AI-based data reduction for next-generation large-scale nuclear physics collider experiment

The modern large-scale nuclear physics (NP) experiments in high-energy particle colliders utilize streaming-readout electronics to digitize detector response at O(100) Tbps bandwidth. Prominent examples at Brookhaven National Lab (BNL) include the sPHENIX experiment at Relativistic Heavy Ion Collider (RHIC), which is under construction, and the experiments proposed for the Electron-Ion Collider (EIC), planned for the 2030s . One of the main challenges for these streaming readout systems is to manage the data rate with sufficient data reduction in real time so the end-data fit persistent storage for offline analysis, which is typically at O(1000) times smaller and O(100) Gbps. Such data reduction traditionally is achieved via real-time high level triggers, which select and save a small subset of collisions of interest. Although triggering is applicable to high energy collider experiments such as those at the Large Hardron Collider at CERN, it is insufficient for these nuclear physics experiments which study diverse collision topologies. And traditional triggering approach is inefficient to preserve the max information harvested from the operation of colliders that costs O(100)M per year to DOE. Meanwhile, in recent years, ML-based high-throughput data reduction has emerged as a promising approach to efficiently preserve max information for a given space of persistent storage, e.g. via AI data compression, feature extraction, and noise filtering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Selecting step sizes in sensitivity analysis by finite differences

This paper deals with methods for obtaining near-optimum step sizes for finite difference approximations to first derivatives with particular application to sensitivity analysis. A technique denoted the finite difference (FD) algorithm, previously described in the literature and applicable to one derivative at a time, is extended to the calculation of several simultaneously. Both the original and extended FD algorithms are applied to sensitivity analysis for a data-fitting problem in which derivatives of the coefficients of an interpolation polynomial are calculated with respect to uncertainties in the data. The methods are also applied to sensitivity analysis of the structural response of a finite-element-modeled swept wing. In a previous study, this sensitivity analysis of the swept wing required a time-consuming trial-and-error effort to obtain a suitable step size, but it proved to be a routine application for the extended FD algorithm herein.

Iott, J.↗

The ISEE-1 ULEWAT flux tape description

The ISEE-1 ULEWAT FLUX tapes of solar proton observations summarized over 1hr intervals are described. Three years of data fit onto one 1600 BPI 9-track magnetic tape. In making the tapes, corrections were made to the ULEWAT basic data tapes in order to, e.g., remove rate spikes and account for changes in instrument response, so that instrument fluxes can be calculated easily from the FLUX tapes without further consideration of instrument performance. Limitations and warnings about the use of the tapes are included.

Mason, G. M.↗

Buckling behavior of Rene 41 tubular panels for a hypersonic aircraft wing

The buckling characteristics of Rene 41 tubular panels for a hypersonic aircraft wing were investigated. The panels were repeatedly tested for buckling characteristics using a hypersonic wing test structure and a universal tension/compression testing machine. The nondestructive buckling tests were carried out under different combined load conditions and in different temperature environments. The force/stiffness technique was used to determine the buckling loads of the panel. In spite of some data scattering, resulting from large extrapolations of the data fitting curve (because of the termination of applied loads at relatively low percentages of the buckling loads), the overall test data correlate fairly well with theoretically predicted buckling interaction curves. Also, the structural efficiency of the tubular panels was found to be slightly higher than that of beaded panels.

Ko, W. L.↗

The ISEE-3 ULEWAT: Flux tape description and heavy ion fluxes 1978-1984

The ISEE ULEWAT FLUX tapes contain ULEWAT and ISEE pool tape data summarized over relatively long time intervals (1hr) in order to compact the data set into an easily usable size. (Roughly 3 years of data fit onto one 1600 BPI 9-track magnetic tape). In making the tapes, corrections were made to the ULEWAT basic data tapes in order to, remove rate spikes and account for changes in instrument response so that to a large extent instrument fluxes can be calculated easily from the FLUX tapes without further consideration of instrument performance.

Mason, G. M.↗

Buckling behavior of Rene 41 tubular panels for a hypersonic aircraft wing

The buckling characteristics of Rene 41 tubular panels for a hypersonic aircraft wing were investigated. The panels were repeatedly tested for buckling characteristics using a hypersonic wing test structure and a universal tension/compression testing machine. The nondestructive buckling tests were carried out under different combined load conditions and in different temperature environments. The force/stiffness technique was used to determine the buckling loads of the panels. In spite of some data scattering resulting from large extrapolations of the data-fitting curve (because of the termination of applied loads at relatively low percentages of the buckling loads), the overall test data correlate fairly well with theoretically predicted buckling interaction curves. Also, the structural efficiency of the tubular panels was found to be slightly higher than that of beaded panels.

Ko, W. L.↗

SANS quantification of bound water in water-soluble polymers across multiple concentration regimes

Contrast-variation small-angle neutron scattering (CV-SANS) is a widely used technique for quantifying hydration water in soft matter systems, but it is predominantly applied in the dilute regime or for systems with a well-defined structure factor. Here, CV-SANS was used to quantify the number of hydration water molecules associating with three water-soluble polymers with different critical solution temperatures and types of water–solute interactions in dilute, semidilute, and concentrated solution through the exploration of novel methods of data fitting and analysis. Multiple SANS fitting workflows with varying levels of model assumptions were evaluated and compared to give insight into SANS model selection. These fitting pathways ranged from general, model-free algorithms to more standard form and structure factor fitting. In addition, Monte Carlo bootstrapping was evaluated as a method to estimate parameter uncertainty through simulation of technical replicates. The most robust fitting workflow for dilute solutions was found to be form factor fitting without CV-SANS ( i.e. polymer in 100% D 2 O). For semidilute and concentrated solutions, while the model-free approach can be mathematically defined for CV-SANS data, the addition of a structure factor imposes physical constraints on the optimization problem, suggesting that the optimal fitting pathway should include appropriate form and structure factor models. The measured hydration numbers were consistent with the number of tightly bound water molecules associated with each monomer unit, and the concentration dependence of the hydration number was largely governed by the chemistry-specific interactions between water and polymer. Polymers with weaker water–polymer interactions ( i.e. those with fewer hydration water molecules) were found to have more bound water at higher concentrations than those with stronger water–polymer interactions due to the increase in the number of forced water–polymer contacts in the concentrated system. This SANS-based method to count hydration water molecules can be applied to polymers in any concentration regime, which will lead to improved understanding of water–polymer interactions and their impact on materials design.

36 MATERIALS SCIENCE↗

Strong Upper-Plate Heterogeneity at the Hikurangi Subduction Margin (North Island, New Zealand) Imaged by Adjoint Tomography

We use earthquake-based adjoint tomography to invert for three-dimensional structure of the North Island, New Zealand, and the adjacent Hikurangi subduction zone. The study area, having a shallow depth to the plate interface below the North Island, offers a rare opportunity for imaging material properties at an active subduction zone using land-based measurements. Starting from an initial model derived using ray tomography, we perform iterative model updates using spectral element and adjoint simulations to fit waveforms with periods ranging from 4–30 s. We perform 28 model updates using an L-BFGS optimization algorithm, improving data fit and introducing P- and S-wave velocity changes of up to ±30%. Resolution analysis using point spread functions show that our measurements are most sensitive to heterogeneities in the upper 30 km. The most striking velocity changes coincide with areas related to the active Hikurangi subduction zone. Lateral velocity structures in the upper 5 km correlate well with New Zealand geology. Here, the inversion reveals increased along-strike heterogeneity on the margin. In Cook Strait we observe a low-velocity zone interpreted as deep sedimentary basins. In the central North Island, low-velocity anomalies are linked to surface geology, and we relate velocity structures at depth to crustal magmatic activity below the Taupō Volcanic Zone.

58 GEOSCIENCES↗

A detailed analysis of the high-resolution X-ray spectra of NGC 3516: variability of the ionized absorbers

The 1.5 Seyfert galaxy NGC 3516 presents a strong time variability in X-rays. We re-analyzed the nine observations performed in 2006 October by XMM-Newton and Chandra in the 0.3 to 10 keV energy band. An acceptable model was found for the XMM-Newton data fitting the EPIC-PN and RGS spectra simultaneously; later, this model was successfully applied to the contemporary Chandra high-resolution data. The model consists of a continuum emission component (power law + blackbody) absorbed by four ionized components (warm absorbers), and 10 narrow emission lines. Three absorbing components are warm, producing features only in the soft X-ray band. The fourth ionization component produces Fe XXV and Fe XXVI in the hard-energy band. We study the time response of the absorbing components to the well-detected changes in the X-ray luminosity of this source and find that the two components with the lower ionization state show clear opacity changes consistent with gas close to photoionization equilibrium. These changes are supported by the models and by differences in the spectral features among the nine observations. On the other hand, the two components with higher ionization state do not seem to respond to continuum variations. The response time of the ionized absorbers allows us to constrain their electron density and location. We find that one component (with intermediate ionization) must be located within the obscuring torus at a distance 2.7 × 10{sup 17} cm from the central engine. This outflowing component likely originated in the accretion disk. The three remaining components are at distances larger than 10{sup 16}-10{sup 17} cm. Two of the absorbing components in the soft X-rays have similar outflow velocities and locations. These components may be in pressure equilibrium, forming a multi-phase medium, if the gas has metallicity larger than the solar one (≳ 5 Z {sub ☉}). We also search for variations in the covering factor of the ionized absorbers (although partial covering is not required in our models). We find no correlation between the change in covering factor and the flux of the source. This, in connection with the observed variability of the ionized absorbers, suggests that the changes in flux are not produced by this material. If the variations are indeed produced by obscuring clumps of gas, these must be located much closer in to the central source.

79 ASTRONOMY AND ASTROPHYSICS↗

Detailed analysis of low energy plasma data under the Voyager Uranus data analysis program

Research effort included the PLS data analysis program where modifications to the data fitting procedure and elimination of possible noise and electron contamination were made. The analysis code corrections were used in checking the Neptune data gathered during the Voyager 2 encounter and for analyzing selected plasma spectra from the warm Io torus. A major task accomplished was the summary of Uranus-related research in the U.S. National Report to the International Union of Geodesy and Geophysics for the 1987 - 1990 quadrennium. A limited amount of work was accomplished on assessing the Pedersen conductivity of the ionosphere and comparing it with inferred values from shielding by the Uranian ring current. Under this grant there has been a great deal of effort expended on identifying and classifying plasma waves and oscillations in the magnetosheath and solar wind downstream from Uranus. Large amplitude oscillations in plasma parameters are found in the magnetosheath, with density changes of up to a factor of ten occurring on times scales of minutes. New algorithms developed for analyzing the inbound bow shock crossing of Neptune will probably be applied to a more detailed analysis of the Uranus shock in the near future.

Mcnutt, Ralph L., Jr.↗

A Gaussian Process Enhancement to Linear Parameter Varying Models

Simulation and analysis for modern engineering systems now routinely requires the merging of multiple disciplines, physical-domains, time-scales, and data sets — all at ever increasing levels. These capabilities are especially needed in the domain of Advanced Air Mobility, where rapidly emerging vehicle designs are significantly more complex, while having to be both cost-effective and safe. To meet these engineering challenges, machine learning methods are an attractive option for merging models and data across multiple areas while providing uncertainty quantification and maintaining computational efficiency. This paper examines the use of Gaussian process machine learning to generalize and enhance the commonly used class of quasi-Linear Parameter Varying models for fast full-envelope simulation while also supporting control system design and analysis with model uncertainty. Gaussian process machine learning is selected because it: can fuse multiple data sets, enables an easy trade-off between data fitting and smoothing, provides model uncertainty quantification, scales well with increasing complexity, and does not generally require starting from a large training data set. To demonstrate the benefits of the approach, a robust stability analysis with Gaussian process uncertainty is shown for a NASA reference design of an electric quad-rotor air-taxi concept vehicle with motor parameter uncertainty.

Gaussian Process↗

Time-Dependent Reliability Analysis of Nuclear Hybrid Energy Systems

This paper focuses on developing a framework for comprehensive time-dependent reliability analysis of a nuclear hybrid energy system (NHES) design. Subsystem interactions of this complex integrated system under stochastic electricity demand and load-following operation capabilities requires dynamic reliability assessment of the NHES in component, subsystem and system level. In order to capture the dynamic operational behavior of critical systems and components, a detailed thermal hydraulic model of the NHES is generated using Modelica and tested under different electricity demand histories generated using RAVEN. The physical data (e.g., valve position, flow rates, inlet outlet water temperatures) gathered from Modelica used to calculate time-dependent failure rates through fitting data into the piecewise Weibull distribution. Optimum maintenance interval is calculated, and maintenance cost estimates based on the calculated reliability metric is presented for the selected component. Component time-dependent failure rates fed into the subsystem reliability model built with non-Markovian Stochastic Petri Nets. The component and subsystem failure rates are calculated and updated in every hour to characterize the system behavior appropriately and to understand operational reliability of the NHES design in a given time.

Guler Yigitoglu, Askin↗

Attitude and translation control of a low-altitude Gravsat

Dynamic simulation of the Gravsat's attitude and translation control system is used to provide an upper bound for the fuel supply and test the feasibility of the preliminary design. A preliminary design is made for the disturbance compensation system (DISCOS) sensor, the thruster control laws, reaction wheel control laws, and the onboard state estimators. The sensor analysis and noise measurements show no problems in scaling the Triad navigation satellite sensor design up to meet the Gravsat requirements, except for proof mass center-of-mass offset. A promising technique is proposed to measure and eliminate this error. The covariance analysis confirms that a sophisticated post-flight data fit will be necessary to reconstruct a scientifically useful proof mass state. The DISCOS sensor will have to be continuously calibrated from the inflight data to achieve this reconstruction.

Ray, J. C.↗

Precision geodesy and geodynamics using Starlette laser ranging

The French Starlette satellite, launched in February 1975, was the first satellite specifically designed to minimize the effects of nongravitational forces and to obtain the highest possible accuracy for laser range measurements. It has been found that Starlette represents a valuable complement to the U.S. geodetic satellites. In an analysis of Starlette laser ranging data conducted by Marsh and Williamson (1978), it was concluded that by tailoring a gravity model to a specific satellite and observation period, substantial improvements in data fits and, therefore, in orbit accuracy could be achieved. In the present analyses, a tailored gravity model has been derived for Starlette from the data acquired for the 4-year period, 1975-1978. Attention is given to the solution for geodetic and geodynamic parameters and polar motion.

Marsh, J. G.↗

Longitudinal Solid P Longitudinal Solid Polarized Target for CL get for CLAS12 and Study of Spin AS12 and Study of Spin Structure of Nucleons

A suite of experiments measuring target-spin observables in electron-nucleon scattering (dubbed Run Group C) was conducted at Jefferson Lab's Hall B in Newport News, VA with a new polarized nuclear target known as 'APOLLO' (Ammonia POLarized LOngitudinally). This innovative target is engineered to seamlessly integrate with the advanced 12GeV CEBAF (Continuous Electron Beam Accelerator Facility) accelerator and the Hall B CLAS12 (12 GeV CEBAF Large Acceptance Spectrometer) detector array. The 'APOLLO' target harnesses the power of Dynamic Nuclear Polarization (DNP) to achieve longitudinal polarization of solid ammonia, thereby creating a net polarization in both protons (NH3) and deuterons (ND3). These samples are subjected to a 5 Tesla magnetic field produced by the CLAS12 spectrometer central solenoid, cooled to 1 K using helium evaporation, and subsequently exposed to microwave radiation at a frequency of 140 GHz. This project was made possible through a collaborative partnership between the Jefferson Lab Target Group, Old Dominion University, University of Virginia, Christopher Newport University, and the CLAS Collaboration. Integrating the APOLLO target into CLAS12 presented unique challenges because of its specific spatial and dimensional constraints. This dissertation will highlight the innovative solutions developed to address these challenges, including the creation of a new target material transport system, the incorporation of superconducting magnetic correction coils, and the development of a customized Nuclear Magnetic Resonance (NMR) system. In addition to a detailed description of the development and operation of APOLLO, and an overview of RGC, this document also describes the global data fitting of proton asymmetries in the Deep Inelastic scattering (DIS) and the Resonance regions, and it offers insights into the comprehensive analysis of target polarization data.

Pandey, Pushpa↗