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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Adjusting a Full Body Model to Mitigate Inverse Kinematics Artifacts in Opensim

BACKGROUND: In support of Vibration Isolation and Stabilization (VIS) system development for Human Health Countermeasures (HHC) exercise systems in space, such as the European Enhanced Exploration Exercise Device (E4D) [1], dynamic quantities required to model the response of a proposed VIS while considering the effect of VIS motion on the forces between the human and VIS platform were obtained using motion capture data [2]. On occasion, large-amplitude oscillatory spikes were found in the subject’s linear and angular momentum derivatives, affecting analyses that depend on forces and moments derived from motion capture. The purpose of this investigation was to identify causes of these artifacts and techniques for their resolution. METHODS AND RESULTS: To obtain the required human dynamic quantities to drive the VIS simulation, motion capture data was collected containing recorded trajectories of passive retroreflective markers on body landmarks of an exercising subject. Since the full body Rajagopal model [3] was originally used to enhance gait analysis, upper body joints did not require large Ranges of Motion (ROM). We thus modified the Rajagopal model [4, 5] to allow it to be used for upper body intensive exercises like those common to the E4D. OpenSim Inverse Kinematics (IK) [6] was performed using these scaled subject models to generate the joint angles throughout the exercise while minimizing marker error. At times, the arms were observed to ‘snap’ from one configuration to another, causing spike artifacts. Following IK, a custom OpenSim plugin [7] was used to determine the required dynamic quantities including the rates of change of the linear and angular momenta of the human. Since motion capture is recorded at a larger time step than required by the VIS simulation, the human center of mass location was fit with splines and a second derivative taken to obtain the momentum derivative, allowing the VIS simulation to maintain conservation of momentum when appropriate. In a few cases during this stage, artifacts much larger than the expected noise of the second derivatives were introduced. Investigation of cases containing artifacts revealed several modifications that could be made to the OpenSim model to improve IK results. Since OpenSim models use Euler angles and rotation sequences, ‘gimbal lock’ would be encountered in the arms when raised 90 degrees to the side (e.g., T-pose, some hang clean exercise, etc.). This was resolved by reorienting the horizontal axes at the shoulder joint by 45 degrees, placing ‘gimbal lock’ outside common arm ROM, with the arm ROMs adjusted following this change. Elbow and wrist ROMs could also be adjusted to allow realistic motion while at the same time limiting the likelihood of unrealistic orientations. On occasion, the arms flipped backwards when raised above the head. This was prevented by using medial elbow markers in scaling and IK. When medial markers were not available, the acromial joint location in the unscaled model was shifted before model scaling to better align the arm with available markers. Lastly, artifacts which became apparent after taking the second derivatives of spline-fit data were found to occur in cases when the pelvis rotation limit prevented the full range of motion of an exercise. These issues were resolved by unclamping the pelvis rotation limit. Through this investigation, an understanding of conditions leading to IK artifacts was acquired allowing the automation of artifact detection. These artifact detection and mitigation techniques can be applied toward modeling of upper body motions in aerospace and other fields for improved IK results.

C A Bell↗
DOE OSTI2025

A mathematical approach to using the forgetting curve to evaluate experience and training factors in human reliability analysis

Traditional human reliability analysis (HRA) methods have difficulty dealing with the dynamic nature of factors such as time and rely on static and expert-judgment-based assessments of performance-shaping factors (PSFs) across limited levels. In this study, we introduce a mathematical approach for dynamically evaluating the experience and training PSF. Our proposed method integrates the psychological concept of the “forgetting curve” to evaluate how PSFs are impacted by the number of trainings and the time elapsed since training. To confirm the validity of the model, we provide experimental data fitted by identifying the quantitative relationship between training and human performance. This research enables dynamic and objective assessments, thus reducing reliance on subjective expert judgment and improving the accuracy of HRA.

99 - GENERAL AND MISCELLANEOUS↗
DOE OSTI2026

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↗
DOE OSTI2024

CO 2 storage site characterization using ensemble-based approaches with deep generative models

Estimating spatially distributed properties such as permeability from available sparse measurements is a great challenge in efficient subsurface CO 2 storage operations. In this paper, a deep generative model that can accurately capture complex subsurface structure is tested with an ensemble-based inversion method for accurate and accelerated characterization of CO 2 storage sites. We chose Wasserstein Generative Adversarial Network with Gradient Penalty (WGAN-GP) for its realistic reservoir property representation and Ensemble Smoother with Multiple Data Assimilation (ES-MDA) for its robust data fitting and uncertainty quantification capability. WGAN-GP are trained to generate high-dimensional permeability fields from a low-dimensional latent space and ES-MDA then updates the latent variables by assimilating available measurements. Several subsurface site characterization examples including Gaussian, channelized, and fractured reservoirs are used to evaluate the accuracy and computational efficiency of the proposed method and the main features of the unknown permeability fields are characterized accurately with reliable uncertainty quantification. Furthermore, the estimation performance is compared with a widely-used variational, i.e., optimization-based, inversion approach, and the proposed approach outperforms the variational inversion method in several benchmark cases. We explain such superior performance by visualizing the objective function in the latent space: because of nonlinear and aggressive dimension reduction via generative modeling, the objective function surface becomes extremely complex while the ensemble approximation can smooth out the multi-modal surface during the minimization. This suggests that the ensemble-based approach works well over the variational approach when combined with deep generative models at the cost of forward model runs unless convergence-ensuring modifications are implemented in the variational inversion.

42 ENGINEERING↗
DOE OSTI2026

The evolution of coal porosity during pyrolysis

Gasification of coal, municipal waste, or other organic materials is a potential hydrogen source that entails complex thermal decomposition and transport processes. This study provides a multiscale analysis of these processes for sub-bituminous (Usibelli, Healy, Alaska) and lignite (Center, North Dakota) coals and provides data useful for process design. The chemistry, mineralogy, and pore structures of pyrolyzed coal and their evolution with thermal decomposition are discussed. Samples pyrolyzed at 200–1000 °C were analyzed by small-angle neutron scattering; ultra-small, small-, and wide-angle X-ray scattering; and other complementary techniques. Scanning electron microscopy showed new pores in the high-temperature-pyrolyzed material. Upon heating, the coals became progressively denser, and the concentration of hydrogen decreased. Changes in pore volume fell into three temperature ranges: an initial, low-temperature range that, for the Usibelli coal, involved an increase in overall porosity; a mid-temperature range associated with pore volume loss; and a high-temperature range associated with significant porosity increase and char formation. This transformation was paralleled by changes in fractal dimension and correlation length. The higher the pyrolysis temperature the greater the small-pore-volume fraction and overall surface area became. Pyrolysis increased the lateral size of coal crystallites, decreased the amorphous fraction, and increased the aromatics fraction and overall coal rank. Comparisons of neutron and X-ray scattering data and subsequent water uptake studies showed that pre-dried coals can re-hydrate relatively rapidly upon exposure to air, which can significantly affect the porosity calculated from small-angle-scattering data. Fits to the cumulative porosity curves provide a method for modeling the physical and chemical transformation of hydrogen-containing feedstock during gasification.

Anovitz, Lawrence {Larry} [ORNL] (ORCID:0000000226↗
DOE OSTI2026

Adsorption of Mixed Micelles of Polysorbate 80 and Oleic Acid to the Air–Water Interface

The solubilization of long chain amphiphiles with limited solubility into mixed micelles composed of highly soluble surfactants plays a crucial role in modulating the stability and functionality of formulations in pharmaceutical and food systems. Subsequently, the mixed micelles adsorb from solution to the air−water surface, defining the transport mechanism of the insoluble amphiphiles to the interface. We use X-ray reflectivity to measure the composition and provide insight into the structure of these mixed monolayers, and, alongside interfacial tension measurements, we provide an understanding of how the bulk composition determines the surface composition and the dynamics of tension reduction. We use a model system consisting of an insoluble fatty acid, oleic acid (OA), and a soluble micelle-forming surfactant, polysorbate 80 (PS80). PS80 forms spherical micelles, and, above the critical micelle concentration (cmc), OA is readily solubilized inside the micelles. We show that the adsorption of PS80/OA mixed micelles to the air−water interface rapidly reduces the tension and lowers the equilibrium tension in proportion to the OA concentration. X-ray reflectivity data, fit using Parratt-slab models, quantitively demonstrates that the monolayers become enriched with OA, and we show how the OA intercalates into the PS80 monolayers.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Coupling to rotational manifolds to improve gas-phase pump–probe spectroscopic models

The physical picture of gas-phase optical transitions is normally presented as an isolated two-level system balanced by upward and downward processes. Isolated models assume a phenomenological treatment of collisional dephasing but do not strictly account for collisional population exchange with the rotational baths. While this assumption is valid under low-intensity conditions, where excitation is rate-limiting, isolated models can deviate from Beer’s Law at sufficient pressures and monochromatic intensities when both collisional broadening and power broadening are comparable to (or greater than) lifetime broadening, which are not uncommon conditions for cavity enhanced spectroscopies in the mid-IR spectral range. Although this problem has been addressed by rate-equation models for linear absorption measurements, a general treatment for multi-level quantum mechanical models suitable for non-linear absorption measurements (two-photon/two-color/pump–probe) is lacking. Isolated models require physical parameter inputs that disagree with expected values by at least an order of magnitude. These non-physical models undermine the ability to predict non-linear signal strengths under untested conditions and thereby limit the potential to optimize the sensitivity of non-linear spectroscopies and to expand their analytical applications (e.g., new analytes and/or buffer gases, changes in cavity free-spectral-range, changes in intracavity powers or wavelengths, and accurate investigation of physical phenomena). In this study, we derive bath-coupled models for gaseous pump–probe spectroscopy by application of the quantum Lindblad equation and detailed balance. Bath-coupled models are shown to fit data consistently across variations in intensity and agree with all physically expected values.

Cavity ring-down spectroscopy↗
DOE OSTI2024

Magnetic tunnel junction random number generators applied to dynamically tuned probability trees driven by spin orbit torque

Abstract Perpendicular magnetic tunnel junction (pMTJ)-based true-random number generators (RNGs) can consume orders of magnitude less energy per bit than CMOS pseudo-RNGs. Here, we numerically investigate with a macrospin Landau–Lifshitz-Gilbert equation solver the use of pMTJs driven by spin–orbit torque to directly sample numbers from arbitrary probability distributions with the help of a tunable probability tree. The tree operates by dynamically biasing sequences of pMTJ relaxation events, called ‘coinflips’, via an additional applied spin-transfer-torque current. Specifically, using a single, ideal pMTJ device we successfully draw integer samples on the interval [0, 255] from an exponential distribution based on p -value distribution analysis. In order to investigate device-to-device variations, the thermal stability of the pMTJs are varied based on manufactured device data. It is found that while repeatedly using a varied device inhibits ability to recover the probability distribution, the device variations average out when considering the entire set of devices as a ‘bucket’ to agnostically draw random numbers from. Further, it is noted that the device variations most significantly impact the highest level of the probability tree, with diminishing errors at lower levels. The devices are then used to draw both uniformly and exponentially distributed numbers for the Monte Carlo computation of a problem from particle transport, showing excellent data fit with the analytical solution. Finally, the devices are benchmarked against CMOS and memristor RNGs, showing faster bit generation and significantly lower energy use.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Error field measurements with rotating RMP fields for DIII-D H-mode

3D magnetic sensors are employed to identify the amplitude and toroidal phase of error fields (EF) by analyzing the torque balance for magnetic islands entrained by rotating resonant magnetic perturbations (RMPs) in DIII-D H-mode plasmas. This technique of torque balance allows for efficient error field identification, offering a valuable tool for scenario-specific and optimized error field compensation (EFC) and requires only magnetic diagnostics. The torque balance used in this work includes the contribution from electromagnetic torque due to error fields, wall response, RMP fields, and viscous. Results show that viscous torque plays a crucial role, particularly during locked modes and H-mode plasmas, ensuring accurate data fits with lower residuals. The torque balance technique reveals that the L- and H-mode plasmas have distinct EF configurations, and consistent with a model-based EF analysis including MHD response in IPEC and the SURFMN EF simulation. This technique shows great robustness in measuring the intrinsic EF amplitude regardless of its amplitude or toroidal phase. Repeated discharges with EFC disparities exhibit consistent results of intrinsic error field within a reasonable range near the “standard” error field compensation. Additionally, the use of a rotating n = 1 resonant magnetic perturbation offers the advantage of reducing disruption risks by entraining saturated magnetic islands. These findings are instrumental for optimizing EF correction in fusion devices, thereby enhancing tearing mode suppression and overall plasma stability.

3D magnetic sensors↗
DOE OSTI2025

Numerical calculation of Coulomb corrections in forward elastic 𝑝↑⁢𝑝 and 𝑝↑⁢𝐴 scattering

The analysis of RHIC hydrogen gas jet target polarimeter measurements of transverse analyzing powers 𝐴 N ⁡(𝑡) in proton-nucleus scattering requires accurate Coulomb corrections to both spin-flip and nonflip amplitudes. These corrections must cover a wide range of nuclear charges 𝑍 and form factor slopes, with flexibility to vary form factors during data fitting. To avoid technically challenging calculations involving a small but finite fictitious photon mass, the Coulomb correction to the nonflip electromagnetic amplitude with an exponential form factor was related to the corresponding correction for the spin-flip amplitude. This approach allows soft photon contributions to all amplitudes, including those with nonexponential form factors, to be calculated in the massless photon limit using only analytical expressions and numerically stable integrals with nonsingular integrands and finite integration limits. In addition, an absorptive correction to the spin-flip electromagnetic amplitude, which plays a critical role in spin effects in forward polarized proton-nucleus scattering, was accurately evaluated.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Semi-Analytical Hierarchical Bayesian Inference of Nonlinear Model Structure in Stochastic Dynamics: Applied to Compartmental Models of Infectious Diseases

A Bayesian computational framework for parsimonious inference in stochastic nonlinear dynamical systems is presented. This framework enables the concurrent estimation of system states, time-varying parameters, time-invariant parameters, and the optimal sparsity structure of the model parameters. Because differential equation-based models are often simplified mechanistic or phenomenological representations, robust inference from noisy measurement data requires explicit treatment of model error and uncertainty. Model error and time-varying parameters can be represented as random processes, enabling inference while making minimal assumptions about the underlying sources of discrepancy and variability. Adopting stochastic differential equation representations affords the model significant flexibility, but can also render it susceptible to overfitting during statistical inversion, where the inferred model may track noise rather than the underlying signal. To alleviate the effects of overfitting and to enable the discovery of the optimal sparse representation of the time-invariant parameters, a Bayesian sparse learning algorithm is embedded within the framework. This sparse learning framework adopts an approximate hierarchical Bayesian setting defined by a series of semi-analytical expressions. The model structure inference framework is validated using a stochastic compartmental model for tracking and forecasting active cases of an infectious disease. Compartmental models describe population-level infectious disease dynamics through interactions among population fractions grouped by disease state. Mathematically, such models consist of a system of coupled ordinary differential equations. This example adopts an expressive compartmental model that includes multiple possible interactions between disease states, motivated by early uncertainty surrounding COVID-19 reinfection dynamics and their implications for long-term epidemic forecasting. The sparse learning exercise permits the inference of a priori unknown epidemiological dynamics from simulated public health data, discovering the nested compartmental model that optimizes the trade-off between average data-fit and model complexity. It is shown that inducing sparsity among the model parameters eliminates redundant interactions between compartments, equivalently revealing the optimal coupling structure between differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Enabling in-situ LIBS measurements of liquids and slurries

The aim of this work was to explore the use of, and further develop, a in-situ and near-real time LIBS (Laser Induced Breakdown Spectroscopy) instrument to measure liquids relevant to the Hanford site. Specifically, we investigated the feasibility of using LIBS to measure the cation composition of HLW (High Level Waste) streams. In addition, we focused on enabling continuous interrogation of those liquids using the LIBS system. For the first part, we investigated simple (sodium nitrate) solutions as well as various Hanford tank waste surrogate slurries. We showed that relevant cations (such as sodium, aluminum, magnesium, potassium) can be measured via LIBS under good conditions but the results – both quantitative and qualitative (aka reproducibility) – strongly depend on the kind and concentration of the material used. The dependence of signal intensity on concentration was also confirmed by other measurements conducted (sodium chloride and sodium nitrate in water). For the second part, we investigated the influence of varying depth and concentration of a liquid sample. We observed that for the concentrations and fill depths investigated, the liquid fill depth did not impact the emission intensity of the analyte, while concentration is confirmed to have significant effects. The results of the sample depth experiments influenced testing of waste surrogates. In addition, we focused on further developing and adapting data fitting tools in order to efficiently, accurately, and reproducibly identify the elements and materials present in a LIBS spectrum.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Unrolled Video Super-Resolution Network with Autoregressive Prior for the Case of Known Motion

Real-time detection and classification of distant objects is necessary for many national security applications. However, when objects are far from the sensor, they occupy only a small number of pixels in the captured video, limiting the amount of visual detail available for recognition. State-of-the-art classification methods typically rely on high-resolution (HR) video streams to capture characteristic object features, but obtaining such detail is challenging for distant objects that occupy only a few pixels. This motivates the development of video super-resolution (VSR) methods that enhance object classification by recovering fine details from low-pixel representations. Current VSR methods rely either on model-based optimization, which is interpretable but computationally expensive, or on learning-based approaches, which are efficient and high-performing but often lack flexibility and interpretability. In this report, we propose an end-to-end trainable unrolled VSR network, UVSRNet, which super-resolves each frame in a video by exploiting sub-pixel motion between neighboring low-resolution (LR) frames as well as incorporating high-frequency detail from previously super-resolved frames. In particular, by unrolling a plug-and-play (PnP) half-quadratic splitting (HQS) algorithm, we leverage a model-based data-fitting module alongside a learning-based autoregressive prior module. This combination yields a method that maintains the flexibility and interpretability of model-based methods while achieving the performance advantages of learning-based methods.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Quadrature Based Neural Network Learning of Stochastic Hamiltonian Systems

Hamiltonian Neural Networks (HNNs) provide structure-preserving learning of Hamiltonian systems. In this paper, we extend HNNs to structure-preserving inversion of stochastic Hamiltonian systems (SHSs) from observational data. We propose the quadrature-based models according to the integral form of the SHSs’ solutions, where we denoise the loss-by-moment calculations of the solutions. The integral pattern of the models transforms the source of the essential learning error from the discrepancy between the modified Hamiltonian and the true Hamiltonian in the classical HNN models into that between the integrals and their quadrature approximations. This transforms the challenging task of deriving the relation between the modified and the true Hamiltonians from the (stochastic) Hamilton–Jacobi PDEs, into the one that only requires invoking results from the numerical quadrature theory. Meanwhile, denoising via moments calculations gives a simpler data fitting method than, e.g., via probability density fitting, which may imply better generalization ability in certain circumstances. Numerical experiments validate the proposed learning strategy on several concrete Hamiltonian systems. The experimental results show that both the learned Hamiltonian function and the predicted solution of our quadrature-based model are more accurate than that of the corrected symplectic HNN method on a harmonic oscillator, and the three-point Gaussian quadrature-based model produces higher accuracy in long-time prediction than the Kramers–Moyal method and the numerics-informed likelihood method on the stochastic Kubo oscillator as well as other two stochastic systems with non-polynomial Hamiltonian functions. Moreover, the Hamiltonian learning error εH arising from the Gaussian quadrature-based model is lower than that from Simpson’s quadrature-based model. These demonstrate the superiority of our approach in learning accuracy and long-time prediction ability compared to certain existing methods and exhibit its potential to improve learning accuracy via applying precise quadrature formulae.

Mathematics↗
NASA NTRS1973

Role of gas-surface interactions in the reduction of OGO 6 neutral-particle mass spectrometer data

Data obtained with the quadrupole mass spectrometer aboard the OGO 6 satellite show the effects of significant surface interaction processes, including nearly complete recombination of incoming atomic oxygen on the walls of the instrument antechamber plus absorption and desorption of oxygen and carbon monoxide. The observed data are fit by solving the time-dependent continuity equations accounting for production and loss of atomic oxygen, molecular oxygen, and (in the case of mass 28) carbon monoxide. The surface parameters that best fit the data are selected and applied to the determination of ambient densities and their estimated errors.

Hedin, A. E.↗
NASA NTRS1973

Role of gas-surface interactions in the reduction of Ogo 6 neutral particle mass spectrometer data.

Data obtained with the quadrupole mass spectrometer aboard the Ogo 6 satellite show the effects of significant surface interaction processes, including nearly complete recombination of incoming atomic oxygen on the walls of the instrument antechamber plus adsorption and desorption of oxygen and carbon monoxide. The observed data are fit by solving the time-dependent continuity equations accounting for production and loss of atomic oxygen, molecular oxygen, and (in the case of mass 28) carbon monoxide. The surface parameters that best fit the data are selected and applied to the determination of ambient densities and their estimated errors.

Hedin, A. E.↗
NASA NTRS1972

Revision of the standard reference data for thermocouples.

Revision of the International Practical Temperature Scale requires that there be changes for all accurately tabulated thermophysical values. Revised reference data for thermocouples have been generated in a program carried out by the National Bureau of Standards. The new reference data reflect not only revisions in the temperature scale, but also slight changes in the materials themselves and improvements in data fitting methods. A new NBS monograph that contains tables, analytic expressions, various approximations, and explanatory text has been prepared. A general discussion of the project and some specific examples are given.

Powell, R. L.↗