Search NASA⌕ Search

SEARCH · Search NASA

Results for “Function Estimation”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

TU-H-207A-08: Estimating Radiation Dose From Low-Dose Lung Cancer Screening CT Exams Using Tube Current Modulation

Purpose: The purpose of this work is to estimate effective and lung doses from a low-dose lung cancer screening CT protocol using Tube Current Modulation (TCM) across patient models of different sizes. Methods: Monte Carlo simulation methods were used to estimate effective and lung doses from a low-dose lung cancer screening protocol for a 64-slice CT (Sensation 64, Siemens Healthcare) that used TCM. Scanning parameters were from the AAPM protocols. Ten GSF voxelized patient models were used and had all radiosensitive organs identified to facilitate estimating both organ and effective doses. Predicted TCM schemes for each patient model were generated using a validated method wherein tissue attenuation characteristics and scanner limitations were used to determine the TCM output as a function of table position and source angle. The water equivalent diameter (WED) was determined by estimating the attenuation at the center of the scan volume for each patient model. Monte Carlo simulations were performed using the unique TCM scheme for each patient model. Lung doses were tallied and effective doses were estimated using ICRP 103 tissue weighting factors. Effective and lung dose values were normalized by scanspecific 32 cm CTDIvol values based upon the average tube current across the entire simulated scan. Absolute and normalized doses were reported as a function of WED for each patient. Results: For all ten patients modeled, the effective dose using TCM protocols was below 1.5 mSv. Smaller sized patient models experienced lower absolute doses compared to larger sized patients. Normalized effective and lung doses showed some dependence on patient size (R2 = 0.77 and 0.78, respectively). Conclusion: Effective doses for a low-dose lung screening protocol using TCM were below 1.5 mSv for all patient models used in this study. Institutional research agreement, Siemens Healthcare; Past recipient, research grant support, Siemens Healthcare; Consultant, Toshiba America Medical Systems; Consultant, Samsung Electronics.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Direct Estimation of Parameters in ODE Models Using WENDy: Weak-Form Estimation of Nonlinear Dynamics

Abstract We introduce the Weak-form Estimation of Nonlinear Dynamics (WENDy) method for estimating model parameters for non-linear systems of ODEs. Without relying on any numerical differential equation solvers, WENDy computes accurate estimates and is robust to large (biologically relevant) levels of measurement noise. For low dimensional systems with modest amounts of data, WENDy is competitive with conventional forward solver-based nonlinear least squares methods in terms of speed and accuracy. For both higher dimensional systems and stiff systems, WENDy is typically both faster (often by orders of magnitude) and more accurate than forward solver-based approaches. The core mathematical idea involves an efficient conversion of the strong form representation of a model to its weak form, and then solving a regression problem to perform parameter inference. The core statistical idea rests on the Errors-In-Variables framework, which necessitates the use of the iteratively reweighted least squares algorithm. Further improvements are obtained by using orthonormal test functions, created from a set of $$C^{\infty }$$ C ∞ bump functions of varying support sizes.We demonstrate the high robustness and computational efficiency by applying WENDy to estimate parameters in some common models from population biology, neuroscience, and biochemistry, including logistic growth, Lotka-Volterra, FitzHugh-Nagumo, Hindmarsh-Rose, and a Protein Transduction Benchmark model. Software and code for reproducing the examples is available at https://github.com/MathBioCU/WENDy .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Advanced measurement techniques in quantum Monte Carlo: The permutation matrix representation approach

In a typical finite temperature quantum Monte Carlo (QMC) simulation, estimators for simple static observables such as specific heat and magnetization are known. With a great deal of system-specific manual labor, one can sometimes also derive more complicated non-local or even dynamic observable estimators. In contrast, we show that arbitrary static observables can be estimated within the permutation matrix representation (PMR) flavor for any Hamiltonian. We then generalize these results to general imaginary-time correlation functions and non-trivial integrated susceptibilities thereof. Finally, we demonstrate the practical versatility of our method by estimating various non-local, random observables for the transverse-field Ising model on a square lattice and a toy random model.

Permutation matrix representation↗
DOE OSTI2017

Soil microbial beta-diversity is linked with compositional variation in aboveground plant biomass in a semi-arid grassland

Exploring biodiversity linkages between above-ground and below-ground biota is a core topic in ecology, and can have implications on our understanding of ecosystem process stability. Yet, this topic still remains underexplored. Here, we explored diversity linkages, in terms of both alpha- and beta- diversity, between plant and top soil microbial communities in a semi-arid grassland ecosystem. Soil microbial community structure was assessed based on both 16S rRNA and functional genes, and plant community composition was evaluated by traditional “species composition” and a newly-defined “biomass composition”, which includes the information on the biomass of each species. The bacterial alpha-diversity, expressed as the richness and Shannon diversity of 16S rRNA genes, was significantly correlated with plant species richness and Shannon diversity, whereas the alpha-diversity of microbial functional genes showed marginal association with total plant biomass. Microbial beta-diversity, evaluated by 16S rRNA genes, showed close relationship with plant beta-diversity estimated by both “species composition” and “biomass composition”, while the microbial beta-diversity based on functional genes was only associated with the compositional variation in above-ground plant biomass. These results showed that the differences in metabolic potential of soil microbial communities, which is closely related with ecosystem functions, can be better predicted by the variation of plant-derived resources returned to soil, than merely by the species composition of the macro-organism communities.

16s rRNA↗
DOE OSTI2022

β decay of Ba 141

The β-decay strength function of nuclides produced in fission is important as it dictates the distribution of decay energy between electrons, neutrinos, and γ rays and so is critical for calculating decay heat in reactors and for estimating the reactor antineutrino spectrum. Several experimental techniques are available to determine this strength function, including electron spectroscopy, γ-ray calorimetry (TAGS spectroscopy), and detailed, high-resolution spectroscopy with modern large high-purity germanium arrays. This work investigates the decay of the well-known and strongly produced fission fragment 141 Ba. A beam of 141 Cs was implanted at the target position of the Gammasphere and the subsequent decay of the daughter 141 Ba was studied. Extensive decay spectroscopy was possible up to the decay Q value of 3.197(7) MeV, including a significant extension of the level scheme and detailed angular correlation measurements for all levels with greater than 0.25% β feeding. The distribution of the β-decay strength was then inferred and compared to previous calorimetric studies. In conclusion, the agreement was excellent and provides a benchmark for comparing strength function methods and data for a more detailed understanding of the structure of 141 La.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Dynamic Modeling, Parameter Estimation, and Data Reconciliation of a Supercritical Pulverized Coal-Fired Boiler

A high-fidelity, spatially distributed dynamic model of a supercritical pulverized coal-fired boiler was developed for calculating transient thermal profiles of the boiler tubes, which are often not measured or unmeasurable in industrial settings due to the harsh operating conditions. In this work, rigorous models of the subcritical and supercritical steam properties, which can be highly nonlinear, were used to capture the transition between two-phase and one-phase flow, respectively, as supercritical boilers frequently transition through the critical point during load-following operation. The transient boiler model was also used for data reconciliation and parameter estimation for the purpose of model validation against industrial data. An approach was developed to reduce the size of large-scale reconciliation problem by focusing on the implementation of bias terms for the reconciled variables, which can be estimated in place of the reconciled variables themselves. In addition, a functional approximation of the reconciliation problem was defined so that it can be more simply optimized. Combining this functional approximation optimization with reformulation of the reconciliation problem considering measurement biases can significantly improve the tractability of many such large-scale problems. The validated boiler model was then used to study transient responses under load-following operation with specific attention on the calculation of unmeasurable variables.

42 ENGINEERING↗
DOE OSTI2015

MODELING ATMOSPHERIC EMISSION FOR CMB GROUND-BASED OBSERVATIONS

Atmosphere is one of the most important noise sources for ground-based cosmic microwave background (CMB) experiments. By increasing optical loading on the detectors, it amplifies their effective noise, while its fluctuations introduce spatial and temporal correlations between detected signals. We present a physically motivated 3D-model of the atmosphere total intensity emission in the millimeter and sub-millimeter wavelengths. We derive a new analytical estimate for the correlation between detectors time-ordered data as a function of the instrument and survey design, as well as several atmospheric parameters such as wind, relative humidity, temperature and turbulence characteristics. Using an original numerical computation, we examine the effect of each physical parameter on the correlations in the time series of a given experiment. We then use a parametric-likelihood approach to validate the modeling and estimate atmosphere parameters from the polarbear-i project first season data set. We derive a new 1.0% upper limit on the linear polarization fraction of atmospheric emission. We also compare our results to previous studies and weather station measurements. The proposed model can be used for realistic simulations of future ground-based CMB observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

CGC for ultra-peripheral Pb+Pb collisions at the Large Hadron Collider: a more realistic calculation

We provide the first calculation of two-gluon production at mid-rapidity in ultra-peripheral collisions in the Color Glass Condensate framework. To estimate systematic uncertainty associated with poor understanding of the wave function of the nearly real photon, we consider two diametrically different models: the dilute quark-antiquark dipole approximation and a vector meson, in which color charge density is approximated by McLerran-Venugopalan model. In the experimentally relevant range, the target nucleus can be faithfully approximated by a highly saturated state. This simplification enables us to perform efficient numerical simulations and extract the two-gluon correlation functions and the associated azimuthal harmonics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Fast Decomposition of Three-Component Spectra of Fluorescence Quenching by White and Grey Methods of Data Modeling

‘White’ and ‘grey’ methods of data modeling have been employed to resolve the heterogeneous fluorescence from a fluorophore mixture of 9-cyanoanthracene (CNA), 10-chloro-9-cyanoanthracene (ClCNA) and 9,10-dicyanoanthracene (DCNA) into component individual fluorescence spectra. The three-component spectra of fluorescence quenching in methanol were recorded for increasing amounts of lithium bromide used as a quencher. The associated intensity decay profiles of differentially quenched fluorescence of single components were modeled on the basis of a linear Stern-Volmer plot. These profiles are necessary to initiate the fitting procedure in both ‘white’ and ‘grey’ modeling of the original data matrices. ‘White’ methods of data modeling, called also ‘hard’ methods, are based on chemical/physical laws expressed in terms of some well-known or generally accepted mathematical equations. The parameters of these models are not known and they are estimated by least squares curve fitting. ‘Grey’ approaches to data modeling, also known as hard-soft modeling techniques, make use of both hard-model and soft-model parts. In practice, the difference between ‘white’ and ‘grey’ methods lies in the way in which the ‘crude’ fluorescence intensity decays of the mixture components are estimated. In the former case they are given in a functional form while in the latter as digitized curves which, in general, can only be obtained by using dedicated techniques of factor analysis. In the paper, the initial values of the Stern-Volmer constants of pure components were evaluated by both ‘point-by-point’ and ‘matrix’ versions of the method making use of the concept of wavelength dependent intensity fractions as well as by the rank annihilation factor analysis applied to the data matrices of the difference fluorescence spectra constructed in two ways: from the spectra recorded for a few excitation lines at the same concentration of a fluorescence quencher or classically from a series of the spectra measured for one selected excitation line but for increasing concentration of the quencher. The results of multiple curve resolution obtained by all types of the applied methods have been scrutinized and compared. In addition, the effect of inadequacy of sample preparation and increasing instrumental noise on the shape of the resolved spectral profiles has been studied on several datasets mimicking the measured data matrices. .

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Reducing the small-heart effect in pediatric gated myocardial perfusion single-photon emission computed tomography

Background: We compared two reconstruction algorisms and two cardiac functional evaluation software programs in terms of their accuracy for estimating ejection fraction (EF) of small hearts (SH). Methods: The study group consisted of 66 pediatric patients. Data were reconstructed using a filtered back projection (FBP) method without the resolution correction (RC) and an iterative method based on an ordered subset expectation maximization (OSEM) algorithm with the RC. EF was evaluated using two software programs of quantitative gated single-photon emission computed tomography (SPECT) (QGS) and cardioREPO. We compared the EF of gated myocardial perfusion SPECT to echocardiographic measurement (Echo). Results: Forty-eight of 66 patients had an end-systolic volume < 20 mL which was used as the criterion for being included in the SH group, and the SH effect led to an overestimation of EF. While significant differences were observed between Echo (66.9 ± 5.0%) and QGS-FBP without RC (76.9 ± 8.4%, P < .0001), QGS-OSEM with RC (76.6 ± 8.6%, P < .0001), and cardioREPO-FBP without RC (72.1 ± 10.0%, P = .0011), no significant difference was observed between Echo and cardioREPO-OSEM with RC (67.4 ± 6.1%) in SH group. Conclusions: In pediatric gated myocardial perfusion SPECT, the SH effect can be significantly reduced when an OSEM algorithm is used with RC in combination with the specific cardioREPO algorithm.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Toward the understanding of hydration phenomena in aqueous electrolytes from the interplay of theory, molecular simulation, and experiment

We confront the microstructural analysis of aqueous electrolytes and present a detailed account of the fundamentals underlying the neutron scattering with isotopic substitution (NDIS) approach for the experimental determination of ion coordination numbers in systems involving both halides anions and oxyanions. We place particular emphasis on the frequently overlooked ion-pairing phenomenon, identify its microstructural signature in the neutron-weighted distribution functions, and suggest novel techniques to deal with either the estimation of the ion-pairing magnitude or the correction of its effects on the experimentally measured coordination numbers. We illustrate the underlying ideas by applying these new developments to the interpretation of four NDIS test-cases via molecular simulation, as convenient dry runs for the actual scattering experiments, for representative aqueous electrolyte solutions at ambient conditions involving metal halides and nitrates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Concomitant Imaging Dose and Cancer Risk in Image Guided Thoracic Radiation Therapy

Purpose: Kilovoltage cone beam computed tomography (CT) (kVCBCT) imaging guidance improves the accuracy of radiation therapy but imposes an extra radiation dose to cancer patients. This study aimed to investigate concomitant imaging dose and associated cancer risk in image guided thoracic radiation therapy. Methods and Materials: The planning CT images and structure sets of 72 patients were converted to CT phantoms whose chest circumferences (C{sub chest}) were calculated retrospectively. A low-dose thorax protocol on a Varian kVCBCT scanner was simulated by a validated Monte Carlo code. Computed doses to organs and cardiac substructures (for 5 selected patients of various dimensions) were regressed as empirical functions of C{sub chest}, and associated cancer risk was calculated using the published models. The exposures to nonthoracic organs in children were also investigated. Results: The structural mean doses decreased monotonically with increasing C{sub chest}. For all 72 patients, the median doses to the heart, spinal cord, breasts, lungs, and involved chest were 1.68, 1.33, 1.64, 1.62, and 1.58 cGy/scan, respectively. Nonthoracic organs in children received 0.6 to 2.8 cGy/scan if they were directly irradiated. The mean doses to the descending aorta (1.43 ± 0.68 cGy), left atrium (1.55 ± 0.75 cGy), left ventricle (1.68 ± 0.81 cGy), and right ventricle (1.85 ± 0.84 cGy) were significantly different (P<.05) from the heart mean dose (1.73 ± 0.82 cGy). The blade shielding alleviated the exposure to nonthoracic organs in children by an order of magnitude. Conclusions: As functions of patient size, a series of models for personalized estimation of kVCBCT doses to thoracic organs and cardiac substructures have been proposed. Pediatric patients received much higher doses than did the adults, and some nonthoracic organs could be irradiated unexpectedly by the default scanning protocol. Increased cancer risks and disease adverse events in the thorax were strongly related to higher imaging doses and smaller chest dimensions.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2019

Study on impact of Cr and Mo on diffusion of H in 2.25Cr1Mo steel using first-principle calculations

Chrome-molybdenum steel (2.25Cr1Mo steel) is one of the main materials of a steam generator (SG) in high-temperature reactor-pebblebed modules (HTR-PM). It is essential to analyze the source term of tritium in this material, because the behavior of tritium in a SG is important for performing source term analysis in normal and accident conditions. In this article, the diffusion behavior of H atom in 2.25Cr1Mo steel was calculated to estimate the diffusivity of tritium using first-principles density functional theory. To develop and simplify the model of hydrogen diffusion in 2.25Cr1Mo steel, the impact of Cr and Mo on the diffusion of hydrogen in bcc-Fe were first calculated, all the possible diffusion paths were considered, and the minimum energy path was obtained. The diffusion activation energy and pre-exponential factor of the diffusion coefficient were obtained from Vienna Ab initio Simulation Package combined with the climbing image-nudged elastic band method. Here the results indicate that the minimum energy path for the impurity H atom is from one tetrahedral interstitial site to an adjacent tetrahedral interstitial site. The function of the diffusion coefficient of H in 2.25Cr1Mo steel with temperature T can be expressed as D = 1.486 X 10 -7 x ($\frac{-16.350kJ/mol}{RT}$) (M2/s). The diffusion coefficient of our calculation and some of the previous experiments have an excellent quantitative agreement, which indicates the reliability of our crystalline model and the practicability of the present theoretical approach. More importantly, the computational results in this work can be treated as a good screening method to collect reasonable experimental data, which will provide a good reference for tritium source term evaluation in the SG of the HTR-PM.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Fundamental flow characterization for LiCl–KCl molten salt natural circulation closed loop

Molten salt flow characterization presents numerous challenges due to the hostile environments that cause corrosion and erosion of instrumentation and structural components. To support future molten salt flow characterization under natural convection conditions, a LiCl–KCl eutectic molten salt loop was designed and constructed. The primary objective of this loop is to facilitate natural convection as the principal cooling mechanism, enabling insights into the physics behind natural convection in closed loops and serving as a baseline for loop dimensional and non-dimensional analyses. Here, this study summarizes the fundamental functions employed for flow characterization and introduces a model to estimate flow rate when intrusive instrumentation is unavailable. The approach relies solely on structural surface measurements to avoid flow disturbances. Furthermore, a comprehensive thermal–hydraulic analysis is included to deepen understanding of the system’s behavior.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Comparing the performance of density functionals in describing the adsorption of atoms and small molecules on Ni(111)

Reaction networks calculated using density functional theory (DFT) play a crucial role in understanding and predicting heterogeneous catalyst performance. However, results obtained using different density functionals can vary significantly. Accordingly, it is important to assess the accuracy of available functionals in capturing adsorption energetics derived from single-crystal adsorption microcalorimetry experiments. Here, we evaluate how different density functionals describe the molecular and dissociative adsorption of CH 3 I, the adsorption of CH 3 , I, and H, and the dissociative adsorption of CH 4 on Ni(111). First, we identify energetically-preferred adsorbed states for each adsorbate at three surface coverages and calculate the corresponding enthalpies of adsorption at 160 K. Then, we compare the results to experimental measurements published in the literature. For each adsorption process we find at least one functional that is quantitatively accurate within experimental error, but no functional is accurate (within experimental error) for all adsorption processes studied. Even when assuming an additional DFT error of ±20 kJ/mol on top of the experimental error, only the PBE-D3 and RPBE-D3 functionals are accurate for all considered adsorption systems. Overall, these results indicate that quantitative agreement between density functional theory calculations and experimental measurements is both system- and functional-dependent. Importantly, single-crystal adsorption microcalorimetry experiments will continue to play an important role for benchmarking density functionals.

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

An Attractive Way to Correct for Missing Singles Excitations in Unitary Coupled Cluster Doubles Theory

Coupled cluster methods based exclusively on double excitations are comparatively “cheap” and interesting model chemistries, as they are typically able to capture the bulk of the dynamic electron correlation effects. The trade-off in such approximations is that the effect of neglected excitations, particularly single excitations, can be considerable. Using standard and electron-pair-restricted T 2 operators to define two flavors of unitary coupled cluster doubles (UCCD) methods, we investigate the extent to which missing single excitations can be recovered from low-order corrections in many-body perturbation theory (MBPT) within the unitary coupled cluster (UCC) formalism. Here, our analysis includes the derivations of finite-order UCC energy functionals, which are used as a basis to define perturbative estimates of missed single excitations. This leads to the novel UCCD[4S] and UCCD[6S] methods, which consider energy corrections for missing single excitations through fourth- and sixth-order in MBPT, respectively. We also apply the same methodology to the electron-pair-restricted ansatz, but the improvements are only marginal. Our findings show that augmenting UCCD with these post hoc perturbative corrections can lead to UCCSD-quality results.

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

Femtosecond X-ray Spectroscopy Directly Quantifies Transient Excited-State Mixed Valency

Quantifying charge delocalization associated with short-lived photoexcited states of molecular complexes in solution remains experimentally challenging, requiring local element specific femtosecond experimental probes of time-evolving electron transfer. In this study, we quantify the evolving valence hole charge distribution in the photoexcited charge transfer state of a prototypical mixed valence bimetallic iron–ruthenium complex, [(CN) 5 Fe II CNRu III (NH 3 ) 5 ]-, in water by combining femtosecond X-ray spectroscopy measurements with time-dependent density functional theory calculations of the excited-state dynamics. We estimate the valence hole charge that accumulated at the Fe atom to be 0.6 ± 0.2, resulting from excited-state metal-to-metal charge transfer, on an ~ 60 fs time scale. Our combined experimental and computational approach provides a spectroscopic ruler for quantifying excited-state valency in solvated complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗