Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Grain Growth”

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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Defect Engineering of Bi 2 SeO 2 Thermoelectrics

Abstract Bi 2 SeO 2 is a promising n ‐type semiconductor to pair with p ‐type BiCuSeO in a thermoelectric (TE) device. The TE figure of merit zT and, therefore, the device efficiency must be optimized by tuning the carrier concentration. However, electron concentrations in self‐doped n ‐type Bi 2 SeO 2 span several orders of magnitude, even in samples with the same nominal compositions. Such unsystematic variations in the electron concentration have a thermodynamic origin related to the variations in native defect concentrations. In this study, first‐principles calculations are used to show that the selenium vacancy, which is the source of n ‐type conductivity in Bi 2 SeO 2 , varies by 1–2 orders of magnitude depending on the thermodynamic conditions. It is predicted that the electron concentration can be enhanced by synthesizing under more Se‐poor conditions and/or at higher solid‐state reaction temperatures ( T SSR ), which promote the formation of selenium vacancies without introducing extrinsic dopants. The computational predictions are validated through solid‐state synthesis of Bi 2 SeO 2 . More than two orders of magnitude increase are observed in the electron concentration simply by adjusting the synthesis conditions. Additionally, a significant effect of grain boundary scattering on the electron mobility in Bi 2 SeO 2 is revealed, which can also be controlled by adjusting T SSR . By simultaneously optimizing the electron concentration and mobility, a zT of ≈0.2 is achieved at 773 K for self‐doped n ‐type Bi 2 SeO 2 . The study highlights the need for careful control of thermodynamic growth conditions and demonstrates TE performance improvement by varying synthesis parameters according to thermodynamic guidelines.

Novitskii, Andrei [MANA, National Institute for Ma↗
NASA NTRS1992

Elevated-temperature fracture resistances of monolithic and composite ceramics using chevron-notched bend tests

The quasi-static fracture behaviors of monolithic ceramics (SiC, Si3N4, MgAl2O4), self-reinforced monoliths (acicular grained Si3N4, acicular grained mullite), and ceramic matrix composites (SiC whisker/Al2O3 matrix, TiB2 particulate/SiC matrix, SiC fiber/CVI SiC matrix, Al2O3 fiber/CVI SiC matrix) were measured over the temperature range of 20 to 1400 C. The chevron notched, bend bar test geometry was essential for characterizing the elevated temperature fracture resistances of this wide range of quasi-brittle materials during stable crack growth. Fractography revealed the differences in the fracture behavior of the different materials at the various temperatures. The fracture resistances of the self-reinforced monoliths were comparable to those of the composites and the fracture mechanisms were found to be similar at room temperature. However at elevated temperatures the differences of the fracture behavior became apparent where the superior fracture resistance of the self-reinforced monoliths were attributed to the minor amounts of glassy, intergranular phases which were often more abundant in the composites and affected the fracture behavior when softened by elevated temperatures.

Ghosh, Asish↗
NASA NTRS2005

Estimation of ZnSe Slow-Crack-Growth Properties for Design of the Flow Enclosure Accommodating Novel Investigations in Combustion of Solids (FEANICS) Windows

This report reviews some of the literature on the fracture strength, fracture toughness, and crack growth properties of chemical-vapor-deposited ZnSe. The literature was reviewed to determine if the existing data on ZnSe is adequate to design windows for the Flow Enclosure Accommodating Novel Investigations in Combustion of Solids (FEANICS) project. Unfortunately, most of the published reports do not give all of the necessary design parameters despite having measured the data to do so. Further, the original data is not available. The data tabulated herein was determined by digitizing plots in original reprints of the publications. Based on the published data, an estimate of the slow-crack-growth parameters for small cracks in 100 percent humidity was made. For 100 percent humidity, the slow-crack-growth parameters n and A for small crack (or single crystal) failure were estimated. Weibull moduli estimated from bending of beams and circular plates ranged from 4 to 9, while fracture strengths ranged from 29 MPa in water to 72 MPa in dry nitrogen. Fracture toughness measurements yielded ranges, with the lower values representing failure from small flaws within grains and the larger values representing macroscopic cracks. Much of the data analyzed exhibited significant scatter, and the standard deviations were very large.

Salem, Jonathan A.↗
DOE OSTI2020

Controlling Pt Crystal Defects on the Surface of Ni–Pt Core–Shell Nanoparticles for Active and Stable Electrocatalysts for Oxygen Reduction

A strategy of direct growth of Pt on Ni was used to create and control Pt crystal defects on the surface of Ni–Pt core–shell nanoparticles. The control over the types of defects was easily achieved by changing the surfactant system. Here, in this work, two types of crystal defects have been introduced into Ni–Pt core–shell nanoparticles: polycrystalline shells with multiple grain boundaries and step-edge shells with undercoordinated atoms at corners and steps. We show that the step-edge shell has a higher specific activity for the oxygen reduction reaction (ORR), while the thinner polycrystalline shell results in a higher activity per mass and stability. Our results suggest that Ni–Pt core–shell nanoparticles with a thin Pt shell that have high density of crystal defect should be targeted for high performance ORR catalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Bioenergy sorghum’s deep roots: A key to sustainable biomass production on annual cropland

Abstract Bioenergy sorghum has high biomass yield potential, drought resilience, good nitrogen use efficiency, and a root system that contributes to the accumulation of soil organic carbon. In this study, field grown bioenergy sorghum root systems were analyzed during the growing season to characterize their depth, biomass, morphology, anatomy, and gene expression profiles. Bioenergy sorghum roots grew continuously during a 155‐day growing season producing ~175 nodal roots, accumulating ~7 Mg of dry biomass per hectare, and reaching >2 m deep in the soil profile. Nodal roots within 20 cm of the stem were 1–5 mm in diameter, whereas roots deeper in soil profiles were enriched in lateral roots with small diameters (~30–500 µm) enabling growth through soil macropores. In older field‐grown plants, roots with intact endodermal, vascular and inner root tissues were surrounded by degraded or aerenchyma‐filled epidermal and cortical cell layers. Transcriptome analysis of nodal, surface, and deep roots identified >2,500 differentially expressed genes involved in root growth, transport, adaptation, defense, and AMF–root interaction. Deep roots (180–240 cm) differentially expressed genes that regulate lateral root growth. Surface roots (0–20 cm) located mid‐row differentially expressed genes involved in nitrate transport, whereas ammonium transport genes were expressed in surface and deep roots and genes involved in phosphate transport were expressed in nodal, surface, and deep roots. Overall, bioenergy sorghum's long growing season enables root systems to grow deeper and accumulate more biomass than annual grain crops such as maize, attributes that could help restore annual cropland soil organic carbon levels and improve soil productivity. Deep roots active in nutrient transport are positioned to take‐up fertilizer leached deep into soil profiles mitigating potential nutrient run‐off. Bioenergy sorghum's large and deep root system is a key to sustainable production of biofuels, biopower, and bioproducts on annual cropland.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Mesoscale Modeling for Restructuring and Fragmentation in High Burnup UO 2

This report summarizes the mesoscale modeling conducted in fiscal year 2025 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, focusing on the microstructural evolution and restructuring in high burnup UO 2 nuclear fuel and its impact on fuel fragmentation. We developed a pioneering phase-field model to simulate restructuring behavior across different regions of high burnup fuel, including the dark zone and rim region. A grand-potential-based phase-field model is employed to concurrently evaluate subgrain formation and the growth of fission gas bubbles within the fuel. An energy-based subgrain formation criterion was introduced to simulate the restructuring process. The effects of temperature and burnup rate were studied to capture how each of these parameters influences the characteristics of the restructured fuel. Subgrain formation was observed to initiate around existing fission gas bubbles and proceed toward triple junctions, grain boundaries, and grain interiors. Under a given subgrain formation rate, the rate of restructuring increases with rising fuel temperature. The restructuring occurs faster with higher burnup rate. A restructuring bias was observed within the microstructure, due to the variation in defect accumulation when comparing different grains. Microstructures corresponding to the dark zone and rim region can be obtained by parameterizing the model with the appropriate defect production rate, as determined based on the burnup rate and temperature. The predicted microstructures are consistent with experimental observations of the restructured regions. Based on the mesoscale simulations, a mechanistic model for restructuring and grain size evolution was implemented in BISON. Thus, this work provides a first-of-its-kind restructuring model for different regions of high-burnup fuel to BISON. It was shown that the model predicts appropriate grain size evolution along fuel radius as those observed in experiments. This model enables BISON to predict the effect of restructuring on the fission gas release. Finally, the phase-field fracture simulation with dark zone specific microstructures were presented to provide the fragmentation criteria for the dark zone. This work introduces a first-of-its-kind restructuring model for high burnup fuel, significantly enhancing BISON's predictive capabilities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1990

Deformation and fracture of thin sheet aluminum-lithium alloys: The effect of cryogenic temperatures

The objective is to characterize the fracture behavior and to define the fracture mechanisms for new Al-Li-Cu alloys, with emphasis on the role of indium additions and cryogenic temperatures. Three alloys were investigated in rolled product form: 2090 baseline and 2090 + indium produced by Reynolds Metals, and commercial AA 2090-T81 produced by Alcoa. The experimental 2090 + In alloy exhibited increases in hardness and ultimate strength, but no change in tensile yield strength, compared to the baseline 2090 composition in the unstretched T6 condition. The reason for this behavior is not understood. Based on hardness and preliminary Kahn Tear fracture experiments, a nominally peak-aged condition was employed for detailed fracture studies. Crack initiation and growth fracture toughness were examined as a function of stress state and microstructure using J(delta a) methods applied to precracked compact tension specimens in the LT orientation. To date, J(delta a) experiments have been limited to 23 C. Alcoa 2090-T81 exhibited the highest toughness regardless of stress state. Fracture was accompanied by extensive delamination associated with high angle grain boundaries normal to the fatigue precrack surface and progressed microscopically by a transgranular shear mechanism. In contrast the two peak-aged Reynolds alloys had lower toughness and fracture was intersubgranular without substantial delamination. The influences of cryogenic temperature, microstructure, boundary precipitate structure, and deformation mode in governing the competing fracture mechanisms will be determined in future experiments. Results contribute to the development of predictive micromechanical models for fracture modes in Al-Li alloys, and to fracture resistant materials.

Wagner, John A.↗
NASA NTRS1972

Development of single crystal membranes

The design and construction of a high pressure crystal growth chamber was accomplished which would allow the growth of crystals under inert gas pressures of 2 MN/sq m (300 psi). A novel crystal growth technique called EFG was used to grow tubes and rods of the hollandite compounds, BaMgTi7O16, K2MgTi7O16, and tubes of sodium beta-alumina, sodium magnesium-alumina, and potassium beta-alumina. Rods and tubes grown are characterized using metallographic and X-ray diffraction techniques. The hollandite compounds are found to be two or three-phase, composed of coarse grained orientated crystallites. Single crystal c-axis tubes of sodium beta-alumina were grown from melts containing excess sodium oxide. Additional experiments demonstrated that crystals of magnesia doped beta-alumina and potassium beta-alumina also can be achieved by this EFG technique.

Stormont, R. W.↗
NASA NTRS1999

Nucleation, Growth, Annealing, and Coagulation of Refractory Oxides and Metals: Recent Experimental Progress and Applications to Astrophysical Systems

Starting with cooling, refractory vapors diluted in significant quantities of H and He there are four processes that most natural systems will undergo: nucleation, growth, annealing, and coagulation. Nucleation is the processes by which the first stable refractory nuclei form in the vapor. These are the seeds onto which the remaining vapors will condense during the growth stage. Solids of any composition will try to arrange themselves into the least energetic configuration, provided that there is sufficient energy available to support such processes as diffusion and the breaking of chemical bonds. There is a significant activation energy associated with the annealing process in refractory solids due to the relatively high energy of the chemical bonds in solids. The grains formed in most cosmochemical systems are extremely small and often tightly coupled to the gas. Because of their small physical cross sections coagulation may be a very slow process unless there is another driving force involved in addition to normal Brownian motion. In what follows we will briefly cover each of these four stages for refractory oxide and metal grains, although in inverse order.

Nuth, J. A.↗
NASA NTRS2004

NanoSIMS Determination of Carbon and Oxygen Isotopic Compositions of Presolar Graphites from the Murchison Meteorite

Graphite from the Murchison density separate KFC1 (2.15 - 2.20 g/cu cm) has previously been studied by combined SEM/EDX and ion microprobe analysis. These studies found several distinct morphological types of graphites and C isotopic compositions that vary over more than 3 orders of magnitude, clearly establishing their presolar origin. Subsequent TEM measurements of a subset of these particles found abundant embedded crystals of metal (Zr, Mo, Ti, Ru) carbides which were incorporated during the growth of the graphites. A new TEM study of a large set of KFC1 graphites led to the discovery of another type of presolar material, Ru-Fe metal. Here we report results of the C and O isotopic measurements in the same graphite sections, which makes it possible for the first time to directly correlate isotopic and TEM data of KFC1 grains.

Stadermann, F. J.↗
DOE OSTI2021

Oxidation Behavior of Welded Fe-Based and Ni-Based Alloys in Supercritical CO 2

Next-generation supercritical CO 2 (sCO 2 ) power cycles will require different classes of alloy throughout the operational temperatures to optimize tradeoff of creep strength, oxidation performance and cost. This will necessitate joining methods such as welding, which might pose compatibility concerns at the joined interfaces. In this study, similar and dissimilar metal welds were generated from a variety of candidate alloys for sCO 2 systems including ferritic/martensitic steels, austenitic steels, and Ni-based superalloys. Samples were extracted from different regions of the welds and exposed to sCO 2 at 550 °C and 20 MPa for 2500 h, then characterized to understand their behavior in this environment. Unsurprisingly, the local oxidation behavior was largely dictated by the Cr content in the underlying metal. High-Cr austenitic steels and Ni alloys formed slow-growing Cr-rich oxide scales with minimal carburization of the underlying metal, while low-Cr ferritic/martensitic steels formed fast-growing Fe-rich oxide scales with significant carburization. Most welds did not show any unusual oxidation behavior at the interfaces, considering the local Cr content. The one exception was the 347H similar metal weld, where a larger grain size and complex grain structure in the fusion zone led to a significantly higher rate of Fe-rich oxide nodule formation compared to the base metal. This suggests that microstructural changes at joined interfaces can play an important role on the oxidation-limited lifetimes in future sCO 2 systems. The composition changes across the interfaces enabled study of the effect of Fe on the growth rate of Cr-rich oxides and of the origins of the sub-surface recrystallization zone that forms beneath them.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Coupled phase field damage and crystal plasticity analysis of intragranular fracture: The role of crystallographic orientation and voids

Damage evolution in engineering metal alloys at the grain scale exhibits significant microstructural heterogeneity and anisotropy. These heterogeneities create local hotspots for stress and strain localization, leading to void nucleation. Crystal orientation influences the active slip systems around voids, affecting lattice rotation and potentially forming discontinuities. At low triaxiality, voids may change shape due to lower stress, rotation, elongation, and coalescence. At high triaxiality, the correlation between crystal orientation and void growth rate becomes stronger, resembling the behavior observed in isolated single crystals. Therefore, understanding the effects of crystal orientation, heterogeneous strain, and defect evolution is crucial for single crystal fracture characterization. Here, in this work, a coupled phase-field damage (PFD) and crystal plasticity (CP) model is implemented within a finite element framework to analyze crystal deformation and failure. The CP method employs a dislocation density-based constitutive model, while intragranular failure is modeled using an anisotropic PFD method. The PFD model considers both the stored energy due to elastic stretching and the energy release due to defect formation and crack formation. A single crystal Al2219 with an intracrystalline spherical void is chosen to analyze fracture. The study finds that fracture propagation is strongly correlated with crystal orientations. This coupled CP-PFD model provides accurate failure prediction in crystalline materials by incorporating the effects of crystal orientations and existing voids. This study demonstrates how the local microstructure and defects influence plastic deformation and failure mechanisms in metal alloys.

Aluminum alloy↗
DOE OSTI2025

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗
NASA NTRS1990

Influence of low-gravity solidification on heterogeneous nucleation in stable iron-carbon alloys

The effect of gravity on the processes of nucleation and growth of Fe-C alloys during directional solidification was investigated in experiments in which inoculated and uninoculated samples of an alloy of the gray iron type were solidified on ground (1 g) and under the low-gravity (low-g) and high-gravity (high-g) conditions obtained by aircraft parabolic flights. It was found that a higher number of grains was obtained during solidification in high-g condition, due to higher convection in high-g. It was demonstrated that grain multiplication due to convection can contribute up to 23 percent of the total number of grains resulting from heterogeneous nucleation of uninoculated samples. For the case of inoculated samples, it was found that the contribution of the convection-induced nucleation can be as high as 30 percent, but can be zero at very low or very high grain numbers.

Tian, Huameng↗
NASA NTRS1992

Sudbury project (University of Muenster-Ontario Geological Survey): Sr-Nd in heterolithic breccias and gabbroic dikes

One major objective of our Sudbury project was to define origin and age of the huge breccia units below and above the Sudbury Igneous Complex (SIC). The heterolithic Footwall Breccia (FB) represents a part of the uplifted crater floor. It contains subrounded fragments up to several meters in size and lithic fragments with shock features (greater than 10 GPa) embedded into a fine- to medium-grained matrix. Epsilon(sub Nd)-epsilon(sub Sr) relationships point to almost exclusively parautochthonous precursor lithologies. The different textures of the matrix reflect the metamorphic history of the breccia layer; thermal annealing by the overlying hot impact melt sheet (SIC) at temperatures greater than 1000 C resulted in melting of the fine crushed material, followed by an episode of metasomatic K-feldspar growth and, finally, formation of low-grade minerals such as actinolite and chlorite. Isotope relationships in the Onaping breccias (Gray and Green Member) are much more complex. All attempts to date the breccia formation failed: Zircons are entirely derived from country rocks and lack the pronounced Pb loss caused by the heat of the slowly cooling impact melt sheet (SIC). Rb-Sr techniques using either lithic fragments of different shock stages or the thin slab method, set time limits for the apparently pervasive alkali mobility in these suevitic breccias. The data array and the intercept in the plots point to a major Rb-Sr fractionation around 1.54 Ga ago. This model age is in the same range as the age obtained for the metasomatic matrix of the FB. Rb-Sr dating of a shock event in impact-related breccias seems to be possible only if their matrix had suffered total melting by the hot melt sheet (FB) or if they contain a high fraction of impact melt (suevitic Onaping breccias), whereas the degree of shock metamorphism in rock or lithic fragments plays a minor role. In the Sudbury case, however, the impact melt in the seuvitic breccias is devitrified and recrystallized, which changed Rb/Sr ratios quite drastically. Therefore, the Onaping breccias give only age limits for alteration and low-grade metamorphism. The Sm-Nd system was not reset during the Sudbury event; clasts as well as the matrix in the FB and in the Onaping breccias show preimpact 'Archean' Nd isotope signatures.

Buhl, D.↗
NASA NTRS2004

Detailed Microstructural Characterization of the Disk Alloy ME3

The advanced powder metallurgy disk alloy ME3 was designed using statistical screening and optimization of composition and processing variables in the NASA/General Electric/Pratt & Whitney HSR/EPM disk program to have extended durability for large disks at maximum temperatures of 600 to 700 C. Scaled-up disks of this alloy were then produced at the conclusion of that program to demonstrate these properties in realistic disk shapes. The objective of the present study was to assess the microstructural characteristics of these ME3 disks at two consistent locations, in order to enable estimation of the variations in microstructure across each disk and across several disks of this advanced alloy. Scaled-up disks processed in the HSR/EPM Compressor/Turbine Disk program had been sectioned, machined into specimens, and tested in tensile, creep, fatigue, and fatigue crack growth tests by NASA Glenn Research Center, in cooperation with General Electric Engine Company and Pratt & Whitney Aircraft Engines. For this study, microstructures of grip sections from tensile specimens in the bore and rim were evaluated from these disks. The major and minor phases were identified and quantified using transmission electron microscopy (TEM). Particular attention was directed to the .' precipitates, which along with grain size can predominantly control the mechanical properties of superalloy disks.

Gabb, Timothy P.↗
NASA NTRS2016

Hydrogen Embrittlement

Hydrogen embrittlement (HE) is a process resulting in a decrease in the fracture toughness or ductility of a metal due to the presence of atomic hydrogen. In addition to pure hydrogen gas as a direct source for the absorption of atomic hydrogen, the damaging effect can manifest itself from other hydrogen-containing gas species such as hydrogen sulfide (H2S), hydrogen chloride (HCl), and hydrogen bromide (HBr) environments. It has been known that H2S environment may result in a much more severe condition of embrittlement than pure hydrogen gas (H2) for certain types of alloys at similar conditions of stress and gas pressure. The reduction of fracture loads can occur at levels well below the yield strength of the material. Hydrogen embrittlement is usually manifest in terms of singular sharp cracks, in contrast to the extensive branching observed for stress corrosion cracking. The initial crack openings and the local deformation associated with crack propagation may be so small that they are difficult to detect except in special nondestructive examinations. Cracks due to HE can grow rapidly with little macroscopic evidence of mechanical deformation in materials that are normally quite ductile. This Technical Memorandum presents a comprehensive review of experimental data for the effects of gaseous Hydrogen Environment Embrittlement (HEE) for several types of metallic materials. Common material screening methods are used to rate the hydrogen degradation of mechanical properties that occur while the material is under an applied stress and exposed to gaseous hydrogen as compared to air or helium, under slow strain rates (SSR) testing. Due to the simplicity and accelerated nature of these tests, the results expressed in terms of HEE index are not intended to necessarily represent true hydrogen service environment for long-term exposure, but rather to provide a practical approach for material screening, which is a useful concept to qualitatively evaluate the severity of hydrogen embrittlement. The effects of hydrogen gas on mechanical properties such as tensile strength, ductility, fracture, low and high cycle fatigue, crack growth rate, and creep rupture are analyzed with respect to the general trends established from the HEE index values. It is observed that the severity of HE effects is also influenced by environmental factors such as pressure, temperature, and hydrogen gas purity. The severity of HE effects is also influenced by material factors such as surface finish, heat treatment, and product forms, compositions, grain direction, and crystal orientations.

Jonathan A. Lee↗
DOE OSTI2021

Contributions of biofilm-induced flow heterogeneities to solute retention and anomalous transport features in porous media

We report that microbial biofilms are ubiquitous within porous media and the dynamics of their growth influence surface and subsurface flow patterns which impacts the physical properties of porous media and large-scale transport of solutes. Here, a two-dimensional pore-scale numerical model was used to evaluate the impact of biofilm-induced flow heterogeneities on conservative transport. Our study integrates experimental biofilm images of Paenibacillus 300A strain in a microfluidic device packed with cylindrical grains in a hexagonal distribution, with mathematical modeling. Biofilm is represented as a synthetic porous structure with locally varying physical properties that honors the impact of biofilm on the porous medium. We find that biofilm plays a major role in shaping the observed conservative transport dynamics by enhancing anomalous characteristics. More specifically, when biofilm is present, the pore structure in our geometry becomes more spatially correlated. We observe intermittent behavior in the Lagrangian velocities that switches between fast transport periods and long trapping events. Our results suggest that intermittency enhances solute spreading in breakthrough curves which exhibit extreme anomalous slope at intermediate times and very marked late solute arrival due to solute retention. The efficiency of solute retention by the biofilm is controlled by a transport regime which can extend the tailing in the concentration breakthrough curves. These results indicate that solute retention by the biofilm exerts a strong control on conservative solute transport at pore-scale, a role that to date has not received enough attention.

54 ENVIRONMENTAL SCIENCES↗