Search NASA⌕ Search

SEARCH · Search NASA

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 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Impacts of control strategies, the Great Recession and weekday variations on NO 2 columns above North American cities

The Ozone Monitoring Instrument (OMI) has been estimating NO2 columns from space for over 10 years, and these have been used to estimate emissions and emission trends for point and area sources all over the world. In this study we evaluate the trends in NO2 columns over 54 cities in the USA and Canada to identify the long term trends due to air quality policies, the impact of the Great Recession, and the weekday-weekend effect. A multiple linear regression model is used to fit annual, seasonal and weekly factors for individual swath retrievals along with the impact of temperature, wind speed and pixel size. For most cities, the correlation coefficients of the model fit ranges from 0.47 to 0.76. There have been strong reductions in NO2 columns, with annual decreases of up to 7% per year in most cities. During the years of the Great Recession, NO2 columns were as much as 30% lower than they would have been had they followed the linear annual trend. The analysis yielded insights into the timing of the reductions, with some cities in the northwest and in the east experiencing reductions in 2008 already, and most areas back to where they would have been based on the uniform trend by 2011. The analysis also finds that reductions in columns during the weekend vary significantly from city to city, with a range in reductions of 10%-30% on Saturdays, and 20%-50% on Sundays.

de Foy, Benjamin (ORCID:0000000341509922)↗
DOE OSTI2025

Information theory optimization of signals from small-angle scattering measurements

Small-angle X-ray scattering (SAXS) of particles in solution informs on the conformational states and assemblies of biological macromolecules (bioSAXS) outside of cryo- and solid-state conditions. In bioSAXS, the SAXS measurement under dilute conditions is resolution limited, and through an inverse Fourier transform, the measured SAXS intensities directly relate to the physical space occupied by the particles via the P (r)-distribution. Yet, this inverse transform of SAXS data has been historically cast as an ill-posed, ill-conditioned problem requiring an indirect approach. Here, we show that through the applications of matrix and information theories, the inverse transform of SAXS intensity data is a well-conditioned problem. The so-called ill-conditioning of the inverse problem is directly related to the Shannon number. By exploiting the oversampling enabled by modern detectors, a direct inverse Fourier transform of the SAXS data is possible, provided the recovered information does not exceed the Shannon number. The Shannon limit corresponds to the maximum number of significant singular values that can be recovered in a SAXS experiment, suggesting this relationship is a fundamental property of band-limited inverse integral transform problems. This correspondence reduces the complexity of the inverse problem to the Shannon limit and maximum dimension. We propose a hybrid scoring function using an information theory framework that assesses both the quality of the model-data fit as well as the quality of the recovered P (r)-distribution. The hybrid score utilizes the Akaike information criteria and Durbin-Watson statistic that considers parameter-model complexity, i.e., degrees of freedom, and the randomness of the model-data residuals. The described tests and findings extend the boundaries for bioSAXS by completing the information theory formalism initiated by Peter B. Moore to enable a quantitative measure of resolution in SAXS, robustly determine maximum dimension, and more precisely define the best parameter model appropriately representing the observed scattering data.

Rambo, Robert P. [Science and Technology Facilitie↗
DOE OSTI2014

A general methodology for inverse estimation of the elastic and anelastic properties of anisotropic open-cell porous materials—with application to a melamine foam

This paper proposes an inverse estimation method for the characterisation of the elastic and anelastic properties of the frame of anisotropic open-cell foams used for sound absorption. A model of viscoelasticity based on a fractional differential constitutive equation is used, leading to an augmented Hooke's law in the frequency domain, where the elastic and anelastic phenomena appear as distinctive terms in the stiffness matrix. The parameters of the model are nine orthotropic elastic moduli, three angles of orientation of the material principal directions and three parameters governing the anelastic frequency dependence. The inverse estimation consists in numerically fitting the model on a set of transfer functions extracted from a sample of material. The setup uses a seismic-mass measurement repeated in the three directions of space and is placed in a vacuum chamber in order to remove the air from the pores of the sample. The method allows to reconstruct the full frequency-dependent complex stiffness matrix of the frame of an anisotropic open-cell foam and in particular it provides the frequency of maximum energy dissipation by viscoelastic effects. The characterisation of a melamine foam sample is performed and the relation between the fractional-derivative model and other types of parameterisations of the augmented Hooke's law is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Helioseismology with long-range dark matter-baryon interactions

Assuming the existence of a primordial asymmetry in the dark sector, we study how long-range dark matter (DM)-baryon interactions, induced by the kinetic mixing of a new U(1) gauge boson and a photon, affect the evolution of the Sun and, in turn, the sound speed the profile obtained from helioseismology. Thanks to the explicit dependence on the exchanged momenta in the differential cross section (Rutherford-like scattering), we find that DM particles with a mass of ∼10 GeV, kinetic mixing parameter of the order of 10{sup –9}, and a mediator with a mass smaller than a few MeV improve the agreement between the best solar model and the helioseismic data without being excluded by direct detection experiments. In particular, the LUX detector will soon be able to either constrain or confirm our best-fit solar model in the presence of a dark sector with long-range interactions that reconcile helioseismology with thermal neutrino results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Sci—Fri PM: Topics — 06: The influence of regional dose sensitivity on salivary loss and recovery in the parotid gland

Purpose: The Quantitative Analyses of Normal Tissue Effects in the Clinic (QUANTEC 2010) survey of radiation dose-volume effects on salivary gland function has called for improved understanding of intragland dose sensitivity and the effectiveness of partial sparing in salivary glands. Regional dose susceptibility of sagittally- and coronally-sub-segmented parotid gland has been studied. Specifically, we examine whether individual consideration of sub-segments leads to improved prediction of xerostomia compared with whole parotid mean dose. Methods: Data from 102 patients treated for head-and-neck cancers at the BC Cancer Agency were used in this study. Whole mouth stimulated saliva was collected before (baseline), three months, and one year after cessation of radiotherapy. Organ volumes were contoured using treatment planning CT images and sub-segmented into regional portions. Both non-parametric (local regression) and parametric (mean dose exponential fitting) methods were employed. A bootstrap technique was used for reliability estimation and cross-comparison. Results: Salivary loss is described well using non-parametric and mean dose models. Parametric fits suggest a significant distinction in dose response between medial-lateral and anterior-posterior aspects of the parotid (p<0.01). Least-squares and least-median squares estimates differ significantly (p<0.00001), indicating fits may be skewed by noise or outliers. Salivary recovery exhibits a weakly arched dose response: the highest recovery is seen at intermediate doses. Conclusions: Salivary function loss is strongly dose dependent. In contrast no useful dose dependence was observed for function recovery. Regional dose dependence was observed, but may have resulted from a bias in dose distributions.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1972

Line-blanketed model stellar atmospheres applied to Sirius

The primary goal of this analysis is to determine whether the effects of atomic bound-bound transitions on stellar atmospheric structure can be represented well in models. The investigation is based on an approach which is called the method of artificial absorption edges. The method is described, developed, tested, and applied to the problem of fitting a model stellar atmosphere to Sirius. It is shown that the main features of the entire observed spectrum of Sirius can be reproduced to within the observational uncertainty by a blanketed flux-constant model with T sub eff = 9700 K and Log g = 4.26. The profile of H sub gamma is reproduced completely within the standard deviations of the measurements except near line center, where non-LTE effects are expected to be significant. The equivalent width of H sub gamma, the Paschen slope, the Balmer jump, and the absolute flux at 5550 A all agree with the observed values.

Fowler, J. W.↗
NASA NTRS1991

Surface scattering properties estimated from modeling ASAS multiple emission angle reflectance data over the lunar crater volcanic field, Nevada

ASAS multispectral data were acquired over Lunar Lake playa at 12 different lighting and viewing geometries. Atmospheric models and field-based data demonstrate that atmospheric contributions to ASAS data were relatively small and were thus ignored in analyses of surface photometric properties. ASAS data of an artificially roughened playa and cobble-strewn surfaces exhibit trends of increasing reflectance as phase angle decreases. The reflectance of a natural playa surface exhibits Lambertian scattering properties that cannot be replicated with the Hapke (1986) photometric function. The Hapke model does fit the rough and cobble-strewn playa data with single-scattering albedo, asymmetry factor, and roughness terms consistent with known properties of these surfaces. No published photometric model explains the reflectance variations of all three surfaces.

Guinness, Edward A.↗
NASA NTRS1992

Correcting x ray spectra obtained from the AXAF VETA-I mirror calibration for pileup, continuum, background and deadtime

The VETA-I mirror was calibrated with the use of a collimated soft X-ray source produced by electron bombardment of various anode materials. The FWHM, effective area and encircled energy were measured with the use of proportional counters that were scanned with a set of circular apertures. The pulsers from the proportional counters were sent through a multichannel analyzer that produced a pulse height spectrum. In order to characterize the properties of the mirror at different discrete photon energies one desires to extract from the pulse height distribution only those photons that originated from the characteristic line emission of the X-ray target source. We have developed a code that fits a modeled spectrum to the observed X-ray data, extracts the counts that originated from the line emission, and estimates the error in these counts. The function that is fitted to the X-ray spectra includes a Prescott function for the resolution of the detector a second Prescott function for a pileup peak and a X-ray continuum function. The continuum component is determined by calculating the absorption of the target Bremsstrahlung through various filters correcting for the reflectivity of the mirror and convolving with the detector response.

Chartas, G.↗
NASA NTRS1993

Correcting X-ray spectra obtained from the AXAF VETA-I mirror calibration for pileup, continuum, background and deadtime

The VETA-I mirror was calibrated with the use of a collimated soft X-ray source produced by electron bombardment of various anode materials. The FWHM, effective area and encircled energy were measured with the use of proportional counters that were scanned with a set of circular apertures. The pulsers from the proportional counters were sent through a multichannel analyzer that produced a pulse height spectrum. In order to characterize the properties of the mirror at different discrete photon energies one desires to extract from the pulse height distribution only those photons that originated from the characteristic line emission of the X-ray target source. We have developed a code that fits a modeled spectrum to the observed X-ray data, extracts the counts that originated from the line emission, and estimates the error in these counts. The function that is fitted to the X-ray spectra includes a Prescott function for the resolution of the detector a second Prescott function for a pileup peak and a X-ray continuum function. The continuum component is determined by calculating the absorption of the target Bremsstrahlung through various filters, correcting for the reflectivity of the mirror and convolving with the detector response.

Chartas, G.↗
DOE OSTI2016

SU-F-T-104: Determining the NTCP Parameters of Pharyngeal Constrictors and Proximal Esophagus for Radiation Induced Swallowing Problems Recorded Six Months After Radiation Therapy for Head and Neck Tumors

Purpose: To estimate the radiobiological parameters of four NTCP models that describe the dose-response relations of pharyngeal constrictors and proximal esophagus regarding the severity of patient reported swallowing problems 6 months post chemo-radiotherapy. To identify the section/structure that best correlates with the manifestation of the clinical endpoints. Finally, to compare the goodness-of-fit of those models. Methods: Forty-three patients were treated on a prospective multi-institutional phase II study for oropharyngeal squamous cell carcinoma. All the patients received 60 Gy IMRT and they reported symptoms using the novel patient reported outcome version of the CTCAE. We derived the individual patient dosimetric data of superior, medium and inferior sections of pharyngeal constrictors (SPC, MPC and IPC), superior and inferior sections of esophagus (SES and IES) as separate structures as well as combinations. The Lyman-Kutcher-Burman (LKB), Relative Seriality (RS), Logit and Relative Logit (RL) NTCP models were used to fit the patient data. The fitting of the different models was assessed through the area under the receiver operating characteristic curve (AUC) and the Odds Ratio methods. Results: The AUC values were highest for the SPC for Grade ≥ 2 (0.719 for the RS and RL models, and 0.716 for LKB and Logit). For Grade ≥ 1, the respective values were 0.699 for RS, LKB and Logit and 0.676 for RL. For MPC the AUC values varied between 0.463–0.477, for IPC between 0.396–0.458, for SES between 0.556–0.613 and for IES between 0.410–0.519. The Odds Ratio for the SPC was 15.6 (1.7–146.4) for RS, LKB and Logit for NTCP of 55%. Conclusion: All the examined NTCP models could fit the clinical data with similar accuracy. The SPC appear to correlate best with the clinical endpoints of swallowing problems. A prospective study could establish the use of NTCP values of SPC as a constraint in treatment planning.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Electron energy partition in the ‘above-the-looptop’ solar hard X-ray sources

The presentation begins by considering the non-thermal fraction of electron energies (R ε ) in 'above-the-looptop' (ALT). Several spectral models area considered: isothermal model, power-law, and the newly derived kappa distribution. Technically, all 3 non-thermal models can fit the data. In the discussion of flare scenario, energy partition, energization mechanism, and collisionality are considered. It is concluded that the kappa distribution works (R ε ≲ 50%), magnetic reconnection scenario works, and Coulomb collisions may reduce the non-thermal fraction of electron energies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Better practices for inferring ecosystem water use strategy from eddy covariance data

Eddy covariance data are critical for inferring ecosystem water use strategies. Yet, such inferences are sensitive to a range of assumptions applied across studies, hindering our understanding of water use strategies within and across eddy covariance sites. A recent analysis across 151 FLUXNET2015 and AmeriFlux-FLUXNET datasets found that poor model performance was the key driver of non-robust inferences of ecosystem water use strategies. Here, we leverage this previous analysis to (i) identify the specific assumptions that improve inference model performance across most sites, (ii) explain the mechanisms behind the performance improvements, and (iii) check whether better performance improves water use inference. We find that the common practice of fitting a model to canopy conductance (G c ) derived from the evapotranspiration (ET) observations, rather than to observed ET itself, artificially amplifies data errors and degrades the model performance. Next, accounting for vegetation dynamics by applying a growing season filter or incorporating satellite LAI data improves performance, but the former practice may remove soil water stress periods. Lastly, using the leaf-to-air vapor pressure deficit (VPD l ) derived from ET observations as a model input may artificially inflate performance. Based on these results, we recommend selecting observed ET (rather than derived G c ) as the response variable, carefully accounting for vegetation dynamics, and avoiding derived VPD l as a model input; these best practices improve model performance by c. 20% and robustness by c. 80% across all eddy covariance sites. Nevertheless, the performance improvements do not always correspond to more robust inference of water use strategies, as model parameter selection and surface energy budget closure corrections still strongly influence the ecosystem water use parameter estimation in a site-specific manner.

AmeriFlux↗
DOE OSTI2022

Bridging the Gap between the X-ray Absorption Spectroscopy and the Computational Catalysis Communities in Heterogeneous Catalysis: A Perspective on the Current and Future Research Directions

X-ray absorption spectroscopy (XAS) [extended X-ray absorption fine structure (EXAFS) and X-ray absorption near-edge structure (XANES)] is a key technique within the heterogeneous catalysis community to probe the structure and properties of the active site(s) for a diverse range of catalytic materials. Furthermore, the interpretation of the raw experimental data to derive an atomistic picture of the catalyst requires modeling and analysis; the EXAFS data are compared to a model, and a goodness of fit parameter is used to judge the best fit. This EXAFS modeling can often be nontrivial and time-consuming; overcoming or improving these limitations remains a central challenge for the community. Considering these limitations, this Perspective highlights how recent developments in analysis software, increased availability of reliable computational models, and application of data science tools can be used to improve the speed, accuracy, and reliability of EXAFS interpretation. In particular, we emphasize the advantages of combining theory and EXAFS as a unified technique that should be treated as a standard (when applicable) to identify catalytic sites and not two separate complementary methods. Building on the recent trends in the computational catalysis community, we also present a community-driven approach to adopt FAIR Guiding Principles for the collection, analysis, dissemination, and storage of XAS data. Written with both the experimental and theory audience in mind, we provide a unified roadmap to foster collaborations between the two communities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Rapid early gas accretion for the inner Galactic disc: A case for short accretion timescale

Recent observations of the Milky Way and galaxies at high redshifts suggest that galaxy discs were already in place soon after the Big Bang. While the gas infall history of the Milky Way in the inner disc has long been assumed to be characterised by a short accretion timescale, this has not been directly constrained using observations. Using data for the inner regions of the Milky Way recently produced by APOGEE and Gaia and of unprecedented quantity and quality, we aim to derive strong constraints on the infall history of the inner (< 6 kpc) Galaxy (with a focus on stars between 4 and 6 kpc, which we show is an appropriate proxy for the entire inner disc). We implemented gas infall into a chemical evolution model of the Galaxy disc, and used a Schmidt–Kennicutt law to connect the infall to the star formation. We explore a number of models, and two different formulations of the infall law. In one formulation, the infall is non-parametric, and in the other the infall has an explicitly exponential form. We fit the model parameters to the time–[Si/Fe] distribution of solar vicinity stars, and the metallicity and [Si/Fe] distribution function of stars with a galactocentric radius of between 4 and 6 kpc from APOGEE. Our results point to a fast, early gas accretion, and an upper limit on the accretion timescale of around 2 Gyr in the inner disc of the Milky Way. This suggests that at least half the baryons were in place within 2–3 Gyr of the Big Bang, and that half the stars of the inner disc formed within the first 5 Gyr, during the thick disc formation phase. This implies that the stellar mass of the inner disc is dominated by the thick disc, supporting our previous work, and that the gas accretion onto the inner disc was rapid and early.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Transferability of data-driven, many-body models for CO 2 simulations in the vapor and liquid phases

Here, extending on the previous work by Riera et al. [J. Chem. Theory Comput. 16, 2246–2257 (2020)], we introduce a second generation family of data-driven many-body MB-nrg models for CO 2 and systematically assess how the strength and anisotropy of the CO 2 –CO 2 interactions affect the models’ ability to predict vapor, liquid, and vapor–liquid equilibrium properties. Building upon the many-body expansion formalism, we construct a series of MB-nrg models by fitting one-body and two-body reference energies calculated at the coupled cluster level of theory for large monomer and dimer training sets. Advancing from the first generation models, we employ the charge model 5 scheme to determine the atomic charges and systematically scale the two-body energies to obtain more accurate descriptions of vapor, liquid, and vapor–liquid equilibrium properties. Challenges in model construction arise due to the anisotropic nature and small magnitude of the interaction energies in CO 2 , calling for the necessity of highly accurate descriptions of the multidimensional energy landscape of liquid CO 2 . These findings emphasize the key role played by the training set quality in the development of transferable, data-driven models, which, accurately representing high-dimensional many-body effects, can enable predictive computer simulations of molecular fluids across the entire phase diagram.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

The Latent Dirichlet Allocation model applied to airborne LiDAR data: A case study on mapping forest degradation associated with fragmentation and fire in the Amazon region

1. LiDAR data are being increasingly used to provide a detailed characterization of the vertical profile of forests. This characterization enables the generation of new insights on the influence of environmental drivers and anthropogenic disturbances on forest structure as well as on how forest structure influences important ecosystem functions and services. Unfortunately, extracting information from LiDAR data in a way that enables the spatial visualization of forest structure, as well as its temporal changes, is challenging due to the high dimensionality of these data. 2. In this study, we show how the Latent Dirichlet Allocation model applied to LiDAR data (LidarLDA) can be used to identify forest structural types and how the relative abundance of these forest types changes throughout the landscape. The code to fit this model is made available through the open-source r package LidarLDA in github. We illustrate the use of LidarLDA both with simulated data and data from a large-scale fire experiment in the Brazilian Amazon region. 3. Using simulated data, we demonstrate that LidarLDA accurately identifies the number of forest types as well as their spatial distribution and absorptance probabilities. For the empirical data, we found that LidarLDA detects both landscape-level patterns in forest structure as well as the strong interacting effect of fire and forest fragmentation on forest structure based on the experimental fire plots. More specifically, LidarLDA reveals that proximity to forest edge exacerbates the impact of fires, and that burned forests remain structurally different from unburned areas for at least 7 years, even when burned only once. Importantly, LidarLDA generates insights on the 3D structure of forest that cannot be obtained using more standard approaches that just focus on top-of-the-canopy information (e.g. canopy height models based on LiDAR data). 4. By enabling the mapping of forest structure and its temporal changes, we believe that LidarLDA will be of broad utility to the ecological research community.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Light output measurements and computational models of microcolumnar CsI scintillators for x-ray imaging

Purpose: The authors report on measurements of light output and spatial resolution of microcolumnar CsI:Tl scintillator detectors for x-ray imaging. In addition, the authors discuss the results of simulations aimed at analyzing the results of synchrotron and sealed-source exposures with respect to the contributions of light transport to the total light output. Methods: The authors measured light output from a 490-μm CsI:Tl scintillator screen using two setups. First, the authors used a photomultiplier tube (PMT) to measure the response of the scintillator to sealed-source exposures. Second, the authors performed imaging experiments with a 27-keV monoenergetic synchrotron beam and a slit to calculate the total signal generated in terms of optical photons per keV. The results of both methods are compared to simulations obtained with hybridMANTIS, a coupled x-ray, electron, and optical photon Monte Carlo transport package. The authors report line response (LR) and light output for a range of linear absorption coefficients and describe a model that fits at the same time the light output and the blur measurements. Comparing the experimental results with the simulations, the authors obtained an estimate of the absorption coefficient for the model that provides good agreement with the experimentally measured LR. Finally, the authors report light output simulation results and their dependence on scintillator thickness and reflectivity of the backing surface. Results: The slit images from the synchrotron were analyzed to obtain a total light output of 48 keV{sup −1} while measurements using the fast PMT instrument setup and sealed-sources reported a light output of 28 keV{sup −1}. The authors attribute the difference in light output estimates between the two methods to the difference in time constants between the camera and PMT measurements. Simulation structures were designed to match the light output measured with the camera while providing good agreement with the measured LR resulting in a bulk absorption coefficient of 5 × 10{sup −5} μm{sup −1}. Conclusions: The combination of experimental measurements for microcolumnar CsI:Tl scintillators using sealed-sources and synchrotron exposures with results obtained via simulation suggests that the time course of the emission might play a role in experimental estimates. The procedure yielded an experimentally derived linear absorption coefficient for microcolumnar Cs:Tl of 5 × 10{sup −5} μm{sup −1}. To the author’s knowledge, this is the first time this parameter has been validated against experimental observations. The measurements also offer insight into the relative role of optical transport on the effective optical yield of the scintillator with microcolumnar structure.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗