Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model 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 757 records · Page 42

Relativistic Iron Emission and Disk Reflection in Galactic Microquasar XTE J1748-288

We report evidence for an Fe K(alpha) fluorescence line feature and disk reflection in the very high, high-, and low-state X-ray spectra of the Galactic microquasar XTE J1748-288 during its 1998 June outburst. Spectral analyses are made on data gathered throughout the outburst by the Rossi X-Ray Timing Explorer Proportional Counter Array. Gaussian line, relativistic disk emission line, and ionized disk reflection models are fitted to the data. In the very high state the line profile appears strongly redshifted, consistent with disk emission from the innermost stable orbits around a maximally rotating Kerr black hole. In the high state the line profile is less redshifted and increasingly prominent. The low-state line profile is very strong (approx. 0.5 keV equivalent width) and centered at 6.7 +/- 0.10 keV; disk line emission model fits indicate that the inner edge of the disk fluctuates between approx. 20Rg and approx. 100Rg in this state. The disk reflection fraction is traced through the outburst; reflection from an ionized disk is preferred in the very high and high states, and reflection from a relatively neutral disk is preferred in the low state. We discuss the implications of our findings for the binary system dynamics and accretion flow geometry in XTE J1748-288.

Miller, J. M.↗

The energy spectra of solar flare electrons

A survey of 50 electron energy spectra from .1 to 100 MeV originating from solar flares was made by the combination of data from two spectrometers onboard the International Sun Earth Explorer-3 spacecraft. The observed spectral shapes of flare events can be divided into two classes through the criteria of fit to an acceleration model. This standard two step acceleration model, which fits the spectral shape of the first class of flares, involves an impulsive step that accelerates particles up to 100 keV and a second step that further accelerates these particles up to 100 MeV by a single shock. This fit fails for the second class of flares that can be characterized as having excessively hard spectra above 1 MeV relative to the predictions of the model. Correlations with soft X-ray and meter radio observations imply that the acceleration of the high energy particles in the second class of flares is dominated by the impulsive phase of the flares.

Evenson, P. A.↗

Training material models using gradient descent algorithms

High temperature design requires accurate constitutive models to describe material inelastic deformation and failure behavior. Oftentimes, calibrating accurate models devolves into the problem of fitting the model parameters against experimental test data. Here, we present the pyopmat package, an open source framework for calibrating constitutive models against experiment data subjected to various loading conditions using machine learning techniques. The package calculates the exact gradient of the model response with respect to the parameters using a combination of automatic differentiation and the adjoint method. Given this exact gradient, we compare the performance of several gradient-based optimization techniques in fitting realistic constitutive models against data. Here, we demonstrate the efficiency and accuracy of our package through example problems using both synthetic data, generated using known parameter sets, under monotonic and cyclic loading conditions and also with an example applying the techniques developed here to actual high temperature creep-fatigue test data.

36 MATERIALS SCIENCE↗

ML Pipeline (Machine Learning Pipeline) [SWR-22-35]

The Machine Learning Pipeline (ML Pipeline) allows the user to fit a model to predict a dependent variable, y, based on a feature set, X. ML Pipeline gives the ability to automatically generate features with its Feature Engineering function, automatically select the most important features with its Feature Selection function, fit the model with its Fit function, and evaluate the model with its Evaluate function. ML Pipeline's functionality can not only help the user predict values for the desired variable it also gives the user a better understanding of which features are important and how effective the model is. ML Pipeline accomplishes all this, while remaining simple to use and interpret.

Schiek, Andrew↗

The Atacama Cosmology Telescope: Cosmological Parameters from the 2008 Power Spectrum

We present cosmological parameters derived from the angular power spectrum of the cosmic microwave background (CMB) radiation observed at 148 GHz and 218 GHz over 296 deg(exp 2) with the Atacama Cosmology Telescope (ACT) during its 2008 season. ACT measures fluctuations at scales 500 < l < 10,000. We fit a model for the lensed CMB, Sunyaev-Zel'dovich (SZ), and foreground contribution to the 148 GHz and 218 GHz power spectra, including thermal and kinetic SZ, Poisson power from radio and infrared point sources, and clustered power from infrared point sources. At l = 3000, about half the power at 148 GHz comes from primary CMB after masking bright radio sources. The power from thermal and kinetic SZ is estimated to be Beta(sub 3000) is identical to 6.8 +/- 2.9 mu K (exp 2), where Beta (sub l) is identical to l(l + 1) C(sub l)/2pi. The IR Poisson power at 148 GHz is Bewta(sub 3000) 7.8 +/- 0.7 muK(exp 2) (C(sub l) = 5.5 +/- 0.5 nK(exp 2)), and a clustered IR component is required with Beta (sub 3000) = 4.6 +/- 0.9 muK(exp 2), assuming an analytic model for its power spectrum shape. At 218 GHz only about 15% of the power, approximately 27 mu K(exp 2), is CMB anisotropy at l = 3000. The remaining 85% is attributed to IR sources (approximately 50% Poisson and 35% clustered), with spectral index alpha = 3.69 +/- 0.14 for flux scaling as S(nu) varies as nu(sup alpha). We estimate primary cosmological parameters from the less contaminated 148 GHz spectrum, marginalizing over SZ and source power. The ACDM cosmological model is a good fit to the data (chi square/dof = 29/46), and ACDM parameters estimated from ACT+Wilkinson Microwave Anisotropy Probe (WMAP) are consistent with the seven-year WMAP limits, with scale invariant n(sub s) = 1 excluded at 99.7% confidence level (CL) (3 sigma). A model with no CMB lensing is disfavored at 2.8 sigma. By measuring the third to seventh acoustic peaks, and probing the Silk damping regime, the ACT data improve limits on cosmological parameters that affect the small-scale CMB power. The ACT data combined with WMAP give a 6 sigma detection of primordial helium, with Y(sub p) = 0.313 +/- 0.044, and a 4 sigma detection of relativistic species, assumed to be neutrinos, with N(sub eff) = 5.3 +/- 1.3 (4.6 +/- 0.8 with BAO+H(sub 0) data). From the CMB alone the running of the spectral index is constrained to be d(sub s) / d ln k = -0,034 +/- 0,018, the limit on the tensor-to-scalar ratio is r < 0,25 (95% CL), and the possible contribution of Nambu cosmic strings to the power spectrum is constrained to string tension G(sub mu) < 1.6 x 10(exp -7) (95% CL),

Dunkley, J.↗

Fourier-like Thermal Relaxation of Nanoscale Explosive Hot Spots

Hot spots are local regions of high temperature that are widely considered to govern explosive initiation. Hot spot dynamics rests on a delicate balance between heat generation due to chemical reactions and heat loss through thermal conduction, making accurate determinations of the conductivity under extreme conditions a key component of predictive explosive models. We develop here an approach to directly determine the thermal transport properties of explosive hot spots with realistic initial structures through a combination of molecular dynamics (MD) and diffusive heat equation (HEq) modeling. Effective thermal conductivity values are determined by fitting HEq models to MD predictions of long timescale hot spot relaxation. The approach is applied to model hot spots in the molecular crystalline explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) for a range of shock strengths and two limiting cases for impact orientation. Isotropic and anisotropic HEq models yield similar results, despite TATB exhibiting some of the largest and most anisotropic thermal conductivity values for explosive near normal conditions. The conductivity is found to be a strong function of density, which parametrically captures dependence on temperature, pressure, and material state. The associated root-mean-square errors of the fitted HEq models are approximately 5% of MD predicted final equilibrium temperatures. Here, the conductivity values determined here for TATB hot spots are considerably larger than those used in a prior hot spot criticality study, which may significantly impact predictions for critical hot spot sizes. The approach provides a convenient foundation for determining the effective thermal conductivity for hot spot problems in other explosives and directly yields information on reasonable approximations that might be taken in higher-level models for those materials.

36 MATERIALS SCIENCE↗

Backus Effect on a Perpendicular Errors in Harmonic Models of Real vs. Synthetic Data

Measurements of geomagnetic scalar intensity on a thin spherical shell alone are not enough to separate internal from external source fields; moreover, such scalar data are not enough for accurate modeling of the vector field from internal sources because of unmodeled fields and small data errors. Spherical harmonic models of the geomagnetic potential fitted to scalar data alone therefore suffer from well-understood Backus effect and perpendicular errors. Curiously, errors in some models of simulated 'data' are very much less than those in models of real data. We analyze select Magsat vector and scalar measurements separately to illustrate Backus effect and perpendicular errors in models of real scalar data. By using a model to synthesize 'data' at the observation points, and by adding various types of 'noise', we illustrate such errors in models of synthetic 'data'. Perpendicular errors prove quite sensitive to the maximum degree in the spherical harmonic expansion of the potential field model fitted to the scalar data. Small errors in models of synthetic 'data' are found to be an artifact of matched truncation levels. For example, consider scalar synthetic 'data' computed from a degree 14 model. A degree 14 model fitted to such synthetic 'data' yields negligible error, but amplifies 4 nT (rmss) added noise into a 60 nT error (rmss); however, a degree 12 model fitted to the noisy 'data' suffers a 492 nT error (rmms through degree 12). Geomagnetic measurements remain unaware of model truncation, so the small errors indicated by some simulations cannot be realized in practice. Errors in models fitted to scalar data alone approach 1000 nT (rmss) and several thousand nT (maximum).

Voorhies, C. V.↗

Glass Property-Composition Models for Support of Hanford WTP LAW Facility Operation

Current plans for the River Protection Project envision starting to vitrifying low-activity waste (LAW) by 2023 using a Direct Feed Low-Activity Waste (DFLAW) approach and subsequently using a full-pretreatment approach. The Hanford Tank Waste Treatment and Immobilization Plant (WTP) LAW Facility will be operated and controlled using a LAW glass formulation algorithm (GFA), which requires several inputs based on research and development results. LAW glass property-composition models for several product quality and processing properties are key inputs for the LAW GFA. It is envisioned that the preliminary LAW GFA discussed by Kim and Vienna (2012) will be used for commissioning and initial radioactive operations of the WTP LAW Facility under Bechtel National, Inc. using the DFLAW approach. Then, an updated LAW GFA will be developed for implementation by the WTP operating contractor that takes over after WTP LAW Facility commissioning. This report documents the enhanced LAW glass property-composition models developed for use in the updated LAW GFA. The properties for which models were developed include Product Consistency Test (PCT) response, Vapor Hydration Test (VHT) response, viscosity at 1150 °C, electrical conductivity at 1150 °C, melter SO 3 tolerance at 1150 °C, and K-3 refractory corrosion at 1208 °C. Table S.1 lists the tables in this report that contain the recommended models for each of these properties. The model types recommended include partial quadratic mixture (PQM) models for viscosity, electrical conductivity, melter SO 3 tolerance and K-3 corrosion, bias corrected PQM model (bcPQM) for PCT, and logistic PQM model for VHT. The fits of model and validation subsets were found to be well predicted by the recommended models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Glass Property-Composition Models for Support of Hanford WTP LAW Facility Operation

Current plans for the River Protection Project envision starting to vitrifying low-activity waste (LAW) by 2023 using a Direct Feed Low-Activity Waste (DFLAW) approach and subsequently using a full-pretreatment approach. The Hanford Tank Waste Treatment and Immobilization Plant (WTP) LAW Facility will be operated and controlled using a LAW glass formulation algorithm (GFA), which requires several inputs based on research and development results. LAW glass property-composition models for several product quality and processing properties are key inputs for the LAW GFA. It is envisioned that the preliminary LAW GFA discussed by Kim and Vienna (2012) will be used for commissioning and initial radioactive operations of the WTP LAW Facility under Bechtel National, Inc. using the DFLAW approach. Then, an updated LAW GFA will be developed for implementation by the WTP operating contractor that takes over after WTP LAW Facility commissioning. This report documents the enhanced LAW glass property-composition models developed for use in the updated LAW GFA. The properties for which models were developed include Product Consistency Test (PCT) response, Vapor Hydration Test (VHT) response, viscosity at 1150 °C, electrical conductivity at 1150 °C, melter SO3 tolerance at 1150 °C, and K-3 refractory corrosion at 1208 °C. Table S.1 lists the tables in this report that contain the recommended models for each of these properties. The model types recommended include partial quadratic mixture (PQM) models for viscosity, electrical conductivity, melter SO 3 tolerance and K-3 corrosion, bias corrected PQM model (bcPQM) for PCT, and logistic PQM model for VHT. The fits of model and validation subsets were found to be well predicted by the recommended models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The magnetic anomaly of the Ivreazone

A magnetic field survey was made in the Ivreazone in 1969/70. The results were: significant anomaly of the vertical intensity is found. It follows the basic main part of the Ivrea-Verbano zone and continues to the south. The width of the anomaly is about 10 km, the maximum measures about +800 gamma. The model interpretation shows that possibly the anomaly belongs to an amphibolitic body, which in connection with the Ivrea-body was found by deep seismic sounding. Therefore, the magnetic anomaly provides further evidence for the conception that the Ivrea-body has to be regarded as a chip of earthmantle material pushed upward by tectonic processes.

Vertical Intensity of Magnetic Field↗

Images of the earth's aurora and geocorona from the Dynamics Explorer mission

Several series of auroral and geocoronal images from th Dynamics Explorer Mission are presented and analyzed. It is found that the motion of the transpolar arc of a theta aurora is correlated with the y-component of the IMF; the arc motion is in the general direction of the y-component. In addition, a sequence of global images of a small auroral substorm shows the initial development of intense luminosities in a relatively small spatial region or 'bright spot' in the premidnight sector of the auroral oval and a subsequent appearance of an expanding area of lesser intensities at low latitudes and contiguous to the midnight boundary of the bright spot. Finally, a series of images of the geocorona in scattered solar Ly-alpha emissions is used to obtain a best-fit spherical model of atomic hydrogen densities in the earth's exosphere; a Chamberlain model provides an adequate fit to radial distances of 4.5 earth-radii, beyond which an exponential fit is used. The geocoronal tail is detected as an asymmetric increase in scattered Ly-alpha intensities in the antisolar direction.

Frank, L. A.↗

Protein crystal nucleation kinetics using relative light-scattering techniques: Center director's discretionary fund

Light-scattering intensity measurements are a sensitive method for following changes in the hydrodynamic radius of particles in solution. The approach used in this report utilizes the light-scattering intensity ratios of a polydisperse to a monodisperse system. By numerically modeling the process, and fitting the model curves to the data, estimates have been obtained for the dimerization equilibrium constant, the dimer + dimer yields tetramer equilibrium constant, and the association rate constant for the dimerization process.

Pusey, Marc Lee↗

Detection of bearing failure in mechanical devices using neural networks

We present a novel time-domain method for the detection of faulty bearings that has direct applicability to monitoring the health of the turbo pumps on the Space Shuttle Main Engine. A feed-forward neural network was trained to detect modelled roller bearing faults on the basis of the periodicity of impact pulse trains. The network's performance was dependent upon the number of pulses in the network's input window and the signal-to-noise ratio of the input signal. To test the model's validity, we fit the model's parameters to an actual vibration signal generated by a faulty roller element bearing and applied the network trained on this model to detect faults in actual vibration data. When this network was tested on the actual vibration data, it correctly identified the vibration signal as a fault condition 76 percent of the time.

Burne, Richard A.↗

Passive Microwave Remote Sensing of Rainfall Considering the Effects of Wind and Nonprecipitating Clouds

It has long been shown both in theory and in observation that emission from rain drops in a raining cloud results in upwelling brightness temperature above that caused by the sea surface alone. High brightness temperatures at microwave frequencies (e.g. 37 and 19 GHz) have usually been quantitatively associated with rainfall using physical or statistical models. By comparing concurrent special sensor microwave/imager and radar data, however, we noticed many cases where there is no appreciable rainfall in a field of view (FOV) which exhibits high brightness temperature (T(sub B)) at 37 and 19 GHz. On the basis of calculations and past literature it is shown that such high brightness temperatures can be caused by nonprecipitating clouds and by wind. The effect of the wind is to create wave and high-emissivity foam on the sea surface. A model is developed to relate T(sub B) to the fractional coverage of rain, f, within a FOV. The parameters of the model are calibrated by fitting the model to the observed brightness temperature and fractional rain coverage data. The critical parameter of the model, T(sub B min.), which is the threshold brightness temperature for the presence of rain, depends on the strength of the storm. The strength of the storm is characterized by the fraction of the FOVs within a large area that have T(sub B) higher than 240 K, which is readily obtainable from satellite data alone. The instantaneous FOV rain rate R can then be obtained through the f approximately R relationship which is empirically derived using radar data. An algorithm has been proposed based on the T(sub B) approximately f and f approximately R relationship. Application of the algorithm to TOGA-COARE and Darwin storms results in reasonable instantaneous FOV rain estimate. When averaged over the entire radar scan, a more accurate and unbiased areal rain estimate can be achieved.

Qihang Li.↗

Gamma-Ray Observations of the Supernova Remnant RX J0852.0-4622 with the Fermi Large Area Telescope

We report on gamma-ray observations of the supernova remnant (SNR) RX J0852.04622 with the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. In the Fermi-LAT data, we find a spatially extended source at the location of the SNR. The extension is consistent with the SNR size seen in other wavelengths such as X-rays and TeV gamma rays, leading to the identification of the gamma-ray source with the SNR. The spectrum is well described as a power law with a photon index of = 1.85 0.06 (stat)+0.18 0.19 (sys), which smoothly connects to the H.E.S.S. spectrum in the TeV energy band. We discuss the gamma-ray emission mechanism based on multiwavelength data. The broadband data can be fit well by a model in which the gamma rays are of hadronic origin. We also consider a scenario with inverse Compton scattering of electrons as the emission mechanism of the gamma rays. Although the leptonic model predicts a harder spectrum in the Fermi-LAT energy range, the model can fit the data considering the statistical and systematic errors.

Tanaka, T.↗

NASA 5.2%-Scale Semispan High Lift Common Research Model Wind Tunnel Test at the 14- by 22-Foot Subsonic Tunnel

A wind tunnel test of the NASA 5.2%-scale Semispan High Lift Common Research Model (CRM-HL) was conducted at the 14- by 22-Foot Subsonic Tunnel (14x22) in the fall of 2022. The main purpose of the test was a model checkout before shipping it overseas. The model was shipped to Germany in the spring of 2023, and the team wanted to be sure that instrumentation in the model was working, all of the model parts fit together well, and model changes could be done smoothly, as designed while installed in a wind tunnel. Force and moment and pressure data were collected during the test, and while these data were a secondary objective, they can be used for for tunnel-to-tunnel verification of the model, and for database supplementation for the NASA CRM-HL ecosystem. The test was run from August 15, 2022 through September 8, 2022. Eight configurations were tested at Mach 0.2 and 0.26, over an alpha range from -4deg to +24deg. All of the configurations that were tested accomplished testing of all model parts. This test campaign was a unique opportunity to test a complicated model design in a wind tunnel that is less expensive to test in and has easier model access than the tunnel that it was designed to run in. The test was invaluable in providing lessons learned about building up and running the model, and the model was returned to the model manufacturer from December 2022 to February 2023 for updates to make model changes easier, to make design fixes on the model, and to do clean up work on the model.

CRM-HL↗

NASA 5.2%-Scale Semispan High Lift Common Research Model Wind Tunnel Test at the 14- by 22-Foot Subsonic Tunnel

A wind tunnel test of the NASA 5.2%-scale Semispan High Lift Common Research Model (CRM-HL) was conducted at the 14- by 22-Foot Subsonic Tunnel (14x22) in the fall of 2022. The main purpose of the test was a model checkout before shipping it overseas. The model was shipped to Germany in the spring of 2023, and the team wanted to be sure that instrumentation in the model was working, all of the model parts fit together well, and model changes could be done smoothly, as designed while installed in a wind tunnel. Force and moment and pressure data were collected during the test, and while these data were a secondary objective, they can be used for for tunnel-to-tunnel verification of the model, and for database supplementation for the NASA CRM-HL ecosystem. The test was run from August 15, 2022 through September 8, 2022. Eight configurations were tested at Mach 0.2 and 0.26, over an alpha range from -4deg to +24deg. All of the configurations that were tested accomplished testing of all model parts. This test campaign was a unique opportunity to test a complicated model design in a wind tunnel that is less expensive to test in and has easier model access than the tunnel that it was designed to run in. The test was invaluable in providing lessons learned about building up and running the model, and the model was returned to the model manufacturer from December 2022 to February 2023 for updates to make model changes easier, to make design fixes on the model, and to do clean up work on the model.

Semispan↗

Suomi NPP VIIRS Solar Diffuser Screen Transmittance Model and Its Applications

The visible infrared imaging radiometer suite on the Suomi National Polar-orbiting Partnership satellite calibrates its reflective solar bands through observations of a sunlit solar diffuser (SD) panel. Sunlight passes through a perforated plate, referred to as the SD screen, before reaching the SD. It is critical to know whether the SD screen transmittance measured prelaunch is accurate. Several factors such as misalignments of the SD panel and the measurement apparatus could lead to errors in the measured transmittance and thus adversely impact on-orbit calibration quality through the SD. We develop a mathematical model to describe the transmittance as a function of the angles that incident light makes with the SD screen, and apply the model to fit the prelaunch measured transmittance. The results reveal that the model does not reproduce the measured transmittance unless the size of the apertures in the SD screen is quite different from the design value. We attribute the difference to the orientation alignment errors for the SD panel and the measurement apparatus. We model the alignment errors and apply our transmittance model to fit the prelaunch transmittance to retrieve the "true" transmittance. To use this model correctly, we also examine the finite source size effect on the transmittance. Furthermore, we compare the product of the retrieved "true" transmittance and the prelaunch SD bidirectional reflectance distribution function (BRDF) value to the value derived from on-orbit data to determine whether the prelaunch SD BRDF value is relatively accurate. The model is significant in that it can evaluate whether the SD screen transmittance measured prelaunch is accurate and help retrieve the true transmittance from the transmittance with measurement errors, consequently resulting in a more accurate sensor data product by the same amount.

sensors↗