Search NASA⌕ Search

SEARCH · Search NASA

Results for “Random systems”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A review of spatial Markov models for predicting pre-asymptotic and anomalous transport in porous and fractured media

Heterogeneity across a broad range of scales in geologic porous media often manifests in observations of non-Fickian or anomalous transport. While traditional anomalous transport models can successfully make predictions in certain geological systems, increasing evidence suggests that assumptions relating to independent and identically distributed increments constrain where and when they can be reliably applied. A relatively novel model, the Spatial Markov model (SMM), relaxes the assumption of independence. The SMM belongs to the family of correlated continuous time random walks and has shown promise across a wide range of transport problems relevant to natural porous media. It has been successfully used to model conservative as well as more recently reactive transport in highly complex flows ranging from pore scales to much larger scales of interest in geology and subsurface hydrology. In this review paper we summarize its original development and provide a comprehensive review of its advances and applications as well as lay out a vision for its future development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A graphics processing unit accelerated sparse direct solver and preconditioner with block low rank compression

We present the GPU implementation efforts and challenges of the sparse solver package STRUMPACK. The code is made publicly available on github with a permissive BSD license. STRUMPACK implements an approximate multifrontal solver, a sparse LU factorization which makes use of compression methods to accelerate time to solution and reduce memory usage. Multiple compression schemes based on rank-structured and hierarchical matrix approximations are supported, including hierarchically semi-separable, hierarchically off-diagonal butterfly, and block low rank. Here, in this paper, we present the GPU implementation of the block low rank (BLR) compression method within a multifrontal solver. Our GPU implementation relies on highly optimized vendor libraries such as cuBLAS and cuSOLVER for NVIDIA GPUs, rocBLAS and rocSOLVER for AMD GPUs and the Intel oneAPI Math Kernel Library (oneMKL) for Intel GPUs. Additionally, we rely on external open source libraries such as SLATE (Software for Linear Algebra Targeting Exascale), MAGMA (Matrix Algebra on GPU and Multi-core Architectures), and KBLAS (KAUST BLAS). SLATE is used as a GPU-capable ScaLAPACK replacement. From MAGMA we use variable sized batched dense linear algebra operations such as GEMM, TRSM and LU with partial pivoting. KBLAS provides efficient (batched) low rank matrix compression for NVIDIA GPUs using an adaptive randomized sampling scheme. The resulting sparse solver and preconditioner runs on NVIDIA, AMD and Intel GPUs. Interfaces are available from PETSc, Trilinos and MFEM, or the solver can be used directly in user code. We report results for a range of benchmark applications, using the Perlmutter system from NERSC, Frontier from ORNL, and Aurora from ALCF. For a high frequency wave equation on a regular mesh, using 32 Perlmutter compute nodes, the factorization phase of the exact GPU solver is about 6.5× faster compared to the CPU-only solver. The BLR-enabled GPU solver is about 13.8× faster than the CPU exact solver. For a collection of SuiteSparse matrices, the STRUMPACK exact factorization on a single GPU is on average 1.9× faster than NVIDIA’s cuDSS solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Novel In-situ Patterning Technique to Fabricate Single Quantum Dots for Quantum Photonics

Photon sources able to emit single or entangled photon pairs are key components in quantum information systems. Semiconductor quantum dots (QDs) are promising candidates due to their high efficiencies and ease of integration with other photonic or electronic components. State-of-the-art QDs, however, are limited to certain emission wavelengths and specific applications due to material choice constraints and their randomness in shape/size. This project is focused on developing a novel in-situ patterning technique to realize QDs with a broad emission range, shape/size control and the ability to emit single/entangled photons. Our approach has two key elements: (1) In-situ patterning via arsenic-induced nanovoid etching on antimonide surfaces and (2) In-filling of nanovoids to form QDs. By closely controlling the experimental conditions, it is shown that this technique can be used to realize III-V QDs in As 2 - etched nanovoids on a GaSb surface. The exposure to As 2 in terms of substrate temperature, time and flux is found to have a significant impact on the process variables such as nanovoid depth, QD dimensions etc. An in-depth analysis of the etch mechanism reveals that by controlling the As 2 exposure, uniform 3-dimensional nanostructures with varying areal densities can be obtained without an in-filling step. Preliminary optical characterization of these nanostructures shows that these QDs may be relevant for realizing emitters in the telecom wavelength range.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Monte Carlo Dropout Uncertainty Quantification of Long Short-Term Memory Autoencoder Anomaly Detection in a Liquid Sodium Cold Trap

Advanced high-temperature fluid reactors, such as sodium-cooled fast reactors (SFRs) and molten salt–cooled reactors (MSCRs), require coolant purification systems to prevent fluid contamination and local freezing that can lead to plugging. Liquid sodium purification can be achieved with a cold trap, where the sodium temperature is reduced to a near-freezing point to precipitate out impurities. Automation of monitoring of the cold trap performance with machine learning algorithms can aid in early detection of incipient anomalies. An efficient approach to loss-of-coolant–type anomaly detection in a cold trap monitored with more than two dozen thermal-hydraulic sensors consists of a long short-term memory (LSTM) autoencoder. This work develops the uncertainty quantification of the LSTM autoencoder performance for cold trap anomaly detection using the Monte Carlo (MC) dropout method. The MC dropout methodology creates a distribution of sister distributions that all slightly differ from each other because of random neurons being turned off for testing. The variances of the sister network distributions are used to make an uncertainty interval. Our analysis shows that the uncertainty in the autoencoder performance is largest near the peak of the anomaly signal. Using the MC dropout method, we investigate the uncertainty in the anomaly detection with missing sensor inputs. This capability allows the reactor operator to evaluate resilience of the anomaly detection system and to make informed decisions about continuity of operation in the event of sensor failure.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Entanglement phase transitions in random stabilizer tensor networks

Here we explore a class of random tensor network models with “stabilizer” local tensors which we name random stabilizer tensor networks (RSTNs). For RSTNs defined on a two-dimensional square lattice, we perform extensive numerical studies of entanglement phase transitions between volume-law and area-law entangled phases of the one-dimensional boundary states. These transitions occur when either (a) the bond dimension D of the constituent tensors is varied or (b) the tensor network is subject to random breaking of bulk bonds, implemented by forced measurements. In the absence of broken bonds, we find that the RSTN supports a volume-law entangled boundary state with bond dimension D ≥ 3 where D is a prime number, and an area-law entangled boundary state for D = 2. Upon breaking bonds at random in the bulk with probability p, there exists a critical measurement rate p c for each D ≥ 3 above which the boundary state becomes area-law entangled. To explore the conformal invariance at these entanglement transitions for different prime D, we consider tensor networks on a finite rectangular geometry with a variety of boundary conditions, and extract universal operator scaling dimensions via extensive numerical calculations of the entanglement entropy, mutual information, and mutual negativity at their respective critical points. Our results at large D approach known universal data of percolation conformal field theory, while showing clear discrepancies at smaller D, suggesting a distinct entanglement transition universality class for each prime D. We further study universal entanglement properties in the volume-law phase and demonstrate quantitative agreement with the recently proposed description in terms of a directed polymer in a random environment.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Flood Susceptibility Mapping Using Machine Learning and Geospatial-Sentinel-1 SAR Integration for Enhanced Early Warning Systems

This study presents a comprehensive framework for flood susceptibility mapping by integrating geospatial factors with both statistical and machine learning models. Thirteen Flood-related factors, including DEM, slope, TWI, NDVI, etc., are extracted as features of models, and historical flood data derived from Sentinel-1 SAR from 2018 to 2023 are used as the target variables of the models. These datasets are analyzed using a frequency-based statistical model and three machine learning models, including Random Forest, XGBoost, and CNN, to generate flood susceptibility maps. The performance of each model is evaluated through AUC; and SHAP scores are separately generated for Machine learning (ML) models to explain each feature contribution in the ML model. The generated susceptibility maps are validated by high-flood-risk locations monitored by flood sensors, BLE inundation models, and flood-prone areas suggested by the Local Community Task Force. The results indicate that the XGBoost model outperforms all other models, with an AUC of 0.92 and demonstrates the highest alignment with recommended high-flood-risk locations, while the frequency-based statistical model showed the weakest performance with an AUC of 0.65. SHAP value graphs highlight the elevation, slope, and TWI as the most influential features across all models. The susceptibility maps generated by the machine learning model show strong agreement with the BLE map and high-flood-risk areas identified by the local Community Task Force.

Google Engine↗
DOE OSTI2021

Geometrical defect detection for additive manufacturing with machine learning models

This study proposed a scheme based on Machine Learning (ML) models to detect geometric defects of additively manufactured objects. The ML models are trained with synthetic 3D point clouds with defects and then applied to detect defects in actual production. Using synthetic 3D point clouds rather than experimental data could save a huge amount of training time and costs associated with many prints for each design. Besides distance differences of individual points between source and target point clouds, this scheme uses a new concept called “patch” to capture macro-level information about nearby points for ML training and implementation. Numerical comparisons of prediction results on experimental data with different shapes showed that the proposed scheme outperformed the existing Z-difference method in the literature. Five ML methods (Bagging of Trees, Gradient Boosting, Random Forest, K-nearest Neighbors and Linear Supported Vector Machine) were compared under various conditions, such as different point cloud densities and defect sizes. Bagging and Random Forest were found the two best models regarding predictability; and the right patch size was found to be at 20. The proposed ML-based scheme is applicable to in-situ defect detection during additive manufacturing with the aid of a proper 3D data acquisition system.

Additive manufacturing↗
DOE OSTI2025

Automation-Accelerated Electrolyte Design Mitigates Solubility Competition between Redox-Active Molecules and Supporting Salts

In nonaqueous redox-flow batteries (NRFBs), redox-active organic molecules (ROMs) and supporting salts compete for solvation sites, limiting achievable energy density. We combine automated high-throughput experimentation (HTE) with camera-based saturation monitoring and quantitative NMR to measure paired (ROM, salt) solubilities across single and mixed organic solvents. Using 2,1,3-benzothiadiazole (BTZ) with lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) as a model system, we find that a binary m-xylene/acetonitrile mixture dissolves ≈3 M of both BTZ and LiTFSI─surpassing the previously reported 2 M ceiling for neat acetonitrile─by leveraging complementary solvation (MX is BTZ-philic and salt-phobic; ACN stabilizes LiTFSI). A random-forest model (RMSE ≈ 0.24) trained on solvent descriptors highlights log P and salt concentration as dominant predictors and predicts MX/ACN ≈0.3/0.7 (v/v) to be near-optimal. These formulations retain practical viscosity and ∼5 mS·cm –1 conductivity at high loading. In conclusion, the workflow provides a reproducible, data-centric route to NRFB electrolyte design and motivates an open, standardized dual-solute solubility resource for accelerated electrolyte discovery.

Electrolytes↗
DOE OSTI2021

Observation and Control of Unidirectional Ballistic Dynamics of Nanoparticles at a Liquid–Gas Interface by 4D Electron Microscopy

Understanding and controlling the dynamics of active Brownian objects far from equilibrium are fundamentally important for emerging technologies such as artificial micro/nanomotors for drug deliveries and noninvasive microsurgery. However, direct observation and control of unidirectional propulsion of individual nanoscale objects are technically challenging due to the required spatiotemporal resolution. Here, we report in situ visualization and manipulation of unidirectional superfast ballistic dynamics of a single-photon-activated gold nanoparticle (NP) along the liquid–gas interface by four-dimensional electron microscopy (4D EM) at nanometer and nanosecond scales. We observed that, upon repetitive femtosecond laser excitation, the NP at the liquid–gas interface exhibits a continuously superfast unidirectional translation with a linear dependence of its root mean squared velocity (νrms) on either the laser fluence or repetition rate. Under a single femtosecond pulse excitation, the NP exhibits a superfast ballistic translation at the nanosecond time scale. Combined experiment and physical modeling reveals that the superfast unidirectional, ballistic translation is driven by unidirectional random impulsive forces arising from the nanobubbles (NBs) induced by enhanced laser heating as a result of plasmonic excitation, which is controllable by tuning the laser characteristics. Furthermore, this directional plasmonic NB-propulsion mechanism sheds light on the design of light-controllable artificially intelligent micro/nanomotor systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Limit cycle oscillations in the zonal-flow-catalyzed interactions of ion-temperature-gradient turbulence

Limit-cycle oscillations are studied for ion temperature gradient turbulence, which, in the absence of large diamagnetic (mean) shear flows, saturates through energy transfer from unstable modes to large-scale stable modes via zonal-flow intermediary modes. Oscillations of zonal flow and turbulence levels are strongly constrained by the reactive, largely non-dissipative character of the zonal flows. Since existing predator–prey models for observed oscillations in experiments do not include energy transfer through zonal flows to stable modes, low-order fluid models with this physics are constructed and investigated. A simple three-wave truncation produces low-amplitude zonal flows that slowly oscillate around a zero mean, with turbulence oscillations between coupled wavenumbers that exceed linear frequencies by orders of magnitude. This inconsistency with experimental observations is caused by the weak non-linear drive of zonal flows in three-wave systems and the lack of multiple-wavenumber turbulent interactions. A more comprehensive model that preserves multiple wavenumber interactions within the context of conservative zonal-flow-mediated energy transfer to stable modes accurately reflects observed dynamics when the phase between stable and unstable modes is occasionally randomized.

Li, P. -Y. (ORCID:0000000295254171)↗
DOE OSTI2021

Substantial reduction of the anisotropy in the critical current densities J c of Ni-doped CaKFe 4 As 4 single crystals by chemical and irradiation-induced disorder

Here, we demonstrate that the anisotropy in the critical current densities, J c , of CaK(Fe 1–x Ni x ) 4 As 4 (x = 0, 0.015, 0.025, and 0.030) single crystals reduces with increasing chemical and irradiation-induced disorder. The magnetic field dependences of J c are analyzed by performing magnetization measurements with H applied parallel and perpendicular to the crystallographic c-axis. The results show that undoped crystals display large anisotropies in Jc due to an enhancement of the vortex pinning with H applied parallel to the crystallographic ab-planes. This anisotropy reduces substantially as Ni addition increases. Moreover, we found that random disorder introduced by proton irradiation enhances mainly the vortex pinning for H parallel to the c-axis. Consequently, using adequate fluencies, the vortex pinning at low temperatures in both undoped and doped samples becomes close to isotropic. These results make the CaKFe 4 As 4 system promising for applications that require isotropic J c under magnetic fields.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Hot-lines topology and the fate of the spin resonance mode in three-dimensional unconventional superconductors

In the quasi-two-dimensional (quasi-2D) copper- and iron-based superconductors, the onset of superconductivity is accompanied by a prominent peak in the magnetic spectrum at momenta close to the wave-vector of the nearby antiferromagnetic state. Such a peak is well described in terms of a spin resonance mode, i.e., a spin-1 exciton theoretically predicted for quasi-2D superconductors with a sign-changing gap. The same theories, however, indicate that such a resonance mode should be absent in a three-dimensional (3D) system with a spherical Fermi surface. This raises the question of the fate of the spin resonance mode in layered unconventional superconductors that are not strongly anisotropic, such as certain heavy-fermion compounds and potentially the newly discovered nickelate superconductor NdNiO 2 . Here, we use the random phase approximation to calculate the dynamical spin susceptibility of 3D superconductors with a d x 2 –y 2 -wave gap symmetry and corrugated cylindrical-like Fermi surfaces. By varying the out-of-plane hopping anisotropy t z /t, we demonstrate that the appearance of a spin resonance mode is determined by the topology of the hot lines, i.e., lines on the Fermi surface that are connected by the magnetic wave vector. For an in-plane antiferromagnetic wave vector, the hot lines undergo a topological transition from open lines to closed loops at a critical t z /t value. The closed hot lines cross the nodal superconducting lines, making the spin resonance mode overdamped and incoherent. In contrast, for an out-of-plane antiferromagnetic wave vector, the hot lines remain open and the spin resonance mode remains sharp. We discuss the experimental implications of our results for the out-of-plane dispersion of the spin resonance mode and, more generally, for inelastic neutron scattering experiments on unconventional superconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The variability structure function of the highest luminosity quasars on short time-scales

ABSTRACT The stochastic photometric variability of quasars is known to follow a random-walk phenomenology on emission time-scales of months to years. Some high-cadence rest-frame optical monitoring in the past has hinted at a suppression of variability amplitudes on shorter time-scales of a few days or weeks, opening the question of what drives the suppression and how it might scale with quasar properties. Here, we study a few thousand of the highest luminosity quasars in the sky, mostly in the luminosity range of $L_{\rm bol}$$=[46.4, 47.3]$ and redshift range of $z=[0.7, 2.4]$. We use a data set from the NASA/Asteroid Terrestrial-impact Last Alert System facility with nightly cadence, weather permitting, which has been used before to quantify strong regularity in longer term rest-frame-UV variability. As we focus on a careful treatment of short time-scales across the sample, we find that a linear function is sufficient to describe the UV variability structure function. Although the result can not rule out the existence of breaks in some groups completely, a simpler model is usually favoured under this circumstance. In conclusion, the data are consistent with a single-slope random walk across rest-frame time-scales of $\Delta t=[10, 250]$ d.

Tang, Ji-Jia (ORCID:0000000218600886)↗
DOE OSTI2024

Destabilizing high-capacity high entropy hydrides via earth abundant substitutions: From predictions to experimental validation

The vast chemical space of high entropy alloys (HEAs) makes trial-and-error experimental approaches for materials discovery intractable and often necessitates data-driven and/or first principles computational insights to successfully target materials with desired properties. In the context of materials discovery for hydrogen storage applications, a theoretical prediction-experimental validation approach can vastly accelerate the search for substitution strategies to destabilize high-capacity hydrides based on benchmark HEAs, e.g. TiVNbCr alloys. Here, in this study, machine learning predictions, corroborated by density functional theory calculations, predict substantial hydride destabilization with increasing substitution of earth-abundant Fe content in the (TiVNb) 75 Cr 25-x Fe x system. The as-prepared alloys crystallize in a single-phase bcc lattice for limited Fe content x < 7, while larger Fe content favors the formation of a secondary C14 Laves phase intermetallic. Short range order for alloys with x < 7 can be well described by a random distribution of atoms within the bcc lattice without lattice distortion. Hydrogen absorption experiments performed on selected alloys validate the predicted thermodynamic destabilization of the corresponding fcc hydrides and demonstrate promising lifecycle performance through reversible absorption/desorption. This demonstrates the potential of computationally expedited hydride discovery and points to further opportunities for optimizing bcc alloy ↔ fcc hydrides for practical hydrogen storage applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Molecular dynamics studies of helium bubble effects on grain boundary fracture vulnerabilities in an Fe 70 Ni 11 Cr 19 –1%H austenitic stainless steel

Comprehensive molecular dynamics tensile test simulations have been performed to study the delamination processes of seven different grain boundaries / cleavage planes (Σ1{111}, Σ3{111}, Σ5{100}, Σ7{111}, Σ9{411}, Σ11{311}, and R{100}/{411}) containing a helium bubble. Combinations of a variety of conditions are explored including different strain rates, system dimensions, bubble density, bubble radius, bubble pressure, and temperature. We found that in general, grain boundaries absorb less energies with decreasing strain rate but increasing bubble areal density, bubble pressure, bubble radius, and temperature. The propensity of grain boundary delamination is sensitive to grain boundary type: The random grain boundary R{100}/{411} is one of the most brittle boundaries whereas the Σ1{111} cleavage plane and the Σ3{111} twin boundary are two of the toughest boundaries. The sorted list of grain boundary fracture vulnerability obtained from our dynamic tensile test simulations differs from the one obtained from our decohesion energy calculations, confirming the important role of plastic deformation during fracture. Finally, detailed mechanistic analyses are performed to interpret the simulated results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

A flexible class of priors for orthonormal matrices with basis function-specific structure

Statistical modeling of high-dimensional matrix-valued data motivates the use of a low-rank representation that simultaneously summarizes key characteristics of the data and enables dimension reduction. Low-rank representations commonly factor the original data into the product of orthonormal basis functions and weights, where each basis function represents an independent feature of the data. However, the basis functions in these factorizations are typically computed using algorithmic methods that cannot quantify uncertainty or account for basis function correlation structure a priori. While there exist Bayesian methods that allow for a common correlation structure across basis functions, empirical examples motivate the need for basis function-specific dependence structure. We propose a prior distribution for orthonormal matrices that can explicitly model basis function-specific structure. The prior is used within a general probabilistic model for singular value decomposition to conduct posterior inference on the basis functions while accounting for measurement error and fixed effects. We discuss how the prior specification can be used for various scenarios and demonstrate favorable model properties through synthetic data examples. Finally, we apply our method to two-meter air temperature data from the Pacific Northwest, enhancing our understanding of the Earth system’s internal variability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Native Architecture of Wheat Straw Cell Walls: A Unified Model from X-ray Scattering and Solid-State NMR

Plant secondary cell walls constitute the dominant reservoir of renewable biomass, comprising tightly packed cellulose, hemicellulose, and lignin at the nanoscale. Recent advances in solid-state NMR spectroscopy and the availability of small-angle X-ray scattering for biomass characterization have led to an accumulation of experimental data on cell wall organization, yet no explicit structure model has simultaneously satisfied both Xray and NMR observations. Using wheat straw as a model system, we propose a structural framework consistent with current knowledge of cellulose biosynthesis, X-ray scattering data, and one- and two-dimensional 13 C solid-state NMR spectra. In this model, 18-chain elementary fibrils align in parallel and populate the cross-section at random. Arabinose-substituted xylan shows no conformational dependence for cellulose-binding in wheat, and only a minor fraction of 2-fold xylan appears in close proximity to cellulose, unlike in Arabidopsis, where xylan is more tightly attached to the cellulose surface. While NMR data cannot unambiguously resolve the internal arrangement of the 18 glucan chains, X-ray scattering profiles uniquely constrain the fibril size and exclude the possibility of tight bundling in the intact walls. The specific interaction between the matrix polymers and the cellulose elementary fibrils must be reconsidered in light of the small interfibril spaces, which bring the matrix components into spatial proximity with cellulose even in the absence of attractive interactions. These findings provide fundamental molecular-level insight into cellulose fibril architecture and matrix−polymer interactions, resolving longstanding discrepancies between spectroscopic and scattering data and advancing our understanding of biopolymer assembly into structurally and functionally versatile lignocellulosic biomaterials.

Carbohydrates↗
DOE OSTI2024

Heisenberg-limited Hamiltonian learning for interacting bosons

We develop a protocol for learning a class of interacting bosonic Hamiltonians from dynamics with Heisenberg-limited scaling. For Hamiltonians with an underlying bounded-degree graph structure, we can learn all parameters with root mean square error ϵ using ${\mathcal{O}}(1/\epsilon )$ total evolution time, which is independent of the system size, in a way that is robust against state-preparation and measurement error. In the protocol, we only use bosonic coherent states, beam splitters, phase shifters, and homodyne measurements, which are easy to implement on many experimental platforms. A key technique we develop is to apply random unitaries to enforce symmetry in the effective Hamiltonian, which may be of independent interest.

computer science↗