Search NASA⌕ Search

SEARCH · Search NASA

Results for “Grain boundary fracture”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Insights into elevated temperature tensile deformation mechanisms and kink banding in additively manufactured tungsten

Additive manufacturing (AM) potentially enables fabrication and repair of plasma facing components (PFCs) of tungsten. However, deformation mechanisms in AM-fabricated tungsten under service-relevant thermomechanical conditions remain unknown majorly due to challenges in achieving crack-free W. This study reveals the deformation mechanisms operative under tension at 800 °C and 1200 °C in the electron beam melting powder bed fusion (EBM-PBF) fabricated W. Textured columnar grains with mixed <001 >/<111 >|| build-direction were observed. Strong anisotropy in tensile properties persisted across test temperatures, and all specimens exhibited extensive deformation-induced banding phenomena. Grain boundary character analysis of the deformed specimens indicated that these deformation bands were kink bands with tilt grain boundaries prominently present at the band-matrix interface. Interestingly, the material within the kink bands rotated toward higher resolved shear stress, providing insights into the origins of kink band formation in body-centered cubic refractory metals. Intragranular misorientation axis analysis revealed that plasticity was dominated by {110} <111 > slip at 800 °C, whereas additional (213)[11⁢̄1] and (112)[11⁢̄1] slip systems activated at 1200 °C. A dense low-angle boundary networks observed near fractured surface indicated the onset of dynamic recrystallization. Results reveal plasticity governing mechanisms in AM-fabricated W under power plant-relevant conditions, and provide insights into kink band attributes in refractory metals.

Mayes, Riley [ORNL] (ORCID:000900089307010X)↗
DOE OSTI2024

Predicted Fracture Tendency of Naturally Occurring Aluminum Surface Coatings under Tensile Loading

Naturally occurring coatings on aluminum metal, such as its oxide or hydroxide, serve to protect the material from corrosion. Understanding the conditions under which these coatings mechanically fail is therefore expected to be an important aspect of predictive models for aluminum component lifetimes. To this end, we develop and apply a molecular dynamics (MD) modeling framework for conducting tension tests that is capable of isolating factors governing the mechanical strength as a function of coating chemistry, defect morphology, and variables associated with the loading path. We consider two representative materials, including γ-Al 2 O 3 and γ-Al(OH) 3 (i.e., oxide and hydroxide), both of which form readily as aluminum surface coatings. Our results indicate that defects have a significant bearing on the strength of aluminum oxide, with grain boundaries serving to reduce the strain at failure from ε zz = 0.300 to 0.219, relative to perfect single crystal. Our simulations also predict that porosity lowers the elastic stiffness and yield strength of the oxide. Relative to perfect crystal, we find porosity factors of 5%, 10% and 20% decrease the yield stress by 26%, 36% and 53%, respectively. MD predicts that perfect hydroxide and oxide single crystal have respective strains at failure of 0.08 and 0.31 under tensile uniaxial strain loading, and that the corresponding yield stresses are respectively 1.6 and 11.1 GPa. These data indicate that the hydroxide is substantially more susceptible to mechanical failure than the oxide. Our results, coupled with literature findings that indicate hot and humid conditions favor formation of hydroxide and defective oxide coatings, indicate the potential for a complicated dependence of aluminum corrosion susceptibility and stress corrosion cracking on aging history.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Microstructural characterization of cold-worked 316 stainless steel flux thimble tubes irradiated up to 100 dpa in a commercial Pressurized Water Reactor

Two flux thimble tubes (FTT) made of 15% cold-worked 316 stainless steel (SS) were harvested from Ringhals Pressurized Water Reactor (PWR) Unit 2, with peak damages of 76 and 100 displacements per atom (dpa) after 29 and 34 years’ service, respectively. Specimens sectioned from parent tubes were comprehensively characterized with nominal damage levels of ~0, ~41, ~74, 76, and 100 dpa at a nominal temperature range of 285-323 °C. Both FTTs contained helium and hydrogen gases as transmutation products. The helium follows a production rate of ~9.8 appm/dpa, while environmental factors complicate hydrogen production obscuring an exact H/dpa ratio. Irradiation-induced dislocation loops, nano-cavities, solute clusters, and microsegregation were all observed. The dislocation loops and nano-cavities indicated saturation at 41 dpa. The solute clusters continued to evolve with Ni-Si clusters formed at 41 dpa, and Ni-Si-Mn-P clusters formed at 74 and 100 dpa, but neither clusters exhibited distinct diffraction patterns at any damage levels. Solute clusters were observed to frequently be co-located with dislocation loops, but fully decorated loops were rarely detected. Significant radiation-induced segregation (RIS) was observed around grain boundaries at all damage levels. The modified inverse Kirkendall (MIK) model captured the RIS behavior of major elements. Large cavities within or around an Mn-S rich region were observed for the first time. Through all the damage levels, void swelling is always below 0.05%, making significant dimensional change unlikely in core internals when used at similar conditions. Meanwhile, the role of overwhelming nanocavities, presumably helium bubbles, should be considered in other potential degradation mechanisms, including irradiation-assisted stress corrosion cracking, embrittlement, and loss of fracture toughness, which remain the concerns for extended operation of nuclear power plants.

Post Irradiation Examination, Isotope Dilution Mas↗
DOE OSTI2026

Effects of temperature and dose rate on ion-irradiated γ-LiAlO 2 pellets

Defect accumulation and microstructural evolution during ion irradiation at elevated temperatures are governed by competing processes of defect production, driven by the dose rate, and defect recovery, controlled by diffusion, interaction, and annihilation. Here, this study investigates the effects of irradiation temperature and the dose rate on microstructural evolution, deuterium retention, and lithium volatilization in γ-LiAlO 2 pellets subjected to sequential He + and D + ion irradiation. Experiments were performed to a total fluence of 3 × 10 17 (He + + D + )/cm 2 at 623, 673, 723, and 773 K with an average He + dose rate of 7.7 × 10 −4 dpa/s, and to 2 × 10 17 (He + + D + )/cm 2 at 773 K with dose rates of 6.8 × 10 −5 , 2.9 × 10 −4 , and 7.3 × 10 −4 dpa/s. At 623 K, the microstructure was dominated by cavities and fractures with no observable precipitate formation, while small precipitates emerged at 673 K. Increasing the irradiation temperature to 723–773 K promoted the formation of larger, faceted LiAl 5 O 8 precipitates, and surface amorphization, accompanied by pronounced lithium depletion and H–D isotopic exchange. At 773 K, medium and high dose rates produced an amorphized surface layer over a crystalline subsurface containing LiAl 5 O 8 precipitates and blisters at the crystalline–amorphous interface, whereas low-dose-rate irradiation preserved surface crystallinity with cavities distributed in the matrix, around precipitates, and along grain boundaries. Precipitate morphology was anisotropic with limited size dependence on the dose rate. These results elucidate the coupled effects of temperature and the dose rate and demonstrate that sequential He + and D 2 + irradiation at 773 K reproduces key microstructural features and H isotope behavior observed in neutron-irradiated γ-LiAlO 2 at 573 K.

dose rate effects↗
DOE OSTI2025

Elasto-viscoplastic fast Fourier transform modeling framework for assessing microstructural effects on stress intensity factors characterizing fracture toughness

A large-strain elasto-viscoplastic fast Fourier transform (LS-EVPFFT) model with non-periodic (NP) velocity-based boundary conditions is adapted to simulate the sensitivity of stress intensity factors on microstructure for 304L stainless steel. The material was characterized via electron backscattered diffraction (EBSD) serial-sectioning to obtain a measured 3-D microstructural cell to perform simulations. The NP-LS-EVPFFT model, including the simulation setup and boundary conditions, was verified using a crystal plasticity finite element (CPFE) model. To this end, the generation of meshes of notched specimens was developed, which involved creating Python scripts for mesh “cutting” in Abaqus, and Sculpt scripts in Cubit for meshing of the measured microstructural cell processed with DREAM.3D. The complexity of the mesh preparation highlighted the advantages of the FFT-based model, which circumvents the mesh generation process. Given the efficiency of the FFT-based model, statistical distribution of stress intensity factors in function of crystal orientation at the crack tip, grain structure, and crystallographic texture surrounding the crack tip were predicted. Further, the distributions reveal about 10% variation of stress intensity factors with microstructure with the most significant sensitivity found to be the crystal orientation at the crack tip. The methodology developed in this work is discussed as a practical simulation tool for predicting the sensitivity of stress intensity factors on microstructural variability in metallic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Gradient nanostructured steel with superior tensile plasticity

Nanostructured metallic materials with abundant high-angle grain boundaries exhibit high strength and good radiation resistance. While the nanoscale grains induce high strength, they also degrade tensile ductility. We show that a gradient nanostructured ferritic steel exhibits simultaneous improvement in yield strength by 36% and uniform elongation by 50% compared to the homogenously structured counterpart. In situ tension studies coupled with electron backscattered diffraction analyses reveal intricate coordinated deformation mechanisms in the gradient structures. The outermost nanolaminate grains sustain a substantial plastic strain via a profound deformation mechanism involving prominent grain reorientation. This synergistic plastic co-deformation process alters the rupture mode in the post-necking regime, thus delaying the onset of fracture. The present discovery highlights the intrinsic plasticity of nanolaminate grains and their significance in simultaneous improvement of strength and tensile ductility of structural metallic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Updimensioning strategy derived from synthetic equiaxed grain structures for approximating 3D grain size distributions from 2D visualizations with 1D parameters

We generated synthetic equiaxed grain structures using computer graphics software to explore the relationship between various grain size determination methods and true three-dimensional (3D) grain diameters. Mirroring grain measurement techniques, the synthetic 3D grain structures are imaged as 2D micrographs which are measured to yield 1D grain size parameters. Synthetic grain structures provide data at a mass scale and permit exploration of both polished and fractured surface micrographs, revealing one-to-one correspondence between exposed 2D grain cross-sections and individual 3D grains. Analysis of this correspondence yielded a procedure to approximate 3D equiaxed grain size and volume distributions based on the mode of the 2D fractograph grain size distribution. The 3D approximation procedure is shown to be less susceptible to different imaging conditions that affect small, undiscernible grains compared to the standard planimetric and linear intercept methods, which by design also tend to underestimate the 3D grain diameter. The procedure requires larger sample sizes to lower variance and a deeper analysis which could become more practical with machine learning (ML) models for grain boundary segmentation, which synthetic grain structures can help train. This work lays the foundation for analyzing other grain distributions such as columnar and composite grains in similar depth.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

A Quantum Mechanical Study of Quartz (101) Interfacial Boundaries

Density functional theory (DFT) calculations were performed on periodic and molecular models to explore the energetics of bonding modes at quartz–quartz interfaces, common grain contacts in sandstones. Four interface types were modeled: H 2 O-mediated, silanol H-bond mediated, O–O peroxide bonds, and siloxane bonds. Each type of interface may exist in the subsurface of the Earth, depending upon water activity, temperature, and pressure, as interfacial structures are stable within DFT molecular dynamics simulations up to 473 K. The results predict interfacial energies SiOSi > SiOOSi > SiOH-6H 2 O > 2SiOH, as expected based on similar bond strengths in other systems. The proposed SiOOSi interface is not generally accounted for in chemomechanical models but could be important in sandstones at depth in the Earth. The implications for intergrain fracturing and H2 generation are discussed.

Hydrogen↗
DOE OSTI2023

Non‐Linear Kinetics of The Lithium Metal Anode on Li 6 PS 5 Cl at High Current Density: Dendrite Growth and the Role of Lithium Microstructure on Creep

Abstract Interfacial instability, viz., pore formation in the lithium metal anode (LMA) during discharge leading to high impedance, current focusing induced solid–electrolyte (SE) fracture during charging, and formation/behaviour of the solid–electrolyte interphase (SEI), at the anode, is one of the major hurdles in the development of solid‐state batteries (SSBs). Also, understanding cell polarization behaviour at high current density is critical to achieving the goal of fast‐charging battery and electric vehicle. Herein, via in situ electrochemical scanning electron microscopy (SEM) measurements, performed with freshly deposited lithium microelectrodes on transgranularly fractured fresh Li6PS5Cl (LPSCl), the LiǀLPSCl interface kinetics are investigated beyond the linear regime. Even at relatively small overvoltages of a few mV, the LiǀLPSCl interface shows non‐linear kinetics. The interface kinetics possibly involve multiple rate‐limiting processes, i.e., ion transport across the SEI and SE|SEI interfaces, as well as charge transfer across the LiǀSEI interface. The total polarization resistance R P of the microelectrode interface is determined to be ≈ 0.8 Ω cm 2 . It is further shown that the nanocrystalline lithium microstructure can lead to a stable LiǀSE interface via Coble creep along with uniform stripping. Also, spatially resolved lithium deposition, i.e., at grain surface flaws, grain boundaries, and flaw‐free surfaces, indicates exceptionally high mechanical endurance of flaw‐free surfaces toward cathodic load (>150 mA cm −2 ). This highlights the prominent role of surface defects in dendrite growth.

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

Friction-based riveting technique for AZ31 magnesium alloy

Rivet materials used in the aerospace and automotive industries are typically alloys of aluminum (Al) or steel. Magnesium (Mg) is an attractive alternative because it is 33% and 75% lighter than Al and steel, respectively. However, a primary barrier for using Mg as rivet material is its poor formability at room temperature, resulting in brittle fracture of the rivet head during cold riveting. Hot riveting is possible but suffers from high cycle time due to the need for pre-heating. A new friction-based riveting technique, Rotating Hammer Riveting (RHR), is reported to fully form AZ31 Mg rivet heads in a mere 0.23 seconds cycle time. Heat and pressure generated through severe plastic deformation during the process was sufficient to form the Mg rivet head without the need for a separate pre-heating operation. Due to dynamic recrystallization, AZ31 Mg grains in the rivet head were significantly refined from about 15 to 2.6 µm during RHR, which favors the formability of Mg rivet by triggering grain boundary sliding and reducing plastic anisotropy of Mg. In addition, RHR joints showed a metallurgical bond between the rivet head and top sheet, which eliminates a significant pathway for corrosion.

36 MATERIALS SCIENCE↗
DOE OSTI2022

From force chains to nonclassical nonlinear dynamics in cemented granular materials

In this letter, we present evidence for a mechanism responsible for the nonclassical nonlinear dynamics observed in many cemented granular materials that are generally classified as mesoscopic nonlinear elastic materials. We demonstrate numerically that force chains are created within the complex grain-pore network of these materials when subjected to dynamic loading. The interface properties between grains along with the sharp and localized increase of the stress occurring at the grain-grain contacts leads to a reversible decrease of the elastic properties at macroscopic scale and peculiar effects on the propagation of elastic waves when grain boundary properties are appropriately considered. These effects are observed for relatively small amplitudes of the elastic waves, i.e., within tens of microstrain, and relatively large wavelengths, i.e., orders of magnitude larger than the material constituents. The mechanics are investigated numerically using the hybrid finite-discrete-element method and match those observed experimentally using nonlinear resonant ultrasound spectroscopy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Optimization of the microstructure and mechanical properties of electron beam welded high-strength medium-entropy alloy (NiCoCr) 94 Al 3 Ti 3

For this work, the weldability and effects of subsequent heat treatment on the microstructure and mechanical properties of the medium-entropy alloy (MEA) (NiCoCr) 94 Al 3 Ti 3 have been explored using electron beam welding (EBW). It was found that EBW joints consist of a supersaturated f.c.c. matrix with coarse, asymmetric columnar grains (17.7–24 μm). The fraction of twins and Σ3 n CSL boundaries in weld joints was only 10% of that in the base metal, a value which increased to 30% after aging at 800 °C. The aging also led to improvements in the hardness, tensile strength, and work hardening rate, a result of the precipitation of a high volume fraction (0.29–0.33) of ultrafine (15–17.4 nm dia.) coherent spherical L1 2 -structured particles. The contributions of the L1 2 nanoparticles and Hall-Petch strengthening to the yield strength were estimated to be 537–573 MPa and 105–123 MPa, respectively. The MEA shows excellent weldability and, after aging, good strength (YS~1098 MPa, UTS~1368 MPa) and elongation to failure (~11%) by ductile fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A strong fracture-resistant high-entropy alloy with nano-bridged honeycomb microstructure intrinsically toughened by 3D-printing

Strengthening materials via conventional “top-down” processes generally involves restricting dislocation movement by precipitation or grain refinement, which invariably restricts the movement of dislocations away from, or towards, a crack tip, thereby severely compromising their fracture resistance. In the present study, a high-entropy alloy Al0.5CrCoFeNi is produced by the laser powder-bed fusion process, a “bottom-up” additive manufacturing process similar to how nature builds structures, with the microstructure resembling a nano-bridged honeycomb structure consisting of a face-centered cubic (fcc) matrix and an interwoven hexagonal net of an ordered body-centered cubic B2 phase. While the B2 phase, combined with high-dislocation density and solid-solution strengthening, provides strength to the material, the nano-bridges of dislocations connecting the fcc cells, i.e., the channels between the B2 phase on the cell boundaries, provide highways for dislocation movement away from the crack tip. Consequently, the nature-inspired microstructure imparts the material with an excellent combination of strength and toughness.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Analysis of time-dependent strain heterogeneity in shales using X-ray microscopy and digital volume correlation

In this paper, we report the results of a 4-dimensional (3 spatial dimensions plus time) in-situ creep experiment on a milli-core shale plug from Barnett formation in Texas. An outstanding challenge in shale mechanics is connecting observable characteristics to predictable mechanical and hydraulic behavior. We employ novel imaging and analysis techniques to statistically link sample mineralogy and structure to observed deformation. In particular, we focus on mechanisms of porosity reduction during transient response to step changes in stress, analogous to the aseismic deformation that occurs after an episode of fracturing. Our main goal in this study is to take the first steps to visualize and capture the small-scale deformation mechanisms involved in creep of shales, using imaging techniques. To do so, we conducted a uniaxial creep experiment inside an Xradia Versa 520 micro-CT system for 12 h at room temperature under a constant uniaxial stress of 25 MPa. We employed a Digital Volume Correlation (DVC) method to measure micro-scale deforspmation of the sample by comparing the resulting volumetric images with a high-resolution scan of the sample obtained prior to the creep experiment. The results indicate that the time-dependent deformation in the compliant clay plus kerogen matrix constituents facilitates rigid grain rotation and compaction of intergranular porosity at the boundary between rigid grains and matrix.

58 GEOSCIENCES↗
DOE OSTI2024

Microstructure and mechanical properties of friction stir weld performed on neutron-irradiated 304L steel with helium

As nuclear power plants (NPPs) approach or exceed 40–60 years of service, it may become necessary to repair damaged neutron-irradiated components to prolong longevity. However, fusion welding repair of irradiated steels and metallic materials is challenging because of helium-induced cracking. This study used friction stir welding (FSW) to address helium-related issues, such as helium bubble formation and grain boundary cracking. The microstructure, helium-induced degradation, and mechanical properties of a friction stir weld produced on neutron-irradiated 304L stainless steel with approximately 5.2 appm of helium were characterized. The analysis focused on variations in grain size, texture, and the morphology of helium-induced damage in the friction stir weld. Mechanical properties were characterized on the irradiated base metal (BM) and metallurgical zones of the friction stir weld: the stir zone (SZ), the thermo-mechanically affected zone (TMAZ), and the heat affected zone (HAZ). Only minor scattered porosity in the SZ and TMAZ and a few short microcracks (below 20 μm in length) in the TMAZ were observed, indicating limited helium-induced degradation. Tensile tests revealed good mechanical properties and fractography analysis demonstrated predominantly ductile fracture. In conclusion, the results highlight the immediate and substantial benefits of the FSW approach for repairing or joining helium-containing irradiated materials in NPPs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

LDRD 22A1059-068FP Tailoring the Properties of Multi-Phase Materials Through the Use of Correlative Microscopy and Machine Learning - Poster

High strength alloys with good ductility, hardness, and toughness are needed to meet stringent design requirements for extreme environments. One complication in this pursuit is the evidence that metals rarely exhibit both high strength and good fracture toughness as the underlying mechanisms work in opposition. An exception to this behavior is found in multiphase alloys that form complex microstructures of mixed phases with variable grain sizes and shapes that provide increased fracture toughness by the arrangement of their constituent elements. We propose to explore this phenomenon using state-of-the-art machine learning (ML) techniques in a new and novel manner to identify and correlate the critical microstructural features in a Titanium-10Vanadium-2Iron-3Aluminum (Ti-10V-2Fe-3Al) alloy that is reported to exhibit high strength and fracture toughness. Additionally, we will employ multiple, complementary characterization techniques such as optical microscopy, electron backscatter diffraction (EBSD), energy dispersive spectroscopy (EDS) and scanning electron microscopy to provide multi-layer, quantitative ground truth measures of the microstructures. This data will be used to train a Convolutional Neural Network (CNN) in a semi-supervised environment to identify key microstructural features such as ? platelet dimensions and locations and ?/? phase boundaries and correlate those features with the strength and toughness. Here the ? and ? nomenclature refers to hexagonal close pack (hcp) and body center cubic (bcc) crystal structures, respectively. Previous work has focused on popular alloys and typically used one characterization technique. This research is focused on a promising titanium alloy, uses multiple complimentary characterization tools to provide precise microstructural information and correlates to improved fracture toughness. The resulting ML tool can be trained for additional microstructural features, different alloy(s), and or target mechanical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Analysis of Deformation and Fracture Mechanisms in the Harvested Low-Dose Baffle Former Bolt via Advanced Mechanical Tests

This report details the production of tensile specimens from in-service irradiated baffle former bolts (BFBs) and presents the results of a pilot in-situ SEM/EBSD test conducted at Oak Ridge National Laboratory’s Low Activation Materials Development and Analysis (LAMDA) facility. The SEM/EDS results confirmed that the material is AISI 316L steel, not AISI 347, and that no Nb addition was detected in several analyzed locations. EBSD analysis of the microstructure revealed in-grain misorientation gradients indicative of cold work, likely caused by material processing, and no retained ferrite was identified. During straining, dislocation channel formation was the primary deformation mechanism, with twinning also observed in favorably oriented grains at local strains as small as ~1.5%. No strain-induced phase instability was observed, likely due to the fact that the Ni and Cr content reached the upper limits for the 316L steel specification. Multiple microfracture events were observed during the tensile test. Fine strain-induced pores formed at the channel-grain boundary intersection points, while coarser microcracks were associated with non-metallic inclusions. Data analysis is ongoing, and additional tests are planned in the near future to provide more statistics and details.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Inertia-induced mixing and reaction maximization in laminar porous media flows

Solute transport and biogeochemical reactions in porous and fractured media flows are controlled by mixing, as are subsurface engineering operations such as contaminant remediation, geothermal energy production, and carbon sequestration. Porous media flows are generally regarded as slow, so the effects of fluid inertia on mixing and reaction are typically ignored. Here, we demonstrate through microfluidic experiments and numerical simulations of mixing-induced reaction that inertial recirculating flows readily emerge in laminar porous media flows and dramatically alter mixing and reaction dynamics. An optimal Reynolds number that maximizes the reaction rate is observed for individual pore throats of different sizes. This reaction maximization is attributed to the effects of recirculation flows on reactant availability, mixing, and reaction completion, which depend on the topology of recirculation relative to the boundary of the reactants or mixing interface. Recirculation enhances mixing and reactant availability, but a further increase in flow velocity reduces the residence time in recirculation, leading to a decrease in reaction rate. The reaction maximization is also confirmed in a flow channel with grain inclusions and randomized porous media. Interestingly, the domain-wide reaction rate shows a dramatic increase with increasing Re in the randomized porous media case. This is because fluid inertia induces complex three-dimensional flows in randomized porous media, which significantly increases transverse spreading and mixing. This study shows how inertial flows control reaction dynamics at the pore scale and beyond, thus having major implications for a wide range of environmental systems.

Chen, Michael A. (ORCID:0000000294614710)↗