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 1,243 records · Page 69

Classifying and analyzing small-angle scattering data using weighted k nearest neighbors machine learning techniques

A consistent challenge for both new and expert practitioners of small-angle scattering (SAS) lies in determining how to analyze the data, given the limited information content of said data and the large number of models that can be employed. Machine learning (ML) methods are powerful tools for classifying data that have found diverse applications in many fields of science. Here, ML methods are applied to the problem of classifying SAS data for the most appropriate model to use for data analysis. The approach employed is built around the method of weighted k nearest neighbors (wKNN), and utilizes a subset of the models implemented in the SasView package (https://www.sasview.org/) for generating a well defined set of training and testing data. The prediction rate of the wKNN method implemented here using a subset of SasView models is reasonably good for many of the models, but has difficulty with others, notably those based on spherical structures. A novel expansion of the wKNN method was also developed, which uses Gaussian processes to produce local surrogate models for the classification, and this significantly improves the classification accuracy. Further, by integrating a stochastic gradient descent method during post-processing, it is possible to leverage the local surrogate model both to classify the SAS data with high accuracy and to predict the structural parameters that best describe the data. The linking of data classification and model fitting has the potential to facilitate the translation of measured data into results for both novice and expert practitioners of SAS.

97 MATHEMATICS AND COMPUTING↗

Fast Gaussian Process Estimation for Large-Scale In Situ Inference using Convolutional Neural Networks

Exascale computing will bring with it significant I/O limitations. One foreseeable consequence of such restrictions is that the user can save only a small fraction of complex simulation data to disk for subsequent analysis. An alternative is to fit statistical models to data in situ, that is, inside the simulation as it runs. This option requires extremely fast statistical estimation to avoid slowing down the simulation. Gaussian processes (GPs) have state-of-the-art predictive performance for modeling spatial data. However, standard estimation methods for GPs scale quite poorly to large data sets as parameter estimation requires inverting a covariance matrix to the size of the data set. In the presented work, we use a convolutional neural network (CNN) to predict the GP parameters for a spatial data set, from a simulation or otherwise, rather than optimize the parameters directly. Here, our presented case study models spatial data from E3SM, the Department of Energy’s Exascale climate model. The CNN is trained on synthetic data simulated from GP models with known parameters and then applied to data from the climate simulation. In the presented examples, the neural network scheme produces parameter estimates that compare well with standard methods such as maximum likelihood estimation in predictive performance but is obtained four orders of magnitude faster.

big data↗

Measurements of the Influence of Integral Length Scale on Stagnation Region Heat Transfer

The purpose was twofold: first, to determine if a length scale existed that would cause the greatest augmentation in stagnation region heat transfer for a given turbulence intensity and second, to develop a prediction tool for stagnation heat transfer in the presence of free stream turbulence. Toward this end, a model with a circular leading edge was fabricated with heat transfer gages in the stagnation region. The model was qualified in a low turbulence wind tunnel by comparing measurements with Frossling's solution for stagnation region heat transfer in a laminar free stream. Five turbulence generating grids were fabricated; four were square mesh, biplane grids made from square bars. Each had identical mesh to bar width ratio but different bar widths. The fifth grid was an array of fine parallel wires that were perpendicular to the axis of the cylindrical leading edge. Turbulence intensity and integral length scale were measured as a function of distance from the grids. Stagnation region heat transfer was measured at various distances downstream of each grid. Data were taken at cylinder Reynolds numbers ranging from 42,000 to 193,000. Turbulence intensities were in the range 1.1 to 15.9 percent while the ratio of integral length scale to cylinder diameter ranged from 0.05 to 0.30. Stagnation region heat transfer augmentation increased with decreasing length scale. An optimum scale was not found. A correlation was developed that fit heat transfer data for the square bar grids to within +4 percent. The data from the array of wires were not predicted by the correlation; augmentation was higher for this case indicating that the degree of isotropy in the turbulent flow field has a large effect on stagnation heat transfer. The data of other researchers are also compared with the correlation.

Vanfossen, G. James↗

Age-specific case data reveal varying dengue transmission intensity in US states and territories

Dengue viruses (DENV) are endemic in the US territories of Puerto Rico, American Samoa, and the US Virgin Islands, with focal outbreaks also reported in the states of Florida and Hawaii. However, little is known about the intensity of dengue virus transmission over time and how dengue viruses have shaped the level of immunity in these populations, despite the importance of understanding how and why levels of immunity against dengue may change over time. These changes need to be considered when responding to future outbreaks and enacting dengue management strategies, such as guiding vaccine deployment. We used catalytic models fitted to case surveillance data stratified by age from the ArboNET national arboviral surveillance system to reconstruct the history of recent dengue virus transmission in Puerto Rico, American Samoa, US Virgin Islands, Florida, Hawaii, and Guam. We estimated average annual transmission intensity (i.e., force of infection) of DENV between 2010 and 2019 and the level of seroprevalence by age group in each population. We compared models and found that assuming all reported cases are secondary infections generally fit the surveillance data better than assuming all cases are primary infections. Using the secondary case model, we found that force of infection was highly heterogeneous between jurisdictions and over time within jurisdictions, ranging from 0.00008 (95% CrI: 0.00002–0.0004) in Florida to 0.08 (95% CrI: 0.044–0.14) in American Samoa during the 2010–2019 period. For early 2020, we estimated that seropositivity in 10 year-olds ranged from 0.09% (0.02%–0.54%) in Florida to 56.3% (43.7%–69.3%) in American Samoa. In the absence of serological data, age-specific case notification data collected through routine surveillance combined with mathematical modeling are powerful tools to monitor arbovirus circulation, estimate the level of population immunity, and design dengue management strategies.

60 APPLIED LIFE SCIENCES↗

Photometry of the 4686 A emission line of gamma(2) Velorum from the South Pole

An automated optical telescope located at the Amundsen-Scott South Pole station on Antarctica, has been used to obtain more than 78 h of photometry of the He II emission line (4686 A) of the spectroscopic binary gamma(2) Velorum. These data were obtained on seven different days during the 1987 austral winter; the longest continuous run spans 19 h. Two independent period search techniques have been used to search for periodic behavior in the strength of the He II emission line of this Wolf-Rayet star. They are: (1) power spectrum analysis and (2) a first-order sine function fit to the data using least squares. Various multiplicities of a period on the order of 1.3 h with amplitudes of a few percent are found in most of these data. According to recent theoretical models of Wolf-Rayet stars, fluctuations in the He II emission line may indicate vibrational instability in gamma(2) Vel. These pulsations may, in turn, give rise to shocks which propagate outward and which may provide the necessary conditions for periodic changes in the state of a given region of the atmosphere to occur.

Taylor, Maryjane↗

The Construction of Curves and Surfaces Using Numerical Optimization Techniques

Numerical optimization techniques are playing an increasing role in curve and surface construction. Often difficult problems in curve and surface construction, especially when some aspect of shape control is involved, can be phrased as a constrained optimization problem. Four such classes of problems are explored: parametric curve fitting with non-linear shape constraints; explicit surface fitting with linear shape constraints; surface fitting to scattered data giving rise to ill-posed problems; finally, variable knot problems. In each of these problems there is a nonlinear aspect: either the shape of the curve or surface is important for manufacturing or engineering reasons or the shape affects the convergence of numerical algorithms which use the curve or surface or the placement of knots affects the accuracy of the fits. In all cases the class of functions used is that of parametric spline curves and tensor or direct product spline surfaces. The reason for choosing this class is that splines provide flexible models that are easily evaluated and stored. Furthermore, the B-spline representation of splines leads to convenient expressions for shape control over regions.

Ferguson, D. R.↗

Amplitude analysis of the ${D}_s^{+}$ → π - π + π + decay

A Dalitz plot analysis of the ${D}_s^{+}$ → π - π + π + decay is presented. The analysis is based on proton-proton collision data recorded by the LHCb experiment at a centre of-mass energy of 8 TeV and corresponding to an integrated luminosity of 1.5 fb -1 . The resonant structure of the decay is obtained using a quasi-model-independent partial-wave analysis, in which the π + π - S-wave amplitude is parameterised as a generic complex function determined by a fit to the data. The S-wave component is found to be dominant, followed by the contribution from spin-2 resonances and a small contribution from spin-1 resonances. The latter includes the first observation of the ${D}_s^{+}$ → ω(782)π + channel in the ${D}_s^{+}$ → π - π + π + decay. The resonant structures of the ${D}_s^{+}$ → π - π + π + and ${D}^{+}$ → π - π + π + decays are compared, providing information about the mechanisms for the hadron formation in these decays.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Sulfur Oxidation and Contrail Precursor Chemistry

Sulfuric acid (H2SO4), formed in commercial aircraft operations via Fuel-S (right arrow) SO2 + SO3 (right arrow) H2SO4, plays an important role in affecting the global climate change through atmospheric chemical reactions and radiative forcing. Measurement of the sulfur oxidation rates is critical to the understanding of the contrail formation. The principle reaction pathway is SO2 + O + M (right arrow) SO3 + M. Although there are many measurements for the rates of this reaction, it has never been measured in the temperature and pressure regime available to aircraft operation. In this investigation, a series of experiments were performed behind the reflected shock waves in a shock tube. OH radicals were produced in lean, shock heated SO2/H2/O2/Ar mixtures. The reaction progress was followed using OH absorption spectroscopy at 310 nm. The data were analyzed with the aid of computer modeling/simulation. The mean value of the rate coefficients of R21 determined is k(sub 21,0)/[M]= 3.9 x 10(exp 15) cm(sup 6) per square mole per second at T = 960 - 1150 K and rho = 16-30 micromole per cubic centimeter with uncertainty limits of plus or minus 30%. A non-Arrhenius fit to our data together with all existing data gives k(sub 21,0)/[M] = 1.3 x 10(exp 24) T (exp -2.5) exp(-2350 K/T) cm(sup 6) per square mole per second at T = 300 - 2500 K with the same uncertainty limits given above. The calculated conversion of S(IV)(SO2) to S(VI) (SO3 + H2SO4) was about 2% in our experimental conditions.

DeWitt, Kenneth↗

Micrometeoroid Impacts and Optical Surface Scatter in Space Environment

This presentation will discuss the results of an attempt to use laboratory test data and empirically derived models to quantify the degree of surface damage and associated light scattering that might be expected from hypervelocity particle impacts in space environment. Micrometeoroid sizes are predicted to be predominantly in the mass range 10(exp -5) g or less, with most having diameters near 1 micrometer, but some larger than 120 micrometers, with velocities near 20 km/s. In a laboratory test, latex (p = 1.1 g/cu cm) and iron (7.9 g/cu cm) particles with sizes ranging from 0.75 micrometer to 1.60 micrometer and with velocities ranging from 2.0 km/s to 18.5 km/s, were shot at a dielectric coated gold mirror. Scanning electron and atomic force microscopy were used to measure crater dimensions that were then associated with particle impact energies. These data were then fitted to empirical models derived from solar cell and other spacecraft surface components returned from orbit, as well as studies of impact craters studied on glassy materials returned from the lunar surface, to establish a link between particle energy and impact crater dimension. Published estimates of the Martian and Saturnian meteoroid environments were used as the sources of particle flux estimates. From these data, an estimate of total expected damaged area was computed and this result produced an estimate of expected surface scatter from the modeled environment.

Heaney, James B.↗

Focus history of the Hubble Space Telescope: Launch to May 1993

Since the launch of the Hubble Space Telescope (HST) the secondary mirror of the telescope has been moved several times in order to collimate the telescope and also to define a position of best focus. In addition to these moves the focus position changes over time because of water desorption by the graphite epoxy in the metering truss. The authors report here the focus history of the telescope based on a knowledge of the mirror moves made and an analysis of desorption monitoring data obtained by the Faint Object Camera (FOC) in the F/96 mode and of the routine calibration data obtained by the Wide Field and Planetary Cameras. Focus values are extracted using two different methods. In the first method the distance between the center of the point spread function (PSF) and the shadows of the pads supporting the HST primary mirror are related to the focus error. In the second method an analytical formula for the PSF with variable aberration coefficients is fitted to the data. Focus positions derived from the two methods show good agreement. The data show that a desorption of about 83 microns has taken place since Aug. 16, 1990. The desorption has clearly not leveled off as expected from the trend of the earlier data. Long term variations of the secondary mirror position of approximately 3-15 microns from the 'best' focus position have been observed. Variations of the order of 2-5 microns over an orbital period have also been noted. Focus changes resulting from secondary mirror movements greater than approximately 5 microns changes the point spread function significantly and makes deconvolution and quantitative measurements difficult.

Hasan, H.↗

Carbon Dioxide Adsorption on a 5A Zeolite Designed for CO2 Removal in Spacecraft Cabins

Carbon dioxide adsorption data were obtained for a 5A zeolite manufactured by AlliedSignal Inc. (Des Plaines, Illinois). The material is planned for use in the Carbon Dioxide Removal Assembly (CDRA) for U.S. elements of the International Space Station. The family of adsorption isotherms covers a temperature range of O to 250 C, and a pressure range of 0.001 to 800 torr. Coefficients of the Toth equation are fit to the data. Isosteric heats of adsorption are derived from the equilibrium loading data.

Mulloth, Lila M.↗

Cosmological constraints using Minkowski functionals from the first year data of the Hyper Suprime-Cam

We use Minkowski functionals to analyse weak lensing convergence maps from the first-year data release of the Subaru Hyper Suprime-Cam (HSC-Y1) survey. Minkowski functionals provide a description of the morphological properties of a field, capturing the non-Gaussian features of the Universe matter-density distribution. Using simulated catalogues that reproduce survey conditions and encode cosmological information, we emulate Minkowski functionals predictions across a range of cosmological parameters to derive the best-fit from the data. By applying multiple scales cuts, we rigorously mitigate systematic effects, including baryonic feedback and intrinsic alignments. From the analysis, combining constraints of the angular power spectrum and Minkowski functionals, we obtain S8≡σ8Ωm/0.3=0.808−0.046+0.033 and Ωm=0.293−0.043+0.157⁠. These results represent a 40 per cent improvement on the S8 constraints compared to using power spectrum only. Minkowski functionals results are consistent with other two-point, and higher order statistics constraints using the same data, being in agreement with CMB results from the Planck S8 measurements. Our study demonstrates the power of Minkowski functionals beyond two-point statistics to constrain and break the degeneracy between Ωm and σ8⁠.

79 ASTRONOMY AND ASTROPHYSICS↗

DESI 2024: reconstructing dark energy using crossing statistics with DESI DR1 BAO data

Here, we implement Crossing Statistics to reconstruct in a model-agnostic manner the expansion history of the universe and properties of dark energy, using DESI Data Release 1 (DR1) BAO data in combination with one of three different supernova compilations (PantheonPlus, Union3, and DES-SN5YR) and Planck CMB observations. Our results hint towards an evolving and emergent dark energy behaviour, with negligible presence of dark energy at z ≳ 1, at varying significance depending on data sets combined. In all these reconstructions, the cosmological constant lies outside the 95% confidence intervals for some redshift ranges. This dark energy behaviour, reconstructed using Crossing Statistics, is in agreement with results from the conventional w 0 –w a dark energy equation of state parametrization reported in the DESI Key cosmology paper. Our results add an extensive class of model-agnostic reconstructions with acceptable fits to the data, including models where cosmic acceleration slows down at low redshifts. We also report constraints on H 0 r d from our model-agnostic analysis, independent of the pre-recombination physics.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurements and Modeling of Soot Formation and Radiation in Microgravity Jet Diffusion Flames

This is a computational and experimental study for soot formation and radiative heat transfer in jet diffusion flames under normal gravity (1-g) and microgravity (0-g) conditions. Instantaneous soot volume fraction maps are measured using a full-field imaging absorption technique developed by the authors. A compact, self-contained drop rig is used for microgravity experiments in the 2.2-second drop tower facility at NASA Lewis Research Center. On modeling, we have coupled flame structure and soot formation models with detailed radiation transfer calculations. Favre-averaged boundary layer equations with a k-e-g turbulence model are used to predict the flow field, and a conserved scalar approach with an assumed Beta-pdf are used to predict gaseous species mole fraction. Scalar transport equations are used to describe soot volume fraction and number density distributions, with formation and oxidation terms modeled by one-step rate equations and thermophoretic effects included. An energy equation is included to couple flame structure and radiation analyses through iterations, neglecting turbulence-radiation interactions. The YIX solution for a finite cylindrical enclosure is used for radiative heat transfer calculations. The spectral absorption coefficient for soot aggregates is calculated from the Rayleigh solution using complex refractive index data from a Drude- Lorentz model. The exponential-wide-band model is used to calculate the spectral absorption coefficient for H20 and C02. It is shown that when compared to results from true spectral integration, the Rosseland mean absorption coefficient can provide reasonably accurate predictions for the type of flames studied. The soot formation model proposed by Moss, Syed, and Stewart seems to produce better fits to experimental data and more physically sound than the simpler model by Khan et al. Predicted soot volume fraction and temperature results agree well with published data for a normal gravity co-flow laminar flames and turbulent jet flames. Predicted soot volume fraction results also agree with our data for 1-g and 0-g laminar jet names as well as 1-g turbulent jet flames.

Ku, Jerry C.↗

Hot-spot model for inertial confinement fusion implosions with an applied magnetic field

Imposing a magnetic field on inertial confinement fusion implosions magnetizes the electrons in the compressed fuel; this suppresses thermal losses, which increases temperature and fusion yield. Indirect-drive experiments at the National Ignition Facility with 12 and 26 T applied magnetic fields demonstrate up to 40% increase in temperature, 3× increase in fusion yield, and indicate that magnetization alters the radial temperature profile [Moody et al., Phys. Rev. Lett. 129, 195002 (2022); Lahmann et al., APS DPP (2022)]. In this work, we develop a semi-analytic hot-spot model, which accounts for the two-dimensional (2D) Braginskii anisotropic heat flow due to an applied axial magnetic field. First, we show that hot-spot magnetization alters the radial temperature profile, increasing the central peakedness, which is most pronounced for moderately magnetized implosions (with 8–14 T applied field), compared to both unmagnetized (with no applied field) and highly magnetized (with 26 T or higher applied field) implosions. This model explains the trend in the experimental data, which finds a similarly altered temperature profile in the 12 T experiment. Next, we derive the hot-spot model for gas-filled (Symcap) implosions, accounting for the effects of magnetization on the thermal conduction and in changing the radial temperature (and density) profiles. Using this model, we compute predicted central temperature amplification and yield enhancement scaling with the applied magnetic field. The central temperature fits the experimental data accurately, and the discrepancy in the yield suggests a systematic (independent of applied field) degradation, such as mix, and additional degradation in the reference unmagnetized shot, such as reduced laser drive, increased implosion asymmetry, or the magnetic field suppressing ablator mixing into the hot-spot.

Alpha particles↗

O2 reduction at the IFC orbiter fuel cell O2 electrode

O2 reduction Tafel data were obtained for the IFC Orbiter fuel cell O2 electrode (Au-10 percent Pt catalyst) at temperatures between 24 and 81 C. BET measurements gave an electrode surface area of about 2040 sq cm per sq cm of geometric area. The Tafel data could be fitted to three straight line regions. For current densities less than 0.001 A/sq cm, the slope was essentially independent of temperature with a value of about 0.032 V/decade. Above 0.001 A/sq cm, the two regions, designated in the present study as the 0.04 and 0.12 V/decate regions, were temperature dependent. The apparent energies of activation for these two regions were about 9.3 and 6.5 kcal/mol, respectively. Tafel data (1 atmosphere O2) were extrapolated to 120 C for predicting changes in overpotential with increasing temperature. A mechanism is presented for O2 reduction.

Fielder, William L.↗

The O2 reduction at the IFC modified O2 fuel cell electrode

The International Fuel Corporation (IFC) state of the art (SOA) O2 electrode (Au-10 percent Pt electrocatalyst by weight) is currently being used in the alkaline H2-O2 fuel cell in the NASA Space Shuttle. Recently, IFC modified O2 electrode, as a possible replacement for the SOA electrode. In the present study, O2 reduction data were obtained for the modified electrode at temperatures between 23.3 and 91.7 C. BET measurements gave an electrode BET surface area of about 2070 sq. cm/sq. cm of geometric surface area. The Tafel data could be fitted to two straight line regions. The slope for the lower region, designated as the 0.04 V/decade region, was temperature dependent, and the transfer coefficient was about 1.5. The 'apparent' energy of activation for this region was about 19 kcal/mol. An O2 reduction mechanism for this 0.04 region is presented. In the upper region, designated as the 0.08 V/decade region, diffusion may be the controlling process. Tafel data are presented to illustrate the increase in performance with increasing temperature.

Fielder, William L.↗

Development of a Generic Creep-Fatigue Life Prediction Model

The objective of this research proposal is to further compile creep-fatigue data of steel alloys and superalloys used in military aircraft engines and/or rocket engines and to develop a statistical multivariate equation. The newly derived model will be a probabilistic fit to all the data compiled from various sources. Attempts will be made to procure the creep-fatigue data from NASA Glenn Research Center and other sources to further develop life prediction models for specific alloy groups. In a previous effort [1-3], a bank of creep-fatigue data has been compiled and tabulated under a range of known test parameters. These test parameters are called independent variables, namely; total strain range, strain rate, hold time, and temperature. The present research attempts to use these variables to develop a multivariate equation, which will be a probabilistic equation fitting a large database. The data predicted by the new model will be analyzed using the normal distribution fits, the closer the predicted lives are with the experimental lives (normal line 1 to 1 fit) the better the prediction. This will be evaluated in terms of a coefficient of correlation, R 2 as well. A multivariate equation developed earlier [3] has the following form, where S, R, T, and H have specific meaning discussed later.

Goswami, Tarun↗