Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model fitting”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29

Neural Posterior Estimation for Scalable and Accurate Inverse Parameter Inference in Li-Ion Batteries

Diagnosing the internal state of Li-ion batteries is critical for battery research, operation of real-world systems, and prognostic evaluation of remaining lifetime. By using physics-based models to perform probabilistic parameter estimation via Bayesian calibration, diagnostics can account for the uncertainty due to model fitness, data noise, and the observability of any given parameter. However, Bayesian calibration in Li-ion batteries using electrochemical data is computationally intensive even when using a fast surrogate in place of physics-based models, requiring many thousands of model evaluations. A fully amortized alternative is neural posterior estimation (NPE). NPE shifts the computational burden from the parameter estimation step to data generation and model training, reducing the parameter estimation time from minutes to milliseconds, enabling real-time applications. The present work shows that NPE can infer parameters equally or more accurately than Bayesian calibration, even if it leads to higher voltage reconstruction errors. We also demonstrate that the higher computational costs for data generation are tractable even in high-dimensional cases (ranging from 6 to 27 estimated parameters). The NPE method also offers several interpretability advantages over Bayesian calibration, such as local parameter sensitivity to specific regions of the voltage curve. The NPE method is demonstrated using an experimental fast charge dataset, with parameter estimates validated against measurements of loss of lithium inventory and loss of active material. The implementation is made available in a companion repository (https://github.com/NatLabRockies/BatFIT).

25 ENERGY STORAGE↗

Kinetic assessment of pulp mill-derived lime mud calcination in high CO 2 atmosphere

The chemical pulping of biomass involves the recycling of calcium through the calcination of lime mud, which is mostly comprised of calcium carbonate (CaCO 3 ). Lime mud decomposes under elevated temperatures to generate calcium oxide (CaO) and carbon dioxide (CO 2 ), the kinetics of which are strongly influenced by the CO 2 partial pressure and temperature. Oxy-fuel combustion and electrified lime kilns for lime mud calcination are intriguing methods to decarbonize this highly polluting operation within the biomass pulping industry. However, the high CO 2 concentration in oxy-fuel and electrified calcination processes alters the kinetics and overall reactivity of lime mud. For the first time, a model-fitting method is used to determine the kinetic parameters for lime mud calcination under a wide range of temperatures (550 °C–1250 °C) and under different concentrations of CO2 (0 %, 15 %, 50 %, and 90 %). A kinetic model is developed that accurately predicts the reaction rates as a function of temperature and CO 2 concentration. The apparent activation of energy for lime mud calcination is elevated under a high CO 2 environment. Relative to inert gas (N 2 , Ar), the temperature window for calcination is much smaller under high CO 2 environments. The presence of Na in lime mud does not seem to affect calcination under a high CO 2 environment. Finally, particle size variation does not have a significant effect on calcination under a high CO 2 environment.

09 BIOMASS FUELS↗

NEOWISE Reactivation Mission Year Three: Asteroid Diameters and Albedos

The Near-Earth Object Wide-field Infrared Survey Explorer (NEOWISE) reactivation mission has completed its third year of surveying the sky in the thermal infrared for near-Earth asteroids and comets. NEOWISE collects simultaneous observations at 3.4 and 4.6 μ m of solar system objects passing through its field of regard. These data allow for the determination of total thermal emission from bodies in the inner solar system, and thus the sizes of these objects. In this paper, we present thermal model fits of asteroid diameters for 170 NEOs and 6110 Main Belt asteroids (MBAs) detected during the third year of the survey, as well as the associated optical geometric albedos. We compare our results with previous thermal model results from NEOWISE for overlapping sample sets, as well as diameters determined through other independent methods, and find that our diameter measurements for NEOs agree to within 26% (1 σ ) of previously measured values. Diameters for the MBAs are within 17% (1 σ ). This brings the total number of unique near-Earth objects characterized by the NEOWISE survey to 541, surpassing the number observed during the fully cryogenic mission in 2010.

79 ASTRONOMY AND ASTROPHYSICS↗

The second-moment climatology of the GATE rain rate data

The first part of this paper presents the description of the GARP (Global Atmospheric Research Program) Atlantic Tropical Experiment (GATE) 1 rain-rate data and its two-dimensional spectral and correlation characteristics, which has made it possible to accomplish the following: to show the concentration of a significant power along the frequency axis in the spatiotemporal spectra; to detect a diurnal cycle (which has a range of variation of about 3.4-5.4 mm/n) as one of the sources of bias in the rain statistics of satellite data; to study the distinction between the north-south and east-west transport of spatial rain-rate field and character of its anisotropy; to evaluate the scales of the distinction between second-moment estimates associated with ground and satellite samples; and to determine the appropriate spatial and temporal scales of simple linear stochastic models fitted to averaged rain-rate fields. The second part of this paper is devoted to an analysis of the diffusion of the rain rate by establishing a relationship between the parameters of the multivariate autoregressive model and the coefficients of a diffusion equation. This analysis led to the use of rain data to estimate the rain advection velocity as well as other coefficients of the diffusion equation of the corresponding field. The results obtained can be used for comparison with corresponding estimates of other sources of data (satellite, Tropical Oceans Global Atmosphere Coupled Ocean - Atmosphere Response Experiment (TOGA, COARE) or simulated by physical models), for generating multiple samples of any size, for solving the inverse problems of some of the hydrodynamic equations, and in some other areas of rain data analysis and modeling.

Polyak, Ilya↗

PAHs as Probes of Photodissociation Regions in M17 and the Orion Bar

We have obtained narrow-band images of M17 SW and the Orion Bar in the PAH emission bands and pedestal (3.29, 3. 36, 3. 42 microns) The emission bands at 3.3 and 3.4 microns arise from the photodissociation regions (PDRs) between ionized gas and molecular clouds. In both M17 SW and the Orion Bar, the PDRs are nearly edge-on, providing excellent opportunities for comparing models of PDRs with observations. We observe an exponential drop in the strength of the 3.3 micron emission with a 1/e distance of 9 arcsec in Orion and 5 arcsec in M17 SW, in good agreement with previous observations. These results show that the two regions are very similar, and they imply that the mean density is 2.4 times higher in the Orion Bar than in the M17 SW PDR. However, we also find that in the Orion Bar, the ratio of the 3.4 micron emission to the 3.3 micron emission is consistent with the 1/e distance of 3 arcsec determined from PDR models fit to the molecular H and CO in the Orion Bar. We are presently investigating how the main band can imply that the UV field is dropping with a 1/e distance of 9 arcsec while the model PDR predicts a 1/e distance of 3 arcsec.

Sloan, G. C.↗

Gamma-ray Burst Afterglows as Probes of Environment and Blastwave Physics: Absorption by Host Galaxy Gas and Dust, Circumburst Media and the Distribution of P

We use a new approach to obtain limits on the absorbing columns towards a sample of 10 Gamma-ray Bursts observed by BeppoSAX from simultaneous fits to X-ray, optical and IR data, in counts space and including the effects of metallicity. For half the afterglows the best-fitting model to the SED includes SMC-like extinction (as opposed to LMC or MW) and in one LMC-like extinction, and in no cases is there a preference for MW-like extinction. Gas-to-dust ratios generally do not match those of the 3 standard and most well-known extinction models of SMC, LMC and MW, but tend to be higher. We compare the results from this method to those of previous works using other methods. We constrain the jet models for a subsample of the bursts by constraining the cooling break position and power law spectral slopes, allowing the injected electron energy index to be measured. We derive secure values of p from our spectral fits and comparison with the temporal optical and X-ray slopes for 4 afterglows. The mean of these single value, suggesting that either external factors such as circumburst medium play a strong role or that the microphysics is not identical for each GRB. For GRB 971214 we find that the circumburst medium has a wind-like density profile and the cooling frequency appears to be moving to higher frequencies.

Starling, R. L. C.↗

Comparison of simulated neutrino emission models with data on Supernova 1987A

Here we compare models of supernova (SN) neutrino emission with the Kamiokande II data on SN 1987A using the Bayesian approach. These models are taken from simulations and are representative of current one-dimensional SN models. We find that models with a brief accretion phase of neutrino emission are the most favored. This result is not affected by varying the overall flux normalization or considering neutrino oscillations. We also check the compatibility of the best-fit models with the data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Single-crystal elasticity of (Al,Fe)-bearing bridgmanite up to 82 GPa

Abstract Thermoelastic properties of mantle candidate minerals are essential to our understanding of geophysical phenomena, geochemistry, and geodynamic evolutions of the silicate Earth. However, the lower-mantle mineralogy remains much debated due to the lack of single-crystal elastic moduli (Cij) and aggregate sound velocities of (Al,Fe)-bearing bridgmanite, the most abundant mineral of the planet, at the lower mantle pressure-temperature (P-T) conditions. Here we report single-crystal Cij of (Al,Fe)-bearing bridgmanite, Mg0.88Fe0.1Al0.14Si0.90O3 (Fe10-Al14-Bgm) with Fe3+/ΣFe = ~0.65, up to ~82 GPa using X-ray diffraction (XRD), Brillouin light scattering (BLS), and impulsive stimulated light scattering (ISLS) measurements in diamond-anvil cells (DACs). Two crystal platelets with orientations of (–0.50, 0.05, –0.86) and (0.65, –0.59, 0.48), that are sensitive to deriving all nine Cij, are used for compressional and shear wave velocity (νP and νS) measurements as a function of azimuthal angles over 200° at each experimental pressure. Our results show that all Cij of singe-crystal Fe10-Al14-Bgm increase monotonically with pressure with small uncertainties of 1–2% (±1σ), except C55 and C23, which have uncertainties of 3–4%. Using the third-order Eulerian finite-strain equations to model the elasticity data yields the aggregate adiabatic bulk and shear moduli and respective pressure derivatives at the reference pressure of 25 GPa: KS = 326 ± 4 GPa, µ = 211 ± 2 GPa, KS′ = 3.32 ± 0.04, and µ′ = 1.66 ± 0.02 GPa. The high-pressure aggregate νS and νP of Fe10-Al14-Bgm are 2.6–3.5% and 3.1–4.7% lower than those of MgSiO3 bridgmanite end-member, respectively. These data are used with literature reports on bridgmanite with different Fe and Al contents to quantitatively evaluate pressure and compositional effects on their elastic properties. Comparing with one-dimensional seismic profiles, our modeled velocity profiles of major lower-mantle mineral assemblages at relevant P-T suggest that the lower mantle could likely consist of about 89 vol% (Al,Fe)-bearing bridgmanite. After considering uncertainties, our best-fit model is still indistinguishable from pyrolitic or chondritic models.

Geochemistry & Geophysics↗

Preliminary Results for Using Uncertainty and Out-of-distribution Detection to Identify Unreliable Predictions.

As machine learning (ML) models are deployed into an ever-diversifying set of application spaces, ranging from self-driving cars to cybersecurity to climate modeling, the need to carefully evaluate model credibility becomes increasingly important. Uncertainty quantification (UQ) provides important information about the ability of a learned model to make sound predictions, often with respect to individual test cases. However, most UQ methods for ML are themselves data-driven and therefore susceptible to the same knowledge gaps as the models themselves. Specifically, UQ helps to identify points near decision boundaries where the models fit the data poorly, yet predictions can score as certain for points that are under-represented by the training data and thus out-of-distribution (OOD). One method for evaluating the quality of both ML models and their associated uncertainty estimates is out-of-distribution detection (OODD). We combine OODD with UQ to provide insights into the reliability of the individual predictions made by an ML model.

97 MATHEMATICS AND COMPUTING↗

Thermal stresses, differential subsidence, and flexure at oceanic fracture zones

Geosat geoid undulations over four Pacific fracture zones have been analyzed. After correcting for the isostatic thermal edge effect, the amplitudes of the residuals are shown to be proportional to the age offset. The shape of the residuals seems to broaden with increasing age. Both geoid anomalies and available ship bathymetry data suggest that slip must sometimes occur on the main fracture zone or secondary faults. Existing models for flexure at fracture zones cannot explain the observed anomalies. A combination model accounting for slip and including flexure from thermal stresses and differential subsidence is presented. This model accounts for lateral variations in flexural rigidity from brittle and ductile yielding due to both thermal and flexural stresses and explains both the amplitudes and the shape of the anomalies along each fracture zone. The best fitting models have mechanical plate thicknesses that are described by the depth to the 600-700 C isotherms.

Wessel, Pal↗

Quantifying mean, variability, and uncertainty in indoor radon exposure in Pennsylvania using random forest and quantile regression forest models

Radon is a naturally occurring radioactive gas that poses a serious health risk as the primary cause of lung cancer in non-smokers. Despite the well-known adverse association with health outcomes, current radon exposure assessments are limited to county-level or average-level estimates, which fail to capture regional variability. This study uses Machine Learning models, including Random Forest (RF) and Quantile Regression Forest (QRF), to estimate the indoor radon concentrations at the ZCTA (Zip code tabulation area)-level and characterize uncertainties in model estimates. Incorporating geological, meteorological, and building-specific data, the models aim to improve radon risk assessment by capturing mean exposure, variability, and extreme concentration levels. Processed radon test data (n = 718,111) were analyzed using average, variability, and quantile prediction methods. Models that estimate the average radon exposure at the ZCTA-level can yield promising model-fit results, but they do not capture the underlying variability of indoor radon exposure within a ZCTA. We utilize volatility analyses to identify characteristics indicative of high variability of indoor radon exposure. We also show that a QRF model can be used to estimate upper quantiles of residential radon exposure, thereby uncovering localized areas of elevated exposure that were not apparent in mean estimates. The results highlighted the need for a deep characterization of exposure risk and show that regions with moderate average exposure levels could still harbor extreme outliers with implications for evaluating health risks. Utilizing multiple radon exposure models allows for a deeper characterization of radon risk within a geographic area and can better identify high-risk areas. The results from this study provide a foundation for developing mitigation strategies and examining associations between radon exposure and health outcomes at fine scales. Future research should extend the geographic scope and incorporate additional environmental risk factors to establish a comprehensive framework for risk assessment.

Lee, Heechan [ORNL]↗

Recovery of atmospheric refractivity profiles from simulated satellite-to-satellite tracking data

Techniques for recovering atmospheric refractivity profiles from simulated satellite-to-satellite tracking data are documented. Examples are given using the geometric configuration of the ATS-6/NIMBUS-6 Tracking Experiment. The underlying refractivity model for the lower atmosphere has the spherically symmetric form N = exp P(s) where P(s) is a polynomial in the normalized height s. For the simulation used, the Herglotz-Wiechert technique recovered values which were 0.4% and 40% different from the input values at the surface and at a height of 33 kilometers, respectively. Using the same input data, the model fitting technique recovered refractivity values 0.05% and 1% different from the input values at the surface and at a height of 50 kilometers, respectively. It is also shown that if ionospheric and water vapor effects can be properly modelled or effectively removed from the data, pressure and temperature distributions can be obtained.

Murray, C. W., Jr.↗

The AU MIC debris disk: Far-infrared and submillimeter resolved imaging

We present far-infrared and submillimeter maps from the Herschel Space Observatory and the James Clerk Maxwell Telescope of the debris disk host star AU Microscopii. Disk emission is detected at 70, 160, 250, 350, 450, 500, and 850 μm. The disk is resolved at 70, 160, and 450 μm. In addition to the planetesimal belt, we detect thermal emission from AU Mic’s halo for the first time. In contrast to the scattered light images, no asymmetries are evident in the disk. The fractional luminosity of the disk is 3.9×10{sup −4} and its milimeter-grain dust mass is 0.01 M{sub ⊕} (±20%). We create a simple spatial model that reconciles the disk spectral energy distribution as a blackbody of 53 ± 2 K (a composite of 39 and 50 K components) and the presence of small (non-blackbody) grains which populate the extended halo. The best-fit model is consistent with the “birth ring” model explored in earlier works, i.e., an edge-on dust belt extending from 8.8 to 40 AU, but with an additional halo component with an r{sup −1.5} surface density profile extending to the limits of sensitivity (140 AU). We confirm that AU Mic does not exert enough radiation force to blow out grains. For stellar mass-loss rates of 10–100 times solar, compact (zero porosity) grains can only be removed if they are very small; consistently with previous work, if the porosity is 0.9, then grains approaching 0.1 μm can be removed via corpuscular forces (i.e., the stellar wind).

79 ASTRONOMY AND ASTROPHYSICS↗

Generating Exploration Mission-3 Trajectories to a 9:2 NRHO using Machine Learning

The purpose of this thesis is to design a machine learning algorithm platform that provides expanded knowledge of mission availability through a launch season by improving trajectory resolution and introducing launch mission forecasting. The specific scenario addressed in this paper is one in which data is provided for four deterministic translational maneuvers through a mission to a Near Rectilinear Halo Orbit (NRHO) with a 9:2 synodic frequency. Current launch availability knowledge under NASA's Orion Orbit Performance Team is established by altering optimization variables associated to given reference launch epochs. This current method can bean abstract task and relies on an orbit analyst to structure a mission based on an established mission design methodology associated to the performance of Orion and NASA's Space Launch System. Introducing a machine learning algorithm trained to construct mission scenarios within the feasible range of known trajectories reduces the required interaction of the orbit analyst by removing the needed step of optimizing the orbit to fit an expected translational response required of the spacecraft. In this study, k-Nearest Neighbor and Bayesian Linear Regression successfully predicted classical orbital elements for the launch windows observed. However both algorithms had limitations due to their approaches to model fitting. Training machine learning algorithms off of classical orbital elements introduced a repetitive approach to reconstructing mission segments for different arrival opportunities through the launch window and can prove to be a viable method of launch window scan generation for future missions.

Guzman, Esteban↗

BatteryPro: A Python Toolkit for Battery Data Analysis and Machine Learning Predictions

Analyzing battery test data for research & development can be time-consuming since battery tests often run on the order of months to years, generating large volumes of data. BatteryPro is a comprehensive Python package and software designed to facilitate advanced analysis and performance predictions for battery test data. Developed for battery researchers, it supports data types from widely used battery testing instruments, including MACCOR and Biologic cycling systems. The software provides a variety of tools for extracting and plotting key battery parameters such as time, voltage, capacity, current, and pressure. In addition to its extensive data analysis capabilities, BatteryPro features a dedicated machine learning module that employs a Bayesian Gaussian Mixture Model (GMM) to predict battery performance and degradation. Users can generate synthetic capacity fade data, calculate fade metrics, and leverage predictive models to forecast long-term battery behavior. The software's graphical user interface (GUI) enhances usability, allowing researchers to upload, merge, and analyze multiple data files with full customizability. The GUI also supports machine learning predictions, enabling users to fit models and make predictions based on selected data and parameters. BatteryPro is built using QtDesigner, scikit-learn, matplotlib, and pandas, ensuring a high level of customization, flexibility, and accuracy in battery data analysis. This tool aims to empower researchers with the ability to perform detailed battery analysis and make informed predictions, ultimately advancing the field of battery research.

25 - ENERGY STORAGE↗

Parameterization of photon beam dosimetry for a linear accelerator

Purpose: In radiation therapy, accurate data acquisition of photon beam dosimetric quantities is important for (1) beam modeling data input into a treatment planning system (TPS), (2) comparing measured and TPS modeled data, (3) the quality assurance process of a linear accelerator’s (Linac) beam characteristics, (4) the establishment of a standard data set for comparison with other data, etcetera. Parameterization of the photon beam dosimetry creates a data set that is portable and easy to implement for different applications such as those previously mentioned. The aim of this study is to develop methods to parameterize photon beam dosimetric quantities, including percentage depth doses (PDDs), profiles, and total scatter output factors (S{sub cp}). Methods: S{sub cp}, PDDs, and profiles for different field sizes, depths, and energies were measured for a Linac using a cylindrical 3D water scanning system. All data were smoothed for the analysis and profile data were also centered, symmetrized, and geometrically scaled. The S{sub cp} data were analyzed using an exponential function. The inverse square factor was removed from the PDD data before modeling and the data were subsequently analyzed using exponential functions. For profile modeling, one halfside of the profile was divided into three regions described by exponential, sigmoid, and Gaussian equations. All of the analytical functions are field size, energy, depth, and, in the case of profiles, scan direction specific. The model’s parameters were determined using the minimal amount of measured data necessary. The model’s accuracy was evaluated via the calculation of absolute differences between the measured (processed) and calculated data in low gradient regions and distance-to-agreement analysis in high gradient regions. Finally, the results of dosimetric quantities obtained by the fitted models for a different machine were also assessed. Results: All of the differences in the PDDs’ buildup and the profiles’ penumbra regions were less than 2 and 0.5 mm, respectively. The differences in the low gradient regions were 0.20% ± 0.20% (<1% for all) and 0.50% ± 0.35% (<1% for all) for PDDs and profiles, respectively. For S{sub cp} data, all of the absolute differences were less than 0.5%. Conclusions: This novel analytical model with minimum measurement requirements was proved to accurately calculate PDDs, profiles, and S{sub cp} for different field sizes, depths, and energies.

61 RADIATION PROTECTION AND DOSIMETRY↗

Prediction Uncertainties beyond the Range of Experience: A Case Study in Inertial Confinement Fusion Implosion Experiments

Scientists often predict physical outcomes, e.g., experimental results, with the assistance of computer codes that, at their best, only coarsely approximate reality. Coarse predictions are challenging in critical part due to the multitude of seemingly arbitrary yet consequential decisions that must be made such as choice of relevant data, calibration of code parameters, and construction of empirical discrepancy forms. In this paper, we present a case study in the context of inertial confinement fusion (ICF) implosion experiments where extrapolative predictions are needed with quantified uncertainties. The purpose of this case study is to reflect relevant statistical methods, as applied to ICF model fitting and prediction, to document the numerous decisions that must be made in the prediction pipeline, to extend a complex example in extrapolation to the uncertainty quantification (UQ) community, and to reflect on the challenges we encountered supporting extrapolations with imperfect models and thereby recommend several future research directions. We conclude with a discussion about the UQ community's role in less than ideal predictive scenarios like our ICF exercise.

97 MATHEMATICS AND COMPUTING↗

SU-E-T-74: A Simple Model of Dose Falloff for Highly Convergent SBRT Beams

Purpose: The purpose of this study was to create a simple model of dose falloff outside the target for a SBRT system employing a highly convergent beam geometry. With an estimate of the rate of reduction in dose away from a target of a given size, dose to tissues at a given distance from the target can be estimated. Our hypothesis is that dose at a given distance away from the edge of a target can be approximated using an inverse square model. Methods: With the GammaPod™ system used for this study, thirty-six rotating Co-60 beams with circular cross section converge at an isocenter and then diverge distal to the isocenter. Dose falloff beyond the target is mainly due to reduction in beam overlap, which in turn depends on projected area of all beams and projected beam area swept out during rotation. Beams also decrease in intensity beyond the isocenter. The GammaPod™ TPS calculates dose via superposition of Monte Carlo-generated kernels, each corresponding to the convergence of beams at an isocenter within a stereotactic breast cup. Results: The factors above yield a simple inverse square-based model that has been compared with TPS predictions and measurements for a variety of targets. Good agreement is seen considering the simple nature of the model. Fitting of the treatment planning and measured dose falloff data with a power law expression yields a negative exponent close to 2, confirming the inverse square nature of dose falloff. Conclusion: It is possible to apply a simple inverse square dose falloff model for SBRT or SRS systems such as the GammaPod which employ a highly-convergent beam geometry. This can be useful in estimating the “worst case scenario” for dose to critical organs at a known distance from the edge of the target volume. Authors are employees of Xcision Medical Systems.

60 APPLIED LIFE SCIENCES↗