Search NASA⌕ Search

SEARCH · Search NASA

Results for “method of averaging”

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 451 records · Page 25

Predicting Li-Ion Battery Capacity Fade Using Early-Life Data and a Hybrid Data-Driven Gaussian Process-Bayesian Regression Approach

Accurately predicting Li-ion battery capacity trajectories using early-life data can dramatically improve battery-life understandings and be used to rapidly evaluate design/cost/performance trade-offs when developing new battery materials. Accurate early-life predictions enable researchers to quickly iterate over cell designs and material precursor properties without consistently cycling cells to failure. To this end, we present a toolbox that uses a combined Gaussian Process and Bayesian regression approach that capitalizes on signals other than just capacity (e.g., dQ/dV, voltage drops) to rapidly predict capacity-fade trajectories. The prediction tool uses Bayesian regression to fit functional forms, e.g., power law, sigmoids, etc., to predict capacity-fade dynamics. By fitting functional forms, the capacity fade can be interrogated at any point in the future, allowing for early cell-failure prediction. Additionally, Bayesian regression allows for accurate uncertainty estimates that account for cell-to-cell variability (aleatoric uncertainty) and the lack of observation data (epistemic uncertainty). By only using early cycle data to predict the capacity fade trajectory, uncertainty bounds at end-of-life can be extremely large. The large uncertainty bounds are further exacerbated because there is no systematic way to define the prior distribution of the functional forms' parameters. We improve our the predicted trajectory confidence interval of our predicted trajectory using two methods. First, we shows that a small amount of held-out cycling data is sufficientuse some train cells, that have been cycled to failure to derive information regarding the appropriate prior distributions for the functional forms' parameters of the functional form, effectively leading to data-driven priors.. We propose constructing the data-driven priors by first running a Bayesian regression starting with uninformed priors to generate intermediate cell-specific posterior parameter distributions. These posterior distributions are combined using a Ggaussian mixture model for each parameter to create the data-driven priors. These mixture models serve as the data-driven prior distributions for the parameters for. Second, we derive multiple features, e.g., C_dchg 0.5 DoD 0.5, log (|mean(dQ/dV_(w_3-w_0 ) (V)|), etc., from the train cellsheld-out cycling data, identify which the features are that best predicting capacity at early/mid-life cycles, and then create Ggaussian process regression models that are used for predicting capacity at early/mid-life cycles for the test cells (see blue dots with error bars in Fig 1b). Finally, these predicted data-points are used in addition to the actual early cycle data capacity fade to construct the Bayesian regression trajectory for the test cell s. Notably. We note that these two methods are complementary and can be combined with each other. We evaluate the performance of our proposed method on an testing open-source dataset from Iowa State University and Iowa Lakes Community College (ISU-ILCC). This dataset comprises of 251 nickel-manganese-cobalt/graphite Lithium-ion cells that are cycled under 63 different conditions. We compute the mean average percentage error (MAPE) and negative log predictive density (NLPD) to quantify the efficacy of our method. Our initial findings suggest that, when only few observations are available, for test cells, when using only Bayesian regression with uninformed priors, a power law functional provides the most accurate predictions. with very few data points. However, asHowever, a the number of data points increases, a twin sigmoidal function becomes more accurate as the number of observations further increases. We also find that using as little as 10% of the data set towards generating data-driven priors can lead to significant improvement in prediction accuracy when using early cycle data. Lastly, we found that augmenting early-cycle data with Gaussian process-predicted capacity data for Bayesian regression greatly improves the prediction accuracy. We will present a comprehensive comparison of our methods to other methods available in the literature and apply this method to additional battery datasets.

42 ENGINEERING↗

Evaluation of the PhaseNet Model Applied to the IMS Seismic Network

Producing a complete and accurate set of signal detections is essential for automatically building and characterizing seismic events of interest for nuclear explosion monitoring. Signal detection algorithms have been an area of research for decades, but still produce large quantities of false detections and misidentify real signals that must be detected to produce a complete global catalog of events of interest. Deep learning methods have shown promising capabilities in effectively characterizing seismic signals for complex tasks such as identifying phase arrival times. We use the PhaseNet model, a UNet-based Neural Network, trained on local distance data from northern California to predict seismic arrivals on data from the International Monitoring System (IMS) global network. We use an analyst-curated bulletin generated from this data set to compare the performance of PhaseNet to that of the Short-Term Average/Long-Term Average (STA/LTA) algorithm. We find that PhaseNet has the potential of outperforming traditional processing methods and recommend the training of a new model with the IMS data to achieve optimal performance.

58 GEOSCIENCES↗

Reduced-Order Modeling Method for Phase-Change Thermal Energy Storage Heat Exchangers

Thermal energy storage can facilitate the effective utilization of renewable energy. To speed up the design process of thermal energy storage devices, it is critical to develop fast and accurate modeling methods for phase change material embedded heat exchangers (PCM HXs). This study developed and compared two approximation-assisted reduced-order PCM HX models for the simulation of thermal storage components and systems, which were verified against a validated finite-volume model. They are a pure black-box model, and a grey-box model based on the Number of Transfer Units (or effectiveness-NTU) approach. We used the reduced-order models to predict the performance of standalone PCM HX for 1000 cases with different fluid inlet conditions and HX designs. We also integrated the models with a vapor compression system to predict the compressor energy consumption and total charging time under various conditions. The results show that overall, the black-box model gave more accurate results than the grey-box model. On average, the mean absolute deviation in the PCM HX fluid outlet temperature was 0.05 K and 0.1 K for the black-box model and grey-box model, respectively. The grey-box model gave larger temperature deviation toward the end of the phase change process, due to the simplified two-node PCM representation. On the system level, the mean absolute deviation in compressor energy consumption was 0.2% and 0.3% for the black-box model and grey-box model, respectively. Their mean absolute deviation in total charging time was 1.1% and 2.6%, respectively. In terms of computation efficiency, the system simulation speedup ratio gained by using the reduced-order models was 11 to 57. System simulation time was decreased from an average of 1465 s to an average of 59 s. This shows that the proposed reduced-order modeling methods can be used to predict system performance of PCM HXs with less than 3% accuracy penalty, and 25 times less computational time than finite-volume models, enabling faster design and evaluation of PCM thermal storage devices.

25 ENERGY STORAGE↗

A deep learning model for automatic analysis of cavities in irradiated materials

Transmission electron microscopy (TEM) is a commonly used technique in materials science for defect investigation. Quantitative analysis of defects is important for understanding the properties of a material, but manual analysis of TEM micrographs can be time-consuming and prone to error, especially when the defects have irregular shapes rather than spherical shapes. Many existing methods or deep learning models do not handle a wide range of sizes for the same object type within a single image. In this work, we present a framework that enables users to train an instance segmentation model called Mask R- CNN on any microstructure dataset, perform multi-detection on the same image at different scales, and obtain properties (e.g., size, area) of the objects based on the desired shape (e.g., circle, ellipse, rectangle). Additionally, we have developed a parallel detection module that uses multiple GPUs to increase the efficiency of the object detection process. We demonstrate the capabilities of our framework using a set of TEM images of cavities with different shapes, size distributions, and background contrasts. Finally, we show that the performance of our model in terms of density, size, and swelling of the cavities is comparable to the human average and that our model achieves the highest recall value compared to existing methods due to the use of image multi-rescaling.

36 MATERIALS SCIENCE↗

Maximizing machine learning interatomic potential transferability for the discovery of the novel stellated octadecagon Bi18-Pt24 cage structure

Achieving true transferability remains the central challenge for Machine Learning Interatomic Potentials (ML-IAPs) in modeling complex bimetallic nanoclusters across their vast potential energy surfaces. We systematically investigate data selection strategies to optimize the Chebyshev Interaction Model for Efficient Simulation (ChIMES) potential for the Bi-Pt nanoclusters by comparing three innovative sampling methods: Principal Component Analysis (PCA)/k-means (structural diversity), t-distributedStochasticNeighborEmbedding (t-SNE)/k-means (force-space diversity), and hierarchical clustering. Quantitatively, the PCA/k-means strategy proved most effective for global accuracy, yielding the lowest force errors and achieving energy root mean square errors (RMSE) values competitive with Density Functional Theory (DFT), demonstrating excellent accuracy (19.16meV/atom). Structural validation on 34 unique DFT-optimized isomers further confirmed the potential’s high fidelity, with the best model PCA/k-means reproducing structures with an average root mean square deviation (RMSD) of 0.10 Å. However, the t-SNE methods, by maximizing diversity in the force space, demonstrated superior extrapolative power, leading to the more precise prediction of a novel stellated octadecagon Bi18⁢Pt24 cage structure, demonstrating the potential for exploring previously unseen morphologies. Our results establish a clear methodology for strategic data sampling that successfully maximizes ML-IAP transferability, providing an accurate and computationally efficient tool that accelerates the theoretical discovery of complex bimetallic architectures.

Vangheluwe, Raphaël [Université Paris-Saclay, CNRS↗

A comparative study of machine learning models for predicting the state of reactive mixing

Mixing phenomena are important mechanisms controlling flow, species transport, and reaction processes in fluids and porous media. Accurate predictions of reactive mixing are critical for many Earth and environmental science problems such as contaminant fate and remediation, macroalgae growth, and plankton biomass evolution. Here, to investigate the evolution of mixing dynamics under different scenarios (e.g., anisotropy, fluctuating velocity fields), a finite-element-based numerical model was built to solve the fast, irreversible bimolecular reaction-diffusion equations to simulate a range of reactive-mixing scenarios. A total of 2,315 simulations were performed using different sets of model input parameters comprising various spatial scales of vortex structures in the velocity field, time-scales associated with velocity oscillations, the perturbation parameter for the vortex-based velocity, anisotropic dispersion contrast (i.e., ratio of longitudinal-to-transverse dispersion), and molecular diffusion. The outputs comprised concentration profiles of reactants and products. The inputs to and outputs from these simulations were concatenated into feature and label matrices, respectively, to train 20 different machine learning (ML) models intended to emulate system behavior. These 20 ML emulators, based on linear methods, Bayesian methods, ensemble learning methods, and multilayer perceptrons (MLPs), were trained to classify the state of mixing and predict three quantities of interest (QoIs) characterizing species production, decay (i.e., average concentration, square of average concentration), and degree of mixing (i.e., variances of species concentration). Unsurprisingly, linear classifiers and regressors failed to reproduce the QoIs; however, ensemble methods (classifiers and regressors) and the MLP model accurately classified the state of reactive mixing and the QoIs. Among ensemble methods, random forest and decision-tree-based AdaBoost faithfully predicted the QoIs. At run time, trained ML emulators produced results times faster than the finite-element simulations. Due to their low computational expense and high accuracy, ensemble and MLP models are excellent emulators for these numerical simulations and great utilities in uncertainty quantification exercises, which can require 1,000s of forward model runs.

97 MATHEMATICS AND COMPUTING↗

Novel pore size-controlled, susceptibility matched, 3D-printed MRI phantoms

We report the design concept and fabrication of MRI phantoms, containing blocks of aligned microcapillaires that can be stacked into larger arrays to construct diameter distribution phantoms or fractured, to create a “powder-averaged” emulsion of randomly oriented blocks for vetting or calibrating advanced MRI methods, that is, diffusion tensor imaging, AxCaliber MRI, MAP-MRI, and multiple pulsed field gradient or double diffusion-encoded microstructure imaging methods. The goal was to create a susceptibility-matched microscopically anisotropic but macroscopically isotropic phantom with a ground truth diameter that could be used to vet advanced diffusion methods for diameter determination in fibrous tissues. Two-photon polymerization, a novel three-dimensional printing method is used to fabricate blocks of capillaries. Double diffusion encoding methods were employed and analyzed to estimate the expected MRI diameter. Susceptibility-matched microcapillary blocks or modules that can be assembled into large-scale MRI phantoms have been fabricated and measured using advanced diffusion methods, resulting in microscopic anisotropy and random orientation. This phantom can vet and calibrate various advanced MRI methods and multiple pulsed field gradient or diffusion-encoded microstructure imaging methods. We demonstrated that two double diffusion encoding methods underestimated the ground truth diameter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Absolute quantitation of propranolol from 200-μm regions of mouse brain and liver thin tissues using laser ablation-dropletProbe-mass spectrometry

Rationale: The ability to quantify drugs and metabolites in tissue with sub-mm resolution is a challenging but much needed capability in pharmaceutical research. To fill this void, a novel surface sampling approach combining laser ablation with the commercial dropletProbe automated liquid surface sampling system (LA-dropletProbe) was developed and presented here.Methods: Parylene C-coated 200 × 200 μm tissue regions of mouse brain and kidney thin tissue sections were analyzed for propranolol by laser ablation of tissue directly into a preformed liquid junction. Propranolol was detected by high performance liquid chromatography with positive ion mode electrospray ionization tandem mass spectrometry. Quantitation was achieved via application of a stable isotope-labeled internal standard and an external calibration curve.Results: Absolute concentration of propranolol determined from 200 × 200 µm tissue regions were compared to propranolol concentrations obtained from 2.3-mm-diameter tissue punches of adjacent, non-coated sections using standard bulk tissue extraction protocols followed by regular HPLC-MS/MS analysis. Average concentration of propranolol in both organs determined by the two employed methods agreed within ±12%. Furthermore, the relative abundance of phase II hydroxypropranolol glucuronide metabolites were recorded and found to be consistent with previous results.Conclusions: This work illustrates that depositing a thin parylene C layer onto thin tissue prior to analysis, which seals the surface and prevents direct liquid extraction of the drug from the tissue, coupled to the novel LA-dropletProbe surface sampling system is a viable approach for sub-mm resolution quantitative drug distribution analysis.

59 BASIC BIOLOGICAL SCIENCES↗

Optimization using pathwise algorithmic derivatives of electromagnetic shower simulations

Among the well-known methods to approximate derivatives of expectancies computed by Monte-Carlo simulations, averages of pathwise derivatives are often the easiest one to apply. Computing them via algorithmic differentiation typically does not require major manual analysis and rewriting of the code, even for very complex programs like simulations of particle-detector interactions in high-energy physics. However, the pathwise derivative estimator can be biased if there are discontinuities in the program, which may diminish its value for applications. This work integrates algorithmic differentiation into the electromagnetic shower simulation code HepEmShow based on G4HepEm, allowing us to study how well pathwise derivatives approximate derivatives of energy depositions in a sampling calorimeter with respect to parameters of the beam and geometry. We found that when multiple scattering is disabled in the simulation, means of pathwise derivatives converge quickly to their expected values, and these are close to the actual derivatives of the energy deposition. Additionally, we demonstrate the applicability of this novel gradient estimator for stochastic gradient-based optimization in a model example.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Langmuir turbulence in suspended kelp farms

This study investigates the influence of suspended kelp farms on ocean mixed layer hydrodynamics in the presence of currents and waves. We use the large eddy simulation method, where the wave effect is incorporated by solving the wave-averaged equations. Distinct Langmuir circulation patterns are generated within various suspended farm configurations, including horizontally uniform kelp blocks and spaced kelp rows. Intensified turbulence arises from the farm-generated Langmuir circulation, as opposed to the standard Langmuir turbulence observed without a farm. The creation of Langmuir circulation within the farm is attributed to two primary factors depending on farm configuration: (i) enhanced vertical shear due to kelp frond area density variability, and (ii) enhanced lateral shear due to canopy discontinuity at lateral edges of spaced rows. Both enhanced vertical and lateral shear of streamwise velocity, representing the lateral and vertical vorticity components, respectively, can be tilted into downstream vorticity to create Langmuir circulation. This vorticity tilting is driven by the Craik–Leibovich vortex force associated with the Stokes drift of surface gravity waves. In addition to the farm-generated Langmuir turbulence, canopy shear layer turbulence is created at the farm bottom edge due to drag discontinuity. The intensity of different types of turbulence depends on both kelp frond area density and the geometric configuration of the farm. The farm-generated turbulence has substantial consequences for nutrient supply and kelp growth. These findings also underscore the significance of the presence of obstacle structures in modifying ocean mixed layer characteristics.

Bo, Tong (ORCID:0000000260300561)↗

Mesoporous Crystalline Niobium Oxide with a High Surface Area: A Solid Acid Catalyst for Alkyne Hydration

A mesoporous crystalline niobium oxide with tunable pore sizes was synthesized via the sol–gel-based inverse micelle method. The material shows a surface area of 127 m 2 /g, which is the highest surface area reported so far for crystalline niobium oxide synthesized by soft template methods. The material also has a monomodal pore size distribution with an average pore diameter of 5.6 nm. A comprehensive characterization of niobium oxide was performed using powder X-ray diffraction, Brunauer–Emmett–Teller, thermogravimetric analysis, scanning electron microscopy, transmission electron microscopy, UV–vis, and X-ray photoelectron spectroscopy. The material acts as an environmentally friendly, solid acid catalyst toward hydration of alkynes under with excellent catalytic activity (99% conversion, 99% selectivity, and 4.39 h –1 TOF). Brønsted acid sites present in the catalyst were found to be responsible for the high catalytic activity. Lastly, the catalyst was reusable up to five cycles without a significant loss of the activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Biopolymer-Templated Titania Film Formation for Nanostructured Coatings Revealed by Machine Learning-Supported Time-Resolved Analysis

This study presents a machine learning approach to derive the film formation of biopolymer-templated titania nanostructures during spray deposition, in combination with in situ grazing-incidence small-angle X-ray scattering (GISAXS). A neural network trained on synthetic GISAXS data directly predicts domain-size distributions from experimental two-dimensional scattering patterns, capturing the full kinetics of nanostructure evolution with high temporal resolution. The predictions reveal hierarchical size distributions and periodic growth features, consistent with layer-by-layer spray deposition and validated by complementary scanning electron microscopy (SEM) imaging. Quantitative comparison with conventional parametric GISAXS fits shows good qualitative agreement, with systematic differences explained by domain-shape assumptions and resolved by applying a geometric scaling factor. Simulated SEM-like surfaces derived from neural network outputs reproduce the porous, foam-like nanoscale morphology observed experimentally, reinforcing the method’s credibility. This integrated approach enables real-time, nondestructive, statistically averaged monitoring of bulk nanostructure development in functional coatings, offering a scalable methodology to accelerate the characterization and process control of sustainably manufactured nanostructured titania films for energy-related applications such as photocatalysis and photovoltaics.

Heger, JulianEliah↗

Machine learning dielectric screening for the simulation of excited state properties of molecules and materials

Accurate and efficient calculations of absorption spectra of molecules and materials are essential for the understanding and rational design of broad classes of systems. Solving the Bethe–Salpeter equation (BSE) for electron–hole pairs usually yields accurate predictions of absorption spectra, but it is computationally expensive, especially if thermal averages of spectra computed for multiple configurations are required. We present a method based on machine learning to evaluate a key quantity entering the definition of absorption spectra: the dielectric screening. We show that our approach yields a model for the screening that is transferable between multiple configurations sampled during first principles molecular dynamics simulations; hence it leads to a substantial improvement in the efficiency of calculations of finite temperature spectra. We obtained computational gains of one to two orders of magnitude for systems with 50 to 500 atoms, including liquids, solids, nanostructures, and solid/liquid interfaces. Importantly, the models of dielectric screening derived here may be used not only in the solution of the BSE but also in developing functionals for time-dependent density functional theory (TDDFT) calculations of homogeneous and heterogeneous systems. Overall, our work provides a strategy to combine machine learning with electronic structure calculations to accelerate first principles simulations of excited-state properties.

36 MATERIALS SCIENCE↗

Explicit environmental and vibronic effects in simulations of linear and nonlinear optical spectroscopy

Accurately simulating the linear and nonlinear electronic spectra of condensed phase systems and accounting for all physical phenomena contributing to spectral line shapes presents a significant challenge. Vibronic transitions can be captured through a harmonic model generated from the normal modes of a chromophore, but it is challenging to also include the effects of specific chromophore–environment interactions within such a model. We work to overcome this limitation by combining approaches to account for both explicit environment interactions and vibronic couplings for simulating both linear and nonlinear optical spectra. We present and show results for three approaches of varying computational cost for combining ensemble sampling of chromophore–environment configurations with Franck–Condon line shapes for simulating linear spectra. We present two analogous approaches for nonlinear spectra. Simulated absorption spectra and two-dimensional electronic spectra (2DES) are presented for the Nile red chromophore in different solvent environments. Furthermore, employing an average Franck–Condon or 2DES line shape appears to be a promising method for simulating linear and nonlinear spectroscopy for a chromophore in the condensed phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Inexact iterative numerical linear algebra for neural network-based spectral estimation and rare-event prediction

Understanding dynamics in complex systems is challenging because there are many degrees of freedom, and those that are most important for describing events of interest are often not obvious. The leading eigenfunctions of the transition operator are useful for visualization, and they can provide an efficient basis for computing statistics, such as the likelihood and average time of events (predictions). Here, we develop inexact iterative linear algebra methods for computing these eigenfunctions (spectral estimation) and making predictions from a dataset of short trajectories sampled at finite intervals. We demonstrate the methods on a low-dimensional model that facilitates visualization and a high-dimensional model of a biomolecular system. Implications for the prediction problem in reinforcement learning are discussed.

Chemistry↗

Ab initio calculations of monopole sum rules: From finite nuclei to infinite nuclear matter

We compute moments of the isoscalar monopole response of 𝑁 = 𝑍 closed-shell nuclei based on chiral nucleon-nucleon plus three-nucleon interactions. We employ the random-phase approximation (RPA) and two ab initio many-body approaches, the in-medium similarity renormalization group (IMSRG) and coupled-cluster theory (CC). In the IMSRG framework, the moments are obtained as ground-state expectation values, whereas in the CC approach, they are evaluated through excited-state calculations. We find good agreement between the IMSRG and CC results across all nuclei studied. RPA provides a reasonable approximation to the correlated methods if the interaction is soft. From the calculated moments, we extract average energies of the monopole response, compute finite-nucleus incompressibilities, and estimate the incompressibility of symmetric nuclear matter by a fit to a leptodermous expansion. Our extrapolated values are lower than those obtained in nuclear-matter calculations with the same interactions, but the values are consistent with phenomenological ranges.

Bonaiti, Francesca [Michigan State Univ., East Lan↗

Lagrangian analysis for turbulent transport in variable-density turbulence

In this report Lagrangian analysis of materially conserved scalars is applied to the problem of turbulent transport in variable-density flows. The consequences of an additional material conserved quantity, the density, is generally not acknowledged and leads to significant and meaningfully different expressions for turbulent transport in the moment equations. The formal Lagrangian analysis produces gradient transport expressions substantially different from those obtained by the physically intuitive “argument by analogy” method used in computational models. Various intuitive arguments, in Favre and Reynolds averaged settings, are contrasted to the formal Lagrangian results. Using expressions from the formal analysis, we derive consistent gradient transport closures for the turbulent transport terms that appear in the first- and second-order Favre moment equations. Results for coupled multispecies turbulent transport are given. The analysis is limited to variable-density turbulence in which the dilatation of the fluctuating velocity is small. The results are applicable to turbulent combustion and to stellar convection problems in which the density fluctuations are on the order of the mean density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Eddy diffusivity operator in homogeneous isotropic turbulence

Here, we use the recently developed macroscopic forcing method to compute the scale-dependent eddy diffusivity characterizing ensemble-averaged scalar and momentum transport in incompressible homogeneous isotropic turbulence. For scales larger than the energy containing eddies, eddy diffusivity is found to be constant and consistent with the Boussinesq approximation. However, for small scales eddy diffusivity is found to vanish inversely proportional to the wave number. Behavior at all scales is reasonably captured by a nonlocal eddy diffusivity operator modeled as D/$\sqrt{I - l^2∇^2}$, where D is the eddy diffusivity in the Boussinesq limit, and l is a constant on the order of the large-eddy length. These results present a direct measurement of eddy diffusivity in turbulence with implications in turbulence modeling.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗