Search NASA⌕ Search

SEARCH · Search NASA

Results for “percolation scaling”

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

Quantifying Single-Ion Transport in Percolated Ionic Aggregates of Polymer Melts

Single-ion conducting polymers such as ionomers are promising battery electrolyte materials, but it is critical to understand how rates and mechanisms of free cation transport depend on the nanoscale aggregation of cations and polymer-bound anions. We perform coarse-grained molecular dynamics simulations of ionomer melts to understand cation mobility as a function of polymer architecture, background relative permittivity, and corresponding ionic aggregate morphology. In systems exhibiting percolated ionic aggregates, cations diffuse via stepping motions along the ionic aggregates. These diffusivities can be quantitatively predicted by calculating the lifetimes of continuous association between oppositely charged ions, which equal the time scales of the stepping (diffusive) motions. In contrast, predicting cation diffusivity for systems with isolated ionic aggregates requires another time scale. Finally, our results suggest that to improve conductivity the Coulombic interaction strength should be strong enough to favor percolated aggregates but weak enough to facilitate ion dissociation.

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

Phase-field simulations of fission gas bubble growth and interconnection in U-(Pu)-Zr nuclear fuel

Abstract The growth and interconnection of fission gas bubbles in the hotter central regions of U-(Pu)-Zr nuclear fuel has been simulated with a phase-field model. The Cahn-Hilliard equation was used to represent the two-phase microstructure, with a single defect species. The volume fraction of the bubble phase and surface area of the bubble-matrix interface were determined during growth and interconnection. Surface area increased rapidly during the initial stages of growth, then slowed and finally decreased as bubble interconnection began and coarsening acted to reduce surface area. The fraction of the bubbles vented to a simulation domain boundary, f V , was quantified as a measure of the microstructure’s interconnectivity and plotted as a function of porosity p . The defect species diffusivity was varied; although changes in diffusivity significantly affected the microstructure, the plots of f V vs. p did not change significantly. The percolation threshold p c was calculated to be approximately 0.26, depending on the assumed diffusivity and using an initial bubble number density based on experimental observations. This is slightly smaller than the percolation threshold for continuum percolation of overlapping 3D spheres. The simulation results were used to parameterize two different engineering-scale swelling models for U-(Pu)-Zr in the nuclear fuel performance code BISON.

Aagesen, Larry K. (ORCID:000000034936676X)↗
DOE OSTI2020

Modeling CO 2 plume migration using an invasion-percolation approach that includes dissolution

Geological carbon dioxide (CO 2 ) sequestration has been proposed as a viable technique to decrease effective emissions of CO 2 into the atmosphere. However, the security of this sequestration is tied to our understanding of the long-term migration of CO 2 in subsurface. The dissolution of CO 2 in the reservoir brine is one of the main long-term trapping mechanisms. However, the assumptions used in large-scale reservoir simulations usually lead to an overestimation of the dissolution volume. We propose a modified approach based on the macroscopic invasion-percolation (MIP) theory that allows the dissolution of CO 2 into brine. We used a high-resolution geological model to compare the Darcy-, modified MIP-, and classic MIP-based simulation results. We observed a significant shrinkage in the nonaqueous plume volume when dissolution is considered during the MIP simulation. In the case of Darcy-based simulation, the plume was completely trapped inside the reservoir with limited migration even after a thousand-year simulation. On the other hand, the majority of the plume migrated out of the simulated reservoir in the case of MIP. Our approach provides more realistic estimation of the dissolution volume and nonaqueous plume extent while leveraging the computational efficiency enjoyed by MIP.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Improved Understanding of Zircaloy Hydrogen Pickup Mechanisms in BWRs. NSUF Report (Final Report)

Using an instrumented autoclave operating at 288°C and 1100 psi, electrical impedance spectroscopy (EIS) measurements were taken from three neutron irradiated fuel channel samples (Zircaloy-2 and Zircaloy-4). The three high priority fuel channel samples were tested at 288°C using a simulated PWR water with the AC frequency range extending down through 10 -4 Hz, with corresponding DC measurements to validate the resistance measurements. EIS measurements were also taken at room temperature using a simulated PWR water from four neutron irradiated fuel channel samples (Zircaloy-2 and Zircaloy-4). and four water rod samples (Zircaloy-2). EIS measurements were also performed on unirradiated Zircaloy-2 tubes that were previously oxidized in autoclaves at various temperatures for different time durations to obtain oxide layers of different thicknesses. These experiments proved to be more challenging than anticipated because of the small size of the samples, which required development of fixtures to allow valid testing of the materials, some of which are radioactive, under pressure and temperature. Equivalent circuits were used to model the data. Some samples resulted in Nyquist plots where the arcs were too close to each other resulting in the inability to separate the different oxide morphological thicknesses. Microstructural characterization in the form of optical microscopy, scanning electron microscopy and (scanning) transmission electron microscopy was performed to document the microstructural features of the oxidized Zr-2 and Zr-4 samples (oxide and metal/oxide interface) and to interpret the EIS results. At the start of this work, it was assumed that despite the laminar cracks and porosity visible in both optical and SEM, the oxide scale is dense near the metal, effectively preventing the water from contacting the metal directly. While it is generally accepted that the large laminar cracks pores allow the water to percolate into the oxide film, it is unknown just how deep it can go, and how interconnected is the porosity, which will dictate how much the porosity is connected electrically to the water within the oxide scale. The EIS results demonstrated that there is in fact a distribution of oxide thicknesses over the various samples, and within each material condition, that are communicating to the liquid. This is derived from the fact that Nyquist plots have an arc radius whose center point is often below the Z’ axis, implying that instead of a single oxide thickness, these oxides have a distribution of thicknesses. This distribution of thicknesses produces arcs that are not clearly separable in most of the samples.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1992

The distribution of nearby rich clusters of galaxies

Redshifts are acquired for a complete sample of 351 Abell clusters with tenth-ranked galaxy magnitudes (m10) less than or equal to 16.5, including 115 entirely new cluster redshifts. Analysis of the spatial distribution of these clusters reveals no clustering on scales larger than 75/h Mpc. The correlation length is 20.0 (+/-4.3)/h Mpc, consistent with the results from other surveys. The frequency of voids with radii of order 60/h Mpc or less is consistent with the form and amplitude of the observed two-point correlation function. There is no significant difference between the clustering properties of clusters with RC = 0 and RC not less than 1. A percolation analysis yields 23 superclusters, 17 of which are new. The superclusters are not significantly elongated in the radial direction; large-scale peculiar motions are of order 1000 km/s or less.

Postman, Marc↗
DOE OSTI2026

Coexistence of Synchronization and Stochasticity in Thermally Coupled Mott Oscillators

Synchronization is conventionally regarded as a mechanism for suppressing variability and enforcing order in coupled systems, from pendula and lasers to neurons and electronic oscillators. Here, we show that synchronization can also embed stochasticity at finer scales. We observe this phenomenon in thermally coupled VO 2 neuristors, where robust in-phase synchronization at the microsecond scale coexists with spike onset fluctuations at the nanosecond scale, with no fixed leader. The coexistence of order and disorder originates from stochastic domain-level physics of the insulator–metal and metal–insulator transitions, where local variations in transition temperature drive cycle-to-cycle randomness in nucleation, percolation, and relaxation. A stochastic domain model reproduces this effect by generating synchronized spike trains with random lead–lag jitter, and experimental interspike interval statistics confirm the persistence of fine-scale variability despite macroscopic phase locking. These findings establish that synchronization and stochasticity can coexist within the same physical platform, revealing hidden disorder within collective order. Furthermore, this insight reframes synchronization as not purely deterministic, but as a universal context where microscopic variability can persist, with implications for electronics, cryptography, and the fundamental physics of order–disorder coexistence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Entanglement Randomness and Gapped Itinerant Carriers in a Frustrated Quantum Magnet

The quantum spin liquid is a state manifesting extraordinary many-body entanglement, and the material NaYbSe 2 is thought to be one of the most promising candidates for its realization. Through low-temperature heat capacity and thermal conductivity measurements, we identify an apparent contradiction familiar to many quantum spin liquid candidates: While entropy is stored by apparently gapless excitations, the itinerant carriers of entropy are gapped. By studying the compositional series NaYb 𝑥⁢ Lu 1−𝑥⁢ Se 2 across a percolation transition of the magnetic lattice, we suggest that this contradiction can be resolved by the presence of entanglement scales of random sizes. Moreover, as we truncate the scale of entanglement by magnetic dilution, we show that the itinerant magnetic entropy carrier in NaYbSe 2 does not arise from a uniform globally entangled spin ground state but rather materializes through the stochastic propagation of boundaries between locally entangled spin objects.

Lyu, Yuanqi 吕源祺 [University of California, Berkele↗
DOE OSTI2021

Calibration of the Diffusivity Predictions of Centipede Using Approximate Bayesian Computation and Applications in Nyx (Engineering Scale) and Xolotl-MARMOT (Meso-Scale) Simulations

Fission gas evolution and release in UO 2 nuclear fuel are important fuel performance metrics and occur in several distinct stages: 1) nucleation, growth and resolution of intra-granular bubbles, 2) diffusion to grain boundaries and 3) nucleation and growth of bubbles at grain boundaries, which eventually form a connected network (percolation) enabling release of gas from grain boundaries through connections to triple junctions, grain edges or free surfaces. The NE-SciDAC project is developing several computational tools to model this problem, which are connected in a hierarchical multi-scale framework. The information transfer in the multi-scale framework is a critical step that, in addition to best-estimates, should include uncertainty quantification. Despite taking a first-principles multi-scale approach, there is a need to perform parameter calibration to ensure consistency with available experimental data. In the present study, uncertainty quantification (UQ) and parameter calibration is demonstrated for one of the lower length scale codes in the multi-scale framework (Centipede) and then the results, including instances of the propagated uncertainties, are used in other codes within the framework, specifically Nyx and Xolotl-MARMOT. We calibrated the model parameters in Centipede, a computer code used to predict diffusivities of uranium (U) and xenon (Xe) in the context of the simulation of fission gas in uranium oxide (UO 2 ) nuclear fuel. The Centipede code depends on 183 parameters, all of which are subject to uncertainty. The three data sets used in our calibration effort are taken from the literature. This data is available as a set of measurements, including measurement errors. Our goal is to calibrate a statistical model that predicts both the value of the measurement and the uncertainty associated with the measurement. We perform a Bayesian calibration of the model parameters using a dedicated approximate Bayesian computation (ABC) likelihood function. To avoid excessive computational costs, we replace the expensive Centipede simulation code by a higher-order surrogate model, constructed using only the 9 most important parameters. These important parameters are identified by a preliminary global sensitivity analysis (GSA) study. Among the important parameters are T0 (the temperature at which UO 2 is perfectly stoichiometric) and Hf_pO2 (the temperature dependence of the oxygen (O) partial pressure) that should be considered as operating conditions to be estimated along with the other parameters. We consider two different cases: one where we define one set of these operating conditions for all data sets, and one where we define distinct operating condition parameters for each data set. The Xe diffusivities predicted by the latter case show distinct features that could not be observed in the former. Next, we use the diffusivity predictions by Centipede as input to Nyx, a reduced order fuel performance code focused on gas behavior alone, in order to estimate quantities associated with inter-granular bubble formation at conditions specified by the experiments. Finally, the diffusivities obtained from the calibrated Centipede runs were used in coupled Xolotl-MARMOT simulations of intra- and inter-granular gas evolution. The results are compared to simulations using the baseline diffusivities from Turnbull et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Trigger-based Incremental Data Processing with Unified Sync and Async Model

In recent years, more and more applications in the cloud have needs to process large-scale on-line datasets, which evolve over time as new entries are added and existing entries are modified. Several programming frameworks, such as Percolator and Oolong, are proposed for such incremental data processing and can achieve efficient processing with an event-driven abstraction. However, these frameworks are inherently asynchronous, leaving the heavy burden of managing synchronization to applications' developers, which further significantly restricts their usabilities. In this study, we propose a trigger-based incremental computing framework in the cloud, called Domino, with both synchronous and asynchronous mechanisms to coordinate parallel triggers. With this new framework, both synchronous and asynchronous applications can be seamlessly developed. Use cases and extensive evaluation results confirm that it can deliver sufficient performance, and also is easy to use for incremental applications in large-scale distributed computing.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Kibble-Zurek mechanism for nonequilibrium phase transitions in driven systems with quenched disorder

Abstract We show that the Kibble-Zurek mechanism applies to nonequilibrium phase transitions found in driven assemblies of superconducting vortices and colloids moving over quenched disorder where a transition occurs from a plastic disordered flowing state to a moving anisotropic crystal. We measure the density of topological defects as a function of quench rate through the nonequilibrium phase transition, and find that on the ordered side of the transition, the topological defect density ρ d scales as a power law, $${\rho }_{{{{{{{{\rm{d}}}}}}}}}\propto 1/{t}_{{{{{{{{\rm{q}}}}}}}}}^{\beta }$$ ρ d ∝ 1 / t q β , where t q is the quench time duration, consistent with the Kibble-Zurek mechanism. We show that scaling with the same exponent holds for varied strengths of quenched disorder and that the exponents fall in the directed percolation universality class. Our results suggest that the Kibble-Zurek mechanism can be applied to the broader class of systems that exhibit absorbing phase transitions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A new phenomenological model to describe root-soil interactions based on percolation theory

In his paper on net primary productivity of terrestrial communities predicted from climatological data, Rosenzweig (1968) argued that variability in productivity is well accounted for by (evapo)-transpiration, and that water from transpiration is, on global scales, the most variable component in the photosynthesis reaction. The goal of this paper is to investigate whether variability in plant growth on local scales and within species is primarily related to transpiration under several scenarios including different terrain curvature, slope aspect, soil characteristics, and climate ranges. Here, we test the hypothesis that this relationship exists because root growth into the surface soil layers (0–2 m) tends to follow paths with minima in resistance, which in turn maximizes water flow and nutrient delivery rates that regulate growth. The set of all connected paths with individual pore-to-pore flow resistances less than a critical, percolating, value forms a cluster with mass fractal dimensionality, d f . We propose that roots follow paths through the 2D percolation cluster, defining the set of all optimal flow paths, making the 2D value of d f from percolation relevant to root fractal dimensionality. The tortuosity of such optimal paths as defined in percolation theory should then relate root length to root radial extent, linking the parameters of root tortuosity and plant productivity. Our analysis of large data sets across species implies that root radial extent and tree height are both proportional to cumulative transpiration until trees approached maximum height, and their growth rates are proportional to the transpiration rate, not to the moisture content. Local variations in tree height as functions of the variables investigated appear generally consistent with deduced variations in transpiration. Here this correlation is investigated more closely in the context of studies addressing individual tree species.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Test of model of equivalence of tree height growth and transpiration rates in percolation-based phenomenology for root-soil interaction

We recently introduced a phenomenological model to represent root-soil interactions that accounts for variability in tree growth on local scales and within species (Hunt et al., 2020). Here, using publicly available datasets, we examine the developed model to provide additional empirical evidence, which also gives novel insight into two particular tree species, Eucalyptus regnans and Eucalyptus sieberri. In particular, we discovered two novel aspects of the model: a power-law time dependence with a different power-law exponent than usually heretofore applied, and apparently a direct equivalence of results for transpiration and growth rates. Also, although our original expectation was that growth and transpiration were merely proportional, comparison of the entire BAAD data set (Falster et al., 2015) for tree height with our model (Hunt et al., 2020; Fig. 1) already suggested approximate equality of these disparate measurements (volume/area, height). In this present manuscript we find further evidence for equivalence of the two measurements and provide additional theoretical basis for a model assumption regarding constancy of sap flow rates.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2000

Effect of Single-Electron Interface Trapping in Decanano MOSFETs: A 3D Atomistic Simulation Study

We study the effect of trapping/detrapping of a single-electron in interface states in the channel of n-type MOSFETs with decanano dimensions using 3D atomistic simulation techniques. In order to highlight the basic dependencies, the simulations are carried out initially assuming continuous doping charge, and discrete localized charge only for the trapped electron. The dependence of the random telegraph signal (RTS) amplitudes on the device dimensions and on the position of the trapped charge in the channel are studied in detail. Later, in full-scale, atomistic simulations assuming discrete charge for both randomly placed dopants and the trapped electron, we highlight the importance of current percolation and of traps with strategic position where the trapped electron blocks a dominant current path.

Asenov, Asen↗
DOE OSTI2025

Stable Symmetry-Protected Topological Phases in Systems with Heralded Noise

Here, we present a family of local quantum channels whose steady-states exhibit stable mixed-state symmetry-protected topological (SPT) order. Motivated by recent experimental progress on "erasure conversion" techniques that allow one to identify ($\textit{herald}$) decoherence processes, we consider open systems with biased erasure noise, which leads to strongly symmetric heralded errors. We utilize this heralding to construct a local correction protocol that effectively confines errors into short-ranged pairs in the steady-state. Using a combination of numerical simulations and mean-field analysis, we show that our protocol stabilizes SPT order against a sufficiently low rate of decoherence. As the rate of heralded noise increases, SPT order is eventually lost through a directed percolation transition. We further find that while introducing unheralded errors destroys SPT order in the limit of long length- and time-scales, the correction protocol is sufficient for ensuring that local SPT order persists, with a correlation length that diverges as $\xi \sim (1-f_e)^{-1/2}$, where $f_e$ is the fraction of errors that are heralded.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1991

Fluctuation driven electroweak phase transition

We examine the dynamics of the electroweak phase transition in the early Universe. For Higgs masses in the range 46 less than or = M sub H less than or = 150 GeV and top quark masses less than 200 GeV, regions of symmetric and asymmetric vacuum coexist to below the critical temperature, with thermal equilibrium between the two phases maintained by fluctuations of both phases. We propose that the transition to the asymmetric vacuum is completed by percolation of these subcritical fluctuations. Our results are relevant to scenarios of baryogenesis that invoke a weakly first-order phase transition at the electroweak scale.

Gleiser, Marcelo↗
NASA NTRS1992

Fluctuation-driven electroweak phase transition

We examine the dynamics of the electroweak phase transition in the early Universe. For Higgs masses in the range 46 less than or = M sub H less than or = 150 GeV and top quark masses less than 200 GeV, regions of symmetric and asymmetric vacuum coexist to below the critical temperature, with thermal equilibrium between the two phases maintained by fluctuations of both phases. We propose that the transition to the asymmetric vacuum is completed by percolation of these subcritical fluctuations. Our results are relevant to scenarios of baryogenesis that invoke a weakly first-order phase transition at the electroweak scale.

Gleiser, Marcelo↗
DOE OSTI2022

Coupled continuum and network model framework to study catalyst layers of polymer electrolyte fuel cells

The nanostructured thin film (NSTF) catalyst layers which have demonstrated high power densities, mass activities, and exceptional metal and support stability can have limited operational robustness due to their thin thickness and the hydrophilicity of the metal-coated nano whiskers. The dispersed nanostructured thin film (dNSTF) catalyst layers have been developed by dispersing the NSTF Pt whiskers with ionomer and carbon support to increase the thickness and hydrophobicity. Continuum and network models (NM) are coupled through boundary conditions to study the polymer electrolyte fuel cell with a dNSTF cathode catalyst layer. The coupled model combines the computational efficiency of the continuum model with the pore-scale information in the dNSTF cathode catalyst layer of the NM. It captures the special morphology of the partially ionomer/water covered cylindrical whiskers, as well as water percolation through the pore structures and their impact on the cell performance. Here we observe optimal ionomer coverage on whiskers to be 0.5, ionomer to carbon ratio to be 0.9 and higher whisker to carbon ratios to be desired.

08 HYDROGEN↗
NASA NTRS2000

HTMT-class Latency Tolerant Parallel Architecture for Petaflops Scale Computation

Computational Aero Sciences and other numeric intensive computation disciplines demand computing throughputs substantially greater than the Teraflops scale systems only now becoming available. The related fields of fluids, structures, thermal, combustion, and dynamic controls are among the interdisciplinary areas that in combination with sufficient resolution and advanced adaptive techniques may force performance requirements towards Petaflops. This will be especially true for compute intensive models such as Navier-Stokes are or when such system models are only part of a larger design optimization computation involving many design points. Yet recent experience with conventional MPP configurations comprising commodity processing and memory components has shown that larger scale frequently results in higher programming difficulty and lower system efficiency. While important advances in system software and algorithms techniques have had some impact on efficiency and programmability for certain classes of problems, in general it is unlikely that software alone will resolve the challenges to higher scalability. As in the past, future generations of high-end computers may require a combination of hardware architecture and system software advances to enable efficient operation at a Petaflops level. The NASA led HTMT project has engaged the talents of a broad interdisciplinary team to develop a new strategy in high-end system architecture to deliver petaflops scale computing in the 2004/5 timeframe. The Hybrid-Technology, MultiThreaded parallel computer architecture incorporates several advanced technologies in combination with an innovative dynamic adaptive scheduling mechanism to provide unprecedented performance and efficiency within practical constraints of cost, complexity, and power consumption. The emerging superconductor Rapid Single Flux Quantum electronics can operate at 100 GHz (the record is 770 GHz) and one percent of the power required by convention semiconductor logic. Wave Division Multiplexing optical communications can approach a peak per fiber bandwidth of 1 Tbps and the new Data Vortex network topology employing this technology can connect tens of thousands of ports providing a bi-section bandwidth on the order of a Petabyte per second with latencies well below 100 nanoseconds, even under heavy loads. Processor-in-Memory (PIM) technology combines logic and memory on the same chip exposing the internal bandwidth of the memory row buffers at low latency. And holographic storage photorefractive storage technologies provide high-density memory with access a thousand times faster than conventional disk technologies. Together these technologies enable a new class of shared memory system architecture with a peak performance in the range of a Petaflops but size and power requirements comparable to today's largest Teraflops scale systems. To achieve high-sustained performance, HTMT combines an advanced multithreading processor architecture with a memory-driven coarse-grained latency management strategy called "percolation", yielding high efficiency while reducing the much of the parallel programming burden. This paper will present the basic system architecture characteristics made possible through this series of advanced technologies and then give a detailed description of the new percolation approach to runtime latency management.

Sterling, Thomas↗