Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flow cell”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Triggering C-N Coupling on Metal Oxide Nanocomposite for the Electrochemical Reduction of CO 2 and NO x ⁻ to Formamide

The co-electroreduction of CO 2 and NO x ⁻ (NO 3 ⁻ /NO 2 ⁻ ) to generate formamide (HCONH 2 ) offers an opportunity for downstream chemical and polymer manufacturing decarbonization; however, significant challenges lie in the C-N coupling and the associated low product selectivity. Herein, p-block metal oxides are incorporated in copper oxides to provide more accessible active sites for reactant adsorption and activation, tuning the reaction selectivity toward the formamide production. Through in situ Raman and synchrotron-based infrared spectroscopy measurements, C-N bond formation is demonstrated in real-time with the CuO x /BiO x catalyst, where the C-N bond is detected via a *CHO and *NH 2 intermediates formation, in agreement with the density functional theory calculations. When tested in a flow electrolyzer, a formamide yield rate of 134 ± 11 mmol h -1 g cat -1 is reported, the first report of co-electroreduction of CO 2 and NO x ⁻ to formamide beyond conventional H-cell measurements. These new insights on the C-N coupling mechanisms and scale-up capability provide directions for further development of electrocatalysts for the formamide production.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fetch-trap Pairs: Exploring definition of carbon storage prospects to increase capacity and flexibility in areas with competing uses

As CCS grows beyond isolated first-mover projects, the question of how to identify potential storage sites (prospects) is key to siting projects amidst the constraints of neighboring storage projects, hydrocarbon fields and other considerations. To date, prospects have been defined either by the extent of a buoyant trap or by the predicted maximum extent of the CO 2 or pressure plume. Both approaches have worked, but they may struggle in places like the Gulf of Mexico where complex geologic structure can invalidate simple assumptions of radial plume spread and buoyant traps are commonly unattractive due to numerous legacy well penetrations. Here, we propose identifying prospective storage sites by fetch area (drainage cell) rather than by buoyant closure alone. Doing so offers 1) greater freedom in siting injectors to avoid surface and subsurface constraints, including legacy wells that tend to cluster on structural high; and 2) a coherent flow regime in which buoyancy drives all injected CO 2 toward a common high. This strategy gives space to dissipate injection pressure and minimize the number of legacy wells needing review. It may also offer 1) better injectivity by bringing synclines and associated channel axes into play; and 2) potentially improved capacity by tapping a larger rock volume and taking advantage of migration losses.

58 GEOSCIENCES↗
DOE OSTI2022

Novel Solver Algorithms for Nearly Singular Linear Systems Arising in Combustion Modelling

Direct Numerical Simulations of realistic combustion devices are extremely challenging due to the wide separation of scales in the simulation, for example an internal combustion (IC) engine chamber, and the flame thickness of a high-pressure flame. The PeleLMeX solver uses adaptive mesh refinement (AMR) to evolve multi-species reacting flows in the low Mach number limit at the Exascale and relies on an embedded boundary (EB) approach to represent complex geometries. In that framework, the EB geometries often give rise to very small cut-cells along the boundary, which translate into extreme ill-conditioning of the pressure-projection, with eigenvalues that span 15-16 orders of magnitude. In this talk, we focus on the case of a typical IC piston bowl geometry for which we present on a novel approach towards solving these nearly singular linear systems with ILU-based, C-AMG smoothers on massively parallel architectures. In particular, we use scaling and equilibration algorithms to handle the non-normality of the upper triangular factors. This enables us to approximate the highly sequential triangular solve algorithm, embedded in the AMG smoothing-solve phase, with Jacobi iterations. This approximation can be written as a convergent Neumann series whose terms are composed of highly parallel sparse matrix vector multiplications. The result is an algorithm that substantially decreases setup and solve time, compared to state-of-the-art, for these challenging linear systems.

combustion modelling↗
DOE OSTI2022

Flow and transport in three-dimensional discrete fracture matrix models using mimetic finite difference on a conforming multi-dimensional mesh

Here, we present a comprehensive workflow to simulate single-phase flow and transport in fractured porous media using the discrete fracture matrix approach. The workflow has three primary parts: (1) a method for conforming mesh generation of and around a three-dimensional fracture network, (2) the discretization of the governing equations using a second-order mimetic finite difference method, and (3) implementation of numerical methods for high-performance computing environments. A method to create a conforming Delaunay tetrahedralization of the volume surrounding the fracture network, where the triangular cells of the fracture mesh are faces in the volume mesh, that addresses pathological cases which commonly arise and degrade mesh quality is also provided. Our open-source subsurface simulator uses a hierarchy of process kernels (one kernel per physical process) that allows for both strong and weak coupling of the fracture and matrix domains. We provide verification tests based on analytic solutions for flow and transport, as well as numerical convergence. We also provide multiple expositions of the method in complex fracture networks. In the first example, we demonstrate that the method is robust by considering two scenarios where the fracture network acts as a barrier to flow, as the primary pathway, or offers the same resistance as the surrounding matrix. In the second test, flow and transport through a three-dimensional stochastically generated network containing 257 fractures is presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Physics-constrained deep neural network method for estimating parameters in a redox flow battery

Here, in this paper, we present a physics-constrained deep neural network (PCDNN) method for parameter estimation in the zero-dimensional (0D) model of the vanadium redox flow battery (VRFB). In this approach, we use deep neural networks to approximate the model parameters as functions of the operating conditions. This method allows the integration of VRFB computational models as the physical constraints in the parameter learning process, leading to enhanced accuracy of parameter estimation and cell voltage prediction. Using an experimental dataset, we demonstrate that the PCDNN method can estimate model parameters for a range of operating conditions and improve the 0D model prediction of voltage compared to the 0D model prediction with constant operation-condition-independent parameters estimated with traditional inverse methods. We also demonstrate that the PCDNN approach has an improved generalization ability for estimating parameter values for operating conditions not used in the training process.

25 ENERGY STORAGE↗
DOE OSTI2024

DECOVALEX-2023: Task F1 Final Report

DECOVALEX-2023 Task F is a comparison of models and methods for post-closure performance assessment (PA) of a deep geologic repository for radioactive waste. The general aims of Task F are to build confidence in the models, methods, and software used for PA and to stimulate additional research and development in PA methodologies. The task objectives are to motivate development of PA modelling skills and capabilities, to examine the influence of model choices on calculated repository performance, and to compare the uncertainties introduced by model choices to other sources of uncertainty. Task F involves no actual experiment or site. It is a PA modelling exercise that requires the conceptual development of hypothetical repository designs and geologic settings. Because three of the teams were interested in salt and the rest of the teams were interested in crystalline rock, Task F was split into two branches: Task F1 for crystalline rock and Task F2 for salt. This report is for Task F1, crystalline rock. Teams from seven countries (Canada, Czech Republic, Germany, Korea, Sweden, Taiwan, and United States) participated in Task F1. The teams worked together to define the features, events, and processes of the reference case repository and established a set of performance measures. In addition, they defined a set of benchmark problems designed to test and compare modelling capabilities for fracture flow and transport at different scales. The repository design and benchmark problems are documented in a Task Specification that evolved over time as the group honed the specifications. The benchmark problems verified that each team can aptly model flow and transport in fractured media in 1-, 2-, and 3-dimensions. Two general approaches were used for the 3-dimensional benchmarks: discrete fracture network (DFN) and equivalent continuous porous medium (ECPM). DFN modelling involves explicit meshing of each fracture while ECPM modelling aims to capture the effective porosity and directional permeability of each cell in a space-filling mesh as affected by intersecting fractures. In some models, a combination of the two is used, i.e., DFN for large known fractures and ECPM for the rest of the domain. Transport is solved by using either the advection-dispersion equation or particle tracking. Although some variation is observed among model breakthrough curves in the benchmark problems, there is strong agreement in breakthrough behaviour up to at least the 75 th percentile for all benchmarks. At the 90 th percentile, breakthrough results show larger differences, suggesting several models retain substantially higher fractions of tracer in regions of slower moving water. In addition to the flow and transport benchmarks, several teams completed the source term benchmark, verifying capabilities for modelling radionuclide decay and ingrowth, waste package breach, instant release fractions, fuel matrix degradation rates, and radionuclide solubility limitations. The reference case is conceptualized as a generic spent fuel repository at a depth of 450 m in fractured crystalline rock. The repository has 50 parallel backfilled drifts, each with 50 deposition holes 6 m apart. Each deposition hole contains a 4-PWR waste package and bentonite buffer. The rock domain is 5 km in length, 2 km in width, and 1 km in depth. It has 6 deterministic fractured deformation zones and a multitude of stochastic fractures. Teams generally used the ECPM approach for the entire rock or a hybrid approach in which the deterministic fracture zones are modelled with a DFN and the rest of the rock is modelled by ECPM. Of the reference case problems specified, only the results of the initial reference case problem are compared in this report. The initial problem focuses on transport from the deposition holes to the surface, i.e., it neglects waste package performance. Tracers are released at all waste package locations at time zero and tracked for their releases to the near field and ground surface. The water fluxes calculated at the ground surface entry and exit regions of the domain are similar for all models except for two that have considerably lower fluxes. For tracer transport, large differences are observed among models in the magnitude of tracer transported. Much of the difference appears to be due to how the repository is implemented and hence the different degrees of repository simplification. Models that exclude the drifts, buffer, and backfill from the domain tend to show greater release of tracers and radionuclides from the repository. The initial study presented here indicates that major differences in modelling important processes within the repository (e.g., diffusion through buffer and backfill) can produce broadly different release and transport results, especially when those processes are excluded. Even for the models that included all specified features, events, and processes, the results show significant differences and demonstrate the importance of examining multiple modelling approaches in performance assessment. The differences in results observed in this study are expected to motivate teams to either increase complexity in future versions of the reference case models or to improve methods to account for the effects of simplified features and processes. Either way, future improvements in these models are expected to produce results that more closely agree.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Numerical Analysis of Regular Material Point Method and its Application to Multiphase Flows

The material point method (MPM) is gaining wide popularity in engineering research to model and simulate complex multiphase flow dynamics. The method relies on solving the governing equations of motion and transport in a Lagrangian framework using particles also known as material points. The fluid and kinematic properties are stored on the material points while the spatial gradient calculation and temporal integration are performed on a background grid. This Lagrangian framework allows for large deformations, easy integration of constitutive models, and direct import of complex geometries as particles. However, despite their increasing popularity, very few studies have addressed the issues of numerical resolution and stability of MPM techniques. The presence of additional factors such as the number of material points-per-cell, the location of the material points, the CFL-like condition used in time update, and the grid shape functions also increase the complexity of the error analysis when compared to other finite element methods. In this presentation, we analyze the various forms of error incurred in the application of MPM to continuum mechanics and multiphase flows. The effect of the previously mentioned factors on the error dynamics is studied. The application of these principles to canonical and industrial problems is also presented.

high pressure reverse osmosis↗
DOE OSTI2020

Thermomechanical response of thickly tamped targets and diamond anvil cells under pulsed hard x-ray irradiation

In the laboratory study of extreme conditions of temperature and density, the exposure of matter to high intensity radiation sources has been of central importance. Here, we interrogate the performance of multi-layered targets in experiments involving high intensity, hard x-ray irradiation, motivated by the advent of extremely high brightness hard x-ray sources, such as free electron lasers and 4th-generation synchrotron facilities. Intense hard x-ray beams can deliver significant energy in targets having thick x-ray transparent layers (tampers) around samples of interest for the study of novel states of matter and materials’ dynamics. Heated-state lifetimes in such targets can approach the microsecond level, regardless of radiation pulse duration, enabling the exploration of conditions of local thermal and thermodynamic equilibrium at extreme temperature in solid density matter. The thermal and mechanical responses of such thick layered targets following x-ray heating, including hydrodynamic relaxation and heat flow on picosecond to millisecond timescales, are modeled using radiation hydrocode simulation, finite element analysis, and thermodynamic calculations. Assessing the potential for target survival over one or more exposures and resistance to damage arising from heating and resulting mechanical stresses, this study doubles as an investigation into the performance of diamond anvil high pressure cells under high x-ray fluences. Long used in conjunction with synchrotron x-ray radiation and high power optical lasers, the strong confinement afforded by such cells suggests novel applications at emerging high intensity x-ray facilities and new routes to studying thermodynamic equilibrium states of warm, very dense matter.

Meza-Galvez, J. (ORCID:000000032699849X)↗
DOE OSTI2020

Cross-field electron diffusion due to the coupling of drift-driven microinstabilities

In this paper, the nonlinear interaction between kinetic instabilities driven by multiple ion beams and magnetized electrons is investigated. Electron diffusion across magnetic field lines is enhanced by the coupling of plasma instabilities. Here, a two-dimensional collisionless particle-in-cell simulation is performed accounting for singly and doubly charged ions in a cross-field configuration. Consistent with prior linear kinetic theory analysis and observations from coherent Thomson scattering experiments, the present simulations identify an ion-ion two-stream instability due to multiply charged ions (flowing in the direction parallel to the applied electric field) which coexists with the electron cyclotron drift instability (propagating perpendicular to the applied electric field and parallel to the ExB drift). Small-scale fluctuations due to the coupling of these naturally driven kinetic modes are found to be a mechanism that can enhance cross-field electron transport and contribute to the broadening of the ion velocity distribution functions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Achieving high hydrogen productivities of 20 L/L-day via microbial electrolysis of corn stover fermentation products

Microbial electrolysis cells have the potential to generate renewable hydrogen from underused waste streams, however current devices have not reached very large productivity targets using real waste products at any scale. This study used a waste from fermented corn stover known as corn stover fermentation product to reach performance metrics that could be commercialized, if adequately scaled. Average current densities in MECs with mature biofilms reached 17.9 +/- 1.6 A/m 2 at an organic loading rate of 30 grams of chemical oxygen demand per liter of anode volume per day (g/L-day), reaching a maximum current density of 27.2 +/- 2.9 A/m 2 . Hydrogen productivities reached 20.2 +/- 2.0 L of H 2 per L anode volume per day (L/L-day). These represent the highest current densities and highest hydrogen productivities using a complex feedstock in a microbial electrolysis cell. Organic acids and sugars present in untreated substrate were converted at high percentages in MECs, with most above 90% conversion, at organic loading rates of 10, 20, and 30 g/L-day. The effect of periodic high liquid flow rates through the anode on performance was assessed. These tests, called pulsing, showed that hydrogen productivities and current densities increased most dramatically as flow was pulsed every hour at a baseline flow rate of 0.3 mL/min. These productivities show promise for high performance systems, if adequate scale up can be achieved.

09 BIOMASS FUELS↗
DOE OSTI2020

Particle interspacing effects on the mechanical behavior of a Fe–TiB 2 metal matrix composite using FFT-based mesoscopic field dislocation mechanics

This paper presents an application to metal matrix composites (MMCs) of an enhanced elasto-viscoplastic Fast Fourier Transform (EVP-FFT) formulation coupled with a phenomenological continuum Mesoscale Field Dislocation Mechanics (MFDM) theory. Contrary to conventional crystal plasticity, which only accounts for plastic flow and hardening induced by statistically stored dislocations (SSDs), MFDM-EVP-FFT also describes the evolution of polarized geometrically necessary dislocation (GND) density and its effect on both plastic flow and hardening. Numerical results for a Fe–TiB 2 MMC made of a ferrite matrix (α-Fe) and elastic ceramic particles (TiB 2 ) are presented. Full-field simulations are performed using synthetic periodic unit cells representative of the MMC, with single-crystalline and polycrystalline matrix, for different particle interspacing distances. A strong dependence of the predicted equivalent stress, cumulated plastic strain and GND density fields with particle interspacing distance is observed, in contrast with conventional crystal plasticity. Correlations between these mechanical fields and microstructural features, and their influence on local and global mechanical behavior are examined for the different MMC microstructures.

42 ENGINEERING↗
DOE OSTI2024

Cost Competitive Process of Battery Grade Oxides Preparation from Aged LIBS

The increasing demand for lithium-ion batteries (LIBs ) is driving development of advanced recycling and refining methods. In this project, new cathode active materials are prepared from recycled lithium battery materials and subsequently electrochemically evaluated in coin cells. In detail, scrap LIBs were mechanically processed into a metal rich black mass, then reductive acid leached into an aqueous metal solution, and finally high-purity metal hydroxide precursor materials were prepared by selective electrochemical flow precipitation. Closed-loop LIB recycling into active electrode materials will enable US manufacturers to break their reliance on foreign sources of critical materials. For example, electroextraction process used here has potential to significantly reduce chemical & water use as compared to traditional metallurgic techniques. To this end, electroextracted materials refined from used batteries were collected and processed to be tested as precursor cathode active material (pCam).During this project ,samples of high-purity recycled NMC hydroxide (~2kg each) were received from Nth Cycle. In the first received batch some impurities like cadmium, calcium, and sodium are present in quantities that would negatively impact battery performance. Nth Cycle indicated that the impurities such as cadmium are present due to poor battery separation and co-shredding of NiCd batteries with with LIB. Cadmium was not present in the two subsequent batches as Nth Cycle processed cadmium free black mass. Thus, effective recycled battery sorting is critical to closed-loop LIB recycling into active electrode materials.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

Cost Competitive Process of Battery Grade Oxides Preparation from Aged LIBS

The increasing demand for lithium-ion batteries (LIBs) is driving development of advanced recycling and refining methods. In this project, new cathode active materials are prepared from recycled lithium battery materials and subsequently electrochemically evaluated in coin cells. In detail, scrap LIBs were mechanically processed into a metal rich black mass, then reductive acid leached into an aqueous metal solution, and finally high-purity metal hydroxide precursor materials were prepared by selective electrochemical flow precipitation. Closed-loop LIB recycling into active electrode materials will enable US manufacturers to break their reliance on foreign sources of critical materials. For example, electroextraction process used here has potential to significantly reduce chemical & water use as compared to traditional metallurgic techniques. To this end, electroextracted materials refined from used batteries were collected and processed to be tested as precursor cathode active material (pCam). During this project , samples of high-purity recycled NMC hydroxide (~2 kg each) were received from Nth Cycle. In the first received batch some impurities like cadmium, calcium, and sodium are present in quantities that would negatively impact battery performance. Nth Cycle indicated that the impurities such as cadmium are present due to poor battery separation and co-shredding of NiCd batteries with with LIB. Cadmium was not present in the two subsequent batches as Nth Cycle processed cadmium free black mass. Thus, effective recycled battery sorting is critical to closed-loop LIB recycling into active electrode materials.

25 ENERGY STORAGE↗
DOE OSTI2026

Cost Competitive Process of Battery Grade Oxides Preparation from Aged LIBs

The increasing demand for lithium-ion batteries (LIBs) is driving development of advanced recycling and refining methods. In this project, new cathode active materials are prepared from recycled lithium battery materials and subsequently electrochemically evaluated in coin cells. In detail, scrap LIBs were mechanically processed into a metal rich black mass, then reductive acid leached into an aqueous metal solution, and finally high-purity metal hydroxide precursor materials were prepared by selective electrochemical flow precipitation. Closed-loop LIB recycling into active electrode materials will enable US manufacturers to break their reliance on foreign sources of critical materials. For example, electroextraction process used here has potential to significantly reduce chemical & water use as well as waste production as compared to traditional metallurgic techniques. To this end, electroextracted materials refined from used batteries were collected and processed to be tested as precursor cathode active material (pCAM). For this project, two samples of high purity recycled NMC hydroxide (~1 kg each) were received from N th Cycle’s Ohio demonstration facility. The NMC compositions of the two materials are similar, and the levels of impurities have been confirmed. Indeed, impurities like copper, boron and sodium are present in quantities that could negatively impact battery performance. However, some studies have demonstrated that the control of the quantity of elements like copper or boron may improve the electrochemistry properties of Lithium-NMC batteries. The principal objective is to determine the electrochemical performance of these 2 NMC hydroxide batches and clearly determine the effect of the impurities on the performance.

25 ENERGY STORAGE↗
DOE OSTI2024

Residual Dose Rates Analyses from Vent Lines Shield Block on the Target Cart at Spallation Neutron Source

The target carriage supports the mercury target vessel in the Spallation Neutron Source (SNS) in a cantilevered fashion and supplies the liquid mercury flowing to and from the target vessel. The target carriage is designed to move the target vessel into its target monolith position to receive the proton beam and then remove it horizontally back into an adjacent service cell, where the target vessel can be exchanged about 2-3 times per calendar year. At the time of the performed studies, the SNS mercury target can accept a proton beam at 1 GeV energy and 1.4 MW power; however, it will be upgraded to receive 1.3 GeV energy at 2 MW by the end of 2024 year [1]. The cart is practically a shield block on wheels to minimize fluxes and dose rates in the service cell downstream of the carriage. There are vent lines located inside a passage on top of the target carriage, and they are covered by a stepping-up steel shield block (vent lines shield block). Vent lines are permanently attached to the shield block. New vent lines with a different design are going to be installed. These lines are shielded by the new steel block to which the vent lines are attached. The estimate of the residual dose rates from the block after it is extracted from the target carriage is performed to estimate the radiation fields and to decide on the type of storage container to be used.

Popova, Irina I↗
DOE OSTI2021

Deletion of the cytochrome bc complex from Heliobacterium modesticaldum results in viable but non-phototrophic cells

The heliobacteria, a family of anoxygenic phototrophs, possess the simplest known photosynthetic apparatus. Although they are photoheterotrophs in the light, the heliobacteria can also grow chemotrophically via pyruvate metabolism in the dark. In the heliobacteria, the cytochrome bc complex is responsible for oxidizing menaquinol and reducing cytochrome c 553 in the electron flow cycle used for phototrophy. However, there is no known electron acceptor for the mobile cytochrome c 553 other than the photochemical reaction center. We have, therefore, hypothesized that the cytochrome bc complex is necessary for phototrophy, but unnecessary for chemotrophic growth in the dark. Here, we used a two-step method for CRISPR-based genome editing in Heliobacterium modesticaldum to delete the genes encoding the four major subunits of the cytochrome bc complex. Genotypic analysis verified the deletion of the petCBDA gene cluster encoding the catalytic components of the complex. Spectroscopic studies revealed that re-reduction of cytochrome c 553 after flash-induced photo-oxidation was over 100 times slower in the petCBDA mutant compared to the wild-type. Steady-state levels of oxidized P 800 (the primary donor of the photochemical reaction center) were much higher in the petCBDA mutant at every light level, consistent with a limitation in electron flow to the reaction center. The petCBDA mutant was unable to grow phototrophically on acetate plus CO 2 but could grow chemotrophically on pyruvate as a carbon source similar to the wild-type strain in the dark. The mutants could be complemented by reintroduction of the petCBDA gene cluster on a plasmid expressed from the clostridial eno promoter.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Coupled Charge and Radiation Transport Processes in Thermophotovoltaic and Thermoradiative Cells

Accurate modeling of charge transport and both thermal and luminescent radiation is crucial to the understanding and design of radiative thermal energy converters. Charge-carrier dynamics in semiconductors are well-described by the Poisson-drift-diffusion equations, and thermal radiation in emitter-absorber structures can be computed using multilayer fluctuational electrodynamics. These two types of energy flows interact through radiation absorption or luminescence and charge-carrier generation or recombination. However, past research has typically only assumed limited interaction, with thermal radiation absorption as an input for charge-carrier models to predict device performance. To examine this assumption, we develop a fully coupled iterative model of charge and radiation transport in semiconductor devices, and we use our model to analyze near-field and far-field GaSb thermophotovoltaic and thermoradiative systems. By comparing our results with past methods that do not consider cross-influences between charge and radiation transport, we find that a fully coupled approach is necessary to accurately model photon recycling and near-field enhancement of external luminescence. As these effects can substantially alter device performance, our modeling approach can aid in the design of efficient thermophotovoltaic and thermoradiative systems.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Impact of artificial topological changes on flow and transport through fractured media due to mesh resolution

Abstract We performed a set of numerical simulations to characterize the interplay of fracture network topology, upscaling, and mesh refinement on flow and transport properties in fractured porous media. We generated a set of generic three-dimensional discrete fracture networks at various densities, where the radii of the fractures were sampled from a truncated power-law distribution, and whose parameters were loosely based on field site characterizations. We also considered five network densities, which were defined using a dimensionless version of density based on percolation theory. Once the networks were generated, we upscaled them into a single continuum model using the upscaled discrete fracture matrix model presented by Sweeney et al. (2019). We considered steady, isothermal pressure-driven flow through each domain and then simulated conservative, decaying, and adsorbing tracers using a pulse injection into the domain. For each simulation, we calculated the effective permeability and solute breakthrough curves as quantities of interest to compare between network realizations. We found that selecting a mesh resolution such that the global topology of the upscaled mesh matches the fracture network is essential. If the upscaled mesh has a connected pathway of fracture (higher permeability) cells but the fracture network does not, then the estimates for effective permeability and solute breakthrough will be incorrect. False connections cannot be eliminated entirely, but they can be managed by choosing appropriate mesh resolution and refinement for a given network. Adopting octree meshing to obtain sufficient levels of refinement leads to fewer computational cells (up to a 90% reduction in overall cell count) when compared to using a uniform resolution grid and can result in a more accurate continuum representation of the true fracture network.

58 GEOSCIENCES↗