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,765 records · Page 98

Analysis of the Trusted Inertial Terrain-Aided Navigation Measurement Function

The trusted inertial terrain-aided navigation (TITAN) algorithm leverages an airborne vertical synthetic aperture radar to measure the range to the closest ground points along several prescribed iso-Doppler contours. These TITAN minimum-range, prescribed-Doppler measurements are the result of a constrained nonlinear optimization problem whose optimization function and constraints both depend on the radar position and velocity. Owing to the complexity of this measurement definition, analysis of the TITAN algorithm is lacking in prior work. This publication offers such an analysis, making the following three contributions: (1) an analytical solution to the TITAN constrained optimization measurement problem, (2) a derivation of the TITAN measurement function Jacobian, and (3) a derivation of the Cramér-Rao lower bound on the estimated position and velocity error covariance. These three contributions are verified via Monte Carlo simulations over synthetic terrain, which further reveal two remarkable properties of the TITAN algorithm: (1) the along-track positioning errors tend to be smaller than the cross-track positioning errors, and (2) the cross-track positioning errors are independent of the terrain roughness.

TITAN↗

Retrieval of Cloud Liquid Water Using Microwave Signals from LEO Satellites: A Feasibility Study through Simulations

A novel method, using low Earth orbit (LEO) satellite microwave communication links for cloud liquid water measurements, is introduced in this paper. The feasibility of this approach is studied through simulations of the retrieval system including a LEO satellite communicating with a group of ground receivers equipped with signal-to-noise ratio (SNR) estimators, a synthetic cloud attenuation field and a tomographic retrieval algorithm. Rectangular and Gaussian basis functions are considered to define the targeted field. Simulation results suggest that the proposed least-squares based retrieval algorithm produces satisfactory outcomes for both types of basis functions. The root-mean-square error of the retrieved field is around 0.2 dB/km, with the range of the reference field as 0 to 2.35 dB/km. It is also confirmed that the partial retrieval of the cloud field is achievable when a limited number of receivers with restricted locations are available. The retrieval outcomes exhibit properties of high resolution and low error, indicating that the proposed approach has great potential for cloud observations.

54 ENVIRONMENTAL SCIENCES↗

An Intermediate-age Alpha-rich Galactic Population in K2

We explore the relationships between the chemistry, ages, and locations of stars in the Galaxy using asteroseismic data from the K2 mission and spectroscopic data from the Apache Point Galactic Evolution Experiment survey. Previous studies have used giant stars in the Kepler field to map the relationship between the chemical composition and the ages of stars at the solar circle. Consistent with prior work, we find that stars with high [α/Fe] have distinct, older ages in comparison to stars with low [α/Fe]. We provide age estimates for red giant branch (RGB) stars in the Kepler field, which support and build upon previous age estimates by taking into account the effect of α-enrichment on opacity. Including this effect for [α/Fe]-rich stars results in up to 10% older ages for low-mass stars relative to corrected solar mixture calculations. This is a significant effect that Galactic archeology studies should take into account. Looking beyond the Kepler field, we estimate ages for 735 RGB stars from the K2 mission, mapping age trends as a function of the line of sight. We find that the age distributions for low- and high-[α/Fe] stars converge with increasing distance from the Galactic plane, in agreement with suggestions from earlier work. We find that K2 stars with high [α/Fe] appear to be younger than their counterparts in the Kepler field, overlapping more significantly with a similarly aged low-[α/Fe] population. This observation may suggest that star formation or radial migration proceeds unevenly in the Galaxy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analysis of Polarimetry Data with Angular Uncertainties

For a track based polarimeter, such as the Imaging X-ray Polarimetry Explorer (IXPE), the sensitivity to polarization depends on the modulation factor, which is a strong function of energy. In previous work, a likelihood method was developed that would account for this variation in order to estimate the minimum detectable polarization (MDP). That method essentially required that the position angles of individual events should be known precisely. In a separate work, however, it was shown that using a machine-learning method for measuring event tracks can generate track angle uncertainties, which can be used in the analysis. Here, the maximum likelihood method is used as a basis for revising the estimate of the MDP in a general way that can include uncertainties in event track position angles. The resultant MDP depends solely upon the distribution of track angle uncertainties present in the input data. Due to the physics of the IXPE detectors, it is possible to derive a simple relationship between these angular uncertainties and the energy-dependent modulation function as a step in the process.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Collapsar R-process Yields Can Reproduce [Eu/Fe] Abundance Scatter in Metal-poor Stars

It is unclear if neutron star mergers can explain the observed r-process abundances of metal-poor stars. Collapsars, defined here as rotating massive stars whose collapse results in a rapidly accreting disk around a black hole that can launch jets, are a promising alternative. We find that we can produce a self-consistent model in which a population of collapsars with stochastic europium yields synthesizes all of the r-process material in metal-poor ([Fe/H] < - 2.5) stars. Our model reproduces the observed scatter and evolution of scatter of [Eu/Fe] abundances. We find that if collapsars are the dominant r-process site for metal-poor stars, r-process synthesis may be linked to supernovae that produce long γ-ray bursts. Our results also allow for the possibility that core-collapse supernovae beyond those that launch γ-ray bursts also produce r-process material (e.g., potentially a subset of Type Ic-BL supernovae). Furthermore, we identify collapsar jet properties (isotropic energy, engine luminosity, or engine time) that may trace r-process yield and verify that the amount of r-process yield produced per collapsar in our model ( ~ 0.07M⊙) is consistent with other independent estimates. In the future, achieving 0.05 dex precision on distribution scatter or a reliable selection function would further constrain our probe of r-process production. Our model would also hold for another prompt r-process site with a power-law yield, and work is needed to determine if, for example, fast-merging neutron stars can also explain abundance scatter.

Astronomy & Astrophysics↗

Detailed Radiation Dose Rate Evaluations of Commercial Spent Nuclear Fuel Canisters

Electronic devices that are being designed for continuous monitoring of spent nuclear fuel (SNF) canister internal integrity may be subjected to very high radiation levels depending on their placement. Components of electronic devices (e.g., silicon semiconductors, optical glasses, insulators) have a limited tolerance to ionizing radiation. For selected locations on the inner and outer surfaces of a representative dry storage pressurized water reactor (PWR) SNF canister and a representative dry storage boiling water reactor (BWR) SNF canister, this paper provides (1) the dose rate in silicon as a function of decay time and (2) the time-integrated dose in silicon. The analyzed decay time interval is from 5 to 70 years after fuel discharge. Each representative canister contained identical SNF assemblies with bounding radiation sources. Surface dose rates varied significantly as a function of location and SNF assembly decay time. At a 5-year decay time, the estimated maximum dose rate values were 480 Gy/h (PWR canister) and ~366 Gy/h (BWR canister), and the estimated minimum dose rate was ~0.2 Gy/h for both canisters. Over a 65-year period, the dose rate decreased by a factor between 13 and 580, depending on location. The total dose over a time interval could be very large (e.g., estimated maximum value over a 5-year interval is 1.5E+07 Gy) depending on location and the time interval over which dose is accumulated.

Radulescu, Georgeta↗

Quarterly Research Performance Progress Report (Q8)

As part of Task 1, we have started by testing our modeling capabilities by reproducing isothermal DFIT simulations presented in the literature. Once satisfied with the results we have started by targeting the modeling of the DFITs at conducted at well 58-32. We have a identified a specific test (cycle 4 in zone 2) as the most interesting to be model with GEOS hydraulic fracturing module. Thus, we have first produced results with an isothermal model and adjusted model parameters to get a satisfying match with field pressure data. The, we have added thermal effects and compared the modeling results with and without thermal effects to estimate how thermal effects may influence test interpretation. Models seem to suggest that, for small volumes of fluid, thermal effects are moderate. In Task 2, we have adapted GEOS phase-field formulation to be able to simulate near-wellbore hydraulic fracture nucleation and propagation. We have devised a novel formulation that, compared to other existing ones, incorporates rock strengths. We have submitted a journal publication about our work. We are currently employing this phase-field formulation to model the experiments taking place at U Pitt and help us understand the effect of various parameters. In Task 3, we have built a model of the region surround well 16A and have started modeling stage 3 stimulation because of its simpler planar geometry. After calibrating simulation parameters using known analytical solutions, we have simulated the stage 3 stimulation using our isothermal hydraulic fracturing module, varying the permeability field, the stress conditions including different physics to get a better understanding of the numerical challenges and of the effects of varying these parameters on the simulation results. In Task 4 laboratory experimentation, a set of specialized drilling and injection tools has been customized and constructed to accommodate an inclined well with an orientation of up to 30 degrees relative to material anisotropy or principal stress axes. These inclined samples have also undergone thermal stress and hydraulic fracturing at a temperature of 190 degrees Celsius. Furthermore, both vertical and deviated sampling testing setups enable an extended analysis of post-peak pressure behaviors, facilitating post-test pressure analyses such as the G-function, step rate, and fracture reopening measurements. Thus, the key components of in-situ stress estimation can be extracted and validated through our experiment, providing a solid foundation for validating existing in-situ stress estimation theories or proposing new ones. Simultaneously, we are integrating computer vision techniques with traditional experimental fracture observation methods such as multi-overcore/slicing and water-penetration fracture observation. This combination will prove beneficial in populating the hydraulic fracture patterns database, generated under challenging EGS conditions. This approach aims to deepen our understanding of the complexities in EGS reservoirs and pave the way for future data-driven investigations. Additionally, PITT has also equipped the ELE International compression machine, which is now prepared for conducting indirect tensile and fracture toughness tests. These tests will aid in characterizing how rock fabrics influence the resulting fracture patterns. Additionally, we have completed the required personnel training and gained access to Scanning Electron Microscopy (SEM) and Energy Dispersive Spectroscopy (EDS) for conducting more detailed characterization and analysis of rock fabrics, as well as the examination of thermal and hydraulically induced fracture patterns. Thus, the PITT team has effectively demonstrated the capabilities of our experimental apparatuses in exploring the thermal effects, well deviation angles, material anisotropy, and operational choices (such as circulation rate, injection fluid viscosity, and injection rate) and their impact on pressure responses and fracture trajectories under the Utah FORGE conditions.

15 GEOTHERMAL ENERGY↗

Bimodal imaging of functional changes in blood flow using optoacoustic and optical coherent angiography

The results of a complementary study of functional changes of the blood flow in human superficial tissues using optoacoustic (OA) microscopy and angiographic optical coherent tomography (OCT-A) during an occlusive test are demonstrated for the first time. OA microscopy data allow one to estimate the increase in blood content during the occlusive test, while OCT-A data demonstrate a decrease in the number of vessels with a functioning blood flow in the area of interest. The combined approach reported in the paper can be employed in experimental and clinical medicine to study the functional changes of microcirculation in diagnostics of vascular pathologies of superficial tissues and to evaluate the microcirculation response to treatment. (paper)

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Analytical solution and parameter estimation for heat of wetting and vapor adsorption during spontaneous imbibition in tuff

An analytical expression is derived for the thermal response observed during spontaneous imbibition of water into a dry core of zeolitic tuff. Sample tortuosity, thermal conductivity, and thermal source strength are estimated from fitting an analytical solution to temperature observations during a single laboratory test. Here, the closed-form analytical solution is derived using Green's functions for heat conduction in the limit of “slow” water movement; that is, when advection of thermal energy with the wetting front is negligible. The solution has four free fitting parameters and is efficient for parameter estimation. Laboratory imbibition data used to constrain the model include a time series of the mass of water imbibed, visual location of the wetting front through time, and temperature time series at six locations. The thermal front reached the end of the core hours before the visible wetting front. Thus, the predominant form of heating during imbibition in this zeolitic tuff is due to vapor adsorption in dry zeolitic rock ahead of the wetting front. The separation of the wetting front and thermal front in this zeolitic tuff is significant, compared to wetting front behavior of most materials reported in the literature. This work is the first interpretation of a thermal imbibition response to estimate transport (tortuosity) and thermal properties (including thermal conductivity) from a single laboratory test.

36 MATERIALS SCIENCE↗

Insights on finite size effects in ab initio study of CO adsorption and dissociation on Fe 110 surface

Adsorption and dissociation of hydrocarbons on metallic surfaces represent crucial steps on the path to carburization, eventually leading to dusting corrosion. While adsorption of CO molecules on Fe surface is a barrier-less exothermic process, this is not the case for the dissociation of CO into C and O adatoms and the diffusion of C beneath the surface that are found to be associated with large energy barriers. In practice, these barriers can be affected by numerous factors that combine to favour the CO-Fe reaction such as the abundance of CO and other hydrocarbons as well as the presence of structural defects. From a numerical point of view, studying these factors is challenging and a step-by-step approach is necessary to assess, in particular, the influence of the finite box size on the reaction parameters for adsorption and dissociation of CO on metal surfaces. Here, we use density functional theory (DFT) total energy calculations with the climbing-image nudged elastic band method to estimate the adsorption energies and dissociation barriers for different CO coverages with surface supercells of different sizes. We further compute the effect of periodic boundary condition for DFT calculations and find that the contribution from van der Waals interaction in the computation of adsorption parameters is important as they contribute to correcting the finite-size error in small systems. The dissociation process involves carbon insertion into the Fe surface causing a lattice deformation that requires a larger surface system for unrestricted relaxation. We show that, in the larger surface systems associated with dilute CO-coverages, C-insertion is energetically more favourable, leading to a significant decrease in the dissociation barrier. This observation suggests that a large surface system with dilute coverage is necessary for all similar metal-hydrocarbon reactions in order to study their fundamental electronic mechanisms, as an isolated phenomenon, free from finite-size effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Spin asymmetries for vector boson production in polarized p + p collisions

We study the cross section for vector boson (W ± /Z 0 /γ $\star$ ) production in polarized nucleon-nucleon collisions for low transverse momentum of the observed vector boson. For the case where one measures the transverse momentum and azimuthal angle of the vector bosons, we present the cross sections and the associated spin asymmetries in terms of transverse momentum dependent parton distribution functions (TMDs) at tree level within the TMD factorization formalism. To assess the feasibility of experimental measurements, we estimate the spin asymmetries forW ± /Z 0 boson production in polarized proton-proton collisions at the Relativistic Heavy Ion Collider by using current knowledge of the relevant TMDs. Here, we find that some of these asymmetries can be sizable if the suppression effect from TMD evolution is not too strong. The W program at RHIC can, thus, test and constrain spin theory by providing unique information on the universality properties of TMDs, TMD evolution, and the nucleon structure. For example, the single transverse spin asymmetries could be used to probe the well-known Sivers function f$⊥q\atop{1T}$, as well as the transversal helicity distribution g$q\atop{1T}$ via the parity-violating nature of W production.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Facile Electron Transfer to CO 2 during Adsorption at the Metal|Solution Interface

We estimate the rate of electron transfer to CO 2 at the Au (211)|water interface during adsorption in an electrochemical environment under reducing potentials. On the basis of density functional theory calculations at the generalized gradient approximation and hybrid levels of theory, we find electron transfer to the adsorbed *CO 2 to be very facile. This high rate of transfer is estimated by the energy distribution of the adsorbate-induced density of states as well as from the interaction between diabatic states representing neutral and negatively charged CO 2 . Up to 0.62 electrons is transferred to CO 2 , and this charge adiabatically increases with the bending angle to a lower limit of 137°. Finally, we conclude that this rate of electron transfer is extremely fast compared to the time scale of the nuclear degrees of freedom, that is, the adsorption process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Uncertainty analysis of signal deconvolution using a measured instrument response function

A common analysis procedure minimizes the ln-likelihood that a set of experimental observables matches a parameterized model of the observation. The model includes a description of the underlying physical process as well as the instrument response function (IRF). Here, we investigate the National Ignition Facility (NIF) neutron time-of-flight (nTOF) spectrometers, the IRF is constructed from measurements and models. IRF measurements have a finite precision that can make significant contributions to the uncertainty estimate of the physical model’s parameters. Finally, we apply a Bayesian analysis to properly account for IRF uncertainties in calculating the ln-likelihood function used to find the optimum physical parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Optimization-based scatter estimation using primary modulation for computed tomography

Purpose: Scatter reduces the image quality in computed tomography (CT), but scatter correction remains a challenge. A previously proposed primary modulation method simultaneously obtains the primary and scatter in a single scan. However, separating the scatter and primary in primary modulation is challenging because it is an underdetermined problem. In this study, an optimization-based scatter estimation (OSE) algorithm is proposed to estimate and correct scatter. Methods: In the concept of primary modulation, the primary is modulated, but the scatter remains smooth by inserting a modulator between the x-ray source and the object. In the proposed algorithm, an objective function is designed for separating the scatter and primary. Prior knowledge is incorporated in the optimization-based framework to improve the accuracy of the estimation: (1) the primary is always positive; (2) the primary is locally smooth and the scatter is smooth; (3) the location of penumbra can be determined; and (4) the scatter-contaminated data provide knowledge about which part is smooth. Results: The simulation study shows that the edge-preserving weighting in OSE improves the estimation accuracy near the object boundary. Simulation study also demonstrates that OSE outperforms the two existing primary modulation algorithms for most regions of interest in terms of the CT number accuracy and noise. The proposed method was tested on a clinical cone beam CT, demonstrating that OSE corrects the scatter even when the modulator is not accurately registered. Conclusions: The proposed OSE algorithm improves the robustness and accuracy in scatter estimation and correction. This method is promising for scatter correction of various kinds of x-ray imaging modalities, such as x-ray radiography, cone beam CT, and the fourth-generation CT.

60 APPLIED LIFE SCIENCES↗

Metagenomic analysis of the microbiota in the highly compartmented hindguts of six wood- or soil-feeding higher termites

Termites are important contributors to carbon and nitrogen cycling in tropical ecosystems. Higher termites digest lignocellulose in various stages of humification with the help of an entirely prokaryotic microbiota housed in their compartmented intestinal tract. Previous studies revealed fundamental differences in community structure between compartments, but the functional roles of individual lineages in symbiotic digestion are mostly unknown. Furthermore, we conducted a highly resolved analysis of the gut microbiota in six species of higher termites that feed on plant material at different levels of humification. Combining amplicon sequencing and metagenomics, we assessed similarities in community structure and functional potential between the major hindgut compartments (P1, P3, and P4). Cluster analysis of the relative abundances of orthologous gene clusters (COGs) revealed high similarities among woodand litter-feeding termites and strong differences to humivorous species. However, abundance estimates of bacterial phyla based on 16S rRNA genes greatly differed from those based on protein-coding genes. In conclusion, the community structure and functional potential of the microbiota in individual gut compartments are clearly driven by the digestive strategy of the host. The metagenomics libraries obtained in this study provide the basis for future studies that elucidate the fundamental differences in the symbiont-mediated breakdown of lignocellulose and humus by termites of different feeding groups. Furthermore, the high proportion of uncultured bacterial lineages in all samples calls for a reference-independent approach for the correct taxonomic assignment of protein-coding genes.

59 BASIC BIOLOGICAL SCIENCES↗

An Extended Group Additivity Method for Polycyclic Thermochemistry Estimation

Automatic kinetic mechanism generation, virtual high-throughput screening, and automatic transition state search are currently trending applications requiring exploration of a large molecule space. Large-scale search requires fast and accurate estimation of molecules’ properties of interest, such as thermochemistry. Existing approaches are not satisfactory for large polycyclic molecules: considering the number of molecules being screened, quantum chemistry (even cheap density functional theory methods) can be computationally expensive, and group additivity, though fast, is not sufficiently accurate. Here, this paper provides a fast and moderately accurate alternative by proposing a polycyclic thermochemistry estimation method that extends the group additivity method with two additional algorithms: similarity match and bicyclic decomposition. It significantly reduces H f (298 K) estimation error from over 60 kcal/mol (group additivity method) to around 5 kcal/mol, C p (298 K) error from 9 to 1 cal/mol/K,and S(298 K) error from 70 to 7 cal/mol/K. This method also works well for heteroatomic polycyclics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Calculation and interpretation of classical turning surfaces in solids

Abstract Classical turning surfaces of Kohn–Sham potentials separate classically allowed regions (CARs) from classically forbidden regions (CFRs). They are useful for understanding many chemical properties of molecules but need not exist in solids, where the density never decays to zero. At equilibrium geometries, we find that CFRs are absent in perfect metals, rare in covalent semiconductors at equilibrium, but common in ionic and molecular crystals. In all materials, CFRs appear or grow as the internuclear distances are uniformly expanded. They can also appear at a monovacancy in a metal. Calculations with several approximate density functionals and codes confirm these behaviors. A classical picture of conduction suggests that CARs should be connected in metals, and disconnected in wide-gap insulators, and is confirmed in the limits of extreme compression and expansion. Surprisingly, many semiconductors have no CFR at equilibrium, a key finding for density functional construction. Nonetheless, a strong correlation with insulating behavior can still be inferred. Moreover, equilibrium bond lengths for all cases can be estimated from the bond type and the sum of the classical turning radii of the free atoms or ions.

Chemistry↗

Physically based correction of systematic errors of Rotating Shadowband Irradiometers

Accurate measurements of direct normal, diffuse horizontal and global horizontal irradiance (DNI, DHI and GHI) are needed for meteorological studies and are essential for the solar resource assessment at potential solar power plant sites. Often, these potential sites are remote and hence require robust sensors that require minimal maintenance that are not affected strongly by soiling. Therefore, Rotating Shadowband Irradiometers (RSI) are widely used for resource assessment. To achieve the required accuracy, corrections for the raw values of RSIs depending on systematic temperature, incidence angle and spectral errors must be used, and a thorough calibration of the sensor head must be applied. The existing correction functions are derived from comparisons of RSIs to thermopile radiometers at selected sites and therefore empirical. Their accuracy is considered to be site dependent. In this work a new correction and calibration method is presented that removes the systematic errors using a physical approach. It is based on information of the sensor properties as well as measurements of its directional response, and incorporates the atmospheric conditions at the measurement site. In this case, no empiric relations obtained from a specific site are required. The method requires estimates of the current DHI and GHI spectra during each measurement of the RSI. Based on these spectra, a spectral correction, which includes a spectrum dependent temperature correction, can be made without employing empirical relationships. The new physical calibration and correction method is tested at three sites and reaches similar results compared to the empirical functions. This is already achieved with rudimentary estimations of the GHI and DHI spectra and we expect that these estimations can be improved in the future. The results indicate that the physical approach reduces the problematic location dependence of the current calibration and correction methods. The physical correction and calibration method show promise for a further improvement of the RSI accuracy.

47 OTHER INSTRUMENTATION↗