Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “slope”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Rapidity distributions of Z = 1 isotopes and the nuclear symmetry energy from Sn+Sn collisions with radioactive beams at 270 MeV/nucleon

The rapidity distributions of hydrogen isotopes emitted from central collisions of neutron-rich 132 Sn+ 124 Sn and neutron-deficient 108 Sn+ 112 Sn systems at 270 MeV/nucleon were investigated at RIKEN-RIBF. The data are compared with antisymmetrized molecular dynamics (AMD) calculations and the rapidity distributions can be reproduced after adjusting the in-medium nucleon-nucleon cross sections. The double ratios between the two reaction systems taken for the relative yields of deuteron to proton (d/p) and triton to proton (t/p) are further examined in the midrapidity domain, where the adjustments in the AMD calculations do not affect much on them. The d/p and t/p double ratios at midrapidity agree well with the ratio of the system neutron numbers and its squared value, respectively, and the rapidity dependence of these double ratios is consistent with a picture of partial mixing of colliding nuclei. By comparing with the AMD model which shows a strong symmetry energy dependence of the t/p double ratio, the experimental result in the midrapidity domain favors the calculation with a symmetry-energy slope parameter around L = 46 MeV rather than L = 108 MeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Light nuclei collectivity from s NN = 3 GeV Au+Au collisions at RHIC

In high-energy heavy-ion collisions, partonic collectivity is evidenced by the constituent quark number scaling of elliptic flow anisotropy for identified hadrons. A breaking of this scaling and dominance of baryonic interactions is found for identified hadron collective flow measurements in sNN = 3 GeV Au+Au collisions. In this paper, we report measurements of the first- and second-order azimuthal anisotropic parameters, v1 and v2, of light nuclei (d, t, 3He, 4He) produced in sNN = 3 GeV Au+Au collisions at the STAR experiment. An atomic mass number scaling is found in the measured v1 slopes of light nuclei at mid-rapidity. For the measured v2 magnitude, a strong rapidity dependence is observed. Unlike v2 at higher collision energies, the v2 values at mid-rapidity for all light nuclei are negative and no scaling is observed with the atomic mass number. Calculations by the Jet AA Microscopic Transport Model (JAM), with baryonic mean-field plus nucleon coalescence, are in good agreement with our observations, implying baryonic interactions dominate the collective dynamics in 3 GeV Au+Au collisions at RHIC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Probing neutron-skin thickness with free spectator neutrons in ultracentral high-energy isobaric collisions

We show that the yield ratio of free spectator neutrons produced in high-energy 96 Zr+ 96 Zr to 96 Ru+ 96 Ru collisions is a clean probe of the neutron-skin thickness of colliding nuclei and the slope parameter L of the symmetry energy. The idea is demonstrated based on the proton and neutron density distributions via a state-of-the-art Skyrme-Hartree-Fock-Bogolyubov calculation. Among spectator nucleons given by the Glauber model, free spectator neutrons include those from direct production that survive from clusterization as well as those from deexcitation of heavy clusters described by the popular GEMINI model. More free neutrons are produced in collisions of 96 Zr nucleus due to its larger neutron skin, compared to those produced in collisions of 96 Ru nucleus with a smaller neutron skin. The difference of the free spectator neutron yield is further increased with the increasing difference of the neutron-skin thickness between 96 Zr and 96 Ru with a larger L value, and the increase in ultracentral collisions is particularly insensitive to model details and experimental uncertainties. Since the production of free spectator neutrons is not affected by the complicated dynamics in the mid-rapidity region, the ratio of their multiplicities in ultracentral isobaric collisions is a robust observable for constraining the neutron skin and L value.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Constraining nucleon effective masses with flow and stopping observables from the S π RIT experiment

Properties of the nuclear equation of state (EoS) can be probed by measuring the dynamical properties of nucleus-nucleus collisions. In this study, we present the directed flow (υ 1 ), elliptic flow (υ 2 ) and stopping (VarXZ) measured in fixed target Sn + Sn collisions at 270 AMeV with the SπRIT Time Projection Chamber. We perform Bayesian analyses in which EoS parameters are varied simultaneously within the Improved Quantum Molecular Dynamics-Skyrme (ImQMD-Sky) transport code to obtain a multivariate correlated constraint. The varied parameters include symmetry energy, S 0 , and slope of the symmetry energy, L, at saturation density, isoscalar effective mass, $m^{⁎}_{S}$/$m_{N}$ , isovector effective mass, $m^{⁎}_{υ}$/$m_{N}$ and the in-medium cross-section enhancement factor η. We find that the flow and VarXZ observables are sensitive to the splitting of proton and neutron effective masses and the in-medium cross-section. Comparisons of ImQMD-Sky predictions to the SπRIT data suggest a narrow range of preferred values for $m^{⁎}_{S}$/$m_{N}$, $m^{⁎}_{υ}$/$m_{N}$ and η.

Astronomy & Astrophysics↗
DOE OSTI2024

How do mirror charge radii constrain density dependence of the symmetry energy?

It has recently been suggested that differences in the charge radii of mirror nuclei (Δ$R^{mirr}_{ch}$) are strongly correlated with the neutron-skin thickness ($R_{skin}$) of neutron-rich nuclei and with the slope of the symmetry energy ($L$). To test this assumption, we present ab initio calculations of $R_{skin}$ in 48 Ca and 208 Pb, Δ$R^{mirr}_{ch}$ in 36 Ca$-$ 36 S, 38 Ca$-$ 38 Ar, 41 Sc$-$ 41 Ca, 48 Ni$-$ 48 Ca, 52 Ni$-$ 52 Cr, and 54 Ni$-$ 54 Fe mirror pairs, and $L$. Employing the recently developed 34 chiral interaction samples, identified by the history matching approach, we conduct rigorous statistical analysis of correlations among Δ$R^{mirr}_{ch}$, $R_{skin}$ and $L$, accounting for quantified uncertainties from low-energy constants of chiral interaction, chiral effective field theory truncation and many-body method approximation. The ab initio results reveal an appreciable Δ$R^{mirr}_{ch}$ $-$ $L$ correlation in $fp$-shell mirror pairs. However, contrary to previous studies, the present calculation finds that the studied $sd$-shell mirror pairs do not exhibit any Δ$R^{mirr}_{ch}$ $-$ $L$ correlation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Quantum entanglement in nuclear fission

Nuclear fission presents a unique example of quantum entanglement in strongly interacting many-body systems. A heavy nucleus can split into hundreds of combinations of two complementary fragments in the fission process. The entanglement of fragment wave functions is persistent even after separation and impacts the partition of particles and energies between fragments. Based on microscopic dynamical calculations of the fission of 240 Pu, this work finds that dynamical quantum entanglement is indispensable in the appearance of sawtooth distributions of average excitation energies of fragments and thus neutron multiplicities, but not in average neutron excess of fragments. Both sawtooth slopes from particle-number projections are found to be steep – a feature which can be alleviated by random fluctuations. The persistent entanglement is mainly due to non-adiabatic dynamics since the final splitting is so fast that the non-localization of wave functions is kept during the separation. These findings may impact the understanding of quantum entanglement more broadly in mesoscopic systems

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Rheological behavior and fiber spinning of polyacrylonitrile (PAN)/Carbon nanotube (CNT) dispersions at high CNT loading

Rheological studies play an important role in polymer processing including fiber spinning. In the current work, rheological behavior has been studied for polyacrylonitrile (PAN)/carbon nanotube (CNT) dispersions in dimethyl formamide (DMF), where CNT loading is as high as 15 wt% of the total solids (polymer + CNT). The presence of CNTs increased the elastic-like and shear thinning behavior of the dispersions. The terminal slope of the log G' versus log G" curve, which was used to characterize the homogeneity of solution or dispersion, is dependent on the CNT concentration but independent of temperature between 0 and 60 °C. PAN/CNT-10 (90 wt% PAN +10 wt% CNT) and PAN/CNT-15 were gel spun at different conditions, with small diameter capillary at high shear rate and large diameter capillary at relatively low shear rate. Additionally, for PAN/CNT-10, fiber spinning with lower molecular weight polymer showed increased continuous jetting time when compared to spinning with higher molecular weight polymer. In addition, the rheological behavior of PAN/cellulose nanocrystal (CNC) dispersion with 20–60 wt% CNC was also compared to that of the PAN/CNT dispersions. PAN/CNC and PAN/CNT share similarities in both rheological behavior and fiber spinning at high filler loading.

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

Hydrokinetic energy harvesting from slow currents using flow-induced oscillations

To harness marine hydrokinetic energy from slow flows, which constitute the majority of currents, tides, and rivers, new Passive Turbulence Control (PTC), consisting of large turbulence stimulators, is tested experimentally on circular cylinders on springs. Here, this study experimentally investigates the effect of PTC on the onset of Flow-Induced Oscillations (FIO) and particularly the relative onset of Vortex-Induced Vibrations (VIV) and galloping. Experiments are conducted in the Low Turbulence Free Surface Water Channel, University of Michigan. Fixed are: mass ratio m* = 1.48, aspect ratio l/D = 10.29, and total damping ratio ζ = 0.04. Parameters are: cylinder diameter D, spring stiffness K, PTC location and height, and flow speed U$\in$[0.36 m/s-1.45 m/s]. Placing the leading edge of PTC at 40–60° induces high amplitude FIO while placement at 10–20° suppresses FIO. As PTC height increases, VIV and galloping initiate earlier and exhibit higher amplitude with a steeper slope. Lower spring stiffness initiates VIV earlier by reducing the oscillator natural frequency in water. Even though large PTC maintained its effectiveness in initiating galloping early, it has no effect on the earlier initiation of VIV, which starts at a nearly fixed reduced velocity. Lower spring stiffness and large PTC enable power generation at low current speed (0.2 m/s).

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Influence of living grass Roots and endophytic fungal hyphae on soil hydraulic properties

Soil hydraulic properties are often estimated based on laboratory data or pedotransfer functions dependent on soil physical properties, which often do not consider potential impacts of soil roots or fungal hyphae. Here, we first review current knowledge of how these soil biotic components affect hydraulic properties, then we conducted laboratory experiments to specifically test if the presence of roots and mycorrhizal fungi had a significant effect on the hydraulic properties of two soils with contrasting textures: Flint sand and Hamblen silt loam. Soil cores were seeded with (Panicum virgatum) and grown in a greenhouse over three separate growth periods. The endophytic fungus Serendipita indica was injected as liquid inoculant into designated mycorrhizal cores. Saturated hydraulic conductivity (K sat ) measurements were made with a constant head permeameter, and soil water retention curves were obtained by the evaporation method, supplemented at the dry end for Hamblen silt loam with water activity meter data. Retention curve parameters were obtained by fitting the van Genuchten equation to the resulting measurements. Mean root volume ratios were higher in the mycorrhizal inoculated treatment than in the uninoculated treatment for both soils. For Flint sand, analysis of variance revealed that K sat was reduced by the presence of roots as compared to bare soil. This was likely due to roots clogging soil pores. Results also indicated the presence of roots changed the shape of the water retention curve for Flint sand by increasing water content at saturation and by reducing the slope of the curve. These changes suggested roots created additional porosity and broadened the pore-size distribution. The presence of mycorrhizal fungi accentuated the root effects. The influence of roots and mycorrhizal fungi on hydraulic properties was less obvious for the Hamblen silt loam, as none of the treatments differed from each other at p < 0.05. The results highlight the necessity to consider the impact of root and fungal structures on models of soil hydraulic properties.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Remote sensing and statistical analysis of the effects of hurricane María on the forests of Puerto Rico

Widely recognized as one of the worst natural disaster in Puerto Rico's history, hurricane María made landfall on September 20, 2017 in southeast Puerto Rico as a high-end category 4 hurricane on the Saffir-Simpson scale causing widespread destruction, fatalities and forest disturbance. This study focused on hurricane María's effect on Puerto Rico's forests as well as the effect of landform and forest characteristics on observed disturbance patterns. In this work, we used Google Earth Engine (GEE) to assess the severity of forest disturbance using a disturbance metric based on Landsat 8 satellite data composites with pre and post-hurricane María. Forest structure, tree phenology characteristics, and landforms were obtained from satellite data products, including digital elevation model and global forest canopy height. Our analyses showed that forest structure, and characteristics such as forest age and forest type affected patterns of forest disturbance. Among forest types, highest disturbance values were found in sierra palm, transitional, and tall cloud forests; seasonal evergreen forests with coconut palm; and mangrove forests. For landforms, greatest disturbance metrics was found at high elevations, steeper slopes, and windward surfaces. As expected, high levels of disturbance were also found close to the hurricane track, with disturbance less severe as hurricane María moved inland. Results demonstrated that forest and landform characteristics accounted for 34% of the variation in spatial forest spectral disturbance patterns. This study demonstrated an informative regional approach, combining remote sensing with statistical analyses to investigate factors that result in variability in hurricane effects on forest ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Analysis of alternative climate datasets and evapotranspiration methods for the Upper Mississippi River Basin using SWAT within HAWQS

This study reports the application of Soil andWater Assessment Tool (SWAT) within the Hydrologic andWater Quality System(HAWQS) on-line platform, for the UpperMississippi River Basin (UMRB). TheUMRB is an important ecosystem located in the north central U.S. that is experiencing a range of ecological stresses. Specifically, testing of SWAT was performed for: (1) Hargreaves (HG) and Penman-Monteith (PM) PET methods, and (2) Livneh, National Climatic Data Center (NCDC) and Parameter-elevation Regressions on Independent Slopes Model (PRISM) climate datasets. The Livneh-PM combination resulted in the highest average annual water yield of 380.6 mm versus the lowest estimated water yield of 193.9 mm for the Livneh-HG combination, in response to 23-year uncalibrated simulations. Higher annual ET and PET values were predicted with HG method versus the PM method for all three weather datasets in response to the uncalibrated simulations, due primarily to higher HG-based estimates during the growing season. Based on these results, it was found that the HG method is the preferred PET option for the UMRB. Initial calibration of SWAT was performed using the Livneh data and HGmethod for threeMississippi Rivermain stem gauge sites, which was followed by spatial validation at 10 other gauge sites located within the UMRB streamnetwork. Overall satisfactory results were found for the calibration and validation gauge sites, with the majority of R2 values ranging between 0.61 and 0.82, Nash- Sutcliffe modeling efficiency (NSE) values ranging between 0.50 and 0.79, and Kling-Gupta efficiency (KGE) values ranging between 0.61 and 0.84. The results of an additional experimental suite of six scenarios,which represented different combinations of climate data sets and calibrated parameters, revealed that suggested statistical criteria were again satisfied by the different scenario combinations. Overall, the PRISM data exhibited the strongest reliability for the UMRB.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Zr 3 Mn 3 Sn 4 Ga: A new hetero-kagome bilayer antiferromagnet

Here, we report the magnetic and electrical transport properties of single-crystalline Zr 3 Mn 3 Sn 4 Ga, featuring two distinct kagome lattices: a non-magnetic breathing Zr 3 Sn 4 lattice and a magnetic intact Mn 3 Ga lattice. The material undergoes an antiferromagnetic phase transition at T N = 87 K, with neutron diffraction confirming commensurate ordering characterized by k = (1/3,1/3,0). Transport measurements show metallic behavior, a resistivity anomaly near T N , and 12% magnetoresistance at 2 K under 9 T. Deviations from the conventional second-order power law, along with negative magnetoresistance and nonlinear Hall slope variations near T N , suggest strong magneto-electronic coupling. Resonant photoemission spectroscopy identifies Zr 4d and Mn 3d orbitals as dominant contributors to the valence band, linking the material’s unique electronic properties to its kagome layers. Zr 3 Mn 3 Sn 4 Ga offers a valuable platform to study interactions between magnetism and topological electronic bands in hetero-kagome systems using the co- existence of magnetic and non-magnetic kagome layers and its tunable electronic structure.

Crystal structure↗
DOE OSTI2023

Land surface dynamics and meteorological forcings modulate land surface temperature characteristics

This study examines the effect of land cover, vegetation health, climatic forcings, elevation heat loads, and terrain characteristics (LVCET) on land surface temperature (LST) distribution over West Africa (WA). We employ fourteen machine-learning models, which preserve nonlinear relationships, to downscale LST and other predictands while preserving the geographical variability of WA. Our results showed that the random forest model performs best in downscaling predictands. This is important for the sub-region since it has limited access to mainframes to power multiplex machine-learning algorithms. In contrast to the northern regions, the southern regions consistently exhibit healthy vegetation. Also, areas with unhealthy vegetation coincide with hot LST clusters. The positive Normalized Difference Vegetation Index (NDVI) trends in the Sahel underscore rainfall recovery and subsequent Sahelian greening. The southwesterly winds cause the upwelling of cold waters, lowering LST in southern WA and highlighting the cooling influence of water bodies on LST. Identifying regions with elevated LST is paramount for prioritizing greening initiatives, and our study underscores the importance of considering LVCET factors in urban planning. Topographic slope-facing angles, heat loads, and diurnal anisotropic heat all contribute to variations in LST, emphasizing the need for a holistic approach when designing resilient and sustainable landscapes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Effect of electrode rinse solutions on the electrodialysis of concentrated salts

The data presented indicate that modest attention to the chemistry of the electrode rinse solution can yield as much as 56% improvement in process efficiency. When operating electrodialysis (ED) with salt concentrations greater than 10,000 mg/l (as NaCl) the amperage can be relatively high even at low voltages. To maximize the overall rate of ion transport from diluate to concentrate, there is a need to minimize resistances in the electrode cells since these can represent as much as 30% of the total resistances in the entire process and in this particular ED unit, represented resistance roughly equivalent to that within the membrane stack. This work describes the performance of a 10-cell pair, 200 cm 2 (0.02 m 2 ) per membrane, pilot ED unit operated in batch mode at 5 V potential. The feedstock to the stack was varied from 0.5% to 6% NaCl. Results from Volt-amp profiles (2.0–15 V) were used as the rationale for choosing the preferred electrode rinse solution, 90 g/l (kg/m 3 ) disodium sulfate at pH 12.5. The recommended solution of 30 g/l (kg/m 3 ) disodium sulfate at neutral pH was tested against stronger solutions (60, 90, and 120 g/l (kg/m 3 ) disodium sulfate) at neutral pH and with 1 g/l (kg/m 3 ) sodium hydroxide added to yield solutions at pH 12.5. There were stark differences in the performance of the ED unit between the solutions at pH 7 and those at pH 12.5. The slopes of the Volt-amp profiles improved indicating less resistance to ion flow at pH 12.5. Most of the difference can be attributed to the greater conductivity of sodium hydroxide compared to disodium sulfate, where even a modest addition of 1 g/l (kg/m 3 ) sodium hydroxide increased the conductivity of the solution between 8 and 20% depending on the concentration of disodium sulfate. Most notably, the initiation voltage was 2.39 ± 0.08 V with solutions at pH 7 and 1.99 V ± 0.04 V when the rinse solutions were around pH 12.5, a difference of around 0.42 V. Since the full batch runs were performed at 5.0 V, this 0.42 V differential represented a loss of 8% in process efficiency. A series of tests showed a shift in initiation voltage occurred between pH 11.5 and 12.5. Nernst equations for water electrolysis are coupled with a simple flux model that estimates the concentration of hydroxide and hydronium ions at the electrode surfaces. This model indicates that the potential demanded by the electrodes is directly linked to the current. Furthermore, there is a steep drop in the potential demanded by the anode when there is sufficient hydroxide available to generate oxygen from hydroxide as opposed to generation from water. The voltage demand is minimal at pH 12.5 and reflective of the reference standard voltage (just over 1.23 V predicted at 0.1 amps). However, the flux model predicts that the voltage demanded by the electrodes at pH 7 is 1.61 V at 0.1 amps (the lowest amperage measurable with the power pack). In conclusion, comparing this result to the reference standard (E 0 = -1.23 V) yields a satisfactory explanation for the observed difference below reference standard of 0.42 V.

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

The effects of warming and soil chemistry on bacterial community structure in Arctic tundra soils

Soil microbial communities in the Arctic play a critical role in regulating the global carbon (C) cycle. Vast amounts of C are stored in northern latitude soils, and rising temperatures in the Arctic threaten to thaw permafrost, making relatively inaccessible C sources more available for mineralization by soil microbes. Few studies have characterized how microbial community structure responds to thawing permafrost in the context of varying soil chemistries associated with contrasting tundra landscapes. Here, we subjected active layer and permafrost soils from upland and lowland tundra sites on the North Slope of Alaska to a soil-warming incubation experiment and compared soil bacterial community profiles (obtained by 16S rRNA amplicon sequencing) before and after incubation. The influence of soil composition (characterized by mid-infrared [MIR] spectroscopy) on bacterial community structure and class abundance was analyzed using redundancy and correlation analyses. We found increased abundances of Alphaproteobacteria, Gammaproteobacteria, and Bacteroidetes [Sphingobacteriia] post incubation, particularly in permafrost soils. The categorical descriptors site and soil layer had the most explanatory power in our predictive models of bacterial community structure, highlighting the close relationship between soil bacteria and the soil environment. Specific soil chemical attributes characterizing the soil environments that were found to be the best predictors include MIR spectral bands associated with inorganic C, silicates, amide II (C=N stretch), and carboxylics (C-O stretch), and MIR peak ratios representing C substrate quality. Overall, these results further characterize soil bacterial community shifts that may occur as a frozen environment with limited access to C sources, such as is found in undisturbed permafrost, transitions to warmer and more C-available environments, such as is predicted in thawing permafrost due to climate change.

16S rRNA sequencing↗
DOE OSTI2021

A non-intrusive optical approach to characterize heliostats in utility-scale power tower plants: Flight path generation/optimization of unmanned aerial systems

A newly developed in situ non-intrusive optical (NIO) approach has been developed to survey various types of heliostat optical errors for a concentrating solar power (CSP) tower plant. To measure mirror surface slope error, facet canting error, and heliostat tracking error at a sub-milliradian accuracy, NIO requires several reflection images scanned over each individual heliostat. For a utility-scale plant that typically includes more than 10,000 heliostats, an unmanned aerial system (UAS) is crucial for efficient implementation of the NIO method. Here, we develop a flight path generation/optimization algorithm to plan more efficient UAS paths to collect NIO data over a utility-scale heliostat field. The algorithm considers NIO data requirements, all potential constraints, optimization within each subfield, and operational flexibility. Case studies are presented to illustrate the feasibility and robustness of the developed flight path algorithm. The path planning algorithm may also find applications elsewhere, such as drone-driven imaging under extreme conditions.

14 SOLAR ENERGY↗
DOE OSTI2020

Phase transformation in boron under shock compression

Using first-principles molecular dynamics, we calculated the equation of state and shock Hugoniot of various boron phases. Here, we find a large mismatch between Hugoniots based on existing knowledge of the equilibrium phase diagram and those measured by shock experiments, which could be reconciled if the α-B 12 /β phases are significantly over-pressurized in boron under shock compression. Our results also indicate that there exist an anomaly and negative Clapeyron slope along the melting curve of boron at 100 GPa and 1500–3000 K. These results enable an in-depth understanding of matter under shock compression, in particular the significance of compression-rate dependence of phase transitions and kinetic effects in experimental measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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↗