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 937 records · Page 52

Expert system shell to reason on large amounts of data

The current data base management systems (DBMS's) do not provide a sophisticated environment to develop rule based expert systems applications. Some of the new DBMS's come with some sort of rule mechanism; these are active and deductive database systems. However, both of these are not featured enough to support full implementation based on rules. On the other hand, current expert system shells do not provide any link with external databases. That is, all the data are kept in the system working memory. Such working memory is maintained in main memory. For some applications the reduced size of the available working memory could represent a constraint for the development. Typically these are applications which require reasoning on huge amounts of data. All these data do not fit into the computer main memory. Moreover, in some cases these data can be already available in some database systems and continuously updated while the expert system is running. This paper proposes an architecture which employs knowledge discovering techniques to reduce the amount of data to be stored in the main memory; in this architecture a standard DBMS is coupled with a rule-based language. The data are stored into the DBMS. An interface between the two systems is responsible for inducing knowledge from the set of relations. Such induced knowledge is then transferred to the rule-based language working memory.

Giuffrida, Gionanni↗

Incorporating Endmember Variability into Spectral Mixture Analysis Through Endmember Bundles

Variation in canopy structure and biochemistry induces a concomitant variation in the top-of-canopy spectral reflectance of a vegetation type. Hence, the use of a single endmember spectrum to track the fractional abundance of a given vegetation cover in a hyperspectral image may result in fractions with considerable error. One solution to the problem of endmember variability is to increase the number of endmembers used in a spectral mixture analysis of the image. For example, there could be several tree endmembers in the analysis because of differences in leaf area index (LAI) and multiple scatterings between leaves and stems. However, it is often difficult in terms of computer or human interaction time to select more than six or seven endmembers and any non-removable noise, as well as the number of uncorrelated bands in the image, limits the number of endmembers that can be discriminated. Moreover, as endmembers proliferate, their interpretation becomes increasingly difficult and often applications simply need the aerial fractions of a few land cover components which comprise most of the scene. In order to incorporate endmember variability into spectral mixture analysis, we propose representing a landscape component type not with one endmember spectrum but with a set or bundle of spectra, each of which is feasible as the spectrum of an instance of the component (e.g., in the case of a tree component, each spectrum could reasonably be the spectral reflectance of a tree canopy). These endmember bundles can be used with nonlinear optimization algorithms to find upper and lower bounds on endmember fractions. This approach to endmember variability naturally evolved from previous work in deriving endmembers from the data itself by fitting a triangle, tetrahedron or, more generally, a simplex to the data cloud reduced in dimension by a principal component analysis. Conceptually, endmember variability could make it difficult to find a simplex that both surrounds the data cloud and has vertices that are realistic endmember spectra with reflectances between 0 and 1. In this paper, we create endmember bundles and bounding fraction images for an AVIRIS subscene simulated with a plant canopy radiative transfer model. The simulated subscene is spatially patterned after a subscene from the AVIRIS image acquired August, 1993 over La Copita, Texas. In addition, for comparison, we performed a traditional unmixing with image endmembers.

Bateson, C. Ann↗

Substructure at High Speed. II. The Local Escape Velocity and Milky Way Mass with Gaia eDR3

Abstract Measuring the escape velocity of the Milky Way is critical in obtaining the mass of the Milky Way, understanding the dark matter velocity distribution, and building the dark matter density profile. In Necib & Lin, we introduced a strategy to robustly measure the escape velocity. Our approach takes into account the presence of kinematic substructures by modeling the tail of the stellar distribution with multiple components, including the stellar halo and the debris flow called the Gaia Sausage (Enceladus). In doing so, we can test the robustness of the escape velocity measurement for different definitions of the “tail” of the velocity distribution and the consistency of the data with different underlying models. In this paper, we apply this method to the Gaia eDR3 data release and find that a model with two components is preferred, although results from a single-component fit are also consistent. Based on a fit to retrograde data with two bound components to account for the relaxed halo and the Gaia Sausage, we find the escape velocity of the Milky Way at the solar position to be v esc = 445 − 8 + 25 km s −1 . A fit with a single component to the same data gives v esc = 472 − 12 + 17 km s −1 . Assuming a Navarro−Frenck−White dark matter profile, we find a Milky Way concentration of c 200 = 19 − 7 + 11 and a mass of M 200 = 4.6 − 0.8 + 1.5 × 10 11 M ⊙ , which is considerably lighter than previous measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

Multimodality in the Search for New Physics in Pulsar Timing Data and the Case of Kination-amplified Gravitational-wave Background from Inflation

We investigate the kination-amplified inflationary gravitational-wave background (GWB) interpretation of the signal recently reported by various pulsar timing array (PTA) experiments. Kination is a post-inflationary phase in the expansion history dominated by the kinetic energy of some scalar field, characterized by a stiff equation of state w = 1. Within the inflationary GWB model, we identify two modes that can fit the current data sets (NANOGrav and EPTA) with equal likelihood: the kination-amplification (KA) mode and the ordinary, no-kination-amplification (no-KA) mode. The multimodality of the likelihood motivates a Bayesian analysis with nested sampling. We analyze the free spectra of current PTA data and mock free spectra constructed with higher signal-to-noise ratios using nested sampling. The analysis of the mock spectrum designed to be consistent with the best fit to the NANOGrav 15 yr (NG15) data successfully reveals the expected bimodal posterior for the first time while excluding the reheating mode that appears in the fit to the current NG15 data, making a case for our correct and comprehensive treatment of potential multimodal posteriors arising from future PTA data sets. The resultant Bayes factor is $\mathcal{B}$ $\equiv$ Z no–KA /Z KA = 2.9 ± 1.9, indicating comparable statistical significance between the two modes. Given the theoretical model-building challenges of producing highly blue-tilted primordial tensor spectra, the KA mode has the advantage of requiring less blue primordial spectra, compared with the no-KA mode. The synergy between future cosmic microwave background polarization, pulsar timing, and laser interferometer measurements of gravitational waves will help resolve the ambiguity implied by the multimodal posterior in PTA-only searches.

Cosmology↗

Representation of the Geosynchronous Plasma Environment for Spacecraft Charging Calculations

Historically, our ability to predict and postdict surface charging has suffered from both a lack of reliable secondary emission and backscattered electron yields and poor characterization of the plasma environment. One difficulty lies in the common practice of fitting the plasma data to a Maxwellian or Double Maxwellian distribution function, which may not represent the data well for charging purposes. For 13 years Los Alamos National Laboratory (LANL) has been accumulating measurements of electron and proton spectra from Magnetospheric Plasma Analyzer (MPA) instruments aboard a series of geosynchronous satellites. These data provide both a plasma characterization and the potential of the instrument ground. We use electron and ion flux spectra measured by the LANL MPA to examine how the use of different spectral representations of the charged particle environment in computations of spacecraft potentials during magnetospheric substorms affects the accuracy of the results. We calculate the spacecraft potential using both the measured fluxes and several different fits to these fluxes. These flux measurements and fits have been corrected for the difference between the measured and calculated potential. The potentials computed using the measured fluxes, the best available material properties of graphite carbon, and a secondary electron escape fraction of 81%, are within a factor of three of the measured potential for nearly all the data. Using a Kappa fit to the electron distribution function and a Maxwellian fit to the ion distribution function gives agreement similar to the calculations using the actual data. Alternative spectral representations, including Maxwellian and double Maxwellian for both species, lead to less satisfactory agreement between predicted and measured potentials.

Davis, V. A.↗

Data Assimilation Into a Coupled Ocean Atmosphere Model: Application to the 1997-1998 El Nino

As part of JPL's ocean data assimilation effort to study ocean circulation and seasonal-interannual climate variability, sea level anomaly observed by TOPEX altimeter, together with sea surface temperature and wind stress data, are assimilated into a simple coupled ocean atmosphere model of the tropical Pacific. Model-data consistency is examined. Impact of the assimilation (as initialization) on El Nino Southern Oscillation (ENSO) forecasts is evaluated. The coupled model consists of a shallow water component with two baroclinic modes, an Ekman shear layer, a simplified mixed-layer temperature equation, and a statistical atmosphere based on dominant correlations between historical surface temperature and wind stress anomaly data. The adjoins method is used to fit the coupled model to the data over various six-month periods from late 1996 to early 1998 by optimally adjusting the initial state, model parameters, and basis functions of the statistical atmosphere. On average, the coupled model can be fitted to the data to approximately within the data and representation errors (5 cm, 0.5 C, and 10 sq m/sq m for sea level, surface temperature, and pseudo wind stress anomalies, respectively). The estimated fields resemble observed spatio-temporal structure reasonably well. Hindcasts/forecasts of the 1997/1998 El Nino initialized from forced estimated ocean states and parameters are much more realistic than those simply initialized from ocean states (see figure below). In particular, the ability of the model to produce significant warming beyond the initial state is dramatically improved. Parameter estimation, which compensates for some model errors, is found to be important to obtaining better fits of the model to data and to improving forecasts.

Lee, Tong↗

OmicsMLMentor: A Web Application for Guided Machine Learning Analysis of Omics Data

Expression-based omics technologies (e.g. proteomics, metabolomics, transcriptomics, etc.) increasingly rely on supervised and unsupervised machine learning (ML) models to find key biomolecules distinguishing conditions, identify natural groupings in biological data, or generate predictions for outcomes of interest. Fitting ML models to omics data presents several challenges, including handling missing data, selecting a normalization method, choosing a valid model, and optimizing hyperparameters, all requiring statistical programming skills to address these challenges. Thus, the open-source web application SLOPE was designed to lower the barrier to ML modeling for omics data. SLOPE supports the fitting of 15 ML models (10 supervised and 5 unsupervised) tailored to omics datasets, such as proteomics, metabolomics, lipidomics, and transcriptomics. SLOPE offers several omics-specific features, including methods for handling missingness (imputation, conversion, removal), normalization tests, ranking of models based on the structure of a user’s data and user input, and optimal hyperparameter selections using cross-validation splits. By streamlining ML workflows for omics analysis, SLOPE address critical gaps in existing online web tools, facilitating a broader adoption of these models for omics research. Here, SLOPE is applied to data from a lignin exposure study to highlight the workflow for fitting both supervised and unsupervised models to data.

lipidomics↗

A smoothing algorithm using cubic spline functions

Two algorithms are presented for smoothing arbitrary sets of data. They are the explicit variable algorithm and the parametric variable algorithm. The former would be used where large gradients are not encountered because of the smaller amount of calculation required. The latter would be used if the data being smoothed were double valued or experienced large gradients. Both algorithms use a least-squares technique to obtain a cubic spline fit to the data. The advantage of the spline fit is that the first and second derivatives are continuous. This method is best used in an interactive graphics environment so that the junction values for the spline curve can be manipulated to improve the fit.

Smith, R. E., Jr.↗

A new equation of state for copper

A new copper equation of state is developed utilizing the available experimental data in addition to recent theoretical calculations. Semi-empirical models are fit to the data and the results are tabulated in the SNL SESAME format. Comparison to other copper EOS tables are given, along with recommendations of which tables provide the best accuracy.

36 MATERIALS SCIENCE↗

Progress in Genetic Algorithm Fitting of X-Ray Fluorescence Data for Absorption Spectroscopy of Liquid Surfaces

X-ray fluorescence (XRF) spectroscopy of liquid surfaces is a promising approach to elemental analysis of forensic samples. Elements in the mixture can be identified from the XRF data because chemical elements are associated with unique spectra. This report describes progress in the development of automated identification of unknown mixtures of elements from XRF spectra using a genetic algorithm (GA). The XRF spectra are measured at the Advanced Photon Source (APS) at the liquid-liquid interface during solvent extraction process. The usage of the GA is demonstrated and future implementation for related problems are then discussed.

47 OTHER INSTRUMENTATION↗

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - Simulated Marine Hydrokinetic Tidal Turbine

The U.S. Department of Energy and National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis, hydrogen compression and storage, and variable hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) initiative. This dataset is part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with other energy technologies. This dataset contains inputs and outputs from simulations of a floating marine hydrokinetic turbine over approximately half a tidal cycle (~6.6 hours). Inflow conditions were derived from field measurements in Alaska’s Cook Inlet and represent a tidal environment in which the current speed ramps from near 0 m/s to a peak of 3 m/s and back. The original acoustic doppler current profiler dataset is publicly available on the Marine and Hydrokinetic Data Repository. In a full tidal cycle, the flow reverses and the rotor would reorient; this reversal was not modeled. In the Cook Inlet campaign , turbulence intensity was similar in both directions. Two inflow cases are included. In the first case, labeled “raw” in the files, the measured current time series was used directly in the InflowWind module of OpenFAST. Speed and direction were applied as a function of time and elevation, uniformly in the horizontal direction. With full spatial coherence, this approach captures high turbulent variability and results in pronounced power fluctuations, so it is considered a conservative, near-worst-case representation of loading. In the second case, labeled “average” in the files, a 30-minute moving average was applied to extract the slowly varying mean speed. The residual fluctuations about this mean were used to generate spatially varying, full-field turbulence inputs with TurbSim, giving a more physically realistic representation of the inflow across the rotor disk. Two random realizations were used to produce distinct inflow conditions for two OpenFAST simulations representing a two-turbine array. The same turbulence intensity is applied across the full time series, producing larger fluctuations at the start and end, where the mean speed is low. The second case is the more appropriate framework for performance and power assessment but overpredicts turbulence at lower flow speeds and underpredicts it at higher speeds. As the floating platform moves and the rotor changes its x-position, Taylor’s frozen turbulence hypothesis used by InflowWind assumes a constant rather than a time-varying mean velocity, introducing some inaccuracy in the velocity plane sampling. The turbine modeled is the 500-kW Reference Model 1, a horizontal-axis two-bladed hydrokinetic turbine on a four-column floating semisubmersible substructure . Simulations were performed using OpenFAST v4.1 with the Reference Open Source Controller (ROSCO) v2.10. All input files required to reproduce the simulations are included. The electrolyzer is a 1.25-MW proton exchange membrane type MC250 system manufactured by Nel . This unit supports up to 2.5 MW, but NLR has only a single 1.25-MW stack. The datasets report hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. The system controls hydrogen production by varying direct current applied to the stack, from a maximum of 3,000 A to a minimum safe operating current of 300 A, or 10%. Because the current–voltage characteristic changes as the stack ages and efficiency degrades, the actual minimum safe operating power changes over time. The simulated tidal turbine time series data was translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz. Each zip file represents a single tidal electrolysis experiment and is named: {technology}_{inflow method}_{number of 500 kW tidal turbines connected} For instance, “tidal-500kW-RM1_average_2.zip” is a 6-hour experiment using the 500-kW tidal reference model, scaled by 2x (1-MW) to better match the electrolyzer maximum of 1.25MW, fed with the 30-minute moving average current case. Each zip folder contains the following files: A .csv file of raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production in kilograms per hour, electrolysis power consumption, and input wave power. A .csv file combines all tidal profiles as "combined_tidal_experiments.csv." A separate experiment, “characterization_200.zip,” shows the MC250 electrolyzer steady-state response with 30-minute load steps over 5 hours and is accessible with this entry.

08 HYDROGEN↗

High-dispersion spectroscopic observations of Venus near superior conjunction. I - The carbon dioxide band at 7820 A

Nine plates of the 7820-A CO2 band were taken in 1971. A curve-of-growth analysis of the CO2 lines indicates a rotational temperature of 241 + or - 2 K, with an average slope to the curve of growth of 0.60 + or - 0.03. The Venus phase angle ranged from 7.2 to 10.7 deg. The equivalent widths of the 1971 data fall on a smooth curve fit through the 1969 data for this band; there does not appear to be any discontinuity in the phase curve at small phase angles.

Young, L. D. G.↗

Analysis of the physical atomic forces between noble gas atoms, alkali ions and halogen ions

The physical forces between atoms and molecules are important in a number of processes of practical importance, including line broadening in radiative processes, gas and crystal properties, adhesion, and thin films. The components of the physical forces between noble gas atoms, alkali ions, and halogen ions are analyzed and a data base for the dispersion forces is developed from the literature based on evaluations with the harmonic oscillator dispersion model for higher order coefficients. The Zener model of the repulsive core is used in the context of the recent asymptotic wave functions of Handler and Smith; and an effective ionization potential within the Handler and Smith wave functions is defined to analyze the two body potential data of Waldman and Gordon, the alkali-halide molecular data, and the noble gas crystal and salt crystal data. A satisfactory global fit to this molecular and crystal data is then reproduced by the model to within several percent. Surface potentials are evaluated for noble gas atoms on noble gas and salt crystal surfaces with surface tension neglected. Within this context, the noble gas surface potentials on noble gas and salt crystals are considered to be accurate to within several percent.

Wilson, J. W.↗

Atomic forces between noble gas atoms, alkali ions, and halogen ions for surface interactions

The components of the physical forces between noble gas atoms, alkali ions, and halogen ions are analyzed and a data base developed from analysis of the two-body potential data, the alkali-halide molecular data, and the noble gas crystal and salt crystal data. A satisfactory global fit to this molecular and crystal data is then reproduced by the model to within several percent. Surface potentials are evaluated for noble gas atoms on noble gas surfaces and salt crystal surfaces with surface tension neglected. Within this context, the noble gas surface potentials on noble gas and salt crystals are considered to be accurate to within several percent.

Wilson, J. W.↗

A high temperature thermal diffusivity determination procedure for solids and liquids

A new method for measuring the thermal diffusivity of materials at high temperatures is presented. The method is applicable to solids on Earth, and to liquids in the reduced gravity environment of space. It is especially suited to levitated liquid metals at elevated temperatures where thermal diffusivity data is not available. The method is applied in two parts, such that lumped analysis is valid in the first part, and Fourier's law of conduction is valid in the second. In both parts, the spherical specimen is assumed to have been heated to a desired temperature and cooled. An inverse conduction problem is then formulated and solved using Laplace transformation techniques. Using this solution, sample sizes, and experimentally obtained surface temperature history, the thermal diffusivity is determined by minimizing a function that satisifies the heat balance at the surface. Minimization is performed using a modified quasilinearization algorithm. Accuracy is very sensitive to error in the temperature data and increases with better curve fits to the temperature data. An error analysis is also performed, and the effect of errors in the various parameters on the evaluated thermal diffusivity is determined. An experimental study for solids on Earth is suggested, before development for implementation in space.

Bayazitoglu, Yildiz↗

Analysis of stellar occultation data for planetary atmospheres. I - Model fitting, with application to Pluto

Consideration is given to an analytic model for a stellar-occultation light curve developed for a small, spherically symmetric planetary atmosphere that includes thermal and molecular weight gradients in a region that overlies an extinction layer. The model incorporates two equivalent sets of parameters. One set specifies the occultation light curve in terms of signal levels, times, and time intervals. The other set specifies physical parameters of the planetary atmosphere. Equations are given for the transforming between the sets of parameters, including their errors and correlation coefficients. Detailed numerical calculations are presented for a benchmark case. The results obtained are consistent with the isothermal prediction of the 'methane-thermostat' model of Pluto's atmosphere.

Elliot, J. L.↗

Preliminary experimental results for a cryogenic brush seal configuration

Preliminary fluid nitrogen flow data are reported for a five-brush, ceramic-coated-rub-runner brush seal system, where the brushes and the rub runner were placed at each end of a centrally pressurized multifunction tester ('back-to-back' set of brushes) and tested at rotor speeds of 0, 10, 18, and 22.5 krpm. After testing, both the brushes and the ceramic-coated rub runner appeared pristine. The coating withstood both the thermomechanical and dynamic loadings with minor wear track scarring. The bristle tips showed some indication of material shearing (smearing) wear. The Ergun porous flow equation was applied to the brush seal data. The Ergun relation, which required heuristic information to characterize the coefficients, fit the gaseous data but was in poor agreement with the fluid results. The brush seal exit conditions were two phase. Two-phase, choked-flow design charts were applied but required one data point at each rotor speed to define the (C(sub f)A x Constant) flow and area coefficients. Reasonable agreement between prediction and data was found, as expected, but such methods are not to be construed as two-phase-flow brush seal analyses.

Carlile, J. A.↗

Rapid soft X-ray variability in NGC 3227

Lightcurves from the Rosat position sensitive proportional counter are presented for the Seyfert galaxy NGC3227 at two different epochs separated by two years. These lightcurves show statistically significant variability, the most rapid being a factor of 2 increase in a period of 20 min. Spectral analysis shows a simple power law featuring only Galactic absorption, can usually characterize the observed spectrum in both the low and high states. If included in the extraction of the high state spectrum, the simple power law model no longer provides an entirely acceptable fit. The addition of either a warm or cold absorber enables the data to be fitted equally well.

Stringfellow, G. S.↗