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 505 records · Page 28

Thermal Output of WK-Type Strain Gauges on Various Materials at Cryogenic and Elevated Temperatures

Strain gage apparent strain (thermal output) is one of the largest sources of error associated with the measurement of strain when temperatures and mechanical loads are varied. In this paper, experimentally determined apparent strains of WK-type strain gages, installed on both metallic and composite-laminate materials of various lay-ups and resin systems for temperatures ranging from -450 F to 230 F are presented. For the composite materials apparent strain in both the 0 ply orientation angle and the 90 ply orientation angle were measured. Metal specimens tested included: aluminum-lithium alloy (Al-LI 2195-T87), aluminum alloy (Al 2219-T87), and titanium alloy. Composite materials tested include: graphite-toughened-epoxy (IM7/997- 2), graphite-bismaleimide (IM7/5260), and graphite-K3 (IM7/K3B). The experimentally determined apparent strain data are curve fit with a fourth-order polynomial for each of the materials studied. The apparent strain data and the polynomials that are fit to the data are compared with those produced by the strain gage manufacturer, and the results and comparisons are presented. Unacceptably high errors between the manufacture's data and the experimentally determined data were observed (especially at temperatures below - 270-F).

Kowalkowski, Matthew K.↗

Global fit of electron and neutrino elastic scattering data to determine the strange quark contribution to the vector and axial form factors of the nucleon

We present a global fit of neutral-current elastic (NCE) neutrino-scattering data and parity-violating electron-scattering (PVES) data with the goal of determining the strange quark contribution to the vector and axial form factors of the proton. Previous fits of this form included data from a variety of PVES experiments (PVA4, HAPPEx, G0, SAMPLE) and the NCE neutrino and anti-neutrino data from BNL E734. These fits did not constrain the strangeness contribution to the axial form factor G A s ( Q 2 ) at low Q 2 very well because there was no NCE data for Q 2 < 0.45 GeV 2 . Our new fit includes for the first time MiniBooNE NCE data from both neutrino and antineutrino scattering; this experiment used a hydrocarbon target and so a model of the neutrino interaction with the carbon nucleus was required. Three different nuclear models have been employed: a relativistic Fermi gas model, the superscaling approximation model, and a spectral function model. We find a tremendous improvement in the constraint of G A s ( Q 2 ) at low Q 2 compared to previous work, although more data is needed from NCE measurements that focus on exclusive single-proton final states, for example from MicroBooNE. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Error enhancement in geomagnetic models derived from scalar data

Models of the main geomagnetic field are generally represented by a scalar potential gamma expanded in a finite number of spherical harmonics. Very accurate observations of F were used, but indications exist that the accuracy of models derived from them is considerably lower. One problem is that F does not always characterize gamma uniquely. It is not clear whether such ambiguity can be encountered in deriving gamma from F in geomagnetic surveys, but there exists a connection, due to the fact that the counterexamples of Backus are related to the dipole field, while the geomagnetic field is dominated by its dipole component. If the models are recovered with a finite error (i.e. they cannot completely fit the data and consequently have a small spurious component), this connection allows the error in certain sequences of harmonic terms in gamma to be enhanced without unduly large effects on the fit of F to the model.

Stern, D. P.↗

Performance characterization of nanofiltration, reverse osmosis, and ion exchange technologies for acetic acid separation

A major obstacle to widespread implementation of bio-based fuels and chemicals is the lack of efficient and costeffective separation methods. To purify acetic acid produced by biochemical conversion of biomass via anaerobic digestion, this work employs two commonly used separation technologies: (1) ion-exchange (IX) resin and (2) pressure-driven membranes. This study tested five commercially available strong- and weak-base anion-exchange resins and five commercially available nanofiltration (NF) and reverse osmosis (RO) membranes. The pH of the feed solution significantly affected adsorption capacity. At pH 6.3, a strong-base IX resin (IRN-78) performed best (95.1% acetate removal). With strong-base IX resins, the Langmuir isotherm model fit well, whereas for weak-base IX resins, the Freundlich isotherm provided a better fit. A pseudo-second-order kinetic model fit well for both IRN-78 and IRA-67. Regarding membrane separation, RO (BW30XFR membrane) achieved the highest rejection (98.6% acetate rejection), whereas an NF membrane (NF*) achieved the best combination of permeate flux (105 L/(h∙m 2 )) and rejection (83.1% acetate rejection). For membrane performance, the experimental data were fit using the solution diffusion model. Increased pH in the feed solution lowered permeate flux but increased acetic acid rejection. Finally, then the acetic acid concentration in the feed solution increased, both permeate flux and acetic acid rejection decreased for membrane NF*.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ring Current Development During Storm Main Phase

The development of the ring current ions in the inner magnetosphere during the main phase of a magnetic storm is studied. The temporal and spatial evolution of the ion phase space densities in a dipole field are calculated using a three dimensional ring current model, considering charge exchange and Coulomb losses along drift paths. The simulation starts with a quiet time distribution. The model is tested by comparing calculated ion fluxes with Active Magnetospheric Particle Tracer Explorers/CCE measurement during the storm main phase on May 2, 1986. Most of the calculated omnidirectional fluxes are in good agreement with the data except on the dayside inner edge (L less than 2.5) of the ring current, where the ion fluxes are underestimated. The model also reproduces the measured pitch angle distributions of ions with energies below 10 keV. At higher energy, an additional diffusion in pitch angle is necessary in order to fit the data. The role of the induced electric field on the ring current dynamics is also examined by simulating a series of substorm activities represented by stretching and collapsing the magnetic field lines. In response to the impulsively changing fields, the calculated ion energy content fluctuates about a mean value that grows steadily with the enhanced quiescent field.

Fok, Mei-Ching↗

PSR J0007+7303 in the CTA1 SNR: New Gamma-ray Results from Two Years of Fermi-LAT Observations

One of the main results of the Fermi Gamma-Ray Space Telescope is the discovery of gamma-ray selected pulsars. The high magnetic field pulsar, PSR J0007+7303 in CTA1, was the first ever to be discovered through its gamma-ray pulsations. Based on analysis of 2 years of LAT survey data, we report on the discovery of I-ray emission in the off-pulse phase interval at the approx. 6sigma level. The flux from this emission in the energy range E > or =::: 100 MeV is F(sub 100) = (1.73+/-0.40) x 10(exp -8) photons/sq cm/s and is best fitted by a power law with a photon index of Gamma = 2.54+/-0.14. The pulsed gamma-ray flux in the same energy range is F(sub 100) = (3.95+/-0.07) x 10(exp -7) photons/sq cm/s and is best fitted by an exponentially-cutoff power-law spectrum with a photon index of Gamma = 1.41+/-0.23 and a cutoff energy E(sub c) = 4.04+/-0.20 GeV. We find no flux variability neither at the 2009 May glitch nor in the long term behavior. We model the gamma-ray light curve with two high-altitude emission models, the outer gap and slot gap, and find that the model that best fits the data depends strongly on the assumed origin of the off-pulse emission. Both models favor a large angle between the magnetic axis and observer line of sight, consistent with the nondetection of radio emission being a geometrical effect. Finally we discuss how the LAT results bear on the understanding of the cooling of this neutron star.

Abdo, A. A.↗

GPS measurements of strain accumulation across the Imperial Valley, California: 1986-1989

The Global Positioning System (GPS) data collected in southern California from 1986 to 1989 indicate considerable strain accumulation across the Imperial Valley. Displacements are computed at 29 stations in and near the valley from 1986 to 1988, and at 11 sites from 1988 to 1989. The earlier measurements indicate 5.9 +/- 1.0 cm/yr right-lateral differential velocity across the valley, although the data are heavily influenced by the 1987 Superstition Hills earthquake sequence. Some measurements, especially the east-trending displacements, are suspects for large errors. The 1988 to 1989 GPS displacements are best modeled by 5.2 +/- 0.9 cm/yr of valley crossing deformation, but rates calculated from conventional geodetic measurements (3.4 to 4.3 cm/yr) fit the data nearly as well. There is evidence from GPS and Very Long Base Interferometry (VLBI) observations that the present slip rate along the southern San Andreas fault is smaller than the long-term geologic estimate, suggesting a lower earthquake potential than is currently assumed. Correspondingly, a higher earthquake potential is indicated for the San Jacinto fault. The Imperial Valley GPS sites form part of a 183 station network in southern California and northern Baja California, which spans a cross-section of the North American-Pacific plate boundary.

Larsen, Shawn↗

Global Positioning System measurements of strain accumulation across the Imperial Valley, California - 1986-1989

The Global Positioning System (GPS) data collected in southern California from 1986 to 1989 indicate considerable strain accumulation across the Imperial Valley. Displacements are computed at 29 stations in and near the valley from 1986 to 1988, and at 11 sites from 1988 to 1989. The earlier measurements indicate 5.9 =/- 1.0 cm/yr right-lateral differential velocity across the valley, although the data are heavily influenced by the 1987 Superstition Hills earthquake sequence. Some measurements, especially the east-trending displacements, are suspects for large errors. The 1988 to 1989 GPS displacements are best modeled by 5.2 =/- 0.9 cm/yr of valley crossing deformation, but rates calculated from conventional geodetic measurements (3.4 to 4.3 cm/yr) fit the data nearly as well. There is evidence from GPS and Very Long Base Interferometry (VLBI) observations that the present slip rate along the southern San Andreas fault is smaller than the long-term geologic estimate, suggesting a lower earthquake potential than is currently assumed. Correspondingly, a higher earthquake potential is indicated for the San Jacinto fault. The Imperial Valley GPS sites form part of a 183 station network in southern California and northern Baja California, which spans a cross-section of the North American-Pacific plate boundary.

Larsen, Shawn↗

Inverting infrasound data for the seismoacoustic source time functions and surface spall at the Source Physics Experiments Phase II: Dry Alluvium Geology

This report presents the infrasound data recorded as part of the Source Physics Experiment - Phase 2, Dry Alluvium Geology. This experiment, also known colloquially as DAG, consisted of four underground chemical explosions at the Nevada National Security Site. We focus our analysis on only the fourth explosion (DAG-4) as we determined that this was the only event that produced clear source-generated infrasound energy as recorded by the DAG sensors. We analyze the data using two inversion methods. The first method is designed to estimate the point-source seismoacoustic source time functions, and the second inversion method is designed to estimate the first-order characteristics (e.g. horizontal dimensions and maximum amplitude) of the actual spall surface. For both analysis methods, we are able to fit the data reasonably well, with various assumptions of the source model. The estimated seismoacoustic source appears to be a combination of a buried, isotropic explosion with a maximum amplitude of ~2 x 10 9 Nm and a vertically oriented force, applied to the Earth's surface with a maximum amplitude of 4 x 10 7 N. We use the vertically oriented force to simulate surface spall. The estimated spall surface has an approximate radius of ~40 m with a maximum acceleration magnitude in the range of 0.8 to 1.5 m/s/s. These estimates are approximately similar to the measured surface acceleration at the site.

58 GEOSCIENCES↗

Computer codes for the evaluation of thermodynamic properties, transport properties, and equilibrium constants of an 11-species air model

The computer codes developed provide data to 30000 K for the thermodynamic and transport properties of individual species and reaction rates for the prominent reactions occurring in an 11-species nonequilibrium air model. These properties and the reaction-rate data are computed through the use of curve-fit relations which are functions of temperature (and number density for the equilibrium constant). The curve fits were made using the most accurate data believed available. A detailed review and discussion of the sources and accuracy of the curve-fitted data used herein are given in NASA RP 1232.

Thompson, Richard A.↗

MB-Fit: Software infrastructure for data-driven many-body potential energy functions

Many-body potential energy functions (MB-PEFs), which integrate data-driven representations of many-body short-range quantum mechanical interactions with physics-based representations of many-body polarization and long-range interactions, have recently been shown to provide high accuracy in the description of molecular interactions from the gas to the condensed phase. Here, we present MB-Fit, a software infrastructure for the automated development of MB-PEFs for generic molecules within the TTM-nrg (Thole-type model energy) and MB-nrg (many-body energy) theoretical frameworks. Besides providing all the necessary computational tools for generating TTM-nrg and MB-nrg PEFs, MB-Fit provides a seamless interface with the MBX software, a many-body energy and force calculator for computer simulations. Given the demonstrated accuracy of the MB-PEFs, particularly within the MB-nrg framework, we believe that MB-Fit will enable routine predictive computer simulations of generic (small) molecules in the gas, liquid, and solid phases, including, but not limited to, the modeling of quantum isomeric equilibria in molecular clusters, solvation processes, molecular crystals, and phase diagrams.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reanalysis of x ray emission from M87 II: Multiphase models

In a previous paper by the authors (hereafter TB), it was shown that the hot intracluster gas around M87 could not be adequately explained in terms of a spherically symmetric, single phase model. It was found that although data from the Einstein satellite High Resolution Imager (HRI), Imaging Proportional Counter (IPC), and Focal Point Crystal Spectrometer (FPCS) could indeed be simultaneously explained by a single phase model, data from the Solid State Spectrometer (SSS) and optically determined mass estimates could not be similarly explained. A qualitative discussion was then given to indicate how the adoption of a specific multiphase model could plausibly improve upon the single phase model, but a detailed discussion was delayed until now. The x ray data is reanalyzed for the gas around M87 in the same spirit as the work in TB. That is, by assuming several cooling flow models, the shortcomings are improved upon. The procedure of TB was adopted and the data reanalyzed in the following manner: a multiphase model is assumed which can compute the surface brightness as seen by HRI and IPC. Also, the line fluxes can be computed of the lines seen by the FPCS and the mass profile after assuming hydrostatic equilibrium. The parameters of the model are then adjusted to fit the data. A check is then made to see whether the resulting model is consistent with the SSS spectrum and the equivalent width of the 7 keV complex of Fe lines as seen by large field of view instruments.

Tsai, John C.↗

Infrared Communications for Small Spacecraft: From a Wireless Bus to Cluster Concepts

Nanosatellites operating singly or in clusters are anticipated for future space science missions. To implement this new communications paradigm, we are approaching cluster communications by first developing an infrared (IR) intra-craft wireless bus capability, following initially the MIL-STD-1553B protocol. Benefits of an IR wireless bus are low mass, size, power, and cost, simplicity of implementation, ease of use, minimum EMI, and efficient and reliable data transfer. Our goals are to maximize the reliable link margin in order to afford greater flexibility in receiver placement, which will ease technology insertion. We have developed a concept demonstration using a high-speed visible-band silicon PIN photodiode and a high-efficiency visible LED operating at a data rate up to 4 Mb/sec. In designing an internal IR wireless bus, we have characterized various candidate materials, emitters, and geometries, assuming a single reflection. Thus, we have measured the bidirectional reflectance distribution function (BRDF) for five different materials characteristic of typical spacecraft structures, which range from nearly Lambertian to highly specular. We have fit our data to empirical BRDF functions and modeled the detected irradiance anywhere in the plane of incidence for a divergent (LED) emitter. We have also determined the angular limits on the link geometry to remain within the required bit error rate by determining the received signal-to-noise ratio (SNR) for minimum values of irradiance received at the detector.

Webb, Suzanne C.↗

Constrained curve fitting for semi-parametric models with radial basis function networks

Common to many analysis pipelines in lattice gauge theory and the broader scientific discipline is the need to fit a semi-parametric model to data. We propose a fit method that utilizes a radial basis function network to approximate the non-parametric component of such models. The approximate parametric model is fit to data using the basin hopping global optimization algorithm. Parameter constraints are enforced through Gaussian priors. The viability of our method is tested by examining its use in a finite-size scaling analysis of the $q$-state Potts model and $p$-state clock model with $q=2,3$ and $p=4,\infty$.

Peterson, Curtis T.↗

Estimating the Properties of Hard X-Ray Solar Flares by Constraining Model Parameters

We wish to better constrain the properties of solar flares by exploring how parameterized models of solar flares interact with uncertainty estimation methods. We compare four different methods of calculating uncertainty estimates in fitting parameterized models to Ramaty High Energy Solar Spectroscopic Imager X-ray spectra, considering only statistical sources of error. Three of the four methods are based on estimating the scale-size of the minimum in a hypersurface formed by the weighted sum of the squares of the differences between the model fit and the data as a function of the fit parameters, and are implemented as commonly practiced. The fourth method is also based on the difference between the data and the model, but instead uses Bayesian data analysis and Markov chain Monte Carlo (MCMC) techniques to calculate an uncertainty estimate. Two flare spectra are modeled: one from the Geostationary Operational Environmental Satellite X1.3 class flare of 2005 January 19, and the other from the X4.8 flare of 2002 July 23.We find that the four methods give approximately the same uncertainty estimates for the 2005 January 19 spectral fit parameters, but lead to very different uncertainty estimates for the 2002 July 23 spectral fit. This is because each method implements different analyses of the hypersurface, yielding method-dependent results that can differ greatly depending on the shape of the hypersurface. The hypersurface arising from the 2005 January 19 analysis is consistent with a normal distribution; therefore, the assumptions behind the three non- Bayesian uncertainty estimation methods are satisfied and similar estimates are found. The 2002 July 23 analysis shows that the hypersurface is not consistent with a normal distribution, indicating that the assumptions behind the three non-Bayesian uncertainty estimation methods are not satisfied, leading to differing estimates of the uncertainty. We find that the shape of the hypersurface is crucial in understanding the output from each uncertainty estimation technique, and that a crucial factor determining the shape of hypersurface is the location of the low-energy cutoff relative to energies where the thermal emission dominates. The Bayesian/MCMC approach also allows us to provide detailed information on probable values of the low-energy cutoff, Ec, a crucial parameter in defining the energy content of the flare-accelerated electrons. We show that for the 2002 July 23 flare data, there is a 95% probability that Ec lies below approximately 40 keV, and a 68% probability that it lies in the range 7-36 keV. Further, the low-energy cutoff is more likely to be in the range 25-35 keV than in any other 10 keV wide energy range. The low-energy cutoff for the 2005 January 19 flare is more tightly constrained to 107 +/- 4 keV with 68% probability.

X-rays↗

Toward particle-resolved accuracy in Euler–Lagrange simulations of multiphase flow using machine learning and pairwise interaction extended point-particle (PIEP) approximation

This study presents two different machine learning approaches for the modeling of hydrodynamic force on particles in a particle-laden multiphase flow. Results from particle-resolved direct numerical simulations (PR-DNS) of flow over a random array of stationary particles for eight combinations of particle Reynolds number (Re) and volume fraction (φ) are used in the development of the models. The first approach follows a two-step process. In the first flow prediction step, the perturbation flow due to a particle is obtained as an axisymmetric superposable wake using linear regression. In the second force prediction step, the force on a particle is evaluated in terms of the perturbation flow induced by all its neighbors using the generalized Faxén form of the force expression. In the second approach, the force data on all the particles from the PR-DNS simulations are used to develop an artificial neural network (ANN) model for direct prediction of force on a particle. Due to the unavoidable limitation on the number of fully resolved particles in the PR-DNS simulations, direct force prediction with the ANN model tends to over-fit the data and performs poorly in the prediction of test data. In contrast, due to the millions of grid points used in the PR-DNS simulations, accurate flow prediction is possible, which then allows accurate prediction of particle force. This hybridization of multiphase physics and machine learning is particularly important, since it blends the strength of each, and the resulting pairwise interaction extended point-particle model cannot be developed by either physics or machine learning alone.

42 ENGINEERING↗

Time-average fluxes of heavy ions at the geostationary orbit

Measurements of average proton, helium, carbon, and oxygen fluxes at 6.6 earth radii are reported. The data represent averages obtained on ATS 6 between June 15 and Oct. 3, 1974. The energy range covered was 0.36-1.1 MeV for protons and approximately 1-4 MeV for the heavier ions. The results indicate that above about 1 MeV the heavier ion fluxes dominate over the proton flux on the energy/ion scale. Using two different spectral dependencies to fit the data, the carbon to oxygen concentration ratio for energies above 1 MeV was found to be 0.43 for the power law spectrum and 0.44 for the exponential spectrum. Thus in either case the abundance ratio is consistent with the solar origin of the particles. Similarly, the helium to oxygen concentration ratio is found to be 0.17 for the power law spectrum and 0.28 for the exponential spectrum. This is inconsistent with the solar wind ratio, which is about 10. The results quoted above are based on a very small portion of the distribution function above 1 MeV.

Konradi, A.↗

Application of HCMM data to regional geologic analysis for mineral and energy resource evaluation

Experimentation with several potentially promising techniques led to the selection of a fairly simple scheme for registration of data from the HCMM using an affine transformation. A method based solely on remote sensing was developed to estimate those meteorological effects which are required for thermal inertia mapping. It assumes that the atmospheric fluxes are spatially invariant and that the solar, sky and sensible heat fluxes can be approximated by a simple mathematical form. Coefficients are determined from least squares method by fitting observational data to the thermal model.

Watson, K.↗