Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Porous flow”

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 523 records · Page 29

Molecular Origin of Wettability Alteration of Subsurface Porous Media upon Gas Pressure Variations

Upon extraction/injection of a large quantity of gas from/into a subsurface system in shale gas production or carbon sequestration, the gas pressure varies remarkably, which may significantly change the wettability of porous media involved. Mechanistic understanding of such changes is critical for designing and optimizing a related subsurface engineering process. Using molecular dynamics simulations, in this work we have calculated the contact angle of a water droplet on various solid surfaces (kerogen, pyrophyllite, calcite, gibbsite, and montmorillonite) as a function of CO 2 or CH 4 gas pressure up to 200 atm at a temperature of 300 K. The calculation reveals a complex behavior of surface wettability alteration by gas pressure variation depending on surface chemistry and structure, and molecular interactions of fluid molecules with surfaces. As the CO 2 gas pressure increases, a partially hydrophilic kerogen surface becomes highly hydrophobic, while a calcite surface becomes more hydrophilic. Considering kerogen and calcite being the major components of a shale formation, we postulate that the wettability alteration of a solid surface induced by a gas pressure change may play an important role in fluid flows in shale gas production and geological carbon sequestration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular-gas-dynamics simulations of turbulent Couette flow over a mean-free-path-scale permeable substrate

Here we report flow statistics and visualizations from molecular-gas-dynamics simulations using the direct simulation Monte Carlo (DSMC) method for turbulent Couette flow in a minimal domain where the lower wall is replaced by an idealized permeable fibrous substrate representative of thermal-protection-system materials for which the Knudsen number is O(10 -1 ). Comparisons are made with smooth-wall DSMC simulations and smooth-wall direct numerical simulations (DNS) of the Navier-Stokes equations for the same conditions. Roughness, permeability, and noncontinuum effects are assessed. In the range of Reynolds numbers considered herein, the scalings of the skin friction on the permeable substrate and of the mean flow within the substrate suggest that they are dominated by viscous effects. While the regenerative cycle characteristic of smooth-wall turbulence remains intact for all cases considered, we observe that the near-wall velocity fluctuations are modulated by the permeable substrate with a wavelength equal to the pore spacing. Additionally, the flow within the substrate shows significant rarefaction effects, resulting in an apparent permeability that is 13% larger than the intrinsic permeability. In contrast, the smooth-wall DSMC and DNS simulations exhibit remarkably good agreement for the statistics examined, despite the Knudsen number based on the viscous length scale being as large as O(10 -1 ). This latter result is at variance with classical estimates for the breakdown of the continuum assumption and calls for further investigations into the interaction of noncontinuum effects and turbulence.

42 ENGINEERING↗

Fluoride-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model

In this report we present work performed in Fiscal Year 2022 that demonstrates the modeling and simulation of a fully coupled neutronics thermal hydraulics reference plant model for a fluoride-cooled high-temperature pebble-bed reactor. The multiphysics model is developed on the Nuclear Regulatory Commission’s Comprehensive Reactor Analysis Bundle (BlueCRAB) available on the Idaho National Laboratory’s high-performance computer, which natively and seamlessly couples Griffin, Pronghorn, and the BISON Multiphysics Object-Oriented Simulation Environment based applications. Griffin provides reactor physics capabilities, including depletion to the equilibrium core, k-eigenvalue, adjoint, and transient. The unique direct equilibrium core capability in Griffin is based on a streamline methodology to spatially deplete the pebbles into burnup groups. Pronghorn solves the porous medium equations for the fluid regions and conduction in the solid regions and incorporates a fluidic diode model to simulate the transition from forced to natural convection during accident scenarios. MOOSE modules solves thermal conduction problems for the pebbles and tristructural isotropic in the pebble-bed core, thus providing the fuel and moderator spatial fields for each pebble burnup group. The neutronics feedback relies primarily on fuel, moderator, and reflector temperatures as as well as the FLiBe salt density. Here, we present results for the uncoupled equilibrium core and perform comparisons to equivalent Monte Carlo models. The power distributions and kinetic parameters obtained with Griffin are consistent with those computed with Griffin. We demonstrate a noticeable improvement with the use of discrete ordinates method (SN) transport. The coupled steady-state equilibrium core provides the initial condition for two time-dependent problems: a control rod withdrawal event and an unprotected loss of flow event. In both cases, the reactor design is self-stabilizing and the solutions are consistent with the expected physics. Although this model is prototypical regarding BlueCRAB’s capabilities, its results are consistent with published work by Kairos Power and other research entities. Significant improvements to the model are planned in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

X-ray imaging of liquid-liquid Mg cavitation

The coupling of dynamic loading platforms with synchrotron light sources have enabled experiments examining void formation in optically opaque materials. Much like spallation in a solid, cavitation occurs when rarefaction waves collide within a liquid. The resulting deformation leads to a complex flow field with numerous voids. It has been suggested that compaction of such a structure may be a significant source of ejecta (i.e. shallow bubble collapse). Experiments imaging supported and unsupported shocks interacting with a free surface in shock melted Mg with synchrotron radiation are presented. Here, the results are used to examine the structure of a cavitation plane that forms in the unsupported case and provide an estimate of the porosity within this region.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

2D-to-3D image translation of complex nanoporous volumes using generative networks

Image-based characterization offers a powerful approach to studying geological porous media at the nanoscale and images are critical to understanding reactive transport mechanisms in reservoirs relevant to energy and sustainability technologies such as carbon sequestration, subsurface hydrogen storage, and natural gas recovery. Nanoimaging presents a trade off, however, between higher-contrast sample-destructive and lower-contrast sample-preserving imaging modalities. Furthermore, high-contrast imaging modalities often acquire only 2D images, while 3D volumes are needed to characterize fully a source rock sample. In this work, we present deep learning image translation models to predict high-contrast focused ion beam-scanning electron microscopy (FIB-SEM) image volumes from transmission X-ray microscopy (TXM) images when only 2D paired training data is available. We introduce a regularization method for improving 3D volume generation from 2D-to-2D deep learning image models and apply this approach to translate 3D TXM volumes to FIB-SEM fidelity. We then segment a predicted FIB-SEM volume into a flow simulation domain and calculate the sample apparent permeability using a lattice Boltzmann method (LBM) technique. Results show that our image translation approach produces simulation domains suitable for flow visualization and allows for accurate characterization of petrophysical properties from non-destructive imaging data.

58 GEOSCIENCES↗

DECOVALEX-2023: Task E Final Report

This is the Task E final report for DECOVALEX-2023. Task E is focused on understanding thermal, two-phase hydrological, and mechanical (TH2M) processes, especially related to predicting brine migration in the excavation damaged zone around a heated excavation in salt. Salt is attractive as a disposal medium for radioactive waste because it is self-healing and is essentially impermeable and essentially non-porous in the far field (away from excavations). Investigation of the short-term (days to years) near-field (centimeters to tens of meters) behavior of salt is important for radioactive waste disposal because this early period strongly controls the amount of brine in a salt repository. Brine leads to corrosion of waste forms and waste packages, and possible dissolution of radionuclides with brine transport being a potential transport vector to the accessible environment. The main test case used in Task E is the ongoing Brine Availability Test in Salt (BATS) heater test located underground at the Waste Isolation Pilot Plant (WIPP) near Carlsbad, New Mexico, USA. The Task was divided into a series of Steps. Step 0 was an introduction to processes in salt, that included matching historical unheated brine inflow data from boreholes at WIPP and matching temperature observations during BATS heater test 1a. Step 1 included validation of models against a thermo-poroelastic analytical solution relevant to heated boreholes in salt, and two-phase flow around an excavation in salt. Step 2 required all the individual components covered in steps 0 and 1 to come together to match observed brine inflow behavior during the BATS 1a heater test. There were a range of approaches from the teams, from mechanistic to prescriptive. Given the uncertainties in the problem, some teams used one- or two-dimensional models of the processes, while other teams included more geometrical complexity in three-dimensional models. The key learning points from Task E have been: • Heat conduction through salt typically requires non-linear thermal conductivity (as a function of temperature), but most models do a good job matching observations, given appropriate adjustments to the applied power and some thermocouple locations. • Thermal pressurization requires coupled thermal-hydrological-mechanical (THM) responses that consider the thermal expansion of the fluid and solid phases. • Initialization of two-phase flow models around a borehole or excavation in salt are more realistically represented as “wetting up”, rather than “drying down” (i.e., the initial state after excavation is mostly dry, rather than mostly wet). • The BATS 1a heater test includes a significant release of brine after the end of heating, which requires a large increase in permeability to recreate. Task E has been a great learning experience for all the teams involved, and feedback from the modeling teams has led to changes in the design of follow-on BATS experiments, which are now ongoing underground at WIPP. There was a balance throughout the task between freedom to model phenomena how each team saw fit, and prescriptiveness in problem design to bring the modeling teams closer together to allow attribution of smaller differences between models to different modeling choices. The modeling approaches seem to go through two phases: an early phase of discovery or testing, and a later phase of refinement and improvement. In future modeling efforts, different field data could be used (e.g., BATS 2) and more time should be included in the processes for teams to make multiple model refinement or even significant changes to their conceptual model or setup, based on lessons learned from the modeling exercise.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Performance Analysis and Limiting Parameters of Cross-flow Membrane-based Liquid-desiccant Air Dehumidifiers

We report that to dehumidify a humid air stream, existing air conditioning (AC) systems substantially overcool the outdoor humid air below its dew point, thereby significantly reducing energy efficiency. Directly capturing humidity, membrane-based liquid-desiccant dehumidification systems separate sensible and latent cooling (SSLC) loads and thus offer a promising pathway for a high-performance AC solution. Design of an energy-efficient SSLC-AC system, however, rests largely on detailed understating of the dehumidification process. While some studies have identified the dehumidification process mainly depends on membrane characteristics, other studies have argued that the process is limited by desiccant liquid or alternatively air thermo-hydraulic physics for typical humid climate conditions. The present study examines performance and physics of the membrane-based liquid-desiccant dehumidification process over a wide range of climate conditions through a novel 3D, two-phase, multi-species CFD model. Decoupling the thermodynamic and hydraulic effects, the study reveals that the dehumidification rate is a linear function of the water vapor pressure potential ($J=α ΔP$) summarizing the system's thermodynamic state. The slope of the curve (i.e., α) depends on hydraulic transport characteristics of the membrane pores, air stream, and desiccant solution. More importantly, it was found that the air dehumidification process is mainly limited by the air-side transport physics for thin liquid-desiccant films and commonly used porous superhydrophobic membranes. Additionally, results show that, depending on ambient/desiccant conditions and physical dehumidifier characteristics, energy effectiveness and dehumidification rate vary from 13 to 34% and from 0.13 to 1.4 g m -2 s -1 , respectively. Therefore, the present study allows to efficiently design future SSLC-based AC systems exhibiting high performance energy metrics.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Modeling Microwave-Enhanced Chemical Vapor Infiltration Process for Preventing Premature Pore Closure

The chemical vapor infiltration (CVI) process involves infiltrating a porous preform with reacting gases that undergo chemical transformation at high temperatures to deposit the ceramic phase within the pores, ultimately leading to a dense composite. The conventional CVI process in composite manufacturing needs to follow an isothermal approach to minimize temperature differences between the external and internal surfaces of the preform, ensuring that reactive gases infiltrate internal pores before external surfaces seal. Here, this study addresses the challenge of premature pore closure in CVI processes through microwave heating. A frequency-domain microwave solver is developed in OpenFOAM to investigate volumetric heating mechanisms within the preform. Through numerical studies, we demonstrate the capability of microwave heating of creating an inside-out temperature inversion. This inversion accelerates reactions proximal to the preform center, effectively mitigating the risk of premature external pore closure and ensuring uniform densification. The results reveal a significant enhancement in temperature inversion when high-permittivity reflectors are incorporated to generate resonant waves. This microwave heating strategy is then coupled with high-fidelity direct numerical simulation (DNS) of reacting flow, enabling the analysis of resulting densification processes. The DNS includes detailed chemistry and realistic diffusion coefficients. The numerical results can be used to estimate the impact of microwave-induced temperature inversion on densification in productions.

42 ENGINEERING↗

Enhanced carbon-transfer and -utilization efficiencies achieved using membrane carbonation with gas sources having a range of CO 2 concentrations

The economic viability of microalgal biofuels relies on increasing productivity in a cost-effective manner. As microalgal biomass contains > 50% carbon, a high rate of CO 2 delivery is required for high productivity, and inefficient CO 2 delivery amplifies operating costs. Membrane carbonation using non-porous hollow fiber membranes can ideally deliver CO 2 without bubble formation and high carbon transfer efficiency. Because CO 2 streams from industrial resources are not 100% CO 2 , the buildup of inert gasses can significantly lower the CO 2 delivery rate when the distal end of the membrane is closed. To overcome the buildup of inert gases, we managed the distal end of the membranes with three different approaches: fully open end, restricted bleed valve, and restricted bleed valve with pH-actuated venting. For all approaches, CO 2 was delivered to membranes ondemand based on a pH set point. Evaluating a wide range of CO 2 concentrations (10% to 100%), we found that all approaches eliminated the buildup of inert gases, could maintain target pH values and gave the same biomass productivities and carbon distributions. However, carbon transfer efficiency depended on the operation of the distal end. Fully open-end operation gave a poor carbon transfer efficiency because of excessive loss of CO 2 from the distal end. However, restricting the exit flow rate to ≤4 cm 3 /min mitigated the problems of excessive CO 2 loss, but without incurring a large loss of CO 2 -delivery flux. For the continuous cultivation, combining a restricted bleed valve with pH-actuated venting improved the carbon-transfer efficiency and -utilization efficiencies up to 85% and 67%, respectively, with a sufficient CO 2 delivery flux.

42 ENGINEERING↗

Modeling Microwave-Enhanced Chemical Vapor Infiltration Process for Preventing Premature Pore Closure

The chemical vapor infiltration (CVI) process involves infiltrating a porous preform with reacting gases that undergo chemical transformation at high temperatures to deposit the ceramic phase within the pores, ultimately leading to a dense composite. The conventional CVI process in composite manufacturing needs to follow an isothermal approach to minimize temperature differences between the external and internal surfaces of the preform, ensuring that reactive gases infiltrate internal pores before external surfaces seal. This study addresses the challenge of premature pore closure in CVI processes through microwave heating. A frequency-domain microwave solver is developed in Open-FOAM to investigate volumetric heating mechanisms within the preform. Through numerical studies, we demonstrate the capability of microwave heating of creating an inside-out temperature inversion. This inversion accelerates reactions proximal to the preform center, effectively mitigating the risk of premature external pore closure and ensuring uniform densification. The results reveal a significant enhancement in temperature inversion when high-permittivity reflectors are incorporated to generate resonant waves. This microwave heating strategy is then coupled with high-fidelity direct numerical simulation (DNS) of reacting flow, enabling the analysis of resulting densification processes. The DNS simulation includes detailed chemistry and realistic diffusion coefficients. The numerical results can be used to estimate the impact of microwave-induced temperature inversion on densification in productions.

Ge, Wenjun↗

Pronghorn Theory Manual

This manual is intended to provide a thorough description of the governing equations; the relations used for material, heat transfer, and fluid flow parameters; and the numerical method utilized in Pronghorn. Detailed derivations of the governing equations are given, beginning from first principles in order to be as explicit as possible about assumptions made along the way and guide any future methods extensions. Additional explanatory or tangential information that is not required for detailed understanding of the governing equations used in Pronghorn is shown in gray boxes in the text. The convention used in this document is to represent all class names (or file names) in typewriter font.

42 ENGINEERING↗

Pore-Scale Oil Connectivity and Displacement by Controlled-Ionic-Composition Waterflooding Using Synchrotron X-Ray Microtomography

Controlled-ionic-composition waterflooding is an economic and effective method to improve oil recovery in carbonate oil reservoirs. Recent studies show controlling the salinity and ionic composition of injection water can alter the wettability of carbonate mineral surfaces. The pore-scale oil connectivity and displacement by controlled-ionic-composition waterflooding in heterogeneous carbonate reservoirs, especially at the early stage, is still unclear. The goal of this study is to examine the role of ion concentrations and types in the oil displacement efficiency and investigate the impact of the waterflooding on the pore-scale oil displacement using the national synchrotron facility. A carbonate rock sample was flooded with synthetic high-salinity water and other water solutions with different sulfate concentrations. The waterflooding processes were visualized with synchrotron X-ray microtomography to follow the evolution of pore-scale oil/brine interactions at typical field flow rates. Experimental results show that the water with lower sulfate concentration and higher salinity did not change the wettability of the pore surfaces. Higher sulfate ion concentrations in the water, in contrast, altered the wettability of carbonate pore surfaces from oil-wet to neutral-wet within the first few minutes of waterflooding. Finally, novel insight was gained on the ability of water with high-sulfate concentration to displace oil in the small pores and through abundant oil channels, which could consequently lead to higher oil recovery from the carbonate rock.

04 OIL SHALES AND TAR SANDS↗

Organic colloid composition in variable-redox porewaters within a mountainous floodplain

Redox gradients, often driven by changes in sediment moisture levels in porous, heterogeneous groundwater systems, create dynamic conditions that may promote the production and transport of colloids within natural waters. While much research has focused on the inorganic composition of colloids, the organic composition remains less well understood. Organic matter (OM) in colloids may associate with minerals, complex metal ions, and serve as an electron donor for microbial respiration; therefore, its composition is of high interest. We examined the composition of porewater OM along a redox gradient in a riparian soil located along the Slate River in Crested Butte, Colorado, USA as a function of depth (90, 130, 200, and 350 cm below ground surface). All depths were oxic to suboxic, except 200 cm, where the products of iron and sulfate reduction were observed concomitant with an increase in dissolved and/or colloidal OM, pH, alkalinity, and conductivity. We investigated the composition of porewater using correlated scanning transmission X-ray microscopy and transmission electron microscopy. We observed a change in the OM chemistry from carboxylate-rich at the 200 cm depth (where it was also enmeshed with non-crystalline iron) to phenol- and substituted-aromatic-rich at other depths. Radiocarbon dating revealed carbon in porewater at 200 cm was younger than depths above and below. Soil porewater can flow down into the underlaying gravel bed during baseflow conditions, thus we speculate whether riparian porewater could transport iron- and carboxylate-rich organic colloids into surrounding surface waters through the gravel bed conduit.

Alpine floodplain↗

Data from Stewart et al. 2026 "Organic Colloid Composition in Variable-Redox Porewaters within a Mountainous Floodplain"

Redox gradients, often driven by changes in sediment moisture levels in porous, heterogeneous groundwater systems, create dynamic conditions that may promote the production and transport of colloids within natural waters. While much research has focused on the inorganic composition of colloids, the organic composition remains less well understood. Organic matter (OM) in colloids may associate with minerals, complex metal ions, and serve as an electron donor for microbial respiration; therefore, its composition is of high interest. We examined the composition of porewater OM along a redox gradient in a riparian soil located along the Slate River in Crested Butte, Colorado, USA as a function of depth (90, 130, 200, and 350 cm below ground surface). All depths were oxic to suboxic, except 200 cm, where the products of iron and sulfate reduction were observed concomitant with an increase in dissolved and/or colloidal OM, pH, alkalinity, and conductivity. We investigated the composition of porewater using correlated scanning transmission X-ray microscopy and transmission electron microscopy. We observed a change in the OM chemistry from carboxylate-rich at the 200 cm depth (where it was also enmeshed with non-crystalline iron) to phenol- and substituted-aromatic-rich at other depths. Radiocarbon dating revealed carbon in porewater at 200 cm was younger than depths above and below. Soil porewater can flow down into the underlaying gravel bed during baseflow conditions, thus we speculate whether riparian porewater could transport iron- and carboxylate-rich organic colloids into surrounding surface waters through the gravel bed conduit. This data package includes geochemical data (.csv) as a function of depth from porewaters in OBJ1. Measurements included are pH, conductivity, TOC, Fe concentration, sulfate concentration, dissolved oxygen. Additionally radiocarbon dating age of carbon in samples is provided. Additionally, this dataset contains metadata files (.csv) describing the methods, file contents, and column headers.

alpine floodplain↗

Thermal Modeling and Simulation of the Packed-bed Thermal Energy Storage Combined with INL Thermal Energy Distribution System

Dynamic Energy Transport and Integration Lab (DETAIL) at Idaho National Laboratory is to support experimental demonstration and validation research on Nuclear-Renewable Hybrid Energy System [1]. The Thermal Energy Distribution System (TEDS) is a thermal-hydraulic flow loop in DETAIL with its own dedicated control system to support the integration of co-located multiple experimental systems, where a packed-bed thermal energy storage (TES) is installed as a thermal buffer and storage unit. Among various TES options, the packed-bed TES is adopted in TEDS because it offers a low-cost single-tank thermal storage option compared to the traditional two-tank TES. However, since the TES tank is filled with granular fillers having different thermophysical properties from those of TES tank wall, there is a potential thermo-mechanical issue to be carefully addressed like thermal ratcheting which may pose a significant design concern for the packed-bed TES tank. Thermal ratcheting is a thermomechanical process caused by the repeated rearrangement of granular filler inside a TES tank during continuous thermal cycling operation of the packed-bed TES system. If the thermally induced stress exceeds the yield strength of a TES tank wall, catastrophic consequences may happen like rupture of the TES tank. Thus, it is important to understand the phenomenon to ensure the robust operation of the packed-bed TES tank. Given that thermal ratcheting is caused by complex interaction of thermal transport in the porous bed and solid mechanics, the accurate prediction of transient thermal behavior of the packed-bed TES tank, which is the focus of this paper, is critical to the reliable thermal ratcheting analysis. This paper discusses the numerical modeling, simulation, and validation studies that are ongoing at INL to investigate the transient thermal behavior of the packed-bed TES. Of particular concern is the transient thermal process occurring in the packed-bed TES unit that is operated in conjunction with the INL TEDS. The main goal of this research is three-fold: (i) provide preliminary insights into the transient thermal behavior of the TEDS TES tank, (ii) support the thermal measurement and validation plan for TEDS experiment, and (iii) provide transient thermal boundary conditions to support the reliable thermal ratcheting analysis of the TEDS TES tank. For the transient thermal modeling and analysis, a CFD model was developed, and the validity of the modeling approach was examined via comparing the numerical simulation results with the experimental data obtained from various design characteristics of packed-bed TES tanks. Then, the present modeling method was applied for the transient thermal analysis of the TEDS TES tank, and the results are discussed along with the potential improvement of data acquisition strategy for the future TEDS experiments for more precise validation study.

25 ENERGY STORAGE↗

Rapid charging made practical in graphite-based lithium batteries: surface-acoustic wave turbulent electrolyte mixing to overcome diffusion limited charging rates [Final Report]

One of the key limits in rapidly recharging a lithium-ion battery is the depletion of lithium ions within the electrolyte adjacent the anode during charging and long diffusion time to overcome this depletion. It also causes dendrite formation, inefficient use of the lithium, and battery degradation over many charge-discharge cycles. Because the liquid electrolyte remains quiescent and unmixed, this depletion layer’s depth rapidly grows to match the anode-cathode separation distance at even modest charge rates. The solution proposed by PI Prof. James Friend and Co-PI Prof. Ping Liu, both from the University of California, San Diego, is to mix the electrolyte and minimize the Li+ ion concentration gradient during charging. Even in the presence of the separator, the charging rate could then be significantly increased. They propose to accomplish this using surface acoustic wave (SAW)-driven acoustic streaming, a technique employing 10-mW fingernail-sized solid-state devices from the telecommunications industry to drive turbulent mixing to submicron length scales in a manner completely compatible with the typical 20700 and 18650 cells used in electric vehicles. 0.1. Turbulent acoustic streaming mixes the electrolyte during charging—even with separator Our proposed 100-MHz SAW device used to recirculate the electrolyte is compatible with lithium- ion battery electrochemistry, as it is made in our lab of single-crystal lithium niobate. Their solution is straightforward to drive from a DC power source alongside the signal provided during battery charging. The device is only needed during charging and does not consume power dur- ing battery discharge. Uniquely, SAW generates extreme accelerations of over 1 billion meters per second squared in the fluid, driving turbulent mixing from centimeter to submicron length scales, even through the porous separator materials that tend to be used in batteries, all while avoiding interfering with the anode’s solid electrolyte interphase layer that forms during use. By employing novel fluid mechanics, the investigators are proposing a new direction for battery research away from direct use of materials science and electrochemistry. The chemistry agnostic solution may be employed in any battery chemistry that makes use of liquid electrolytes, providing a broader transformative benefit to the battery research community. They are to produce a series of prismatic and 20700-cell sized 2 Ah batteries capable of being charged and discharged at least 500 times without more than 20% loss in battery capacity, and to provide analysis tools useful to the battery research and development community for adopting this approach for other battery chemistries and configurations.

25 ENERGY STORAGE↗

A Geo‐Structurally Based Correction Factor for Apparent Dissolution Rates in Fractured Media

Abstract Field measurements of apparent geochemical weathering reaction rates in subsurface fractured porous media are known to deviate from laboratory measurements by multiple orders of magnitude. To date, there is no geologically based explanation for this discrepancy that can be used to predict reaction rates in field systems. Proposed correction factors are typically based on ad hoc characterizations related to geochemical kinetic models. Through a series of high‐fidelity reactive transport simulations of mineral dissolution within explicit 3D discrete fracture networks, we are able to link the geo‐structural attributes with reactive transport observations. We develop a correction factor to linear transition state theory for the prediction of the apparent dissolution rate based on measurable geological properties. The modified rate law shows excellent agreement with numerical simulations, indicating that geological structure could be a primary reason for the discrepancy between laboratory and field observations of apparent dissolution rates in fractured media.

58 GEOSCIENCES↗

A mass-transfer particle-tracking method for simulating transport with discontinuous diffusion coefficients

The problem of a spatially discontinuous diffusion coefficient (D(x)) is one that may be encountered in hydrogeologic systems due to natural geological features or as a consequence of numerical discretization of flow properties. To date, mass-transfer particle-tracking (MTPT) methods, a family of Lagrangian methods in which diffusion is jointly simulated by random walk and diffusive mass transfers, have been unable to solve this problem. This manuscript presents a new mass-transfer (MT) algorithm that enables MTPT methods to accurately solve the problem of discontinuous D(x). To achieve this, we derive a semi-analytical solution to the discontinuous D(x) problem by employing a predictor-corrector approach, and we use this semi-analytical solution as the weighting function in a reformulated MT algorithm. Furthermore, this semi-analytical solution is generalized for cases with multiple 1D interfaces as well as for 2D cases, including a 2 × 2 tiling of 4 subdomains that corresponds to a numerically-generated diffusion field. The solutions generated by this new mass-transfer algorithm closely agree with an analytical 1D solution or, in more complicated cases, trusted numerical results, demonstrating the success of our proposed approach.

42 ENGINEERING↗