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 1,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
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 ↗
DOE OSTI2022

People who inject drugs in metropolitan Chicago: A meta-analysis of data from 1997-2017 to inform interventions and computational modeling toward hepatitis C microelimination

Progress toward hepatitis C virus (HCV) elimination in the United States is not on track to meet targets set by the World Health Organization, as the opioid crisis continues to drive both injection drug use and increasing HCV incidence. A pragmatic approach to achieving this is using a microelimination approach of focusing on high-risk populations such as people who inject drugs (PWID). Computational models are useful in understanding the complex interplay of individual, social, and structural level factors that might alter HCV incidence, prevalence, transmission, and treatment uptake to achieve HCV microelimination. However, these models need to be informed with realistic sociodemographic, risk behavior and network estimates on PWID. We conducted a meta-analysis of research studies spanning 20 years of research and interventions with PWID in metropolitan Chicago to produce parameters for a synthetic population for realistic computational models (e.g., agent-based models). We then fit an exponential random graph model (ERGM) using the network estimates from the meta-analysis in order to develop the network component of the synthetic population.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

A measurement of the proton’s spin structure function g2 at low Q2

JLab E08-027, a measurement of g2p and the longitudinal-transverse (LT) spin polarizability, successfully collected data from March to May, 2012. Nucleon spin structure study has been an active research area, which has attracted a very large effort from both experimentalists and theorists. The spin structure study for the last 2 decades has provided us with many exciting and often surprising results. Recently, new precision results in the low-to-intermediate momentum transfer Q2 region from JLab have provided extensive information on the nucleon structure in the confinement region and the transition region between asymptotic free to confinement. In particular, the extensive comparisons of experimental results with Chiral Perturbation Theory (the effective theory of QCD at low energy) calculations show general good agreements, but strong disagreement in the case of the neutron LT spin polarizability. This experiment completed the measurements of gp2 and the LT spin polarizability on the proton in the low-to-intermediate Q2 region. The experiment used a polarized proton (NH3) target for the first time in Hall A. Scattered electrons were detected by a pair of Hall A high resolution spectrometer (HRS) with a pair of septum magnets. To avoid too much depolarization of the target, beam current was limited to 50-100 nA during the experiment. Since the existing beam current monitors (BCMs), beam position monitors (BPMs) and calibration methods did not work at such a low current range, new BPM and BCM receivers were designed and used for current condition. A pair of super-harps and a tungsten calorimeter were installed to calibrate the BPMs and BCMs. To compensate for the effect of the 2.5/5T transverse magnet field, two chicane dipole magnets were installed. A pair of slow rasters were installed for the first time in Hall A, combining with a pair of fast raster. The standard Hall A DAQ system and the improved high resolution DAQ system were used to record the detector information and the helicity dependent beam information, respectively. In order to achieve the required accuracy of the beam position and angle and reconstruct them event by event at the target location, the data of the BPMs and harps were carefully analyzed. The final uncertainty of BCM after the calibration using the tungsten calorimeter is below 1%, which is important for the asymmetry extraction. Before the acceptance and the dilution factor are available, the models from fitting to the world data were used instead to extract the unpolarized cross section and the dilution factor. The longitudinal and transverse physics asymmetries were extracted and compared with the model, with the study of the radiative correction. The preliminary g1 and g2 results were then extracted using the measured asymmetry and the model.

Zhu, Pengjia↗
NASA NTRS1984

Impact experiments in viscous fluid media

Available phase and group velocity data are inverted by a procedure which includes the effects of transverse anisotropy, anelastic dispersion, sphericity, and gravity. The resulting models, for average Earth, average ocean, and oceanic regions divided according to the age of the ocean floor, are quite different from previous results which ignore the above effects. The models show a low-veocity zone with age dependent anisotropy and velocities higher than derived in previous surface wave studied. The correspondence between the anisotropy variation with age and a physical model based on flow aligned olivine is suggestive. For most of the Earth SHSV in the vicinity of the low-velocity zone. Near the East Pacific Rise, however, SVSH at depth, consistent with ascending flow. Anisotropy is as important as temperature in causing radial and lateral variations in velocity. The models have a high velocity nearly isotropic layer at the top of the mantle that thickens with age. This layer defines the LID, or seismic lithosphere. In the Pacific, the LID thickens with age to a maximum thickness of about 50 km. This thickness is comparable to the thickness of the elastic lithosphere. The LID thickness is thinner than derived using isotropic or pseudo-isotropic procedures A new model for Average Earth is obtained which includes a thin LID. This model extends the fit of a P.R.E.M. type model to shorter period surface waves.

Greeley, R.↗
NASA NTRS1990

Starbursts in interacting galaxies: Observations and models

Starbursts have been a puzzling field of research for more than a decade. It is evident that they played a significant role in the evolution of many galaxies but still quite little is known about the starburst mechanisms. A way towards a better interpretation of the available data is the comparison with evolution models of starburst. The modelling of starbursts and the fitting of such model starbursts to observed data is discussed. The models were applied to a subset of starburst and post-starburst galaxies in a sample of 30 interacting systems. These galaxies are not ultraluminous far infrared (FIR) galaxies but rather ordinary starburst galaxies with FIR luminosities of a few 10(exp 10) to a few 10(exp 11) solar luminosities.

Bernloehr, Konrad↗
NASA NTRS1993

Comparison of distributed reacceleration and leaky-box models of cosmic-ray abundances (Z = 3-28)

A large collection of elemental and isotopic cosmic-ray data has been analyzed using the leaky-box transport model with and without reacceleration in the interstellar medium. Abundances of isotopes and elements with charges Z = 3-28 and energies E = 10 MeV/nucleon-1 TeV/nucleon were explored. Our results demonstrate that reacceleration models make detailed and accurate predictions with the same number of parameters or fewer as standard leaky-box models. Ad hoc fitting parameters in the standard model are replaced by astrophysically significant reacceleration parameters. Distributed reacceleration models explain the peak in secondary-to-primary ratios around 1 GeV/nucleon. They diminish the discrepancy between rigidity-dependent leakage and energy-independent anisotropy. They also offer the possibility of understanding isotopic anomalies at low energy.

Letaw, John R.↗