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.

SourceResearch connectionFrom yearThrough year

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

At least 1,603 records · Page 89

All sources NASA NTRS DOE OSTI

DOE OSTI2022

TDCOSMO VIII. A key test of systematics in the hierarchical method of time-delay cosmography

The largest source of systematic errors in the time-delay cosmography method likely arises from the lens model mass distribution, where an inaccurate choice of model could in principle bias the value of $H$ 0 . A Bayesian hierarchical framework has been proposed which combines lens systems with kinematic data, constraining the mass profile shape at a population level. The framework has been previously validated using a small sample of lensing galaxies drawn from hydro-simulations. The goal of this work is to expand the validation to a more general set of lenses consistent with observed systems, as well as confirm the capacity of the method to combine two lens populations: one which has time delay information and one which lacks time delays and has systematically different image radii. For this purpose, we generated samples of analytic lens mass distributions made of baryons+dark matter and fit the subsequent mock images with standard power-law models. Corresponding kinematics data were also emulated. The hierarchical framework applied to an ensemble of time-delay lenses allowed us to correct the $H$ 0 bias associated with model choice to find $H$ 0 within 1.5σ of the fiducial value. We then combined this set with a sample of corresponding lens systems which have no time delays and have a source at lower $z$, resulting in a systematically smaller image radius relative to their effective radius. The hierarchical framework has successfully accounted for this effect, recovering a value of $H$ 0 which is both more precise (σ ~ 2%) and more accurate (0.7% median offset) than the time-delay set alone. This result confirms that non-time-delay lenses can nonetheless contribute valuable constraining power to the determination of $H$ 0 via their kinematic constraints, assuming they come from the same global population as the time-delay set.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Analytical gradient-based optimization of CALPHAD model parameters

The calibration of CALPHAD (CALculation of PHAse Diagrams) models involves the solution of a very challenging high-dimensional multiobjective optimization problem. Traditional approaches to parameter fitting predominantly rely on gradient-free methods, which while robust, are computationally inefficient and often scale poorly with model complexity. In this work, we introduce and demonstrate a generalizable framework for analytic gradient-based optimization of the parameters of the CALPHAD model enabled by the recently formalized Jansson derivative technique. This method allows for efficient evaluation of gradients of thermodynamic properties at equilibrium with respect to model parameters, even in the presence of arbitrarily complex internal degrees of freedom. Leveraging these semi-analytic gradients, we employ the conjugate gradient (CG) method to optimize thermodynamic model parameters for four binary alloy systems: Cu-Mg, Fe-Ni, Cr-Ni, and Cr-Fe. Across all systems, CG achieves comparable or superior optimality relative to Bayesian ensemble Markov Chain Monte Carlo (MCMC) with improvements in computational efficiency ranging from one to three orders of magnitude. Furthermore, our results establish a new paradigm for CALPHAD assessments in which high fidelity data-rich model calibration becomes tractable using deterministic gradient-informed algorithms.

CALPHAD↗

DOE OSTI2024

Simulating Discharge Curves of an All-Aqueous TRAB to Identify Pathways for Improving System Performance

Thermally regenerative ammonia batteries (TRABs) are an emerging technology that use low temperature heat (T < 150 °C) to recharge a flow battery that produces electrical power on demand. The all-aqueous copper TRAB can provide high power densities and thermal energy efficiencies relative to other devices that harvest energy from waste heat, but its performance is adversely impacted by the crossover of undesired species through the membrane and lower cell voltages compared to conventional batteries. In this work, we developed a numerical model to simulate discharge curves while accounting for crossover inefficiencies without tracking all electrolyte species through the membrane. The model was able to successfully reproduce discharge curves across a diverse range of battery conditions using a single fitting parameter to account for decay of electrode standard potential due to species crossover with minimal error (< 5%). The model was then used to simulate different design scenarios to estimate changes in energy output from alterations to the aspects of the battery electrolyte chemistry. Results from this study are used to identify pathways for improving future TRAB designs with respect to energy capacity and cost-effectiveness of the technology.

25 ENERGY STORAGE↗

DOE OSTI2025

Improved loss functions for machine-learned atomic potentials

Machine learning (ML) has become an invaluable tool across a wide array of domains in science as researchers find new ways to leverage its predictive power. This is especially true in chemistry, where ML is used to fit chemical properties or desirable attributes to the local structure of molecules and materials. In the pursuit of greater accuracy, it is relatively simple to increase the size or complexity of such models, although this often requires simultaneously seeking larger datasets in order to both fit and interpret the larger number of parameters. However, it is equally important to assess the quality and relative importance of the data and how these factors impact the training process. We, therefore, investigate the impact of using different loss functions for training neural network potentials (NNPs), as the loss function defines the error and parameter gradients used to train the NNP. In particular, we test the mean-squared error and Huber loss functions and, using insight from these functions, derive a new loss function based on the Asinh function, which yields significant improvement in the accuracy and generality of NNPs. We show that by discounting/minimizing errors and anomalies in the optimization process, both the Huber and Asinh loss functions improve the training of NNPs, leading to a final potential with a greater effective dimensionality.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

A Principled Framework to Assess the Information-Theoretic Fitness of Brain Functional Sub-Circuits

In systems and network neuroscience, many common practices in brain connectomic analysis are often not properly scrutinized. One such practice is mapping a predetermined set of sub-circuits, like functional networks (FNs), onto subjects’ functional connectomes (FCs) without adequately assessing the information-theoretic appropriateness of the partition. Another practice that goes unchallenged is thresholding weighted FCs to remove spurious connections without justifying the chosen threshold. This paper leverages recent theoretical advances in Stochastic Block Models (SBMs) to formally define and quantify the information-theoretic fitness (e.g., prominence) of a predetermined set of FNs when mapped to individual FCs under different fMRI task conditions. Our framework allows for evaluating any combination of FC granularity, FN partition, and thresholding strategy, thereby optimizing these choices to preserve the important topological features of the human brain connectomes. By applying to the Human Connectome Project with Schaefer parcellations at multiple levels of granularity, the framework showed that the common thresholding value of 0.25 was indeed information-theoretically valid for group-average FCs, despite its previous lack of justification. Our results pave the way for the proper use of FNs and thresholding methods, and provide insights for future research in individualized parcellations.

Duong-Tran, Duy (ORCID:0009000944967575)↗

DOE OSTI2020

Chemical freeze-out parameters of net-kaons in heavy-ion collisions

In this work, we study chemical freeze-out parameters for heavy-ion collisions by performing two different thermal analyses. We analyze results from thermal fits for particle yields, as well as, net-charge fluctuations in order to characterize the chemical freeze-out. The Hadron Resonance Gas (HRG) model is employed for both methods. By separating the light hadrons from the strange hadrons in thermal fits, we study the proposed flavor hierarchy. For the net-charge fluctuations, we calculate the mean-over-variance ratio of the net-kaon fluctuations in the HRG model at the five highest energies of the RHIC Beam Energy Scan (BES) for different particle data lists. We compare these results with recent experimental data from the STAR collaboration in order to extract sets of chemical freeze-out parameters for each list. We focused on particle lists which differ largely in the number of resonant states. By doing so, our analysis determines the effect of the amount of resonances included in the HRG model on the freeze-out conditions. Our findings have potential impact on various other models in the field of relativistic heavy ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Modeling of Thermal Decomposition of TATB-Based Explosive for Safety Analysis

We investigate and model the cook-off behavior of LX-17 to understand the response of explosive systems in abnormal thermal environments. Decomposition has been explored via conventional ODTX (One-Dimensional Time-to-eXplosion), PODTX (ODTX with pressure-measurement), TGA (Thermo-Gravimetric Analysis), and DSC (Differential Scanning Calorimetry) experiments under isothermal and ramped temperature profiles. The data were used to fit reaction rate parameters for proposed schemes in an ALE3D computational model. This model includes chemical reactions, thermo- and hydro-dynamics, and material properties, including thermal expansion, compressibility, and strength. These parameterizations were carried out utilizing a Python evolutionary optimization method on LLNL’s high-performance computing clusters. Additional experiments are being developed to further characterize and monitor decomposition intermediates to improve the model. Once experimentally validated, this model will be scalable to several applications involving LX-17. Furthermore, the optimization methodology developed herein should be applicable to other high explosive materials.

Chemistry - Chemical explosives↗

DOE OSTI2025

Predicting the viscoplastic response of a crystallizing fluoropolymer using transient network theory

We employ a molecular theory of dynamic polymer networks to describe the viscoplastic response of rubbery FK-800, a thermoplastic copolymer of chlorotrifluoroethylene and vinylidene fluoride, over a broad range of thermal histories. The kinetics of crystallization at different annealing temperatures was modeled using a modified Avrami equation, whose parameters were found to evolve through simple relationships over the full temperature range of the rubbery state. By fitting experimental compression data, we discovered predictable trends for the physical parameters in our mechanical model over its full range of crystallinities (up to ≈20%) and provided insights based on molecular-level physics to justify them. Using this, an end-to-end model was developed to predict the yielding and post-yield behavior of rubbery FK-800 for arbitrary thermal histories. The model successfully predicted the highly nonlinear evolution of characteristic mechanical signatures (stiffness, yield point, post-yield drop) throughout the crystallization process. A statistical analysis of variance test was employed to determine that the measured variations in the mechanical behavior of rubbery FK-800 are primarily dictated by its fractional crystallinity, regardless of its exact thermal history.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Dark Energy Survey Year 3: Blue Shear

Modeling the intrinsic alignment (IA) of galaxies poses a challenge to weak lensing analyses. The Dark Energy Survey is expected to be less impacted by IA when limited to blue, star-forming galaxies. The cosmological parameter constraints from this blue cosmic shear sample are stable to IA model choice, unlike passive galaxies in the full DES Y3 sample, the goodness-of-fit is improved and the $\Omega_{m}$ and $S_8$ better agree with the cosmic microwave background. Mitigating IA with sample selection, instead of flexible model choices, can reduce uncertainty in $S_8$ by a factor of 1.5.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Fluid-like shock compression of dilute polymer nanocomposites

The shock Hugoniot of heterogenous mixtures of discrete particles has been experimentally investigated for porous agglomerates and fluids, while the study of full-density solids has been primarily limited to compressed powders and high-dimensional composites. By dispersing ceria nanoparticles in a polyethylene matrix, we are able to examine the hydrodynamic behavior of a nonporous, heterogenous solid in thermal equilibrium during weak shock compression. Phase-driven discontinuities in the Hugoniot particle velocity–shock velocity (u−D) relationship of pure polyethylene are replicated in the nanocomposites but are shifted to lower velocity and to higher pressure with higher particle concentration. The results are explained using an isothermal, two-velocity fluid model under the hydrodynamic approximation. The model, which assumes a theoretical equation-of-state for ceria and either a low-order or high-order fit to the measured polyethylene Hugoniot, reasonably predicts the Hugoniot for two different polyethylene/ceria nanocomposites. Using the model, the mixture Hugoniot is shown to be insensitive to the Hugoniot of the stiffer constituent when the moduli are sufficiently disparate, while dependence on particle density and volume fraction is preserved through fluid-like motion.

Moore, Nathan W. [Sandia National Laboratories (SN↗

DOE OSTI2024

Constraining models for the origin of ultra-high-energy cosmic rays with a novel combined analysis of arrival directions, spectrum, and composition data measured at the Pierre Auger Observatory

The combined fit of the measured energy spectrum and shower maximum depth distributions of ultra-high-energy cosmic rays is known to constrain the parameters of astrophysical models with homogeneous source distributions. Studies of the distribution of the cosmic-ray arrival directions show a better agreement with models in which a fraction of the flux is non-isotropic and associated with the nearby radio galaxy Centaurus A or with catalogs such as that of starburst galaxies. Here, we present a novel combination of both analyses by a simultaneous fit of arrival directions, energy spectrum, and composition data measured at the Pierre Auger Observatory. The model takes into account a rigidity-dependent magnetic field blurring and an energy-dependent evolution of the catalog contribution shaped by interactions during propagation. We find that a model containing a flux contribution from the starburst galaxy catalog of around 20% at 40 EeV with a magnetic field blurring of around 20° for a rigidity of 10 EV provides a fair simultaneous description of all three observables. The starburst galaxy model is favored with a significance of 4.5σ (considering experimental systematic effects) compared to a reference model with only homogeneously distributed background sources. By investigating a scenario with Centaurus A as a single source in combination with the homogeneous background, we confirm that this region of the sky provides the dominant contribution to the observed anisotropy signal. Models containing a catalog of jetted active galactic nuclei whose flux scales with the γ-ray emission are, however, disfavored as they cannot adequately describe the measured arrival directions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Rapid Coal-Ash Characterization using Geophysical Methods & Machine Learning

Coal combustion products (CCP) are challenging to delineate in heterogeneous field settings. Conventional methods (test pits, coring, and laboratory analyses) are labor-intensive, slow, invasive, and provide sparse spatial coverage. This study evaluates whether rapid non-invasive geophysical screening methods—induced polarization (IP), magnetic susceptibility, and nuclear magnetic resonance (NMR) —combined with surface colorimetry (RGB_24), can discriminate CCP-soil mixtures and provide reliable estimates of CCP content. Laboratory measurements were collected on five CCP-soil mixtures (series) and modeled using (i) a linear baseline, (ii) a calibrated non-linear (power-mean) model, and (iii) a machine-learning (ML) Random Forest approach, with validation via leave-one-series-out and site-specific tests. Across the five series, individual signals—particularly IP and magnetic susceptibility—were strongly predictive of ash content but were consistently outperformed by combined models. The pooled calibrated non-linear and ML models captured the observed non-linearity and achieved high accuracy and precision, improving on linear fits. Colorimetry showed the weakest direct relationship with ash content for the tested samples but improved performance when included in multi-signal models. At pre-selected 3.5% decision threshold, calibrated and ML approaches yielded near-perfect classification (Matthews correlation coefficient ˜ 1), suggesting strong practical operability for field screening. Additionally, field-analog tests highlighted the role of endmembers—accuracy declined without access to end-member measurements but was largely recovered by collecting a minimal labeled pair for local recalibration. With end members, accuracy remained high. Globally trained models performed well on three operational unknowns; however, series-specific refits provided the most accurate predictions. Overall, these results highlight the potential of combining rapid geophysics and minimal local calibration for improved coal-ash delineation.

Peshtani, Klaudio↗

DOE OSTI2023

Model validation and selection in metabolic flux analysis and flux balance analysis

13C-Metabolic Flux Analysis (13C-MFA) and Flux Balance Analysis (FBA) are widely used to investigate the operation of biochemical networks in both biological and biotechnological research. Both methods use metabolic reaction network models of metabolism operating at steady state so that reaction rates (fluxes) and the levels of metabolic intermediates are constrained to be invariant. They provide estimated (MFA) or predicted (FBA) values of the fluxes through the network in vivo, which cannot be measured directly. These fluxes can shed light on basic biology and have been successfully used to inform metabolic engineering strategies. Several approaches have been taken to test the reliability of estimates and predictions from constraint-based methods and to compare alternative model architectures. Despite advances in other areas of the statistical evaluation of metabolic models, such as the quantification of flux estimate uncertainty, validation and model selection methods have been underappreciated and underexplored. We review the history and state-of-the-art in constraint-based metabolic model validation and model selection. Applications and limitations of the χ 2 -test of goodness-of-fit, the most widely used quantitative validation and selection approach in 13C-MFA, are discussed, and complementary and alternative forms of validation and selection are proposed. A combined model validation and selection framework for 13C-MFA incorporating metabolite pool size information that leverages new developments in the field is presented and advocated for. Finally, we discuss how adopting robust validation and selection procedures can enhance confidence in constraint-based modeling as a whole and ultimately facilitate more widespread use of FBA in biotechnology.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

A Novel Measurement of the Anomalous Muon Spin Precession Frequency in the Muon $g-2$ Experiment at Fermilab

The Muon $g-2$ Experiment operated at Fermi National Accelerator Laboratory (FNAL, or Fermilab) between 2018 and 2023 to produce the world's most precise measurement of the muon's \textit{anomalous magnetic moment}, $a_\mu = \frac{g_\mu - 2}{2}$, which expresses the relative deviation in the muon's $g$-factor from a baseline theoretical expectation that $g_\mu = 2$. In the Standard Model of particle physics, $g_\mu > 2$ and hence $a_\mu > 0$ by a calculable amount that depends on all possible interactions between the muon and all other fundamental particles, including any potentially undiscovered interactions beyond the Standard Model. For this reason, measurements of the electron anomaly $a_e$ and later the muon anomaly $a_\mu$ have helped guide the development of the Standard Model since the inception of quantum field theory, and the measured value of $a_\mu$ provides a valuable constraint for new hypotheses that extend the Standard Model. As of 2006, the leading measurement and Standard Model prediction for $a_\mu$ exhibited tension at the level of about three standard deviations, motivating an improved measurement at Fermilab that could test the tension more precisely. The experiment functions by storing a polarized beam of $\mu^+$ in a uniform magnetic field, which simultaneously induces circular motion and spin precession. As the stored muons undergo the Michel decay $\mu^+ \to e^+ + \nu_e + \bar{\nu}_\mu$, mediated by the parity-violating weak interaction, the rest-frame $e^+$ emission direction is correlated with the parent $\mu^+$ spin orientation. Boosting into the laboratory frame encodes this correlation in the decay $e^+$ energy, which is higher when the emission (i.e. $\mu^+$ spin direction) is aligned with the $\mu^+$ momentum, and lower when opposite. Detectors then count the rate of high-energy decay $e^+$, which modulates at the difference between the $\mu^+$ revolution and spin precession frequencies. This observed frequency, called the \textit{anomalous spin precession frequency} $\omega_a$, is directly proportional to $a_\mu$. The extraction of $\omega_a$ proceeds by fitting the time spectrum of detected $e^+$, which requires precise modeling of the $\omega_a$ oscillation as well as any perturbations from beam dynamics and detector acceptance. Using the $\omega_a$ analysis presented in this work, based on Runs 4 -- 6 of the Muon $g-2$ Experiment at Fermilab, we find that $a_\mu = 0.001\,165\,920\,738(162)$ with a relative uncertainty of 139 parts per billion.

Barrett, Tyler [Cornell U.]↗

DOE OSTI2021

A best-practice guide to predicting plant traits from leaf-level hyperspectral data using partial least squares regression

Partial least squares regression (PLSR) modelling is a statistical technique for correlating datasets, and involves the fitting of a linear regression between two matrices. One application of PLSR enables leaf traits to be estimated from hyperspectral optical reflectance data, facilitating rapid, high-throughput, non-destructive plant phenotyping. This technique is of interest and importance in a wide range of contexts including crop breeding and ecosystem monitoring. The lack of a consensus in the literature on how to perform PLSR means that interpreting model results can be challenging, applying existing models to novel datasets can be impossible, and unknown or undisclosed assumptions can lead to incorrect or spurious predictions. We address this lack of consensus by proposing best practices for using PLSR to predict plant traits from leaf-level hyperspectral data, including a discussion of when PLSR is applicable, and recommendations for data collection. Further, we provide a tutorial to demonstrate how to develop a PLSR model, in the form of an R script accompanying this manuscript. This practical guide will assist all those interpreting and using PLSR models to predict leaf traits from spectral data, and advocates for a unified approach to using PLSR for predicting traits from spectra in the plant sciences.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

The shock physics of giant impacts: Key requirements for the equations of state

We discuss major challenges in modeling giant impacts between planetary bodies, focusing on the equations of state (EOS). During the giant impact stage of planet formation, rocky planets are melted and partially vaporized. However, most EOS models fail to reproduce experimental constraints on the thermodynamic properties of the major minerals over the required phase space. Here, we present an updated version of the widely-used ANEOS model that includes a user-defined heat capacity limit in the thermal free energy term. Our revised model for forsterite (Mg 2 SiO 4 ), a common proxy for the mantles of rocky planets, provides a better fit to material data over most of the phase space of giant impacts. We discuss the limitations of this model and the Tillotson equation of state, a commonly used alternative model.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Predicting the viability of beta-lactamase: How folding and binding free energies correlate with beta-lactamase fitness

One of the long-standing holy grails of molecular evolution has been the ability to predict an organism’s fitness directly from its genotype. With such predictive abilities in hand, researchers would be able to more accurately forecast how organisms will evolve and how proteins with novel functions could be engineered, leading to revolutionary advances in medicine and biotechnology. In this work, we assemble the largest reported set of experimental TEM-1 β-lactamase folding free energies and use this data in conjunction with previously acquired fitness data and computational free energy predictions to determine how much of the fitness of β-lactamase can be directly predicted by thermodynamic folding and binding free energies. We focus upon β-lactamase because of its long history as a model enzyme and its central role in antibiotic resistance. Based upon a set of 21 β-lactamase single and double mutants expressly designed to influence protein folding, we first demonstrate that modeling software designed to compute folding free energies such as FoldX and PyRosetta can meaningfully, although not perfectly, predict the experimental folding free energies of single mutants. Interestingly, while these techniques also yield sensible double mutant free energies, we show that they do so for the wrong physical reasons. We then go on to assess how well both experimental and computational folding free energies explain single mutant fitness. We find that folding free energies account for, at most, 24% of the variance in β-lactamase fitness values according to linear models and, somewhat surprisingly, complementing folding free energies with computationally-predicted binding free energies of residues near the active site only increases the folding-only figure by a few percent. This strongly suggests that the majority of β-lactamase’s fitness is controlled by factors other than free energies. Overall, our results shed a bright light on to what extent the community is justified in using thermodynamic measures to infer protein fitness as well as how applicable modern computational techniques for predicting free energies will be to the large data sets of multiply-mutated proteins forthcoming.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Fusion Model for Metagenomics

This work highlights the use of an embeddings approach that can encode multiple features and create efficient contextualization of profiled metagenomes derived from microbiome samples using computer vision models and image representations of the abundance profiles. The model's embeddings can be used to cluster existing samples based on multiple conditions and interpretations, and new embeddings can be quickly created for new samples and fitted to existing clusters to characterize them. This has practical applications for unknown, unlabeled microbiome samples. The model's embeddings can be used to cluster existing samples based on multiple conditions and interpretations, and new embeddings can be quickly created for new samples and fitted to existing clusters to characterize them. This has practical applications for unknown, unlabeled microbiome samples.

Valdes, CamiloA [Lawrence Livermore National Labor↗