Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Power Law”

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 343 records · Page 19

Nonlinear Shear Rheology of Entangled Polymer Rings

Steady-state shear viscosity η($\bar{γ}$) of unconcatenated ring polymer melts as a function of the shear rate $\bar{γ}$ is studied by a combination of experiments, simulations, and theory. Experiments using polystyrenes with Z ≈ 5 and Z ≈ 11 entanglements indicate weaker shear thinning for rings compared to linear polymers exhibiting power law scaling of shear viscosity η ~ $\bar{γ}$ –0.56 ± 0.02 , independent of chain length, for Weissenberg numbers up to about 10 2 . Nonequilibrium molecular dynamics simulations using the bead-spring model reveal a similar behavior with η ~ $\bar{γ}$ -0.57 ± 0.08 for 4 ≤ Z ≤ 57. Viscosity decreases with chain length for high $\bar{γ}$. In our experiments, we see the onset of this regime, and in simulations, which we extended to Wi ~ 10 4 , the nonuniversality is fully developed. In addition to a naive scaling theory yielding for the universal regime η ~ $\bar{γ}$ –0.57 , we developed a novel shear slit model explaining many details of observed conformations and dynamics as well as the chain length-dependent behavior of viscosity at large $\bar{γ}$. The signature feature of the model is the presence of two distinct length scales: the size of tension blobs and much larger thickness of a shear slit in which rings are self-consistently confined in the velocity gradient direction and which is dictated by the size of a chain section with relaxation time 1/$\bar{γ}$. These two length scales control the two normal stress differences. In this model, the chain length-dependent onset of nonuniversal behavior is set by tension blobs becoming as small as about one Kuhn segment. This model explains the approximate applicability of the Cox–Merz rule for ring polymers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Yielding Behavior of Bottlebrush and Linear Block Copolymers

Block copolymers can exhibit a pronounced yield stress, but the impact of molecular architecture, chemistry, and self-assembly on macroscopic rheology remains poorly understood. Here, we study the linear-viscoelastic and yield-stress fluid behavior of two architectures—bottlebrush copolymers (with statistical or blocky sequences) and linear diblocks—that self-assemble into body-centered cubic (BCC) spheres and hexagonally close-packed cylinders (HEX). The dynamic properties of these polymers were probed by oscillatory frequency and amplitude sweeps at temperatures well below the order–disorder transition (TODT) to furnish insights into the yielding transition. All BCC-forming polymers have a similar signature of yielding: smaller yield strains (γ yBCC ≈ 0.053 < γ y,HEX ≈ 0.18), sharper solid–liquid transitions, and better reversibility than HEX. Statistical bottlebrushes show the most frequency-independent structural modulus (G 0 ) and no signs of defect relaxation. A simple power-law relationship captures the dependence of the normalized structural modulus (G 0 /RT) on the inter-micelle distance (d) across different architectures and morphologies [G 0 /(RT) = 1.31 × 10 4 (nm 2.6 mol/m 3 ) d –2.6 ]. These studies establish quantitative structure–property relationships that are relevant in contemporary applications, for example, extrusion-based 3D printing.

36 MATERIALS SCIENCE↗

Visualization of Multifractal Superconductivity in a Two-Dimensional Transition Metal Dichalcogenide in the Weak-Disorder Regime

Eigenstate multifractality is a distinctive feature of noninteracting disordered metals close to a metal–insulator transition, whose properties are expected to extend to superconductivity. While multifractality in three dimensions (3D) only develops near the critical point for specific strong-disorder strengths, multifractality in 2D systems is expected to be observable even for weak disorder. Here we provide evidence for multifractal features in the superconducting state of an intrinsic, weakly disordered single-layer NbSe 2 by means of low-temperature scanning tunneling microscopy/spectroscopy. The superconducting gap, characterized by its width, depth, and coherence peaks’ amplitude, shows a characteristic spatial modulation coincident with the periodicity of the quasiparticle interference pattern. The strong spatial inhomogeneity of the superconducting gap width, proportional to the local order parameter in the weak-disorder regime, follows a log-normal statistical distribution as well as a power-law decay of the two-point correlation function, in agreement with our theoretical model. Furthermore, the experimental singularity spectrum f(α) shows anomalous scaling behavior typical from 2D weakly disordered systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Buckling Spectra of Nanoparticle Surfactant Assemblies

Fine control over the mechanical properties of thin sheets underpins transcytosis, cell shape, and morphogenesis. Applying these principles to artificial, liquid-based systems has led to reconfigurable materials for soft robotics, actuation, and chemical synthesis. However, progress is limited by a lack of synthetic two-dimensional membranes that exhibit tunable mechanical properties over a comparable range to that seen in nature. In this study, we show that the bending modulus, B , of thin assemblies of nanoparticle surfactants (NPSs) at the oil–water interface can be varied continuously from sub- k B T to 10 6 k B T , by varying the ligands and particles that comprise the NPS. We find extensive departure from continuum behavior, including enormous mechanical anisotropy and a power law relation between B and the buckling spectrum width. Our findings provide a platform for shape-changing liquid devices and motivate new theories for the description of thin-film wrinkling.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Understanding kinetics of defect annihilation in chemoepitaxy-directed self-assembly.

Directed self-assembly (DSA) of block copolymers (BCP) has attracted considerable interest from the semiconductor industry because it can achieve semiconductor-relevant structures with a relatively simple process and low cost. However, the self-assembling structures can become kinetically trapped into defective states, which greatly impedes the implementation of DSA in high-volume manufacturing. Understanding the kinetics of defect annihilation is crucial to optimizing the process and eventually eliminating defects in DSA. Such kinetic experiments, however, are not commonly available in academic laboratories. To address this challenge, we perform a kinetic study of chemoepitaxy DSA in a 300 mm wafer fab, where the complete defectivity information at various annealing conditions can be readily captured. Through extensive statistical analysis, we reveal the statistical model of defect annihilation in DSA for the first time. The annihilation kinetics can be well described by a power law model, indicating that all dislocations can be removed by sufficiently long annealing time. We further develop image analysis algorithms to analyze the distribution of dislocation size and configurations and discover that the distribution stays relatively constant over time. The defect distribution is determined by the role of the guiding stripe, which is found to stabilize the defects. Although this study is based on polystyrene-b-poly(methyl methacrylate) (PS-b-PMMA), we anticipate that these findings can be readily applied to other BCP platforms as well.

block copolymer↗

Enhanced Interfacial Bonding of Graft Copolymers

To understand how thermoplastic welding strength can be tuned through chemical modifications and macromolecular topology, we combined coarse-grained molecular dynamics (MD) simulations with experimental validation. Our simulations examined the diffusion dynamics of both linear and graft polymers across representative interfaces, revealing that diffusion-controlled interdigitation follows a power law, with the exponent decreasing from 0.34 to 0.11 as grafting density increases from 7.5 to 196% (with side chains grafted to both sides of a monomer unit). The addition of side chains enhances welding efficiency, as dense bottlebrush polymers with high grafting density reach maximum rupture strength faster than linear polymers. However, their saturated rupture strength is lower. This observation is subsequently corroborated by experimental lap-shear tests comparing linear polyethylene with octene grafted polyethylene elastomers. Our MD simulations show that unlike linear polymers, where backbone entanglements dominate, the grafted side chains introduce mechanisms in addition to entanglement dilution. The rapid interdigitation of side chains creates a dense mesh of entropic van der Waals contacts, which can also enhance the film welding. Furthermore, our MD simulations reveal a brittle rupture behavior in linear and comb-like (mildly grafted) polymers, while bottlebrush (densely grafted) polymers display elastomeric behavior with a pronounced stress plateau prior to fracture. Our simulations deconvolute the influence of polymer topology on deformation behavior. The rate of polymer deformation becomes lower than the applied strain rate prior to rupture, and the onset of this deviation is progressively delayed from linear to bottlebrush polymers. This trend highlights the critical role of molecular architecture in governing the mechanical response. In conclusion, these results provide deeper insight into the underlying welding mechanisms of topological polymers and present a potential approach for mitigating the interface anisotropy that is inherent in advanced manufacturing techniques such as fused filament fabrication.

graft copolymers↗

Structural Dynamics in UV Curable Resins Resolved by In Situ 3D Printing X-ray Photon Correlation Spectroscopy

Additive manufacturing (AM) is a promising technique to rapidly produce polymeric materials into complex 3-dimensional (3D) geometries. While AM is widespread and relevant for a range of applications, implementation in industry has outpaced our fundamental understanding of polymer dynamics and structure development during the printing process. Characterization and quantification of such dynamics is necessary to optimize final material properties and design future materials and processes for 3D printing. Here, we utilize X-ray photon correlation spectroscopy (XPCS) to measure spatial and time-resolved, out-of-equilibrium dynamics during direct ink write (DIW) 3D printing. Specifically, we investigate the progression of structural dynamics in a dual cure (UV/thermal) nanocomposite during and directly after printing. As the filament is printed and cured in situ, the relaxation processes of the cross-linking network are measured through the dynamics of inorganic filler particles. The characteristic relaxation time of the dynamics is calculated through the intensity–intensity autocorrelation function g 2 and directly correlated to the printing process parameters, such as printhead velocity and UV light intensity. The time-resolved evolution of nanoscale dynamics follows a power-law dependence as the filament is cured. Additionally, bulk rheological characterizations reveal the macroscopic solidification of the resin, providing correlation of material properties across a wide range of length and time scales. Overall, the measurement of multiscale, out-of-equilibrium dynamics provides insight into the development of structure in polymer nanocomposite filaments during 3D printing and is used to further understand the influence of such parameters on the AM process.

36 MATERIALS SCIENCE↗

Emergent Above-Gap Photoluminescence in Molecularly Engineered Hybrid Bilayer Crystals

Bilayer crystals, built by stacking two-dimensional (2D) covalent monolayers, give rise to coupled excitonic states whose properties are constrained by fixed lattice symmetry and orientation. Replacing one covalent monolayer with a 2D molecular crystal––held together by non-covalent forces––overcomes this limitation, as molecular functional groups afford tunable in-plane lattice geometry, intermolecular spacing, and interlayer coupling, providing a powerful knob for exciton engineering. Here, in this study, we report four-atom-thick hybrid bilayer crystals (HBCs) synthesized by directly growing single-crystalline PDI molecular crystal atop WS 2 monolayers, which exhibit a robust photoluminescence (PL) peak 120 meV above the WS 2 optical band gap alongside a belowgap emission. Both peaks display strong polarization anisotropy—nearing unity for the above-gap emission—and maintain a perfectly linear power-law dependence up to an excitation density of ~10 7 mW/cm 2 , indicative of coexisting localized and delocalized excitonic states. Substituting PDI with a PTCDA monolayer on WS 2 fully quenches PL, demonstrating molecular control over excitonic emission. Lattice scale ab initio GW and GW-BSE calculations reveal a significantly hybridized bilayer band structure in PDI/WS 2 that supports interlayer excitonic species both above and below the WS 2 gap with strong polarization anisotropy, in excellent agreement with experiment. Our work introduces a molecule-based bilayer platform for the bottom-up design and control of excitonic phenomena in atomically thin optoelectronic and quantum materials.

2D materials↗

Distribution Cutoff for Clusters near the Gel Point

The mechanical and dynamic properties of developing networks near the gel point are susceptible to the distribution of clusters coexisting with percolating networks. The distribution of cluster numbers follows a broad power law, wrapped by a cutoff function that rapidly decays at a characteristic size. The form of the cutoff function has been speculated based on known results from lattice percolation and, in certain cases, solved. We obtained this cutoff function from simulated dynamic clusters of polymeric precursor chains using a hybrid Monte Carlo algorithm. The results obtained from three different precursor chain lengths are consistent with each other and are consistent with the expectation from lattice percolation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of Nanocrystalline Graphite from Lignin Sources

Carbon composites are attractive to a variety of high-impact applications, such as carbon fibers, batteries, and vehicle parts, due to their multifunctional properties. The properties of carbon are highly dependent on the allotrope the carbon takes and the functionality, impurities, and defects contained within the structure. The increase in demand for sustainable carbon sources in energy storage devices motivates interest in understanding synthesis parameters of lignin value-added products. Also, as the dependence on oil for fuel decreases, alternative sources for carbon in many applications will be needed. In this work, the thermochemical conversion of lignin powders from different feedstocks was evaluated via small and wide-angle X-ray scattering techniques to resolve the amorphous, disordered, and crystalline domains present in the lignin carbons. Scattering analyses indicated an evolution of hierarchical structures along with an increase in ordered domains as a function of carbonization temperature. Qualitative and quantitative methods were used to describe isotropic scattering intensity profiles at multiple length scales. The use of power law models in the mesoscopic region served as the basis to describe morphological changes related to structural features, for example, graphene stacking, degree of roughness, and surface fractals. Kraft softwood and switchgrass produced carbon powder with the most crystalline domains and the least surface roughness. Softwoods reached the highest degree of crystallinity followed by switchgrass samples and had less variability in particle sizes. These results suggest lignin carbons extracted from softwoods and switchgrass are viable substitutes for graphite. Interpretation of X-ray scattering data from lignin carbon powders elucidates feedstock- and processing-dependent morphological features across multiple length scales providing a straightforward framework to evaluate the feasibility of leveraging lignin carbons for producing tunable application-specific materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evidence for Lignin–Carbohydrate Complexes from Studies of Transgenic Switchgrass and a Model Lignin–Pectin Composite

Lignin–carbohydrate complexes (LCCs) form through interactions of lignin with plant cell wall polysaccharides and are thought to be a significant source of biomass recalcitrance. In this work, we investigated LCCs formed between lignin and pectin homogalacturonan (HG). The structural changes in HG deficient transgenic switchgrass (GAUT4-knockdown, GAUT4-KD) after hot water pretreatment were compared to wild-type plants using small-angle neutron scattering (SANS), which showed that there were ~2.2-fold more lignin aggregates in GAUT4-KD biomass compared to the wild type. This demonstrated that decreased pectin resulted in more lignin redistribution and suggested that interactions between lignin and HG restrict lignin mobility in plant cell walls. To better understand the types of interactions between lignin and pectin, a model composite was prepared by polymerizing either protiated or partially deuterated coniferyl alcohol to form a dehydrogenation polymer (DHP) in the presence of HG. Small-angle X-ray scattering (SAXS) showed that the DHP and HG form a highly interconnected network structure that is not observed in a physical mixture of the individual polymers. Contrast matching SANS revealed the structure of DHP and HG in the composite and showed that the HG forms a swollen interconnected polymer network (power-law exponent, P = 1.5) interspersed with DHP particles (radius of gyration, Rg, 264 Å) that are composed of solvent-accessible DHP polymers (P = 2.3). Fourier transform infrared spectroscopy showed a unique ester absorption band in the DHP/HG composites. Solid-state nuclear magnetic resonance (NMR) analysis also supports interactions between DHP and HG. Overall, this study provides new insights into the relationship between primary and secondary cell wall polymers during cell wall synthesis and how LCCs formed between pectin and lignin could represent a previously unrecognized source of biomass recalcitrance. This knowledge may help develop new approaches to modulate cell wall properties to improve biofuel and bioproduct production.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ultralow-Loading Pt/Zn Hybrid Cluster in Zeolite HZSM-5 for Efficient Dehydroaromatization

Minimizing Pt loading without sacrificing catalytic performance is critical, particularly for designing cost-efficient hydrocarbon transformation catalysts. Here, we show that ultralow-loading (0.001–0.05 wt %) Pt- and Zn-functionalized HZSM-5 catalysts, prepared through simple ion exchange and impregnation, are highly active and stable for light alkane dehydroaromatization (DHA). The specific activity of benzene, toluene, and xylene is up to 8.2 mol/g Pt /min (or 1592 min –1 ) over the 0.001 wt % Pt–Zn 2 /HZSM-5 catalyst during ethane DHA at 550 °C under atmospheric pressure. Additionally, such bimetallic Pt x –Zn y /HZSM-5 catalysts are highly stable in contrast to the monometallic Pt/HZSM-5 catalysts. The rate constant of deactivation (k deactiv ), according to the first-order generalized power law equation model, for the bimetallic catalysts is up to 120 times lower than that of the monometallic counterparts, depending on the Pt loading. In conclusion, this breakthrough is achieved through the formation of the [Pt 1 –Zn n ] δ+ hybrid cluster, instead of Pt 0 cluster–proton adducts, in the micropores of the ZSM-5 zeolite.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Model for Turbulence Spectra in the Equilibrium Range of the Stable Atmospheric Boundary Layer

Stratification can cause turbulence spectra to deviate from Kolmogorov's isotropic -5/3 power law scaling in the universal equilibrium range at high Reynolds numbers. However, a consensus has not been reached with regard to the exact shape of the spectra. Here we propose a shape of the turbulent kinetic energy and temperature spectra in horizontal wavenumber for the equilibrium range that consists of three regimes at small Froude number: the buoyancy subrange, a transition region, and the isotropic inertial subrange through dimensional analysis and substantial revision of previous theoretical approximation. These spectral regimes are confirmed by various observations in the atmospheric boundary layer. The representation of the transition region in direct numerical simulations will require large-scale separation between the Dougherty-Ozmidov scale and the Kolmogorov scale for strongly stratified turbulence at high Reynolds numbers, which is still challenging computationally. In addition, we suggest that the failure of Monin-Obukhov similarity theory in the very stable atmospheric boundary layer is due to the fact that it does not consider the buoyancy scale that characterizes the transition region.

54 ENVIRONMENTAL SCIENCES↗

Adaptive Multirate Mass Transfer (aMMT) Model: A New Approach to Upscale Regional‐Scale Transport Under Transient Flow Conditions

Abstract The long‐term evaluation of regional‐scale groundwater quality needs efficient upscaling methods for transient flow. Upscaling techniques, such as the Multirate Mass Transfer (MRMT) method with constant upscaling parameters, have been used for transport with steady‐state flow, yet the upscaling parameters (i.e., rate coefficients) may be time dependent. This study proposed and validated an adaptive MRMT (aMMT) method by allowing the mass transfer coefficients in MRMT to change with the flow field. Advective‐dispersive contaminant transport simulated in a 3‐D heterogeneous medium was used as a reference solution. Equivalent transport under homogeneous flow conditions was evaluated by applying the MRMT and aMMT models for upscaling. The relationship between mass transfer coefficients and flow rates was fitted under steady‐state flow driven by various hydraulic gradients. A power law relationship was obtained, which was then used to update the mass transfer coefficients in each stress period under transient flow conditions in the aMMT method. Results indicated that for advection‐dominated transport, both the MRMT and aMMT methods can upscale the anomalous transport dynamics affected by subgrid heterogeneity under transient flow conditions. Whereas for diffusion‐dominated systems, the MRMT model failed to capture the tails of tracer breakthrough curves after the boundary condition changed, but the results from the aMMT model were significantly improved. However, if the overall flow direction changed, both MRMT and aMMT failed to represent the breakthrough curve tail generated by the heterogeneous system. The results point toward a promising path for upscaling transport in complex aquifers with transient flow.

Guo, Zhilin↗

River Dynamics Control Transit Time Distributions and Biogeochemical Reactions in a Dam-Regulated River Corridor

Residence Time Distributions (RTDs) exert an important control on biogeochemical translation in the hydrological systems. Previous tracer studies have revealed that RTDs often followed time-invariant exponential, lognormal, or heavy-tailed RTDs that have power-law behavior for long tails in headwater or low-order streams. However, there is a recognition that RTDs can be more complicated and time-variable in response to dynamic hydrological forcing. Here, we use particle tracking to estimate RTDs along the Hanford Reach of the Columbia River and to quantify the influences of river stage fluctuations. Particle tracking is conducted using the velocity field from high-resolution 3D groundwater flow simulations. The effects of hydrological forcing on the residence time distribution were evaluated by varying river flow boundary conditions and releasing particles in different time windows. Our results revealed that dynamic stage fluctuations created rapidly changing losing-gaining conditions in the river and led to highly transient RTDs, which contributes to multiple modes of RTDs. Dam-induced high-frequency (sub-daily) flow variation contributes to the short-time (sub-daily) component of the RTDs. Deviation of the reactant consumption under the single-mode assumption compared to the multimodal RTDs is relatively small (~5%) and appears when the Damköhler number was close to one.More specifically, high-frequency (daily) dam-induced stage variations have higher impacts on biogeochemical reactions with faster reaction rate, and can potentially increase the reactant consumption by 26.68%. Based on these findings, we suggest that current river basin models could be improved by including bank storage and more complex RTDs influenced by both short and long term river stage fluctuations.

54 ENVIRONMENTAL SCIENCES↗

Hyporheic Exchange Driven by Submerged Rigid Vegetation: A Modeling Study

Flow-vegetation interaction affects fluid flow hydraulics and associated material transport in river corridors. Concomitant changes in pressure within the flow field due to the presence of vegetation may act as a driver for the formation of hyporheic flow across the sediment-water interface. This potentially important process, however, has yet to be studied. In order to investigate vegetation-induced hyporheic exchange, a series of numerical models of interlinked surface-subsurface flow modified by plant stems was conducted. Periodically staggered plant stem arrays on a flat sediment bed were considered within a coupled multiphysics computational fluid dynamics approach. Plants were idealized as rigid cylinders and arranged in different streamwise and spanwise spacing distances. Each vegetation array was then subjected to a broad range of flow Reynolds Numbers (Re). The results showed that hyporheic flow occurs in all conditions with the presence of vegetation. The vegetation-induced hyporheic flux is found to be a function of Re via a power law. The flux increases with interstem space until the space reaches the distance that rigid stems no longer affect the flow structures in the vicinity of each other. Larger intervegetation distances lead to a larger hyporheic zone. A direct comparison with bedforminduced hyporheic flow showed that vegetation can induce higher hyporheic flux through relatively shallower exchange zones. The results of all the simulations were synthesized into predictive models for hyporheic flux, bulk residence time and exchange depth based on drag coefficient, vegetation density, and Reynolds Number.

54 ENVIRONMENTAL SCIENCES↗

Understanding the Spatial Organization of Simultaneous Heavy Precipitation Events Over the Conterminous United States

We introduce the idea of simultaneous heavy precipitation events (SHPEs) to understand whether extreme precipitation has a spatial organization manifested as specified tracks or contiguous fields with inherent scaling relationships. For this purpose, we created a database of SHPEs using ground-based precipitation observations recorded by the daily Global Historical Climatology Network across the conterminous United States during 1900–2014. SHPEs are examined for their seasonality, spatial manifestation, orientation, and areal extent. We quantified the spatial distribution of the centroids and principal axes of SHPEs and their quasi-elliptical manifestations, azimuthal orientations, and areal extents on the ground. Four seasons, December-January-February (DJF), March-April-May (MAM),June-July-August (JJA), and September-October-November (SON) are considered to examine the spatial patterns and associated large-scale atmospheric circulations. Results indicate that there are 54, 58, 103, and204 SHPEs in DJF, MAM, JJA, and SON seasons, and their longest stretches range on average between 650and 1,600, between 850 and 1,500, between 950 and 1,550, and between 750 and 1,450 km, respectively. SHPEs in the DJF, MAM, and JJA seasons occur mostly over the Pacific coast and central and midwestern United States, respectively. The atmospheric circulation patterns and mechanisms of precipitable water vapor and moisture transport in the atmosphere are also discussed in relation to these SHPEs. Power laws explain SHPEs' underlying area scaling behavior in all the four seasons, with stronger evidence in DJF and MAM. A seasonal spatial risk model is developed to predict the likelihood of SHPEs. Quantifying the characteristics of SHPEs and modeling their footprints can improve the projections of flood risk and understanding of damages to interconnected infrastructure systems.

54 ENVIRONMENTAL SCIENCES↗

Flow Channeling in Fracture Networks: Characterizing the Effect of Density on Preferential Flow Path Formation

Flow channelization is a commonly observed phenomenon in fractured subsurface media where the flow of fluids is restricted primarily to highly transmissive fracture networks surrounded by a low-permeability rock matrix. The multiscale structural heterogeneity of these networks results in multiscale flow channelization where preferential flow paths form at length scales ranging from the entire system down to the subfracture size. We present an analysis of how one of the largest scales in fractured media, the network density, influences the degree of flow channeling that occurs using an ensemble of semigeneric three-dimensional discrete fracture network (DFN) simulations. We construct 10 DFNs, whose fracture lengths follow a power law distribution, at four densities for a total of 40 networks. We characterize their structure in terms of the network topology and geometry. Eulerian and Lagrangian observations of the steady-state flow fields obtained within the networks are used to quantify the degree of flow channelization at the network scale. We introduce a measure for the importance-ranking/hierarchy of different flow paths in the network using graph-based analysis of Lagrangian transport by which the degree of flow channeling between networks is compared. These flow observations are then linked to the structural properties of the networks. In general, network-scale flow channeling decreases as the network density increases. However, at low densities, there is more uniform flow within the entire connected network than in high-density networks. We also demonstrate how standard transport observables can be used to infer the degree of flow channelization occurring within a fracture network.

58 GEOSCIENCES↗