Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model fitting”

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

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

At least 1,333 records · Page 74

P-V-T equation of state of boron carbide

We report the P-V-T equation of state measurements of B 4 C to 50 GPa and approximately 2500 K in laser-heated diamond anvil cells. We obtain an ambient temperature, third-order Birch–Murnaghan fit to the P-V data that yields a bulk modulus K 0 of 221(2) GPa and derivative, (dK/dP) 0 of 3.3(1). These were used in fits with both a Mie–Grüneisen–Debye model and a temperature-dependent, Birch–Murnaghan equation of state that includes thermal pressure estimated by thermal expansion (α) and a temperature-dependent bulk modulus (dK 0 /dT). The ambient pressure thermal expansion coefficient (α 0 + α 1 T), Grüneisen γ(V) = γ 0 (V/V 0 ) q and volume-dependent Debye temperature, were used as input parameters for these fits and found to be sufficient to describe the data in the whole P-T range of this study.

36 MATERIALS SCIENCE↗

Star cluster classification using deep transfer learning with PHANGS- HST

Currently available star cluster catalogues from the Hubble Space Telescope (HST) imaging of nearby galaxies heavily rely on visual inspection and classification of candidate clusters. Here, the time-consuming nature of this process has limited the production of reliable catalogues and thus also post-observation analysis. To address this problem, deep transfer learning has recently been used to create neural network models that accurately classify star cluster morphologies at production scale for nearby spiral galaxies (D ≲ 20 Mpc). Here, we use HST ultraviolet (UV)–optical imaging of over 20 000 sources in 23 galaxies from the Physics at High Angular resolution in Nearby GalaxieS (PHANGS) survey to train and evaluate two new sets of models: (i) distance-dependent models, based on cluster candidates binned by galaxy distance (9–12, 14–18, and 18–24 Mpc), and (ii) distance-independent models, based on the combined sample of candidates from all galaxies. We find that the overall accuracy of both sets of models is comparable to previous automated star cluster classification studies (~60–80 per cent) and shows improvement by a factor of 2 in classifying asymmetric and multipeaked clusters from PHANGS-HST. Somewhat surprisingly, while we observe a weak negative correlation between model accuracy and galactic distance, we find that training separate models for the three distance bins does not significantly improve classification accuracy. We also evaluate model accuracy as a function of cluster properties such as brightness, colour, and spectral energy distribution (SED)-fit age. Based on the success of these experiments, our models will provide classifications for the full set of PHANGS-HST candidate clusters (N ~ 200 000) for public release.

79 ASTRONOMY AND ASTROPHYSICS↗

Parameter Sensitivity Analysis of the SparTen High Performance Sparse Tensor Decomposition Software (Extended Analysis)

Tensor decomposition models play an increasingly important role in modern data science applications. One problem of particular interest is fitting a low-rank Canonical Polyadic (CP) tensor decomposition model when the tensor has sparse structure and the tensor elements are nonnegative count data. SparTen is a high-performance C++ library which computes a low-rank decomposition using different solvers: a first-order quasi-Newton or a second-order damped Newton method, along with the appropriate choice of runtime parameters. Since default parameters in SparTen are tuned to experimental results in prior published work on a single real-world dataset conducted using MATLAB implementations of these methods, it remains unclear if the parameter defaults in SparTen are appropriate for general tensor data. Furthermore, it is unknown how sensitive algorithm convergence is to changes in the input parameter values. This report addresses these unresolved issues with large-scale experimentation on three benchmark tensor data sets. Experiments were conducted on several different CPU architectures and replicated with many initial states to establish generalized profiles of algorithm convergence behavior.

97 MATHEMATICS AND COMPUTING↗

Re-evaluation of ortho-para- dependence of self pressure-broadening in the ν 1 + ν 3 band of acetylene

Optical frequency comb-referenced measurements of self pressure-broadened line profiles of the R(8) to R(13) lines in the ν 1 + ν 3 combination band of acetylene near 1.52 µ m are reported. The analysis of the data found no evidence for a previously reported systematic alternation in self pressure-broadened line widths with the nuclear spin state of the molecule. This work brought out the need for the use of an accurate line profile model and careful accounting for weak background absorptions due to hot band and lower abundance isotopomer lines. The data were adequately fit using the quadratic speed-dependent Voigt profile model, neglecting the small speed-dependent shift. Parameters describing the most probable and speed-dependent pressure-broadening, most probable shift, and the line strength were determined for each line. Detailed modeling of the results of Iwakuni et al. showed that their neglect of collisional narrowing due to the speed-dependent broadening term combined with the strongly absorbing data recorded and analyzed in transmission mode were the reasons for their results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Designing ITER motional Stark effect line shift (MSE-LS) spectrometers

As a part of ITER beam aided diagnostics, the design of Motional Stark Effect (MSE) diagnostic observing the emission from the Balmer-α line is underway. The physics of Stark splitting shows that the Stark manifold is polarization dependent, and the energy splitting results in a line shift proportional to the electric field. Due to the challenges of maintaining the calibration of the plasma facing mirrors in ITER, the conventional MSE polarimetry measurement technique is replaced with a spectral approach that is deemed more favorable in the ITER environment. The MSE line shift (LS) diagnostic is designed to quantify the Lorentz electric field magnitude by measuring the Stark manifold using visible spectroscopy. In the presence of large magnetic fields and high energy heating beams of 1 MeV, the expected Stark splitting is much larger than in typical devices. The MSE-LS design has unique challenges requiring careful consideration and modeling of its viewing geometry and photon budget. The MSE-LS approach on ITER is promising but has stringent demands on the allowable errors for the statistical and systematic fitting uncertainties. In this study, a full system model and numerical simulations of data for each sightline are completed. For a range of optical transmission fractions, photon noise analysis is conducted to determine the statistical uncertainties. This provides guidance on the spectrometer throughput, dispersion at the detector, optics, and other design choices. A conceptual design of a high throughput spectrometer with a volume phase transmission grating is presented.

Instruments & Instrumentation↗

A Primer on Dose-Response Data Modeling in Radiation Therapy

An overview of common approaches used to assess a dose response for radiation therapy–associated endpoints is presented, using lung toxicity data sets analyzed as a part of the High Dose per Fraction, Hypofractionated Treatment Effects in the Clinic effort as an example. Each component presented (eg, data-driven analysis, dose-response analysis, and calculating uncertainties on model prediction) is addressed using established approaches. Specifically, the maximum likelihood method was used to calculate best parameter values of the commonly used logistic model, the profile-likelihood to calculate confidence intervals on model parameters, and the likelihood ratio to determine whether the observed data fit is statistically significant. The bootstrap method was used to calculate confidence intervals for model predictions. Correlated behavior of model parameters and implication for interpreting dose response are discussed.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Lithium-Ion Battery Diagnostics Using Electrochemical Impedance via Machine-Learning

Diagnosing battery states such as health, state-of-charge, or temperature is crucial for ensuring the safety and reliability of electrochemical energy storage systems. While some states, such as temperature, may be measured using cheap sensors, accurate diagnosis of battery health metrics usually requires time-consuming performance measurements, making them infeasible for use in real-world operation. These health metrics can be measured during lab-testing and then estimated on-line using predictive life models or via state observer algorithms such as Kalman filters, but these predictive methods should be supplemented by actual measurement of battery health whenever possible to ensure reliability. Rapid measurement of battery health may be done by various types of fast diagnostic techniques such as electrochemical impedance spectroscopy (EIS), which can be performed in only a few minutes and require only a fraction of the energy and power needed for a full charge and discharge measurement. But there is a substantial challenge for estimating battery health using EIS data, as EIS is sensitive to cell temperature, state-of-charge, current, and resting time in addition to health. Thus, utilizing EIS data to predict battery capacity requires correcting for all these additional variables, a task that is extremely difficult to handle analytically. This talk utilizes machine-learning methods to estimate the effectiveness of battery capacity prediction from EIS data, leveraging a data set of hundreds of EIS measurements recorded at varying temperature and state-of-charge throughout a 500-day aging study of 32 commercial, large-format NMC-Graphite lithium-ion batteries. Using EIS as input to machine-learning models is complicated by the nonlinear response of impedance to battery health, temperature, and state-of-charge, as well as the collinearity between the impedance response at neighboring frequencies, which can easily lead to overfit models. To train robust models, features from EIS data need to be extracted from the data or some subset of critical frequencies selected. Many approaches for extracting and selecting features from EIS data from electrochemical analysis and machine-learning fields were identified for analysis: using the entire raw spectra; selection of one, two, or many frequencies from the entire spectra; selecting interesting points from the EIS measurement using domain knowledge; fitting EIS with an equivalent-circuit model; calculating statistics on the raw impedance values; and reducing the dimensionality of the data using unsupervised linear (principal component analysis) and non-linear (uniform manifold approximation and projection) methods. These approaches were rigorously compared using a machine-learning pipeline approach, training linear, Gaussian process, and random forest regression models and quantifying performance using cross-validation as well as a held-out test set. An artificial neural network model trained on the raw spectra was also tested. Promising pipelines were fine-tuned via Bayesian hyperparameter optimization using cross-validation loss and training with class-specific weights to counter data set imbalance. The most reliable method for utilizing impedance in this work was the selection of two optimal frequencies through an exhaustive search, resulting in about 2% mean absolute error on test data for both Gaussian process and random forest model architectures. Interrogation of a variety of models reveals critical frequencies of 100 Hz and 103 Hz for this data set, though the optimal set of frequencies is not necessarily intuitive, i.e., the best performing models are not simply those that use impedance at frequencies that have the highest correlation to the relative discharge capacity. The best performing model is an ensemble model, which is able to predict battery capacity with 1.9% mean absolute error for unseen cells using impedance recorded at a variety of temperatures and states-of-charge.

battery↗

HOD-dependent systematics for luminous red galaxies in the DESI 2024 BAO analysis

In this paper, we present the estimation of systematics related to the halo occupation distribution (HOD) modeling in the baryon acoustic oscillations (BAO) distance measurement of the Dark Energy Spectroscopic Instrument (DESI) 2024 analysis. This paper focuses on the study of HOD systematics for luminous red galaxies (LRG). We consider three different HOD models for LRGs, including the base 5-parameter vanilla model and two extensions to it, that we refer to as baseline and extended models. The baseline model is described by the 5 vanilla HOD parameters, an incompleteness factor and a velocity bias parameter, whereas the extended one also includes a galaxy assembly bias and a satellite profile parameter. We utilize the 25 dark matter simulations available in the AbacusSummit simulation suite at z=0.8 and generate mock catalogs for our different HOD models. To test the impact of the HOD modeling in the position of the BAO peak, we run BAO fits for all these sets of simulations and compare the best-fit BAO-scaling parameters α iso and α AP between every pair of HOD models. Furthermore, we do this for both Fourier and configuration spaces independently, using post-reconstruction measurements. We find a 3.3σ detection of HOD systematic for α AP in configuration space with an amplitude of 0.19%. For the other cases, we did not find a 3σ detection, and we decided to compute a conservative estimation of the systematic using the ensemble of shifts between all pairs of HOD models. By doing this, we quote a systematic with an amplitude of 0.07% in α iso for both Fourier and configuration spaces; and of 0.09% in α AP for Fourier space.

79 ASTRONOMY AND ASTROPHYSICS↗

Utilizing Temperature and Brine Inflow Measurements to Constrain Reservoir Parameters During a Salt Heater Test

Brine availability in salt has multiple implications for the safety and design of a nuclear waste storage facility. Brine availability includes both the distribution and transport of brine through a damaged zone around boreholes or drifts excavated into the salt. Coupled thermal, hydrological, mechanical, and chemical processes taking place within heated bedded salt are complex; as part of DECOVALEX 2023 Task E this study takes a parsimonious modeling approach utilizing analytical and numerical one-dimensional simulations to match field measurements of temperature and brine inflow around a heater. The one-dimensional modeling results presented arrive at best-fit thermal conductivity of intact salt, and the permeability and porosity of damaged salt of 5.74 W/m∙K, 10 –17 m 2 , and ≈ 0.02, respectively.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A model-independent data assimilation (MIDA) module and its applications in ecology

Abstract. Models are an important tool to predict Earth system dynamics. An accurate prediction of future states of ecosystems depends on not only model structures but also parameterizations. Model parameters can be constrained by data assimilation. However, applications of data assimilation to ecology are restricted by highly technical requirements such as model-dependent coding. To alleviate this technical burden, we developed a model-independent data assimilation (MIDA) module. MIDA works in three steps including data preparation, execution of data assimilation, and visualization. The first step prepares prior ranges of parameter values, a defined number of iterations, and directory paths to access files of observations and models. The execution step calibrates parameter values to best fit the observations and estimates the parameter posterior distributions. The final step automatically visualizes the calibration performance and posterior distributions. MIDA is model independent, and modelers can use MIDA for an accurate and efficient data assimilation in a simple and interactive way without modification of their original models. We applied MIDA to four types of ecological models: the data assimilation linked ecosystem carbon (DALEC) model, a surrogate-based energy exascale earth system model: the land component (ELM), nine phenological models and a stand-alone biome ecological strategy simulator (BiomeE). The applications indicate that MIDA can effectively solve data assimilation problems for different ecological models. Additionally, the easy implementation and model-independent feature of MIDA breaks the technical barrier of applications of data–model fusion in ecology. MIDA facilitates the assimilation of various observations into models for uncertainty reduction in ecological modeling and forecasting.

58 GEOSCIENCES↗

Self-Supervised T-GCN for Detection of Disturbance and Propagation in Power Grid

Urban power systems increasingly rely on dense sensing to monitor grid reliability, yet disturbance labels are scarce and events are rare. We present a self-supervised spatio-temporal method that detects, localizes, and characterizes grid frequency disturbances across urban areas using only unlabeled data. Our approach trains a tiny Temporal Graph Convolutional Network (T-GCN) to forecast per-site frequency residuals (deviation from 60 Hz). The sensor graph is constructed directly from signals using pre-event Pearson correlation with a cross-correlation lag penalty without geocoding. At inference, node-level anomalies are the model's forecast errors; region-level alarms arise from connected components of high-score nodes. We estimate disturbance propagation by computing per-node arrival times (first persistent exceedance), then fit a planar or time-of-arrival model to obtain direction, speed, and an epicenter proxy. With only three real events collected at decisecond resolution across U.S. cities, we evaluate the T-GCN and report time-to-detect, footprint size, and propagation consistency. We further show that short-window embeddings from the T-GCN's hidden states enable few-shot event-vs-background recognition via a simple prototypical classifier. Despite minimal data and no labels, our system yields fast, spatially coherent detection and interpretable propagation maps, offering a practical, lightweight pathway to city-scale grid resilience analytics.

Niu, Haoran [ORNL] (ORCID:0000000155228297)↗

Layers Can Be Deceiving: A Hopping Model for Small Molecule Diffusion in TATB Crystal

Sorption of small molecule gases in materials can play a significant role in their long‐term stability and compatibility within multi‐material assemblies. While many material properties of the insensitive high explosive TATB (1,3,5‐triamino‐2,4,6‐trinitrobenzene) are well understood, very little is known regarding its permeability to gases. TATB crystal exhibits a graphitic‐like layered packing structure that evokes a mental schema in which the layers form nanoscopic channels, but it is unclear whether this structure promotes gas transport. Here, we use molecular dynamics (MD) simulations to predict transport of small molecules through TATB single crystal. An approach to fit classical force fields is developed to model TATB interactions with H 2 O, He, Ne, and Ar, which is then combined with steered MD to probe gas transport along selected directions in the crystal. We find that small molecule transport occurs via a hopping mechanism that exhibits distinct jumps between interstitial sites and is substantially faster normal to the layers as compared to through them. This result stems from the finding that intralayer junctions between adjacent TATB molecules are the most stable interstitial sites and that energetic barriers are lower for hopping between adjacent layers. An empirical model for diffusion rate based on the MD data shows that the rate decays exponentially with increasing molecular radius and is negligibly small for all molecules larger than He, including common atmospheric gases. These findings have implications for the interpretation of experiments that measure surface area, material response to extreme conditions, and are expected to help constrain models for material aging.

36 MATERIALS SCIENCE↗

Physically interpretable approximations of many-body spectral functions

The rational function approximation provides a natural and interpretable representation of response functions such as the many-body spectral functions. We apply the vector fitting (VFIT) algorithm to fit a variety of spectral functions calculated from the Holstein model of electron-phonon interactions. We show that the resulting rational functions are highly efficient in their fitting of sharp features in the spectral functions, and could provide a means to infer physically relevant information from a spectral data set. The position of the peaks in the approximated spectral function are determined by the location of poles in the complex plane. Additionally, we developed a variant of VFIT that incorporates regularization to improve the quality of fits. With this procedure, we demonstrate it is possible to achieve accurate spectral function fits that vary smoothly as a function of physical conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The metabolic origins of non-photorespiratory CO 2 release during photosynthesis: a metabolic flux analysis

Respiration in the light (R L ) releases CO 2 in photosynthesizing leaves and is a phenomenon that occurs independently from photorespiration. Since R L lowers net carbon fixation, understanding R L could help improve plant carbon-use efficiency and models of crop photosynthesis. Although R L was identified more than 75 years ago, its biochemical mechanisms remain unclear. To identify reactions contributing to R L , we mapped metabolic fluxes in photosynthesizing source leaves of the oilseed crop and model plant camelina (Camelina sativa). We performed a flux analysis using isotopic labeling patterns of central metabolites during 13 CO 2 labeling time course, gas exchange, and carbohydrate production rate experiments. To quantify the contributions of multiple potential CO 2 sources with statistical and biological confidence, we increased the number of metabolites measured and reduced biological and technical heterogeneity by using single mature source leaves and quickly quenching metabolism by directly injecting liquid N 2 ; we then compared the goodness-of-fit between these data and data from models with alternative metabolic network structures and constraints. Our analysis predicted that R L releases 5.2 μmol CO2 g –1 FW h –1 of CO 2 , which is relatively consistent with a value of 9.3 μmol CO 2 g –1 FW h –1 measured by CO 2 gas exchange. The results indicated that ≤10% of R L results from TCA cycle reactions, which are widely considered to dominate R L . Further analysis of the results indicated that oxidation of glucose-6-phosphate to pentose phosphate via 6-phosphogluconate (the G6P/OPP shunt) can account for >% of CO 2 released by R L .

59 BASIC BIOLOGICAL SCIENCES↗

A new method for obtaining model-free viscoelastic material properties from atomic force microscopy experiments using discrete integral transform techniques

Viscoelastic characterization of materials at the micro- and the nanoscale is commonly performed with the aid of force–distance relationships acquired using atomic force microscopy (AFM). The general strategy for existing methods is to fit the observed material behavior to specific viscoelastic models, such as generalized viscoelastic models or power-law rheology models, among others. Here we propose a new method to invert and obtain the viscoelastic properties of a material through the use of the Z-transform, without using a model. We present the rheological viscoelastic relations in their classical derivation and their z-domain correspondence. We illustrate the proposed technique on a model experiment involving a traditional ramp-shaped force–distance AFM curve, demonstrating good agreement between the viscoelastic characteristics extracted from the simulated experiment and the theoretical expectations. We also provide a path for calculating standard viscoelastic responses from the extracted material characteristics. The new technique based on the Z-transform is complementary to previous model-based viscoelastic analyses and can be advantageous with respect to Fourier techniques due to its generality. Additionally, it can handle the unbounded inputs traditionally used to acquire force–distance relationships in AFM, such as ramp functions, in which the cantilever position is displaced linearly with time for a finite period of time.

36 MATERIALS SCIENCE↗

Assessing equation of state-independent relations for neutron stars with nonparametric models

Relations between neutron star properties that do not depend on the nuclear equation of state offer insights on neutron star physics and have practical applications in data analysis. Such relations are obtained by fitting to a range of phenomenological or nuclear physics equation of state models, each of which may have varying degrees of accuracy. In this study we revisit commonly used relations and reassess them with a very flexible set of phenomenological nonparametric equation of state models that are based on Gaussian processes. Our models correspond to two sets: equations of state which mimic hadronic models, and equations of state with rapidly changing behavior that resemble phase transitions. Here we quantify the accuracy of relations under both sets and discuss their applicability with respect to expected upcoming statistical uncertainties of astrophysical observations. We further propose a goodness-of-fit metric which provides an estimate for the systematic error introduced by using the relation to model a certain equation-of-state set. Overall, the nonparametric distribution is more poorly fit with existing relations, with the I–Love–Q relations retaining the highest degree of universality. Fits degrade for relations involving the tidal deformability, such as the binary-Love and compactness-Love relations, and when introducing phase transition phenomenology. For most relations, systematic errors are comparable to current statistical uncertainties under the nonparametric equation of state distributions.

79 ASTRONOMY AND ASTROPHYSICS↗

Determining circuit model parameters from operation data for PV system degradation analysis: $\mathrm{PVPRO}$

Physics-based circuit parameters like series and shunt resistance are essential to provide insights into the degradation status of photovoltaic (PV) arrays. However, calculating these parameters typically requires a full current-voltage characteristic (I-V curve), the acquisition of which involves specific measurement devices and costly methods. Thus, I-V curves of the PV system level are often not available. Here this paper proposes a methodology (PVPRO) to estimate these I-V curve parameters using only operation (string-level DC voltage and current) and weather data (irradiance and temperature). PVPRO first performs multi-stage data pre-processing to remove noisy data. Next, the time-series DC data are used to fit an equivalent circuit single-diode model (SDM) to estimate the circuit parameters by minimizing the differences between the measured and estimated values. In this way, the time evolutions of the SDM parameters are obtained. We evaluate PVPRO on synthetic datasets and find an excellent estimation of both SDM and the key I-V parameters (e.g., open-circuit voltage, short-circuit current, maximum power, etc.) with an average relative error of 0.55%. The performance, especially the extracted degradation rate of parameters, is robust to various measurement noises and the presence of faults. In addition, PVPRO is applied to a 271 kW PV field system. The relative error between the real and estimated operation voltage and current is less than 1%, suggesting that degradation trends are well captured. PVPRO represents a promising open-source tool to extract the time-series degradation trends of key PV parameters from routine operation data.

14 SOLAR ENERGY↗

Multiwavelength Observations of 2HWC J1928+177: Dark Accelerator or New TeV Gamma-Ray Binary?

2HWC J1928+177 is a Galactic TeV gamma-ray source detected by the High Altitude Water Cherenkov (HAWC) Observatory up to ~56 TeV. The HAWC source, later confirmed by the High Energy Stereoscopic System, still remains unidentified as a dark accelerator since there is no apparent supernova remnant or pulsar wind nebula detected in the lower-energy bands. The radio pulsar PSR J1928+1746, coinciding with the HAWC source position, has no X-ray counterpart. In this work, our SED modeling shows that inverse Compton scattering in the putative pulsar wind nebula can account for the TeV emission only if the unseen nebula is extended beyond $r\sim 4^{\prime} $. Alternatively, TeV gamma-rays may be produced by hadronic interactions between relativistic protons from an undetected supernova remnant associated with the radio pulsar and a nearby molecular cloud G52.9+0.1. NuSTAR and Chandra observations detected a variable X-ray point source within the HAWC error circle, potentially associated with a bright infrared (IR) source. The X-ray spectra can be fitted with an absorbed power-law model with N H = (9 ± 3) × 10 22 cm –2 and Γ X = 1.6 ± 0.3 and exhibit long-term X-ray flux variability over the last decade. If the X-ray source, possibly associated with the IR source (likely an O star), is the counterpart of the HAWC source, it may be a new TeV gamma-ray binary powered by collisions between the pulsar wind and stellar wind. Follow-up X-ray observations are warranted to search for diffuse X-ray emission and determine the nature of the HAWC source.

79 ASTRONOMY AND ASTROPHYSICS↗