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↗

The Fermi paradox: An approach based on percolation theory

If even a very small fraction of the hundred billion stars in the galaxy are home to technological civilizations which colonize over interstellar distances, the entire galaxy could be completely colonized in a few million years. The absence of such extraterrestrial civilizations visiting Earth is the Fermi paradox. A model for interstellar colonization is proposed using the assumption that there is a maximum distance over which direct interstellar colonization is feasible. Due to the time lag involved in interstellar communications, it is assumed that an interstellar colony will rapidly develop a culture independent of the civilization that originally settled it. Any given colony will have a probability P of developing a colonizing civilization, and a probability (1-P) that it will develop a non-colonizing civilization. These assumptions lead to the colonization of the galaxy occuring as a percolation problem. In a percolation problem, there will be a critical value of percolation probability, P(sub c). For P less than P(sub c), colonization will always terminate after a finite number of colonies. Growth will occur in 'clusters', with the outside of each cluster consisting of non-colonizing civilizations. For P greater than P(sub c), small uncolonized voids will exist, bounded by non-colonizing civilizations. For P approximately = to P(sub c), arbitrarily large filled regions exist, and also arbitrarily large empty regions.

Landis, Geoffrey A.↗

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↗

Multiphase flow in porous media

A development history and current status evaluation are presented for the theory of permeability and percolation. The microscale phenomena treated in this field have proven difficult to analyze due both to their tortuous geometry and the influence of capilarity. Capilary effects may be not only important but predominant, and are differentiated into those at the fluid-fluid interface, and those involving the existence of a contact line between the solid substrate and this interface. Percolation theory has been borrowed from physics and adapted to the two-phase engineering context.

Adler, Pierre M.↗

Understanding Thermal Transport in Polymer –Silver Nanowire Composites

Understanding thermal transport across filler-polymer interfaces and filler-filler contacts within polymer composites is of great importance for better thermal design of the composites that are widely used in high-performance heat exchangers, energy storage devices, and flexible electronics. Over the past several decades, considerable progress has been made in improving the thermal conductivity of polymer composites, but several key questions concerning the influence of interfacial thermal resistance, or Kapitza resistance, still remain. Firstly, the thermal properties of these composites are highly dependent on thermal transport through the filler network and its contacts. For metallic nanofillers, the thermal conductivity is often estimated using the Wiedemann-Franz law based on electrical conductivity; however, it remains a question whether the Wiedemann-Franz law still holds at nanoscale contacts. Through investigation of silver nanowires of varying sizes, we were able to demonstrate that the Lorenz number for silver nanowire increases with decreasing nanowire diameter. Examination of the corresponding electrical and thermal conductivities indicate that these changes are due to that the relative contribution of phonons becomes more significant as a result of elastic stiffening. Furthermore, we show that for silver nanowires, the contact thermal resistance is ~8 times lower than that of multi-walled carbon nanotube (MWCNT) of similar diameters. Additionally, through systematic studies of electrospun polymer-silver nanowire composite nanofibers, we investigated the impact of interface morphology on the thermal conductivity enhancement of the composite system and probed the value of the Kapitza resistance for individual polymer-filler interfaces. For polymer nanofibers containing continuous, single silver nanowires, the thermal conductivity increases linearly with increasing volume fraction of silver, which is consistent with the prediction of percolation theory for samples above the percolation limit. By comparing this linear trend to the measured thermal conductivities of composite nanofibers with more complex structures, we were able to determine the resistance associated with any additional polymer-filler boundaries. In doing so we find that the thermal boundary resistance for polyvinylpyrrolidone (PVP)–silver interfaces is significantly lower than that of comparable polymer-MWCNT composite systems.Together our studies provide new insights into thermal transport in polymer nanocomposites and should help facilitate the design of high performance polymeric thermal interface material.

Thermal transport↗

Understanding Thermal Transport in Polymer - Silver Nanowire Composites

Understanding thermal transport across filler-polymer interfaces and filler-filler contacts within polymer composites is of great importance for better thermal design of the composites that are widely used in high-performance heat exchangers, energy storage devices, and flexible electronics. Over the past several decades, considerable progress has been made in improving the thermal conductivity of polymer composites, but several key questions concerning the influence of interfacial thermal resistance, or Kapitza resistance, still remain. Firstly, the thermal properties of these composites are highly dependent on thermal transport through the filler network and its contacts. For metallic nanofillers, the thermal conductivity is often estimated using the Wiedemann-Franz law based on electrical conductivity; however, it remains a question whether the Wiedemann-Franz law still holds at nanoscale contacts. Through investigation of silver nanowires of varying sizes, we were able to demonstrate that the Lorenz number for silver nanowire increases with decreasing nanowire diameter. Examination of the corresponding electrical and thermal conductivities indicate that these changes are due to that the relative contribution of phonons becomes more significant as a result of elastic stiffening. Furthermore, we show that for silver nanowires, the contact thermal resistance is ~8 times lower than that of multi-walled carbon nanotube (MWCNT) of similar diameters. Additionally, through systematic studies of electrospun polymer-silver nanowire composite nanofibers, we investigated the impact of interface morphology on the thermal conductivity enhancement of the composite system and probed the value of the Kapitza resistance for individual polymer-filler interfaces. For polymer nanofibers containing continuous, single silver nanowires, the thermal conductivity increases linearly with increasing volume fraction of silver, which is consistent with the prediction of percolation theory for samples above the percolation limit. By comparing this linear trend to the measured thermal conductivities of composite nanofibers with more complex structures, we were able to determine the resistance associated with any additional polymer-filler boundaries. In doing so we find that the thermal boundary resistance for polyvinylpyrrolidone (PVP)–silver interfaces is significantly lower than that of comparable polymer-MWCNT composite systems. Together our studies provide new insights into thermal transport in polymer nanocomposites and should help facilitate the design of high performance polymeric thermal interface material.

Matthew L. Fitzgerald↗

Composite-Based High Performance Electroactive Polymers For Remotely Controlled Mechanical Manipulations in NASA Applications

This program supported investigation of an all-polymer percolative composite which exhibits very high dielectric constant (less than 7,000). The experimental results show that the dielectric behavior of this new class of percolative composites follows the prediction of the percolation theory and the analysis of the conductive percolation phenomena. The very high dielectric constant of the all-polymer composites which are also very flexible and possess elastic modulus not very much different from that of the insulation polymer matrix makes it possible to induce a high electromechanical response under a much reduced electric field (a strain of 2.65% with an elastic energy density of 0.18 J/cu cm can be achieved under a field of 16 MV/m). Data analysis also suggests that in these composites, the non-uniform local field distribution as well as interface effects can significantly enhance the strain responses. Furthermore, the experimental data as well as the data analysis indicate that the conduction loss in these composites will not affect the strain hysteresis.

Zhang, Q. M.↗

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↗