Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mass diffusivity”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A low-mass hub–filament with double centre revealed in NGC 2071-North

We present the first analysis in NGC 2071-North as a resolved hub–filament featuring a double centre. This ∼1.5 × 1.5 pc 2 scale filament hub contains ∼500 M ⊙ . Seen from Planck, magnetic field lines may have facilitated the gathering of material at this isolated location. The energy balance analysis, supported by infalling gas signatures, reveals that these filaments are currently forming stars. Herschel 100 $\mu$m emission concentrates in the hub, at IRAS 05451+0037 and LkHα 316, and presents diffuse lobes and loops around them. We suggest that such a double centre could be formed, because the converging locations of filament pairs are offset, by 2.3 arcmin (0.27 pc). This distance also matches the diameter of a hub ring, seen in column density and molecular tracers, such as HCO + (1–0) and HCN (1–0), that may indicate a transition and the connection between the hub and the radiating filaments. We argue that all of the three components of the emission star LkHα 316 are in physical association. We find that a ∼0.06-pc-sized gas loop, attached to IRAS 05451+0037, can be seen at wavelengths all the way from the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS)-i to Herschel 100 $\mu$m. These observations suggest that both protostars at the double hub centre are interacting with the cloud material. In our 13 CO data, we do not seem to find the outflow of this region that was identified in the 1980s with much lower resolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Non-universal stellar initial mass functions: large uncertainties in star formation rates at z ≈ 2–4 and other astrophysical probes

ABSTRACT We explore the assumption, widely used in many astrophysical calculations, that the stellar initial mass function (IMF) is universal across all galaxies. By considering both a canonical broken-power-law IMF and a non-universal IMF, we are able to compare the effect of different IMFs on multiple observables and derived quantities in astrophysics. Specifically, we consider a non-universal IMF that varies as a function of the local star formation rate, and explore the effects on the star formation rate density (SFRD), the extragalactic background light, the supernova (both core-collapse and thermonuclear) rates, and the diffuse supernova neutrino background. Our most interesting result is that our adopted varying IMF leads to much greater uncertainty on the SFRD at $z \approx 2-4$ than is usually assumed. Indeed, we find an SFRD (inferred using observed galaxy luminosity distributions) that is a factor of $\gtrsim 3$ lower than canonical results obtained using a universal IMF. Secondly, the non-universal IMF we explore implies a reduction in the supernova core-collapse rate of a factor of $\sim 2$, compared against a universal IMF. The other potential tracers are only slightly affected by changes to the properties of the IMF. We find that currently available data do not provide a clear preference for universal or non-universal IMF. However, improvements to measurements of the star formation rate and core-collapse supernova rate at redshifts $z \gtrsim 2$ may offer the best prospects for discernment.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Integrated simulations of premagnetized and self-magnetizing dynamic screw pinch-driven MagLIF

Magnetically driven implosions such as in magnetized liner inertial fusion (MagLIF) on the Z accelerator suffer from magneto-Rayleigh–Taylor instabilities (MRTI) that dynamically redistribute liner mass during implosion, limiting fusion fuel compression and confinement, which ultimately degrades performance. Driving the implosion with an initially helical drive field that dynamically shifts the direction of the magnetic field surrounding the liner (i.e., a dynamic screw pinch, DSP) is a method proposed to mitigate MRTI in-flight and improve target performance. In DSPs, the axial drive magnetic field component implodes the liner and diffuses through the shocked, melted liner material into the fuel throughout the implosion. Liners can be designed to enable enough axial magnetic flux to diffuse through the liner material to effectively magnetize the fuel region without the need of an initial axial magnetic field (i.e., from external field coils). We present results from three-dimensional radiation-magnetohydrodynamic simulations of MagLIF implosions employing drive magnetic fields composed of axial and azimuthal components (a helical drive field). These simulated DSP-driven MagLIF targets demonstrate improved fuel conditions and thermonuclear yield compared to a traditional MagLIF target implosion. Synthetic x-ray radiography of the imploding liner material and x-ray emission images of the fuel region at the time of peak neutron yield rate indicate superior implosion morphology for DSP-MagLIF implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Photodissociation and Theory to Investigate Uranium Oxide Cluster Cations

Uranium oxide cluster cations of the form U n O m + are produced by laser vaporization of a depleted uranium rod in a pulsed supersonic expansion. Ions are mass-analyzed and mass-selected with a time-of-flight spectrometer and studied with UV laser multiphoton dissociation. Cations of the stoichiometry UO 2 (UO 3 ) n + were observed as photofragments from all photodissociated cluster cations. (UO 3 ) n + clusters were also observed to result from dissociation of larger (UO 3 ) n + clusters, with UO 3 neutral as a common leaving group. Electronic structure calculations were used to investigate the stability of the prominent uranium oxide cluster cations using density functional theory (DFT) with the hybrid B3LYP exchange-correlation functional and at the CCSD(T) level and cc-pVnZ-PP basis sets (n = D,T), including diffuse orbitals as computational expense and availability permitted. Clustering energies, relative energies and dissociation energies of the cations are reported. The lowest energy neutral (UO 3 ) n clusters up to n = 3 are rings, n = 4 and 5 are chains with very low energy rings, and n = 6 is 3D. Finally, the lowest energy structures for UO 2 (UO 3 ) n + are composed of uranyl-like UO 2 + units bound by bridging oxygens to other UO 2 2+ units for n = 2 and 3, and for n = 4 a more complex 3D structure is predicted.

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

Air Dehumidification using Ionic Liquid-Based Fiber Bundle Membrane Contactor

Air dehumidification is essential since excess moisture in the buildings causes discomfort to the occupants, encourages the production of air pathogens such as mold or mildew, and causes corrosion and rotting that degrade building materials. Existing moisture removal processes are mainly focused on condensation and desiccant (liquid or solid) techniques with direct contact between air and desiccant. However, these methods are energy-intensive, or desiccant might be lost or cause corrosion in the process. The main objective of this study is to investigate an ionic liquid-based liquid desiccant absorber based on a membrane fiber bundle. The novel system consists of membrane contactors containing bundles of 10,000 polypropylene fibers. The fiber bundles have 0.3 micron diameter, with ionic liquid flowing inside, and air flowing outside. The fibers provide a high contact area among phases: 1.4 m2 contact surface area in a 0.00015 m3 volume. The ionic liquid as a sorbent stipulates the selectivity for air, prevents the loss of solvent in the operation due to negligible volatility, provides fast diffusion due to low viscosity compared to common ionic liquids, and has high affinity and solubility in water. The dehumidification capacity of the membrane technology was experimentally investigated. The experimental results show that the ionic-liquid based membrane system can effectively remove excess moisture from the air. The novel fiber bundle dehumidification system has a system volume of 0.00798 m3 (7.98 L) and active heat and mass transfer surface area of 8.4 m2. It achieved an average dehumidification of 320 ± 25 W with a volumetric air flowrate of 3.1 m3/min (108 ft3/min).

Turnaoglu, Tugba↗
DOE OSTI2022

Emergence of a nocturnal low-level jet from a broad baroclinic zone

An analytical model is presented for the generation of a Blackadar-like nocturnal low-level jet in a broad baroclinic zone. The flow is forced from below (flat ground) by a surface buoyancy gradient and from above (free atmosphere) by a constant pressure gradient force. Diurnally-varying mixing coefficients are specified to increase abruptly at sunrise and decrease abruptly at sunset. With attention restricted to a surface buoyancy that varies linearly with a horizontal coordinate, the Boussinesq-approximated equations of motion, thermal energy, and mass conservation reduce to a system of one-dimensional equations that can be solved analytically. Sensitivity tests with southerly jets suggest that (i) stronger jets are associated with larger decreases of the eddy viscosity at sunset (as in Blackadar theory), (ii) the nighttime surface buoyancy gradient has little impact on jet strength, and (iii) for pure baroclinic forcing (no free-atmosphere geostrophic wind), the nighttime eddy diffusivity has little impact on jet strength, but the daytime eddy diffusivity is very important and has a larger impact than the daytime eddy viscosity. The model was applied to a jet that developed in fair weather conditions over the Great Plains from southern Texas to northern South Dakota on 1 May 2020. The ECMWF Reanalysis v5 (ERA5) for the afternoon prior to jet formation showed that a broad north-south-oriented baroclinic zone covered much of the region. The peak model-predicted winds were in good agreement with ERA5 winds and lidar data from the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) central facility in north-central Oklahoma.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

CFD modeling of high-flux plate-and-frame membrane modules for industrial carbon capture

In this work, we study the application of membrane-based separation systems for carbon capture, considering plate-and-frame membrane modules. The successful deployment of membrane CO2 capture system relies on high-performing membranes as well as effective membrane modules that can fully exploit the developed membranes. A plate-and-frame membrane module is especially attractive for CO2 capture from industrial flue gas due to its lower pressure drop compared to its counterparts such as spiral wound modules and hollow fiber modules. To design better plate-and-frame modules, we investigate their basic unit - a single membrane stack through a combination of computational modeling and experimental investigations. The modeling approach is based on Computational Fluid Dynamics (CFD) to represent a multiphysics problem, including the fluid flow and diffusion processes within a membrane module. We use experimental data collected under different operating conditions to validate the CFD model. Numerical results suggest a good agreement between experiments and model outputs for the CO2 recovery, CO2 mole fraction in the retentate and permeate, and stage-cut. The CFD model is able to predict accurately the flow behavior, providing valuable insights on the effects of fluid dynamics on mass transfer of CO2. We also carry out a sensitivity analysis to identify the effect of key parameters on the CO2 recovery and the CO2 purity of the outlet streams.

Dosso, Cheick↗
DOE OSTI2022

SDSS-IV MaNGA: Understanding Ionized Gas Turbulence Using Integral Field Spectroscopy of 4500 Star-forming Disk Galaxies

The Sloan Digital Sky Survey MaNGA program has now obtained integral field spectroscopy for over 10,000 galaxies in the nearby universe. We use the final MaNGA data release DR17 to study the correlation between ionized gas velocity dispersion and galactic star formation rate, finding a tight correlation in which σ Hα from galactic H ιι regions increases significantly from ~18–30 km s –1 , broadly in keeping with previous studies. In contrast, σ Hα from diffuse ionized gas increases more rapidly from 20–60 km s –1 . Using the statistical power of MaNGA, we investigate these correlations in greater detail using multiple emission lines and determine that the observed correlation of σ Hα with local star formation rate surface density is driven primarily by the global relation of increasing velocity dispersion at higher total star formation rate, as are apparent correlations with stellar mass. Assuming H ιι region models consistent with our finding that σ [O III] < σ Hα < σ [O I] , we estimate the velocity dispersion of the molecular gas in which the individual H ιι regions are embedded, finding values σ Mol = 5–30 km s –1 consistent with ALMA observations in a similar mass range. Finally, we use variations in the relation with inclination and disk azimuthal angle to constrain the velocity dispersion ellipsoid of the ionized gas σ z /σ r = 0.84 ± 0.03 and σ Φ /σ r = 0.91 ± 0.03, similar to that of young stars in the Galactic disk. Our results are most consistent with the theoretical models in which turbulence in modern galactic disks is driven primarily by star formation feedback.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Self-interacting Dark Matter Interpretation of Crater II

Abstract The satellite galaxy Crater II of the Milky Way is extremely cold and exceptionally diffuse. These unusual properties are challenging to understand in the standard model of cold dark matter. We use controlled N -body simulations to investigate the formation of Crater II in self-interacting dark matter (SIDM), where dark matter particles can scatter and thermalize. Taking the orbit motivated by the measurements from Gaia Early Data Release 3, we show a strong self-interacting cross section per particle mass of 60 cm 2 g −1 is favored for Crater II. The simulated SIDM halo, with a 1 kpc core, leads to both a low stellar velocity dispersion and a large half-light radius for Crater II. These characteristics remain robust regardless of the initial stellar distribution.

Astronomy & Astrophysics↗
DOE OSTI2023

Atomic Dynamics of Multi‐Interfacial Migration and Transformations

Redox-induced interconversions of metal oxidation states typically result in multiple phase boundaries that separate chemically and structurally distinct oxides and suboxides. Directly probing such multi-interfacial reactions is challenging because of the difficulty in simultaneously resolving the multiple reaction fronts at the atomic scale. Using the example of CuO reduction in H 2 gas, a reaction pathway of CuO → monoclinic m-Cu 4 O 3 → Cu 2 O is demonstrated and identifies interfacial reaction fronts at the atomic scale, where the Cu 2 O/m-Cu 4 O 3 interface shows a diffuse-type interfacial transformation; while the lateral flow of interfacial ledges appears to control the m-Cu 4 O 3 /CuO transformation. Together with atomistic modeling, it is shown that such a multi-interface transformation results from the surface-reaction-induced formation of oxygen vacancies that diffuse into deeper atomic layers, thereby resulting in the formation of the lower oxides of Cu 2 O and m-Cu 4 O 3 , and activate the interfacial transformations. In conclusion, these results demonstrate the lively dynamics at the reaction fronts of the multiple interfaces and have substantial implications for controlling the microstructure and interphase boundaries by coupling the interplay between the surface reaction dynamics and the resulting mass transport and phase evolution in the subsurface and bulk.

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

Investigating intermolecular interactions among CO 2 , water and PEEK-ionene membrane using cryo ToF-SIMS and isotopic labeling

Cryogenic time-of-flight secondary ion mass spectrometry (cryo ToF-SIMS) has emerged as a powerful tool for investigating molecular interactions, speciation, and dynamics in materials for CO 2 capture. In this study, we apply cryo ToF-SIMS to probe interactions between CO 2 , water, and PEEK-ionene membranes—a promising material for direct CO 2 capture due to its selectivity, durability, and efficiency. Despite this potential, the mechanisms governing CO 2 diffusion and the influence of water vapor on CO 2 behavior remain unclear. To address this, we loaded PEEK-ionene membranes with 13 CO 2 and D 2 O and employed cryo ToF-SIMS to visualize the 3D distribution of CO 2 and water within the membrane. While prior studies suggest that 13 CO 2 is absorbed under ambient conditions, our cryo ToF-SIMS analysis revealed no enhancement of the 13 C/ 12 C ratio, suggesting weak CO 2 -membrane interactions. As a result, CO 2 vaporizes even at low temperatures (−140°C) under vacuum conditions. In contrast, D 2 O displayed a relatively homogeneous distribution in the membrane, suggesting stronger water-membrane interactions via hydrogen bonding (18–20 kJ/mol). Interestingly, CO 2 was not detected in D 2 O-loaded membranes, indicating minimal interference from water vapor on CO 2 diffusion. As a comparison, the cryo ToF-SIMS data show that 13 CO 2 can readily react with a basic Na 2 CO 3 aqueous solution to form NaH 13 CO 3 . These findings demonstrate cryo ToF-SIMS as a critical technique for understanding gas-water-membrane interactions, offering insights for membrane functionalization to improve CO 2 capture efficiency.

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

Hot Water Extraction: Short Rotation Willow, Mixed Hardwoods, and Process Considerations

Short rotation woody crops (SRWC) like shrub willow are highly productive biomass resources of interest for energy and fuel applications. Hot water extraction (HWE) as an upgrading tool to enable the use of willow biomass in pellet applications has been proposed, and is of increasing interest. This study treats willow and mixes of willow and conventional mixed hardwood feedstock with HWE in a tumbling laboratory reactor to elucidate the effects of time, temperature, feedstock mixes, and other process considerations (water:biomass ratio, presteaming, counter-current processing) on mass removals and other extraction outcomes (e.g., sugar, acetate, and furan yields). Results demonstrated alignment of extraction outcomes with P-factor from 155 °C to 175 °C, with a good compromise of removed mass and co-product potential in the range from 575–800 P-factor. The preferred condition was chosen as 575 P-factor. HWE of mixes of willow and hardwood feedstocks showed a linear response of extraction outcomes to willow:hardwood ratios. Testing of water:biomass ratios demonstrated that this is a significant consideration, with each outcome being affected somewhat differently, and indicating that HWE is more diffusion dependent than expected. Presteaming shows little to no effect on extraction outcomes, while multi-stage cooks simulating counter-current operation indicate a significant potential value in counter-current extraction.

09 BIOMASS FUELS↗
DOE OSTI2022

Mass Transport in Membrane Systems: Flow Regime Identification by Fourier Analysis

The numerical calculation of local mass distributions in membrane systems by computational fluid dynamics (CFD) offers indispensable benefits. However, the concept to calculate such distributions in response to separate variations of operation conditions (OCs) makes it difficult to address overall, flow-physics-related questions, which require the consideration of the collective interaction of OCs. It is shown that such understanding-related relationships can be obtained by the analytical solution of the advection–diffusion equation considered. A Fourier series model (FSM) is presented, which provides exact solutions of an advection–diffusion equation for a wide range of OCs. On this basis, a new zeroth-order model is developed, which is very simple and as accurate as the complete FSM for all conditions of practical relevance. Advection-dominated blocked and diffusion-dominated unblocked flow regimes are identified (depending on a Péclet number which compares the flow geometry with a length scale imposed by the flow), which implies relevant requirements for the use of lab results for pilot- and full-scale applications. Analyses reveal the equivalence of variations of OCs, which offers a variety of options to accomplish desired flow regime changes.

Heinz, Stefan (ORCID:0000000248712416)↗
DOE OSTI2020

3D‐Printed Structure Boosts the Kinetics and Intrinsic Capacitance of Pseudocapacitive Graphene Aerogels

Abstract The performance of pseudocapacitive electrodes at fast charging rates are typically limited by the slow kinetics of Faradaic reactions and sluggish ion diffusion in the bulk structure. This is particularly problematic for thick electrodes and electrodes highly loaded with active materials. Here, a surface‐functionalized 3D‐printed graphene aerogel (SF‐3D GA) is presented that achieves not only a benchmark areal capacitance of 2195 mF cm −2 at a high current density of 100 mA cm −2 but also an ultrahigh intrinsic capacitance of 309.1 µF cm −2 even at a high mass loading of 12.8 mg cm −2 . Importantly, the kinetic analysis reveals that the capacitance of SF‐3D GA electrode is primarily (93.3%) contributed from fast kinetic processes. This is because the 3D‐printed electrode has an open structure that ensures excellent coverage of functional groups on carbon surface and facilitates the ion accessibility of these surface functional groups even at high current densities and large mass loading/electrode thickness. An asymmetric device assembled with SF‐3D GA as anode and 3D‐printed GA decorated with MnO 2 as cathode achieves a remarkable energy density of 0.65 mWh cm −2 at an ultrahigh power density of 164.5 mW cm −2 , outperforming carbon‐based supercapacitors operated at the same power density.

asymmetric supercapacitors↗
DOE OSTI2025

A computational investigation of high-flux, plate-and-frame membrane modules for industrial carbon capture

In this work, we study the application of membrane-based separation systems for carbon capture, considering plate-and-frame membrane modules. The successful deployment of membrane CO 2 capture system relies on high-performing membranes as well as effective membrane modules that can fully exploit the developed membranes. A plate-and-frame membrane module is especially attractive for CO 2 capture from industrial flue gas due to its lower pressure drop compared to its counterparts such as spiral wound modules and hollow fiber modules. To design better plate-and-frame modules, we investigate their basic unit - a single membrane stack through a combination of computational modeling and experimental investigations. The modeling approach is based on Computational Fluid Dynamics (CFD) to represent a multiphysics problem, including the fluid flow and diffusion processes within a membrane module. We use experimental data collected under different operating conditions to validate the CFD model. Numerical results suggest a good agreement between experiments and model outputs for the CO 2 recovery, CO 2 mole fraction in the retentate and permeate, and stage-cut. The CFD model is able to predict accurately the flow behavior, providing valuable insights on the effects of fluid dynamics on mass transfer of CO 2 . We also carry out a sensitivity analysis to identify the effect of key parameters on the CO 2 recovery and the CO 2 purity of the outlet streams.

CFD simulation↗
DOE OSTI2022

Diffuse back-illumination temperature imaging (DBI-TI), a novel soot thermometry technique

To meet stringent emissions regulations on soot emissions, it is critical to further advance the fundamental understanding of in-cylinder soot formation and oxidation processes. Among several optical techniques for soot quantification, diffuse back-illumination extinction imaging (DBI-EI) has recently gained traction mainly due to its ability to compensate for beam steering, which if not addressed, can cause unacceptably high measurement uncertainty. Until now, DBI-EI has only been used to measure the amount of soot along the line of sight, and here, we extend the capabilities of a DBI-EI setup to also measure in-cylinder soot temperature. This proof of concept of diffuse back-illumination temperature imaging (DBI-TI) as a soot thermometry technique is presented by implementing DBI-TI in a single cylinder, heavy-duty, optical diesel engine to provide 2-D line-of-sight integrated soot temperature maps. The potential of DBI-TI to be an accurate thermometry technique for use in optical engines is analyzed. The achievable accuracy is due in part to simultaneous measurement of the soot extinction, which circumvents the uncertainty in dispersion coefficients that depend on the optical properties of soot and the wavelength of light utilized. Analysis shows that DBI-TI provides temperature estimates that are closer to the mass-averaged soot temperature when compared to other thermometry techniques that are more sensitive to soot temperature closer to the detector. Furthermore, uncertainty analysis and Monte Carlo (MC) simulations provide estimates of the temperature measurement errors associated with this technique. The MC simulations reveal that for the light intensities and optical densities encountered in these experiments, the accuracy of the DBI-TI technique is comparable or even better than other established optical thermometry techniques. Thus, DBI-TI promises to be an easily implementable extension to the existing DBI-EI technique, thereby extending its ability to provide comprehensive line-of-sight integrated information on soot.

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

Core-shell and egg-shell zeolite catalysts for enhanced hydrocarbon processing

Developing structure-performance relationships with the underlying goal of optimizing known zeolite catalysts involves the manipulation of their physicochemical properties. Here, we systematically assessed the impact of mesoscopic gradients in acid site concentration, which has generally received little attention in the design of zeolite catalysts for hydrocarbon upgrading. A series of core–shell MEL-type zeolites were synthesized with catalytically active ZSM-11 cores and passivated silicalite-2 shells of varying thickness. Our findings revealed that ZSM-11@silicalite-2 particles with ultrathin shells (<10 nm) have enhanced mass transport, characteristic of relatively smaller particles, compared to the corresponding ZSM-11 core. Additionally, catalytic testing using the methanol-to-hydrocarbon (MTH) reaction showed that core–shell zeolites exhibit longer lifetimes, higher total turnovers, and an unexpected promotion of the aromatic cycle in the hydrocarbon pool mechanism. Time-resolved acid titration of core and core–shell catalysts confirmed that the siliceous shell introduces a hydrophobic exterior that impacts molecular diffusion. In comparison, prepared MFI core-shells (ZSM-5@silicalite-1) showed similar enhancement in catalyst performance. Moreover, we prepared egg-shell configurations of each zeolite, silicalite-2@ZSM-11 and silicalite-1@ZSM-5, comprised of an inert core and catalytically active shell. This inverse design of the egg-shell created pseudo nanosheets with total turnovers that were markedly higher than their homogeneous counterparts. Collectively, this study demonstrated that mesoscopic gradients in acid concentration via the design of core–shell and egg-shell zeolites significantly improve catalyst performance over conventional analogues for hydrocarbon upgrading.

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

Characterizing the impact of finite matrix block size on conservative particle transport through three-dimensional fracture networks

Mass transfer of solutes between fractures and the surrounding rock matrix exerts a noticeable signature on the tail of travel time distributions. When the width of the matrix is assumed to be infinite and advective transport through the fracture is sufficiently fast, the tails of the travel time distributions exhibit a classically expected slope of ψ(t) ∝ t -3/2 . However, studies have yet to characterize how solute transfer between fractures via diffusion through finite matrix blocks influences the tail’s slope in three-dimensional fractured media. Here, in this study, we assess the impact of finite matrix block size on breakthrough curve shape at different spatio-temporal scales by con ducting particle tracking simulations in three-dimensional discrete fracture networks. We consider a variety of hydrodynamic and geostructural proper ties to determine their relative impact on the resulting travel time distributions. We observe that the impact of matrix diffusion through a finite block on travel time distributions is similar to that of an infinite matrix block when the fracture spacing is sufficiently large, matrix diffusion is relatively weak, or transport is considered at an early control plane distance. We observe that the converse of these conditions, results in deviations from the classical ψ(t) ∝ t -3/2 scaling. These results provide a first step toward developing a metric to assess when finite block size effects are expected to significantly influence transport.

58 GEOSCIENCES↗