SEARCH · Search NASA

Results for “randomized sampling”

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 847 records · Page 47

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Mapping Rare Earths and Toxics in E-Waste via Hyperspectral Imaging and Machine Learning

Electronic waste (e-waste) presents a mounting challenge to environmental sustainability due to its complex composition, which includes high-value rare earth elements, hazardous organic compounds, and non-recyclable plastics. Accurate and scalable material classification is essential for enabling efficient resource recovery and safe recycling practices. This study introduces a confidence-aware classification pipeline that combines mid-infrared hyperspectral imaging (HSI), spectral angle mapping (SAM), and iterative machine learning to perform pixel-level material identification across e-waste devices. A curated spectral library encompassing artificial materials (e.g., plastic iron oxide, galvanized metals), minerals (e.g., allanite, hematite), and organic compounds (e.g., benzanthracene, toluene) was used to generate pseudo-labels, each assigned a confidence score based on SAM-derived spectral similarity. High-confidence samples from seven consumer electronics—digital cameras, keyboards, laptop fans, modems, motherboards, TV remotes, and speakers—were iteratively expanded and classified using models such as Support Vector Machine (SVM), Random Forest, Gradient Boosting Classifier, Partial Least Squares Discriminant Analysis (PLSDA) and Logistic Regression. The best-performing classifiers achieved macro F1 scores approaching 1.0. Results revealed widespread plastic content (dominated by plastic iron oxide), the presence of rare earth-bearing minerals like cerium-containing allanite, and pervasive detection of hazardous organics such as benzanthracene. Principal Component Analysis (PCA) visualizations and confusion matrices confirmed high separability and robust classification performance. This methodology enables precise, non-destructive, and scalable classification of heterogeneous e-waste streams. It supports automated, hazard-aware sorting in recycling workflows, facilitating selective recovery of critical materials and compliance with circular economy goals. The confidence-aware framework provides a foundation for real-time deployment in industrial settings, offering significant implications for smart e-recycling infrastructure and policy-driven material stewardship.

Circular economy↗

DOE OSTI2019

Multi-scale structural analysis of proteins by deep semantic segmentation

Abstract Motivation Recent advances in computational methods have facilitated large-scale sampling of protein structures, leading to breakthroughs in protein structural prediction and enabling de novo protein design. Establishing methods to identify candidate structures that can lead to native folds or designable structures remains a challenge, since few existing metrics capture high-level structural features such as architectures, folds and conformity to conserved structural motifs. Convolutional Neural Networks (CNNs) have been successfully used in semantic segmentation—a subfield of image classification in which a class label is predicted for every pixel. Here, we apply semantic segmentation to protein structures as a novel strategy for fold identification and structure quality assessment. Results We train a CNN that assigns each residue in a multi-domain protein to one of 38 architecture classes designated by the CATH database. Our model achieves a high per-residue accuracy of 90.8% on the test set (95.0% average per-class accuracy; 87.8% average per-structure accuracy). We demonstrate that individual class probabilities can be used as a metric that indicates the degree to which a randomly generated structure assumes a specific fold, as well as a metric that highlights non-conformative regions of a protein belonging to a known class. These capabilities yield a powerful tool for guiding structural sampling for both structural prediction and design. Availability and implementation The trained classifier network, parser network, and entropy calculation scripts are available for download at https://git.io/fp6bd, with detailed usage instructions provided at the download page. A step-by-step tutorial for setup is provided at https://goo.gl/e8GB2S. All Rosetta commands, RosettaRemodel blueprints, and predictions for all datasets used in the study are available in the Supplementary Information. Supplementary information Supplementary data are available at Bioinformatics online.

Eguchi, Raphael R.↗

DOE OSTI2016

Mapping growth windows in quaternary perovskite oxide systems by hybrid molecular beam epitaxy

Requisite to growing stoichiometric perovskite thin films of the solid-solution A ' 1-x A x B O 3 by hybrid molecular beam epitaxy is understanding how the growth conditions interpolate between the end members A'B O 3 and AB O 3 , which can be grown in a self-regulated fashion, but under different conditions. Using the example of La 1-x Sr x VO 3 , the two-dimensional growth parameter space that is spanned by the flux of the metal-organic precursor vanadium oxytriisopropoxide and composition, x , was mapped out. The evolution of the adsorption-controlled growth window was obtained using a combination of X-ray diffraction, atomic force microscopy, reflection high-energy electron-diffraction (RHEED), and Rutherford backscattering spectroscopy. It is found that the stoichiometric growth conditions can be mapped out quickly with a single calibration sample using RHEED. Once stoichiometric conditions have been identified, the out-of-plane lattice parameter can be utilized to precisely determine the composition x . This strategy enables the identification of growth conditions that allow the deposition of stoichiometric perovskite oxide films with random A-site cation mixing, which is relevant to a large number of perovskite materials with interesting properties, e.g., high-temperature superconductivity and colossal magnetoresistance, that emerge in solid solution A ' 1-x A x B O 3 .

36 MATERIALS SCIENCE↗

DOE OSTI2017

INFLUENCE OF THE GALACTIC GRAVITATIONAL FIELD ON THE POSITIONAL ACCURACY OF EXTRAGALACTIC SOURCES

We investigate the influence of random variations of the Galactic gravitational field on the apparent celestial positions of extragalactic sources. The basic statistical characteristics of a stochastic process (first-order moments, an autocorrelation function and a power spectral density) are used to describe a light ray deflection in a gravitational field of randomly moving point masses as a function of the source coordinates. We map a 2D distribution of the standard deviation of the angular shifts in positions of distant sources (including reference sources of the International Celestial Reference Frame) with respect to their true positions. For different Galactic matter distributions the standard deviation of the offset angle can reach several tens of μ as (microarcsecond) toward the Galactic center, decreasing down to 4–6 μ as at high galactic latitudes. The conditional standard deviation (“jitter”) of 2.5 μ as is reached within 10 years at high galactic latitudes and within a few months toward the inner part of the Galaxy. The photometric microlensing events are not expected to be disturbed by astrometric random variations anywhere except the inner part of the Galaxy as the Einstein–Chvolson times are typically much shorter than the jittering timescale. While a jitter of a single reference source can be up to dozens of μ as over some reasonable observational time, using a sample of reference sources would reduce the error in relative astrometry. The obtained results can be used for estimating the physical upper limits on the time-dependent accuracy of astrometric measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Analyzing the Validity of Brazilian Testing Using Digital Image Correlation and Numerical Simulation Techniques

Characterizing the mechanical behavior of rocks plays a crucial role to optimize the fracturing process in unconventional reservoirs. However, due to the intrinsic anisotropy and heterogeneity in unconventional resources, fracture process prediction remains the most significant challenge for sustainable and economic hydrocarbon production. During the deformation tracking under compression, deploying conventional methods (strain gauge, extensometer, etc.) is insufficient to measure the deformation since the physical attachment of the device is restricted to the size of the sample, monitoring limited point-wise deformation, producing difficulties in data retrieval, and a tendency to lose track in failure points, etc. Where conventional methods are limited, the application of digital image correlation (DIC) provides detailed and additional information of strain evolution and fracture patterns under loading. DIC is an image-based optical method that records an object with a camera and monitors the random contrast speckle pattern painted on the facing surface of the specimen. To overcome the existing limitations, this paper presents numerical modeling of Brazilian disc tests under quasi-static conditions to understand the full-field deformation behaviors and finally, it is validated by DIC. As the direct tensile test has limitations in sample preparation and test execution, the Brazilian testing principle is commonly used to evaluate indirectly the tensile strength of rocks. The two-dimensional numerical model was built to predict the stress distribution and full-field deformation on Brazilian disc under compression based on the assumptions of a homogenous, isotropic and linear elastic material. The uniaxial compression test was conducted using the DIC technique to determine the elastic properties of Spider Berea sandstone, which were used as inputs for the simulation model. The model was verified by the analytical solution and compared with the digital image correlation. The numerical simulation results showed that the solutions matched reasonably with the analytical solutions where the maximum deviation of stress distribution was obtained as 14.59%. The strain evolution (normal and shear strains) and displacements along the central horizontal and vertical planes were investigated in three distinguishable percentages of peak loads (20%, 40%, and 90%) to understand the deformation behaviors in rock. The simulation results demonstrated that the strain evolution contours consistently matched with DIC generated contours with a reasonable agreement. The changes in displacement along the central horizontal and vertical planes showed that numerical simulation and DIC generated experimental results were repeatable and matched closely. In terms of validation, Brazilian testing to measure the indirect tensile strength of rocks is still an issue of debate. The numerical model of fracture propagation supported by digital image correlation from this study can be used to explain the fracturing process in the homogeneous material and can be extended to non-homogeneous cases by incorporating heterogeneity, which is essential for rock mechanics field applications.

58 GEOSCIENCES↗

DOE OSTI2017

A two-stage adaptive stochastic collocation method on nested sparse grids for multiphase flow in randomly heterogeneous porous media

A new computational method is proposed for efficient uncertainty quantification of multiphase flow in porous media with stochastic permeability. For pressure estimation, it combines the dimension-adaptive stochastic collocation method on Smolyak sparse grids and the Kronrod–Patterson–Hermite nested quadrature formulas. For saturation estimation, an additional stage is developed, in which the pressure and velocity samples are first generated by the sparse grid interpolation and then substituted into the transport equation to solve for the saturation samples, to address the low regularity problem of the saturation. Numerical examples are presented for multiphase flow with stochastic permeability fields to demonstrate accuracy and efficiency of the proposed two-stage adaptive stochastic collocation method on nested sparse grids.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Decayheatml

This code is designed to predict and analyze the decay heat generated in molten salt reactors (MSRs) using a hybrid approach that combines machine learning and segmented polynomial fitting. The accurate prediction of decay heat is essential for reactor safety and the optimization of spent fuel storage. The code operates through several key components: 1) Data Architecture: It incorporates a modular data architecture that handles various MSR-specific operational parameters such as power density, humidity content, and air ingress. These parameters are sampled using Sobol sequences to ensure comprehensive coverage of operational uncertainties. 2) Machine Learning Framework: The code employs a diverse set of machine learning models, including polynomial regression, decision trees, random forests, gradient boosting, support vector regression, k-nearest neighbors, multi-layer perceptrons, and symbolic regression. These models are trained to predict decay heat over a wide temporal range, from immediate shutdown up to 10,000 years. 3) Region-Optimized Training: The temporal domain is divided into multiple regions, each modeled separately to capture distinct decay heat characteristics across different time scales. This approach significantly improves the accuracy and interpretability of predictions. 4) Segmented Polynomial Interpretation (SPI): The SPI method translates machine learning predictions into piecewise polynomial equations. These equations are physically interpretable and can be directly integrated into existing engineering workflows and safety analyses. 5) Front-End Interfaces: The code includes both a Jupyter notebook interface for research development and a Streamlit web application for operational deployment. These interfaces allow users to interactively explore decay heat predictions, adjust operational parameters, and visualize results in real-time. 6) Applications: The framework supports various applications, including safety system validation and spent fuel container optimization. It enables real-time evaluation of worst-case decay heat scenarios, informing the design of passive safety systems and optimizing container designs for long-term storage. Overall, this code provides a robust, accurate, and user-friendly tool for predicting decay heat in MSRs, enhancing reactor safety, and optimizing spent fuel management.

Retamales, Mauricio Eduardo Tano [Idaho National L↗

DOE OSTI2020

Risk Assessment of Rare Events in Probabilistic Power Flow via Hybrid Multi-Surrogate Method

The risks associated with rare events threatening the security of power system operation are of paramount importance to power system planners and operators. To analyze the risks caused by high-impact, low-frequency rare events, an immensely large number of samples are typically required for the Monte-Carlo (MC) method on the high-fidelity power system model to achieve a sufficient accuracy, thereby rendering this approach computationally prohibitive. To handle this problem efficiently, it is desirable to construct a surrogate model for the power system response. However, the straightforward MC sampling of the low-fidelity surrogate can lead to biased results in the low-probability tail regions that are vital to risk assessment. Moreover, a single surrogate is unable to handle the topology uncertainties caused by random branch outages. To overcome these issues, we propose a hybrid multi-surrogate (HMS) method based on the polynomial chaos expansion (PCE) with low-probability tail events reevaluated by the high-fidelity model through a probabilistic analysis. This method improves the computational efficiency of the MC method for rare-event risk assessment by leveraging multi-fidelity models while retaining the desired accuracy. Simulations conducted in three test systems verify the excellent performances of the HMS method.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2018

Extended live-tracking and quantitative characterization of wound healing and cell migration with SiR-Hoechst

Highlights: • SiR-Hoechst is a nuclear dye that enables cell tracking without affecting migration. • Single cell analysis provides multiple migration metrics during wound healing. • Cell density-based report of wound healing is less susceptible to subjective bias. • Cell directionality is an important contributor of early wound recovery. Cell migration is essential to many life processes, including immune response, tissue repair, and cancer progression. A reliable quantitative characterization of the cell migration can therefore aid in the high throughput screening of drug efficacy in wound healing and cancer treatments. In this work, we report what we believe is the first use of SiR-Hoechst for extended live tracking and automated analysis of cell migration and wound healing. We showed through rigorous statistical comparisons that this far-red label does not affect migratory behavior. We observed excellent automated tracking of random cell migration, in which the motility parameters (speed, displacement, path length, directionality ratio, persistence time, and direction autocorrelation) obtained closely match those obtained from manual tracking. We also present an analysis framework to characterize the healing of a scratch wound from the perspective of single cells. The use of SiR-Hoechst is advantageous for the crowded environments in wound healing assays because as long as cell nuclei do not overlap, continuous tracking can be maintained even if there is cell-cell contact. In this paper, we report wound recovery based on the number of cells migrating into the wound over time, normalized by the initial cell count prior to the infliction of the wound. This normalized cell count approach is impervious to operator bias during the arbitration of wound edges and is also robust against variability that arises due to differences in the cell density of different samples. Additional wound healing characteristics were also defined based on the evolution of cell speed and directionality during healing. Not unexpected, the wound healing cells exhibited much higher tendency to maintain the same migratory direction in comparison to the randomly migrating cells. The use of SiR-Hoechst thus greatly simplified the automation of single cell and whole population analysis with high spatial and temporal resolution over extended periods of time.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Effect of hatch spacing and laser power on microstructure, texture, and thermomechanical properties of laser powder bed fusion (L-PBF) additively manufactured NiTi

This study systematically evaluates the effects of laser powder bed fusion additive manufacturing (L-PBF-AM) parameters (hatch spacing and laser power) on the thermomechanical behavior and microstructure of Ni 50.8 Ti 49.2 shape memory alloy. The samples were fabricated with hatch spacings from 40 to 240 µm and laser powers of 50 and 100 W at a constant scanning speed of 125 mm/s, resulting in parts with volumetric energy density levels from 55 to 666 J/mm 3 and two sets of linear energy densities of 0.4 and 0.8 J/mm. The results showed a reduced melt pool size and discontinuity of scan tracks with decreased laser power. Additionally, the porosity level was increased with larger hatch spacing and lower laser power. More notably, the transformation temperatures increased, and the critical stress, recoverable strain, and functional stability of samples improved with lower hatch spacing, where the recovery ratio of up to 90% was observed, regardless of the employed laser power. This study also discussed the relationship between the fabrication process and texture formation in the L-PBF-AM process. In conclusion, the advantage of L-PBF-AM was revealed in tailoring the microstructure from highly textured samples in [1 1 1] or [0 0 1] direction when hatch spacing lower than laser beam focused was employed, to the appearance of equiaxed solidification front with island grains and random orientations.

36 MATERIALS SCIENCE↗

DOE OSTI2023

NMR study of Ni 50+x ⁢Ti 50-x strain glasses

Here, we studied Ni 50+x ⁢Ti 50-x with compositions up to x = 2, performing 47 Ti and 49 Ti nuclear magnetic resonance (NMR) measurements from 4 to 400 K. For large x in this system, a strain glass appears in which frozen ferroelastic nanodomains replace the displacive martensite structural transition. Here, we demonstrate that NMR can provide an extremely effective probe of the strain-glass freezing process, with large changes in NMR line shape due to the effects of random strains which become motionally narrowed at high temperatures. At the same time with high-resolution x-ray diffraction we confirm the lack of structural changes in x ≥ 1.2 samples, while we show that there is little change in the electronic behavior across the strain-glass freezing temperature. NMR spin-lattice relaxation time (T 1 ) measurements provide a further measure of the dynamics of the freezing process, and indicate a predominantly thermally activated behavior both above and below the strain-glass freezing temperature. We show that the strain-glass results are consistent with a very small density of critically divergent domains undergoing a Vogel-Fulcher-type freezing process, coexisting with domains exhibiting faster dynamics and stronger pinning.

36 MATERIALS SCIENCE↗

DOE OSTI2022

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called â€¢ BCC_opt â€¢ FCC_opt â€¢ HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string â€œA_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters â€œA_B_Câ€, three files are provided: â€¢ A_B_C_opt.xyz: The optimized geometries in xyz format â€¢ A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. â€¢ A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. â€¢ A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Integration of LIBS with Machine Learning for Real-Time Monitoring of Feedstock in H 2 Gasification Applications

This project, funded by the U.S. Department of Energy (DOE) – Office of Fossil Energy under Award Number DE-FE0032177, aimed to assess the feasibility of an integrated Laser-Induced Breakdown Spectroscopy (LIBS) system with advanced machine learning (ML) models for real-time characterization and potential control of hydrogen gasifiers running on waste materials as feedstocks. This was a multidisciplinary effort that encompassed the acquisition and standardized analysis of individual and blended feedstocks—comprising biomass, coal waste, and plastic waste, followed by the development of a dynamic LIBS bench system for material sample analysis and development of predictive ML models. Comprehensive laboratory testing enabled the creation of a robust elemental dataset that served as the foundation for ML model training. Techniques such as Random Forest, Gradient Boosting, Support Vector Regression, and Neural Networks were employed to predict key feedstock properties, including higher heating value (HHV), moisture content, thermal conductivity, and ash composition with high accuracy. The results were validated against experimental data and demonstrated strong potential for real-time application in gasifier control systems. The project concluded with a study on the integration of the LIBS+ML approach for gasifier control and a techno-economic analysis of the implementation of the approach into hydrogen (H 2 ) gasification systems. Dissemination of results was carried out at a DOE meeting. This work establishes a scalable framework for automated, in-line feedstock quality assessment, offering significant implications for process optimization and emissions reduction in hydrogen production.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2017

Methods of Reverberation Mapping. I. Time-lag Determination by Measures of Randomness

A class of methods for measuring time delays between astronomical time series is introduced in the context of quasar reverberation mapping, which is based on measures of randomness or complexity of the data. Several distinct statistical estimators are considered that do not rely on polynomial interpolations of the light curves nor on their stochastic modeling, and do not require binning in correlation space. Methods based on von Neumann’s mean-square successive-difference estimator are found to be superior to those using other estimators. An optimized von Neumann scheme is formulated, which better handles sparsely sampled data and outperforms current implementations of discrete correlation function methods. This scheme is applied to existing reverberation data of varying quality, and consistency with previously reported time delays is found. In particular, the size–luminosity relation of the broad-line region in quasars is recovered with a scatter comparable to that obtained by other works, yet with fewer assumptions made concerning the process underlying the variability. The proposed method for time-lag determination is particularly relevant for irregularly sampled time series, and in cases where the process underlying the variability cannot be adequately modeled.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Representative elementary volume (REV) of cementitious materials from three-dimensional pore structure analysis

The representative elementary volume (REV) is a fundamental property of a material, but no direct measurements exist for cementitious materials. In this paper, the REV of cement pastes with supplementary cementitious materials (GGBS, PFA, SF) was determined by analysing the three-dimensional pore structure (> 0.2 μm) using laser scanning confocal microscopy (LSCM). The effect of axial distortion inherent to LSCM on 3D pore structure was also investigated. A range of 3D pore parameters was measured using skeletonisation, maximal ball and random walker algorithms. Results show that axial distortion has insignificant effects on most parameters except Euler connectivity, average pore and throat volumes and directional diffusion tortuosities. Most pore parameters become independent of sampling volume at ≈ 60{sup 3} μm{sup 3} except diffusion tortuosities and formation factor. The REV for porosity calculated based on a statistical approach at eight realisations and 5% relative error was found to be ≈ 100{sup 3} μm{sup 3}.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Ensemble variational Fokker-Planck methods for data assimilation

Particle flow filters solve Bayesian inference problems by smoothly transforming a set of particles into samples from the posterior distribution. Particles move in state space under the flow of an McKean-Vlasov-Itˆo process. This work introduces the Variational Fokker-Planck (VFP) framework for data assimilation, a general approach that includes previously known particle flow filters as special cases. The McKean-Vlasov-Itˆo process that transforms particles is defined via an optimal drift that depends on the selected diffusion term. It is established that the underlying probability density - sampled by the ensemble of particles - converges to the Bayesian posterior probability density. For a finite number of particles the optimal drift contains a regularization term that nudges particles toward becoming independent random variables. Based on this analysis, we derive computationally-feasible approximate regularization approaches that penalize the mutual information between pairs of particles, and avoid particle collapse. Moreover, the diffusion plays a role akin to a particle rejuvenation approach that aims to alleviate particle collapse. The VFP framework is very flexible. Different assumptions on prior and intermediate probability distributions can be used to implement the optimal drift, and localization and covariance shrinkage can be applied to alleviate the curse of dimensionality. A robust implicit-explicit method is discussed for the efficient integration of stiff McKean- Vlasov-Itˆo processes. Here, the effectiveness of the VFP framework is demonstrated on three progressively more challenging test problems, namely the Lorenz ’63, Lorenz ’96 and the quasi-geostrophic equations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Dynamics of Shape-Persistent Giant Molecules: Zimm-like Melt, Elastic Plateau, and Cooperative Glass-like

Giant molecules are a new class of soft matter having threedimensional (3D) shapes and composed of chemically linked rigid molecular nanoparticles. Structurally, a 3D cluster of molecular nanoparticles can be one giant molecule or a few giant molecules associated together via specific interactions. The dynamics of clusters that are smaller than a critical diameter (~5 nm) presents a power law relaxation exponent of 0.7 at the high frequency region corresponding to segmental dynamics. Such scaling is similar to the result of the Zimm model although those clusters are neither chain-like nor in solution. Clusters that are larger than this critical diameter and formed by the association of giant molecules exhibit an elastic plateau due to caging of individual giant molecules. Herein, we hypothesize that clusters of such a large size cannot move as a whole, even above the glass transition temperature of the sample. They thus are “cooperative glass-like”. A structural cluster of giant molecules could be abstracted as a dynamical cluster consisting of unlinked but cooperatively mobile beads. As derived in the random first-order transition theory, the cluster loses its mobility and reaches the glassy state when the diameter of the cluster is 6 times larger than the bead diameter. In our cases, we estimate that the critical diameter for these clusters is also approximately 6 times the bead diameter based on the glassy shear modulus of giant molecules. Thus, shape-persistent giant molecules may serve as a bridge between polymers and colloids and a platform to mimic cooperative rearrangements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Optimal Energy Management for Microgrid with Stationary and Mobile Storage

This paper studies energy management in a Microgrid (MG) with solar generation, Battery Energy Management System (BESS) and gridable (V2G) EVs. A two-stage stochastic optimization method is proposed to capture the intermittent solar generation and random EV user behaviors. It is subsequently formulated as a Mixed Integer Linear Programming (MILP) problem. To evalutate the proposed method, real solar generation, loads, BESS and EV data is used in Sample Average Approximation (SAA). Computational results show the correctness of the proposed method as well as steady and tightly bounded optimality gap. Comparisons demonstrate that the proposed stochastic method outperforms its deterministic counterpart at the expense of higher computational cost. It is also observed that moderate number of EVs helps to reduce the overall operational cost of the MG, which sheds light on future EV integration to the smart grid.

Battery Energy Management System (BESS)↗