Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “percolation theory”

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 19 records

A percolation theory for designing corrosion-resistant alloys

Iron-chromium and nickel-chromium binary alloys containing sufficient quantities of chromium serve as the prototypical corrosion-resistant metals owing to the presence of a nanometre-thick protective passive oxide film. Should this film be compromised by a scratch or abrasive wear, it reforms with little accompanying metal dissolution, a key criterion for good passive behaviour. This is a principal reason that stainless steels and other chromium-containing alloys are used in critical applications ranging from biomedical implants to nuclear reactor components. Unravelling the compositional dependence of this electrochemical behaviour is a long-standing unanswered question in corrosion science. Herein, we develop a percolation theory of alloy passivation based on two-dimensional to three-dimensional crossover effects that accounts for selective dissolution and the quantity of metal dissolved during the initial stage of passive film formation. We validate this theory both experimentally and by kinetic Monte Carlo simulation. Our results reveal a path forward for the design of corrosion-resistant metallic alloys.

36 MATERIALS SCIENCE↗

Theoretical power-law relationship between permeability and formation factor

Permeability k and formation factor F are important petrophysical properties that have broad applications to reservoir engineering. Although well-log analysis has advanced reservoir characterizations, direct measurement of k is still not available. Accordingly, resistivity logs and more particularly formation factor have been frequently used to estimate k. In this study, we present a theoretical framework to relate k to F. We use power-law scaling from continuum percolation theory to develop a theoretic relationship between k and F, which is similar in form to a published power-law equation with empirical coefficients. Our result gives a physical interpretation to this empirical relationship. We also compare power-law scaling laws from percolation theory with data derived from experimental measurements and numerical simulations in the literature. Finally, using the Curve Fitting toolbox of MATLAB, the parameters of the power-law equations are determined, and the significance of the non-universal power laws are discussed.

42 ENGINEERING↗

SOC Microstructural Property Estimator

This pre-trained ML model is a tool that uses basic compositional parameters for porous solid oxide cell (SOC) electrodes - the phase fractions and mean particle/pore diameters – as inputs and uses them to estimate additional electrochemical performance parameters: active (i.e., connected) TPB density, all tortuosity factors, and phase pair specific interfacial areas. The electrode is assumed to be composed of two solid phases and a pore phase. The property calculations are performed using neural network regression models trained on a large bank of synthetic electrode microstructural data that NETL has generated using the program DREAM3D (that bank is also hosted on EDX: https://edx.netl.doe.gov/dataset/soc-synthetic-microstructure-bank). This means the generated parameters are based on training from actual measured properties from 3D microstructures, not estimated from geometric simplifications. This tool was developed and is intended to replace percolation theory calculations in models that use hypothetical electrode properties. An example use case would be running SOC performance simulations across a parametric sweep of electrode designs (e.g., varying phase fractions and particle sizes) and assessing how it impacts the electrochemical performance of the SOC. Within the parameter space of the training data (statistics of that parameter space is provided in the readme file), this model achieves sub-5% mean absolute percent errors, an order of magnitude less error than percolation theory across the same parameter space. However, be aware that this tool was developed with parametric simulations in mind, and users are encouraged to assess accuracy for their own specific use case rather than taking accuracy metrics at face value. More info, including a usage guide, is in the included readme file. This tool should be cited with the DOI number provided.

Electrode Microstructure↗

Microstructure and water absorption of ancient concrete from Pompeii: An integrated synchrotron microtomography and neutron radiography characterization

There is renewed interest in using advanced techniques to characterize ancient Roman concrete due to its exceptional durability and low-carbon footprint. In the present work, samples were drilled from the “Hospitium” in Pompeii and were analyzed by synchrotron microtomography (μCT) and neutron radiography to study how the microstructure, including the presence of induced cracks, affects their water adsorption. The water distribution and absorptivity were quantified by neutron radiography. The 3D crack propagation, pore size distribution and orientation, tortuosity, and connectivity were analyzed from μCT results using advanced imaging methods. Porosity was also measured by mercury intrusion porosimetry (MIP) as a reference. Ductile fracture patterns were observed once cracks were introduced. Compared to Portland cement mortar/concrete, the Pompeii samples had relatively high porosity, low connectivity, and a similar coefficient of capillary penetration. In addition, permeability was predicted from models based on percolation theory and pore structure data to evaluate the fluid transport properties. Understanding the microstructure of ancient Pompeii concrete is important because it could inspire the development of modern concrete with high durability.

36 MATERIALS SCIENCE↗

Deciphering the Conflict Between Ion and Electron Percolating Networks in Solid-State Battery Cathodes

High-energy- and power-density solid state batteries require an optimal cathode composition and microstructural arrangement of cathode active material, solid-state electrolyte, conductive carbon, and binder to simultaneously support lithium-ion transport, electron conduction, and storage capacity. The ion and electron conducting phases in solid-state cathodes counteract each other's percolating networks as their mass ratios increase or decrease relative to each other. Here, we investigate targeted mass ratio variations of argyrodite solid electrolyte and two different types of conductive carbon (particles and fibers) in composite LiNi0.8Mn0.1Co0.1O2 (NMC811) solid-state cathodes to ascertain the ionic-electronic tradeoffs in cathode performance. Through ionic and electronic conductivity measurements on composite cathodes, as well as rate-testing and cycling performance in full cells, it is shown that the conductive carbon fibers form a percolative electronic network within the composite at a lower mass ratio (3-5 wt%) than particulate carbon (>5 wt%). The threshold to achieve electronic percolation coincides with higher accessible capacity in the cathode as the active material particles become electronically connected. However, carbon loadings beyond this percolation threshold lead to increased ion transport resistance, arising from disruptions to ionic conduction pathways and degraded contact at the interface between the electrolyte and active materials. Imaging, spectroscopy, and physics-based models quantitatively describe the relationship between carbon and electrolyte compositions and the cell's capacity and rate performance through percolation theory. This work demonstrates the importance of quantitatively understanding percolating networks in solid-state cells and that strategic engineering of conductive carbon morphologies can further increase the energy- and power-density of solid-state cells.

25 ENERGY STORAGE↗

Insulator-to-metal phase transition in a few-layered MoSe 2 field effect transistor

The metal-to-insulator phase transition (MIT) in low-dimensional materials and particularly two-dimensional layered semiconductors is exciting to explore due to the fact that it challenges the prediction that a two-dimensional system must be insulating at low temperatures. Thus, the exploration of MITs in 2D layered semiconductors expands the understanding of the underlying physics. Here we report the MIT of a few-layered MoSe 2 field effect transistor under a gate bias (electric field) applied perpendicular to the MoSe 2 layers. With low applied gate voltage, the conductivity as a function of temperature from 150 K to 4 K shows typical semiconducting to insulating character. Above a critical applied gate voltage, V c , the conductivity becomes metallic (i.e., the conductivity increases continuously as a function of decreasing temperature). Evidence of a metallic state was observed using an applied gate voltage or, equivalently, increasing the density of charge carriers within the 2D channel. We analyzed the nature of the phase transition using percolation theory, where conductivity scales with the density of charge carriers as σ ∝ (n - n c ) δ . The critical exponent for a percolative phase transition, δ(T), has values ranging from 1.34 (at T = 150 K) to 2 (T = 20 K), which is close to the theoretical value of 1.33 for percolation to occur. Thus we conclude that the MIT in few-layered MoSe 2 is driven by charge carrier percolation. Finally, the conductivity does not scale with temperature, which is a hallmark of a quantum critical phase transition.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

A Percolating Path to Green Iron

About 1.9 gigatonnes of steel is produced every year emitting 7% (2.7 gigatonnes) of global CO 2 in the process. More than 50% of the CO 2 emissions come from a single step of steelmaking, known as ironmaking. Hydrogen based direct reduction (HyDR) of iron oxide to iron has emerged as an emissions free ironmaking alternative. However multi-scale phenomena ranging from nanometers to meters inside HyDR reactors exhibit detrimental microstructure evolution which resists gaseous transport of H 2 /H 2 O, slows reaction rates and disrupts continuous reactor operation. To resolve the conundrum between atomic and reactor scales, we devise a percolation-theory model to reconcile nanoscale porosity with macroscopic properties relevant to reactor design models. Using synchrotron nano X-ray computed-tomography, we quantify the evolution of pores in iron oxide pellets, and demonstrate how nano-scale pore networks influence micro and macro-scale flow properties such as permeability, diffusivity and tortuosity. Our new modeling framework bridges the gap between scales and offers the criteria to accelerate HyDR by at least 5x via feedstock-reactor synergies based on percolation.

Paul, Subhechchha↗

A percolating path to green iron

About 1.9 gigatons of steel is produced every year, emitting 8% (3.6 gigatons) of global CO 2 in the process. More than 50% of the CO 2 emissions come from a single step of steel production, known as ironmaking. Hydrogen- based direct reduction (HyDR) of iron oxide to iron has emerged as an emission-free ironmaking alternative. However, multiple physical and chemical phenomena ranging from nanometers to meters inside HyDR reactors alter the microstructure and pore networks in iron oxide pellets, in ways that resist gaseous transport of H 2 /H 2 O, slow reaction rates, and disrupt continuous reactor operation. Using synchrotron nano X-ray computed tomography and percolation theory, we quantify the evolution of pores in iron oxide pellets and demonstrate how nanoscale pore connectivity influences micro- and macroscale flow properties such as permeability, diffusivity, and tortuosity. Our modeling framework connects disparate scales and offers opportunities to accelerate HyDR.

hydrogen↗

Impact of artificial topological changes on flow and transport through fractured media due to mesh resolution

Abstract We performed a set of numerical simulations to characterize the interplay of fracture network topology, upscaling, and mesh refinement on flow and transport properties in fractured porous media. We generated a set of generic three-dimensional discrete fracture networks at various densities, where the radii of the fractures were sampled from a truncated power-law distribution, and whose parameters were loosely based on field site characterizations. We also considered five network densities, which were defined using a dimensionless version of density based on percolation theory. Once the networks were generated, we upscaled them into a single continuum model using the upscaled discrete fracture matrix model presented by Sweeney et al. (2019). We considered steady, isothermal pressure-driven flow through each domain and then simulated conservative, decaying, and adsorbing tracers using a pulse injection into the domain. For each simulation, we calculated the effective permeability and solute breakthrough curves as quantities of interest to compare between network realizations. We found that selecting a mesh resolution such that the global topology of the upscaled mesh matches the fracture network is essential. If the upscaled mesh has a connected pathway of fracture (higher permeability) cells but the fracture network does not, then the estimates for effective permeability and solute breakthrough will be incorrect. False connections cannot be eliminated entirely, but they can be managed by choosing appropriate mesh resolution and refinement for a given network. Adopting octree meshing to obtain sufficient levels of refinement leads to fewer computational cells (up to a 90% reduction in overall cell count) when compared to using a uniform resolution grid and can result in a more accurate continuum representation of the true fracture network.

58 GEOSCIENCES↗

Magnetism in Mixed Valence, Defect, Cubic Perovskites: BaIn 1– x Fe x O 2.5+δ , x = 0.25, 0.50, and 0.75. Local and Average Structures

The series BaIn 1–x Fe x O 2.5+δ , x = 0.25, 0.50, and 0.75, has been prepared under air-fired and argon-fired conditions and studied using X-ray diffraction, d.c. and a.c. susceptibility, Mössbauer spectroscopy, neutron diffraction, X-ray near edge absorption spectroscopy (XANES), and X-ray pair distribution (PDF) methods. While Ba 2 In 2 O 5 (BaInO 2.5 ) crystallizes in an ordered brownmillerite structure, Ibm2, and Ba 2 Fe 2 O 5 (BaFeO 2.5 ) crystallizes in a complex monoclinic structure, P2 1 /c, showing seven Fe 3+ sites with tetrahedral, square planar, and octahedral environments, all phases studied here crystallize in the cubic perovskite structure, Pm$\bar{}3}$m, with long-range disorder on the small cation and oxygen sites. 57 Fe Mössbauer studies indicate a mixed valency, Fe 4+ /Fe 3+ , for both the air-fired and argon-fired samples. The increased Fe 3+ content for the argon-fired samples is reflected in increased cubic cell constants and in the increased Mössbauer fraction. It appears that the Pm$\bar{}3}$m phases are only metastable when fired in argon. From a slightly modified percolation theory for a primitive cubic lattice (taking into account the presence of random O atom vacancies), long-range spin order is permitted for the x = 0.50 and 0.75 phases. Instead, the d.c. susceptibility shows only zero-field-cooled (ZFC) and field-cooled (FC) divergences at ~6 K [5 K] for x = 0.50 and at ~22 K [21 K] for x = 0.75, with values for the argon-fired samples in [ ]. Neutron diffraction data for the air-fired samples confirm the absence of long-range magnetic order at any studied temperature. For the air-fired x = 0.50, a.c. susceptibility data show a frequency-dependent χ'(max) and spin glass behavior, while for x = 0.75, χ'(max) is invariant with frequency, ruling out either a spin glass or a superparamagnetic ground state. These behaviors are discussed in terms of competing Fe 3+ –Fe 3+ antiferromagnetic exchange and ferromagnetic Fe 3+ –Fe 4+ exchange. The PDF and 57 Fe Mössbauer data indicate a local structure at short interatomic distances, which deviates strongly from the average Pm$\bar{}3}$m model. Fe Mössbauer, PDF, and XANES data show a systematic dependence on x and indicate that the Fe 3+ sites are largely fourfold-coordinated and Fe 4+ sites are fivefold- or sixfold-coordinated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Designing polar textures with ultrafast neuromorphic features from atomistic simulations

This review summarizes recent works, all using a specific atomistic approach, that predict and explain the occurrence of key features for neuromorphic computing in three archetypical dipolar materials, when they are subject to THz excitations. The main ideas behind such atomistic approach are provided, and illustration of model relaxor ferroelectrics, antiferroelectrics, and normal ferroelectrics are given, highlighting the important potential of polar materials as candidates for neuromorphic computing. Some peculiar emphases are made in this Review, such as the connection between neuromorphic features and percolation theory, local minima in energy path, topological transitions and/or anharmonic oscillator model, depending on the material under investigation. By considering three different and main polar material families, this work provides a complete and innovative toolbox for designing polar-based neuromorphic systems.

THz pulses↗

Delocalization Transition of a Disordered Axion Insulator

The axion insulator is a higher-order topological insulator protected by inversion symmetry. We show that, under quenched disorder respecting inversion symmetry on average, the topology of the axion insulator stays robust, and an intermediate metallic phase in which states are delocalized is unavoidable at the transition from an axion insulator to a trivial insulator. We derive this conclusion from general arguments, from classical percolation theory, and from the numerical study of a 3D quantum network model simulating a disordered axion insulator through a layer construction. We find the localization length critical exponent near the delocalization transition to be ν = 1.42 ± 0.12. We further show that this delocalization transition is stable even to weak breaking of the average inversion symmetry, up to a critical strength. Finally, we also quantitatively map our quantum network model to an effective Hamiltonian and we find its low-energy k ∙ p expansion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

On the Source of Conductivity in Alkaline Zn Anodes: Zn Percolation and ZnO Bridging

Alkaline Zn anodes are fundamental to commercial Zn-MnO 2 batteries as well as emerging rechargeable designs. In these electrodes, Zn particles are both the active material and source of electronic conduction. However, there are known cases in which electronic connectivity between Zn particles or the current collecting pin is lost even though the battery continues to function. Here we use X-ray computed tomography (CT) of AA batteries to demonstrate several examples of Zn particle-to-particle connectivity breakdown, which is observed even in cases at relatively high discharge rate. This indicates maintenance of the electronic network through the less-conductive ZnO discharge product. We introduce a new equation for electronic conductivity maintained through bridges of ZnO formed between percolating Zn particles. This necessitates discarding the Bruggeman correlation and instead redefining effective electronic conductivity using percolation theory. We demonstrate that such a model for conduction enables prediction of an inverted reaction zone, which is an experimentally observed case in which significant Zn dissolution and ZnO formation occurs heavily near the current collecting pin. Current computational Zn-MnO 2 models never predict an inverted reaction zone, and thus the updated conductivity enables models to accurately explain a wider range of experimental conditions.

Guida, Dominick P. [Northeastern Univ., Boston, MA↗

A New Model for Simulating the Imbibition of a Wetting-Phase Fluid in a Matrix-Fracture Dual Connectivity System

The imbibition experiment is an effective approach for measuring petrophysical properties of porous media, with many such experiments performed over the past decade. Quite some empirical, analytical, and numerical models have been developed to simulate spontaneous imbibition of the wetting phase fluid into porous media, but limitations still exist. In previous studies, the imbibition process has been considered to give a piston-like displacement or the porous medium modeled as multiply-sized pores linked with bonds; both approaches fail to yield comprehensive results due to their neglect of the presence of irregular fractures or nonuniform flow paths through the matrix. By building a numerical model for simulating laboratory-scale experimental data, we performed imbibition tests on several fractured Barnett Shale samples having fractures either parallel ( P ) or transverse ( T ) to the bedding plane and used MATLAB to build a new numerical model by combining the imbibition process in fractures and the matrix using concepts from percolation theory. The experimental data show that the rocks with P -direction fractures have a more steady increase of imbibition rates than the case of T -direction one. As the shale matrix with low pore connectivity hampers the upward water movement, the imbibition rate of shales with T -direction fractures will decrease suddenly after the bottom layer in contact with water is saturated during the initial period. This wetting phase movement (WPM) model can simulate 3D porous media with 2D fractures. The rate of imbibition by fractured porous media is associated with physical parameters such as porosity and fracture distribution (e.g., the number and angle of fractures). Using Monte Carlo methods, we examined fracture parameters and predicted elapsed time and cumulative water imbibition, for the Barnett Shale samples. The results show that the rate of imbibed water mass is sensitive to the number of fractures directly connected to water source, and the connectivity between two neighboring grid cells is a key parameter for the wetting-front progression. The findings of this study can help to better understand the imbibition process with multiple influencing processes and factors in fractured-matrix rocks. Although the experiments, data simulation, and prediction results are based only on Barnett Shale samples, the model is readily applicable to imbibition tests of other fractured rocks to show the spatial and temporal behavior during a dynamic imbibition process that are not easily captured experimentally.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Tailoring Carbide Dispersed Steels: A Path to Increased Strength and Hydrogen Tolerance

The use of transition metal carbides is reported for use as a hydrogen trapping mechanism for ferritic and austenitic steel materials. The program combined computational modeling and simulations to guide experiments towards candidate metal carbide traps, both for interfacial and interior trapping. It was found that interfacial trapping is less effective than interior trapping, with the group IVB transition metal carbides being the most effect internal traps with a loss of carbon. The sub-stoichiometric rocksalt structure accommodate the hydrogen atoms in its octahedral interstices. Using percolation theory, carbon loss of approximately 25% or more was sufficient to ensure an interconnected network of vacancies for such trapping from the surface to the internal sites within the carbide. Using this as a guide, the program developed a means to provide a uniform dispersion of ZrC nanoparticles with either Fe or 304L micron-scale powders which was then consolidated by direct current sintering. Electrolytic hydrogen diffusivity studies confirmed the reduction of hydrogen diffusion in the matrix with increasing ZrC content, which was a linear response over the sample range studied (0.01 to 1.0 wt.%). The consolidated material was micro-tensile tested in either a non-hydrogen or hydrogen charge condition and compared to a control with no carbides. Additions up to 0.05 wt.% ZrC increased the yield strength with no loss in ductility in either the non-hydrogen or hydrogen tested condition. ZrC concentrations above this amount further increased the yield strength at the expense of ductility. While these samples had a lower absolute ductility value prior to failure, the relative change in ductility between the non-hydrogen and hydrogen charge states was less for the carbides than that of the control. Metal-rich ZrC nanoparticles were fabricated through a conformal coating process yielding ZrC0.66 particles that were then incorporated into a metal matrix. Notch fatigue testing in a hydrogen environment was conducted where the number of cycles to failure was found to be less in the control than that of the carbide addition. However, the spread in experimental data and the number of samples tested limits a conclusive outcome based on defects noticed in the gauge section of all the powder processed samples. The collective outcomes of this report provide further insight into the mechanisms by which carbides act as hydrogen traps; a means to process such carbides through powder metallurgy; and their associated mechanical performance in either a non-hydrogen or hydrogen-charged condition.

08 HYDROGEN↗

A generalized Flory-Stockmayer kinetic theory of connectivity percolation and rigidity percolation of cytoskeletal networks

Actin networks are essential for living cells to move, reproduce, and sense their environments. The dynamic and rheological behavior of actin networks is modulated by actin-binding proteins such as α-actinin, Arp2/3, and myosin. There is experimental evidence that actin-binding proteins modulate the cooperation of myosin motors by connecting the actin network. In this work, we present an analytical mean field model, using the Flory-Stockmayer theory of gelation, to understand how different actin-binding proteins change the connectivity of the actin filaments as the networks are formed. We follow the kinetics of the networks and estimate the concentrations of actin-binding proteins that are needed to reach connectivity percolation as well as to reach rigidity percolation. We find that Arp2/3 increases the actomyosin connectivity in the network in a non-monotonic way. We also describe how changing the connectivity of actomyosin networks modulates the ability of motors to exert forces, leading to three possible phases of the networks with distinctive dynamical characteristics: a sol phase, a gel phase, and an active phase. Thus, changes in the concentration and activity of actin-binding proteins in cells lead to a phase transition of the actin network, allowing the cells to perform active contraction and change their rheological properties.

59 BASIC BIOLOGICAL SCIENCES↗

Percolative phase transition in few-layered MoSe 2 field-effect transistors using Co and Cr contacts

The metal-to-insulator phase transition (MIT) in two-dimensional (2D) materials under the influence of a gating electric field has revealed interesting electronic behavior and the need for a deeper fundamental understanding of electron transport processes, while attracting much interest in the development of next-generation electronic and optoelectronic devices. Although the mechanism of the MIT in 2D semiconductors is a topic under debate in condensed matter physics, our work demonstrates the tunable percolative phase transition in few-layered MoSe 2 field-effect transistors (FETs) using different metallic contact materials. Here, we attempted to understand the MIT through temperature-dependent electronic transport measurements by tuning the carrier density in a MoSe 2 channel under the influence of an applied gate voltage. In particular, we have examined this phenomenon using the conventional chromium (Cr) and ferromagnetic cobalt (Co) as two metal contacts. For both Cr and Co, our devices demonstrated n-type behavior with a room-temperature field-effect mobility of 16 cm 2 V −1 s −1 for the device with Cr-contacts and 92 cm 2 V −1 s −1 for the device with Co-contacts, respectively. Further, with low temperature measurements at 50 K, the mobilities increased significantly to 65 cm 2 V −1 s −1 for the device with Cr and 394 cm 2 V −1 s −1 for the device with Co-contacts. By fitting our experimental data to the percolative phase transition theory, the temperature-dependent conductivity data show a transition from an insulating-to-metallic behavior at a bias of ∼28 V for Cr-contacts and ∼20 V for Co-contacts. This cross-over of the conductivity can be attributed to an increase in carrier density as a function of the gate bias in temperature-dependent transfer characteristics. By extracting the critical exponents, we find that the transport behavior in the device with Co-contacts aligns closely with the 2D percolation theory. In contrast, the devices with Cr-contacts deviate significantly from the 2D limit at low temperatures.

36 MATERIALS SCIENCE↗

Universality class for loopless invasion percolation models and a percolation avalanche burst model for hydraulic fracturing

Invasion percolation is a model that was originally proposed to describe growing networks of fractures. In this work, we describe a loopless algorithm on random lattices, coupled with an avalanche-based model for bursts. The model reproduces the characteristic b-value seismicity and spatial distribution of bursts consistent with earthquakes resulting from hydraulic fracturing (“fracking”). We test models for both site invasion percolation and bond invasion percolation. These have differences on the scale of site and bond lengths l. But since the networks are characterized by their large-scale behavior, l << L, we find small differences between scaling exponents. Though data may not differentiate between models, our results suggest that both models belong to different universality classes.

58 GEOSCIENCES↗