Search NASA⌕ Search

SEARCH · Search NASA

Results for “Data 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 73 records · Page 4

Relationships between PEMFC Cathode Kinetic Losses and Contaminants’ Dipole Moment and Adsorption Energy on Pt

A database summarizing the effects of 21 contaminants on the performance of proton exchange membrane fuel cells (PEMFCs) was used to examine relationships between cathode kinetic losses and contaminant physicochemical parameters. Impedance spectroscopy data were employed to obtain oxygen reduction kinetic resistances by fitting data in the 10-158 Hz range to a simplified equivalent circuit. The contaminant dipole moment and the adsorption energy of the contaminant on a Pt surface were chosen as parameters. Dipole moments did not correlate with dimensionless cathode kinetic resistances. In contrast, adsorption energies were quantitatively and linearly correlated with minimum dimensionless cathode kinetic resistances. Contaminants influence the oxygen reduction for contaminant adsorption energies smaller than -24.5 kJ mol -1 , a value near the high limit of the adsorption energy of O 2 on Pt. Dimensionless cathode kinetic resistances linearly increase with decreasing O 2 adsorption energies below -24.5 kJ mol -1 . Measured total cell voltage losses are mostly larger than the cathode kinetic losses calculated from kinetic resistance changes, which indicates the existence of other sources of performance degradation. Modifications to the experimental procedure are proposed to ensure that data are comparable on a similar basis and improve the correlation between contaminant adsorption energy and kinetic cell voltage losses.

25 ENERGY STORAGE↗

Quantifying uncertainties in crystal electric field Hamiltonian fits to neutron data

We systematically examine uncertainties from fitting rare earth single-ion crystal electric field (CEF) Hamiltonians to inelastic neutron scattering data. Using pyrochlore and delafossite structures as test cases, we find that uncertainty in CEF parameters can be large despite visually excellent fits. These results show Yb ^{3+} 3 + compounds have particularly large g g -tensor uncertainty because of the few available peaks. In such cases, additional constraints are necessary for meaningful fits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Empirical and mechanistic transient fission gas release model for high-burnup LOCA conditions

In response to the nuclear industry desire to extend burnup beyond current licensing practices, the US Nuclear Regulatory Commission (NRC) released a research information letter (RIL) that provides a basis for analyzing fuel fragmentation, relocation, and dispersal (FFRD) in light-water reactors. Of the five elements discussed, the most ambiguous is the significance of transient fission gas (FG) release (FGR) (tFGR) and its effects on fuel performance under loss-of-coolant accident conditions. In fresh fuel, FG migration and eventual release is primarily governed by diffusion-based mechanisms at higher temperatures (>1,000 °C). However, the mechanisms governing FGR changes as burnup increases. More recent research indicates that FGR increases as burnup increases, specifically under temperature transient conditions, and this release occurs at lower temperatures with a new release mechanism. This behavior has been attributed to microcracking in the fuel and is likely related to microstructure embrittlement with the presence of over pressurized FG bubbles. The NRC RIL outlines the complexity of the phenomenon and a need for a deeper understanding to adequately address FFRD for regulatory application. Here, therefore, this manuscript intends to summarize the publicly available tFGR data and discuss the observed dependencies (e.g., burnup, heating rate, sample geometry, terminal temperature). An empirical model has been developed and benchmarked against recently published experimental data. However, this empirical model is limited to conditions for which fitting data exist and, therefore, a high-level discussion is included to provide a roadmap for atomistically-informed multiscale modeling in conjunction with experimental data collection to develop a mechanistic tFGR model widely applicable to a broad range of nuclear fuel conditions at high burnup.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Quantitative analysis of total reflection X-ray fluorescence from finely layered structures using XeRay

Total reflection x-ray fluorescence (TXRF) is a widely applicable experimental technique for studying chemical element distributions across finely layered structures at extremely high sensitivity. To promote and facilitate scientific discovery using TXRF, we developed a MATLAB-based software package with a graphical user interface, named XeRay, for quick, accurate, and intuitive data analysis. XeRay lets the user model any layered system, each layer with its independent chemical composition and thickness, and enables fine-tuned data fitting. The accuracy of XeRay has been tested in the analysis of TXRF data from both air/liquid interface and liquid/liquid interfacial studies and has been compared to literature results. In an air/liquid interface study, Ca 2+ sequestration was measured at a Langmuir monolayer of 1-stearoyl-2-oleoyl-sn-glycero-3-phosphatidic acid (SOPA) on a buffer solution of 1 mM CaCl2 at pH 7.5. Data analysis with XeRay reveals that each 1 nm2 of interfacial area contains 2.38 ± 0.06 Ca 2+ ions, which corresponds to a 1:1 ratio between SOPA headgroups and Ca 2+ ions, consistent with several earlier reports. For the liquid/liquid interface study of Sr 2+ enrichment at the dodecane/surfactant/water interface, analysis using XeRay gives a surface enrichment of Sr 2+ at $68^{+6}_{-5}$ Å 2 per ion, consistent with the result published for the same dataset.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A consistent orbital stability analysis for the GJ 581 system

We apply a combination of N-body modeling techniques and automated data fitting with Monte Carlo Markov Chain uncertainty analysis of Keplerian orbital models to RV data to determine long-term stability of the planetary system GJ 581. We find that while there are stability concerns with the four-planet model as published by Forveille et al., when uncertainties in the system are accounted for, particularly stellar jitter, the hypothesis that the four-planet model is gravitationally unstable is not statistically significant. Additionally, the system including proposed planet g by Vogt et al. also shows some stability concerns when eccentricities are allowed to float in the orbital fit, yet when uncertainties are included in the analysis, the system including planet g also cannot be proven to be unstable. We present revised reduced χ{sup 2} values for Keplerian astrocentric orbital fits assuming four-planet and five-planet models for GJ 581 under the condition that best fits must be stable, and we find no distinguishable difference by including planet g in the model. Additionally, we present revised orbital element estimates for each, assuming uncertainties due to stellar jitter under the constraint of the system being gravitationally stable.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluation of models determined by neutron diffraction and proposed improvements to their validation and deposition

The Protein Data Bank (PDB) contains a growing number of models that have been determined using neutron diffraction or a hybrid method that combines X-ray and neutron diffraction. The advantage of neutron diffraction experiments is that the positions of all atoms can be determined, including H atoms, which are hardly detectable by X-ray diffraction. This allows the determination of protonation states and the assignment of H atoms to water molecules. Because neutrons are scattered differently by hydrogen and its isotope deuterium, neutron diffraction in combination with H/D exchange can provide information on accessibility, dynamics and chemical lability. In this study, the deposited data, models and model-to-data fit for all PDB entries that used neutron diffraction as the source of experimental data have been analysed. In many cases, the reported R work and R free values were not reproducible. In such cases, the model and data files were analysed to identify the reasons for this mismatch. The issues responsible for the discrepancies are summarized and explained. The analysis unveiled limitations to the annotation, deposition and validation of models and data, and a lack of community-wide accepted standards for the description of neutron models and data, as well as deficiencies in current model refinement tools. Most of the issues identified concern the handling of H atoms. Since the primary use of neutron macromolecular crystallography is to locate and directly visualize H atoms, it is important to address these issues, so that the deposited neutron models allow the retrieval of the maximum amount of information with the smallest effort of manual intervention. A path forward to improving the annotation, validation and deposition of neutron models and hybrid X-ray and neutron models is suggested.

59 BASIC BIOLOGICAL SCIENCES↗

Studies on Thorium Adsorption Characteristics upon Activated Titanium Hydroxide Prepared from Rosetta Ilmenite Concentrate

The titanium hydroxide prepared from Rosetta ilmenite concentrate has been applied for Th (IV) adsorption from its acid aqueous solutions. The prepared hydroxide is first characterized by both Fourier transform infrared (FT-IR) spectrum and thermogravimetric analysis. The relevant factors affecting the adsorption process have been studied. The obtained equilibrium data fits well with the Langmuir isotherm rather than Freundlich isotherm, while the adsorption kinetic data follow the pseudo-second order model. The different thermodynamic parameters have also been calculated and indicate that the adsorption process is spontaneous.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An Edge Alignment-Based Orientation Selection Method for Neutron Tomography

Neutron computed tomography (nCT) is a 3D char-acterization technique used to image the internal morphology or chemical composition of samples in biology and materials sciences. A typical workflow involves placing the sample in the path of a neutron beam, acquiring projection data at a predefined set of orientations, and processing the resulting data using an analytic reconstruction algorithm. Typical nCT scans require hours to days to complete and are then processed using conventional filtered back-projection (FBP), which performs poorly with sparse views or noisy data. Hence, the main methods in order to reduce overall acquisition time are the use of an improved sampling strategy combined with the use of advanced reconstruction methods such as model-based iterative reconstruction (MBIR). In this paper, we propose an adaptive orientation selection method in which an MBIR reconstruction on previously-acquired measurements is used to define an objective function on orientations that balances a data-fitting term promoting edge alignment and a regularization term promoting orientation diversity. Using simulated and experimental data, we demonstrate that our method produces high-quality reconstructions using significantly fewer total measurements than the conventional approach.

Yang, Diyu↗

Downstream Separation of Formic Acid with Anion-Exchange Resin from Electrocatalytic Carbon Dioxide (CO 2 ) Conversion: Adsorption, Kinetics, and Equilibrium Modeling

The electrocatalytic carbon dioxide reduction reaction (CO 2 RR) that produces liquid formic acid (FA) is one potential route to minimize the CO 2 concentration in the atmosphere by utilizing CO 2 as a feedstock. However, one of the main disadvantages of this route is the high energy demand required for purification through conventional distillation. This method of purification has been shown to result in a large carbon footprint in life-cycle analysis (LCA). Therefore, finding an alternative purification technique that can integrate with an electrochemical CO 2 RR process with the potential to lower the carbon footprint would be of great interest. Under the experimental conditions of the CO 2 RR, FA remains in its basic form as the formate anion in the catholyte. Anion exchange resins with their basic properties can separate the formate anion from the catholyte mixture as an alternative to energy intensive distillation. In this work, three anion exchange resins with increasing basicity order: Amberlite IRA-96 < Amberlite IRA-910 < Ambersep 900 are tested to separate FA in the presence of aqueous solution of potassium bicarbonate catholyte with varying levels of resin and FA. Kinetics and equilibrium studies data for the FA adsorption are interpreted using several kinetics and isotherm models. The kinetics data fit well with a pseudo-first-order model at high initial FA concentrations and a pseudo-second-order model at low initial FA concentrations. The experimental data can be best explained with the Freundlich isotherm model. The fitted results show that strongly basic Ambersep 900 is more effective at separating FA, with a maximum adsorption capacity of 336.7 mg/g compared to Amberlite IRA-96 (275.2 mg/g) and Amberlite IRA-910 (209.2 mg/g) in neutral water. However, weakly basic Amberlite IRA-96 (110.8 mg/g) can separate FA more efficiently from potassium bicarbonate catholyte than Amberlite IRA-910 (99.9 mg/g) and Ambersep 900 (40.0 mg/g). As a result, competitive adsorption occurs in the presence of potassium bicarbonate catholyte, but the formic acid adsorption capacity is still sufficient for adequate separation.

20 FOSSIL-FUELED POWER PLANTS↗

An improved thermodynamic model for the complexation of trivalent actinides and lanthanide with oxalic acid valid to high ionic strength.

The dissociation constants of oxalic acid (Ox), and the stability constants of Am 3+ , Cm 3+ and Eu 3+ with Ox 2– have been determined at 25 °C, over a range of concentration varying from 0.1 to 6.60 m NaClO4 using potentiometric titration and extraction techniques, respectively. The experimental data support the formation of complexes, M(Ox) n 3 – 2n , where (M = Am 3+ , Cm 3+ and Eu 3+ and n = 1 and 2). The dissociation constant and the stability constant values measured as a function of NaClO 4 concentration were used to estimate the Pitzer parameters for the respective interactions of Am 3+ , Cm 3+ and Eu 3+ with Ox. Furthermore, the stability constants data of Am 3+ –Ox measured in NaClO 4 and in NaCl solutions from the literature were simultaneously fitted in order to refine the existing actinide–oxalate complexation model that can be used universally in the safety assessment of radioactive waste disposal. The thermodynamic stability constant: log β 0 101 = 6.30 ± 0.06 and log β 0 102 = 10.84 ± 0.06 for Am 3+ was obtained by simultaneously fitting data in NaCl and NaClO 4 media. Additionally, log β 0 101 = 6.72 ± 0.08 and log β 0 102 = 11.05 ± 0.09 for the Cm 3+ and log β 0 101 = 6.67 ± 0.08 and log β 0 102 = 11.15 ± 0.09 for the Eu 3+ were calculated by extrapolation of data to zero ionic strength in NaClO 4 medium only. For all stability constants, the Pitzer model gives an excellent representation of the data using interaction parameters β (0) , β (1) , and CΦ determined in this work. The thermodynamic model developed in this work will be useful in accurately modeling the potential solubility of trivalent actinides and early lanthanides to ionic strength of 6.60 m in low temperature environments in the presence of Ox. Furthermore, the work is also applicable to the accurate modeling transport of rare earth elements in various environments under the surface conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Paired autoencoders for likelihood-free estimation in inverse problems

Abstract We consider the solution of nonlinear inverse problems where the forward problem is a discretization of a partial differential equation. Such problems are notoriously difficult to solve in practice and require minimizing a combination of a data-fit term and a regularization term. The main computational bottleneck of typical algorithms is the direct estimation of the data misfit. Therefore, likelihood-free approaches have become appealing alternatives. Nonetheless, difficulties in generalization and limitations in accuracy have hindered their broader utility and applicability. In this work, we use a paired autoencoder framework as a likelihood-free estimator (LFE) for inverse problems. We show that the use of such an architecture allows us to construct a solution efficiently and to overcome some known open problems when using LFEs. In particular, our framework can assess the quality of the solution and improve on it if needed. We demonstrate the viability of our approach using examples from full waveform inversion and inverse electromagnetic imaging.

Chung, Matthias (ORCID:0000000178224539)↗

Acclimation of leaf respiration consistent with optimal photosynthetic capacity

Abstract Plant respiration is an important contributor to the proposed positive global carbon‐cycle feedback to climate change. However, as a major component, leaf mitochondrial (‘dark’) respiration ( R d ) differs among species adapted to contrasting environments and is known to acclimate to sustained changes in temperature. No accepted theory explains these phenomena or predicts its magnitude. Here we propose that the acclimation of R d follows an optimal behaviour related to the need to maintain long‐term average photosynthetic capacity ( V cmax ) so that available environmental resources can be most efficiently used for photosynthesis. To test this hypothesis, we extend photosynthetic co‐ordination theory to predict the acclimation of R d to growth temperature via a link to V cmax , and compare predictions to a global set of measurements from 112 sites spanning all terrestrial biomes. This extended co‐ordination theory predicts that field‐measured R d and V cmax accessed at growth temperature ( R d,tg and V cmax,tg ) should increase by 3.7% and 5.5% per degree increase in growth temperature. These acclimated responses to growth temperature are less steep than the corresponding instantaneous responses, which increase 8.1% and 9.9% per degree of measurement temperature for R d and V cmax respectively. Data‐fitted responses proof indistinguishable from the values predicted by our theory, and smaller than the instantaneous responses. Theory and data are also shown to agree that the basal rates of both R d and V cmax assessed at 25°C ( R d,25 and V cmax,25 ) decline by ~4.4% per degree increase in growth temperature. These results provide a parsimonious general theory for R d acclimation to temperature that is simpler—and potentially more reliable—than the plant functional type‐based leaf respiration schemes currently employed in most ecosystem and land‐surface models.

Wang, Han↗

Improved Parametric Models for Explosion Pressure Signals Derived From Large Datasets

Accurate recording and characterization of explosion-induced pressure signals are key components of the forensic analysis of explosion events in the atmosphere. Parametric overpressure models based on several key waveform features (peak overpressure, positive pulse duration, and impulse) are widely used to estimate explosion energy in terms of trinitrotoluene equivalent yield. However, those models are often developed by a limited dataset, including only a few events or recordings at relatively short propagation distances. Here, we develop empirical waveform-parameter models based on a regression analysis of a large set of data curated from four chemical explosion experiments including 16 detonations. We measured peak overpressure and impulse for positive and negative phases from 1000 pressure signals recorded at local ranges (<20 km) with scaled distance up to 8000 m/kg 1/3 . Additionally, the measured waveform parameters showed large variation with respect to observing distances indicating the effects of atmospheric propagation. In this study, a second-order polynomial model was used in a least-squares regression to account for those propagation effects and to improve data fitting. In addition to model parameters for waveform features, we also determined range-dependent model uncertainties based on data variance. The model uncertainty represents the prediction error of our models and can be critical to evaluating the uncertainty of yield estimate.

58 GEOSCIENCES↗

Real-time data reduction at 100 Tbps: Challenge and opportunity for AI-based data reduction for next-generation large-scale nuclear physics collider experiment

The modern large-scale nuclear physics (NP) experiments in high-energy particle colliders utilize streaming-readout electronics to digitize detector response at O(100) Tbps bandwidth. Prominent examples at Brookhaven National Lab (BNL) include the sPHENIX experiment at Relativistic Heavy Ion Collider (RHIC), which is under construction, and the experiments proposed for the Electron-Ion Collider (EIC), planned for the 2030s . One of the main challenges for these streaming readout systems is to manage the data rate with sufficient data reduction in real time so the end-data fit persistent storage for offline analysis, which is typically at O(1000) times smaller and O(100) Gbps. Such data reduction traditionally is achieved via real-time high level triggers, which select and save a small subset of collisions of interest. Although triggering is applicable to high energy collider experiments such as those at the Large Hardron Collider at CERN, it is insufficient for these nuclear physics experiments which study diverse collision topologies. And traditional triggering approach is inefficient to preserve the max information harvested from the operation of colliders that costs O(100)M per year to DOE. Meanwhile, in recent years, ML-based high-throughput data reduction has emerged as a promising approach to efficiently preserve max information for a given space of persistent storage, e.g. via AI data compression, feature extraction, and noise filtering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

SANS quantification of bound water in water-soluble polymers across multiple concentration regimes

Contrast-variation small-angle neutron scattering (CV-SANS) is a widely used technique for quantifying hydration water in soft matter systems, but it is predominantly applied in the dilute regime or for systems with a well-defined structure factor. Here, CV-SANS was used to quantify the number of hydration water molecules associating with three water-soluble polymers with different critical solution temperatures and types of water–solute interactions in dilute, semidilute, and concentrated solution through the exploration of novel methods of data fitting and analysis. Multiple SANS fitting workflows with varying levels of model assumptions were evaluated and compared to give insight into SANS model selection. These fitting pathways ranged from general, model-free algorithms to more standard form and structure factor fitting. In addition, Monte Carlo bootstrapping was evaluated as a method to estimate parameter uncertainty through simulation of technical replicates. The most robust fitting workflow for dilute solutions was found to be form factor fitting without CV-SANS ( i.e. polymer in 100% D 2 O). For semidilute and concentrated solutions, while the model-free approach can be mathematically defined for CV-SANS data, the addition of a structure factor imposes physical constraints on the optimization problem, suggesting that the optimal fitting pathway should include appropriate form and structure factor models. The measured hydration numbers were consistent with the number of tightly bound water molecules associated with each monomer unit, and the concentration dependence of the hydration number was largely governed by the chemistry-specific interactions between water and polymer. Polymers with weaker water–polymer interactions ( i.e. those with fewer hydration water molecules) were found to have more bound water at higher concentrations than those with stronger water–polymer interactions due to the increase in the number of forced water–polymer contacts in the concentrated system. This SANS-based method to count hydration water molecules can be applied to polymers in any concentration regime, which will lead to improved understanding of water–polymer interactions and their impact on materials design.

36 MATERIALS SCIENCE↗

Strong Upper-Plate Heterogeneity at the Hikurangi Subduction Margin (North Island, New Zealand) Imaged by Adjoint Tomography

We use earthquake-based adjoint tomography to invert for three-dimensional structure of the North Island, New Zealand, and the adjacent Hikurangi subduction zone. The study area, having a shallow depth to the plate interface below the North Island, offers a rare opportunity for imaging material properties at an active subduction zone using land-based measurements. Starting from an initial model derived using ray tomography, we perform iterative model updates using spectral element and adjoint simulations to fit waveforms with periods ranging from 4–30 s. We perform 28 model updates using an L-BFGS optimization algorithm, improving data fit and introducing P- and S-wave velocity changes of up to ±30%. Resolution analysis using point spread functions show that our measurements are most sensitive to heterogeneities in the upper 30 km. The most striking velocity changes coincide with areas related to the active Hikurangi subduction zone. Lateral velocity structures in the upper 5 km correlate well with New Zealand geology. Here, the inversion reveals increased along-strike heterogeneity on the margin. In Cook Strait we observe a low-velocity zone interpreted as deep sedimentary basins. In the central North Island, low-velocity anomalies are linked to surface geology, and we relate velocity structures at depth to crustal magmatic activity below the Taupō Volcanic Zone.

58 GEOSCIENCES↗

A detailed analysis of the high-resolution X-ray spectra of NGC 3516: variability of the ionized absorbers

The 1.5 Seyfert galaxy NGC 3516 presents a strong time variability in X-rays. We re-analyzed the nine observations performed in 2006 October by XMM-Newton and Chandra in the 0.3 to 10 keV energy band. An acceptable model was found for the XMM-Newton data fitting the EPIC-PN and RGS spectra simultaneously; later, this model was successfully applied to the contemporary Chandra high-resolution data. The model consists of a continuum emission component (power law + blackbody) absorbed by four ionized components (warm absorbers), and 10 narrow emission lines. Three absorbing components are warm, producing features only in the soft X-ray band. The fourth ionization component produces Fe XXV and Fe XXVI in the hard-energy band. We study the time response of the absorbing components to the well-detected changes in the X-ray luminosity of this source and find that the two components with the lower ionization state show clear opacity changes consistent with gas close to photoionization equilibrium. These changes are supported by the models and by differences in the spectral features among the nine observations. On the other hand, the two components with higher ionization state do not seem to respond to continuum variations. The response time of the ionized absorbers allows us to constrain their electron density and location. We find that one component (with intermediate ionization) must be located within the obscuring torus at a distance 2.7 × 10{sup 17} cm from the central engine. This outflowing component likely originated in the accretion disk. The three remaining components are at distances larger than 10{sup 16}-10{sup 17} cm. Two of the absorbing components in the soft X-rays have similar outflow velocities and locations. These components may be in pressure equilibrium, forming a multi-phase medium, if the gas has metallicity larger than the solar one (≳ 5 Z {sub ☉}). We also search for variations in the covering factor of the ionized absorbers (although partial covering is not required in our models). We find no correlation between the change in covering factor and the flux of the source. This, in connection with the observed variability of the ionized absorbers, suggests that the changes in flux are not produced by this material. If the variations are indeed produced by obscuring clumps of gas, these must be located much closer in to the central source.

79 ASTRONOMY AND ASTROPHYSICS↗