Search NASA⌕ Search

SEARCH · Search NASA

Results for “Separated flows”

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 217 records · Page 12

Intermittent Electron-Only Reconnection at Lunar Mini-Magnetospheres

Lunar crustal magnetic anomalies (LCMA) are sub-ion-gyroradius structures that have been shown to stand off the solar wind (SW) plasma from the Moon’s surface, forming shock-like discontinuities and reflecting incident SW protons. In this Letter, the results of high-resolution, 2D fully kinetic simulations show a bursty electron-only magnetic reconnection in the SW-LCMA interaction region, characterized by the quasiperiodic formation and ejection of magnetic islands and strong parallel electron flows along the X-point separator lines. The islands are observed to modify the magnetic pressure pileup and Hall electric field above the LCMA, leading to sharp increases in reflected protons that drive electromagnetic fluctuations downstream and short distances upstream in the SW.

79 ASTRONOMY AND ASTROPHYSICS↗

Accelerating transients with NekRS: GPU overlapping domain implementation and multi-rate timestepping

The simulation of nuclear transients using Computational Fluid Dynamics (CFD) presents significant computational challenges due to the inherent complexity and the wide separation in temporal scales between various flow physical phenomena. These disparities lead to high computational costs, often making the simulation of transients impractical without advanced techniques. Consequently, multiple research initiatives are being pursued by the NEAMS thermal-hydraulic area, some driven by academic institutions and some by national laboratories. Overall, they are exploring novel methods to make transient simulations more feasible and efficient. This report delves into recent advancements within the CFD code NekRS, specifically those achieved in Fiscal Year 2024 under the CONNECT effort, aimed at improving the performance and feasibility of transient simulations. The first major advancement involves the porting of NekRS to Aurora, one of the Department of Energy’s (DOE) most powerful supercomputers. Additionally, the report discusses the implementation of an overlapping domain capability within NekRS. This novel GPU-accelerated capability allows different spatial regions of the domain to be solved independently, enhancing the code’s efficiency, particularly when running large-scale simulations in complex domains. The scalability of this approach is demonstrated, highlighting its potential to transform how transients are approached in CFD simulations. Lastly, the report focuses on how this overlapping domain capability specifically accelerates transient simulations through multi-rate timestepping. By decoupling different regions and facilitating faster computations, this method offers a promising pathway to making nuclear transient simulations more computationally feasible, addressing one of the critical bottlenecks in the field. Together, these advancements represent a significant leap forward in transient simulation technology, bringing closer the possibility of handling highly complex nuclear scenarios with greater efficiency and accuracy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hypersonic Fluid–Structure Interaction on a Cone–Slice–Ramp Geometry

In this study, fluid–structure interactions were measured between a representative control surface and the hypersonic flow deflected by it. The control surface is simplified as a spanwise finite ramp placed on a longitudinal slice of a cone. The front surface of the ramp contains a thin panel designed to respond to the unsteady fluid loading arising from the shock-wave/boundary-layer interactions. Experiments were conducted at Mach 5 and Mach 8 with ramps of different angles. High-speed schlieren captured the unsteady flow dynamics and accelerometers behind the thin panel measured its structural response. Panel vibrations were dominated by natural modes that were excited by the broadband aerodynamic fluctuations arising in the flowfield. However, increased structural response was observed in two distinct flow regimes: 1) attached or small separation interactions, where the transitional regime induced the strongest panel fluctuations. This was in agreement with the observation of increased convective undulations or bulges in the separation shock generated by the passage of turbulent spots, and 2) large separated interactions, where shear layer flapping in the laminar regime produced strong panel response at the flapping frequency. In addition, panel heating during the experiment caused a downward shift in its natural mode frequencies.

33 ADVANCED PROPULSION SYSTEMS↗

Intermittent Electron-only Reconnection at Lunar Mini-magnetospheres

Abstract Lunar crustal magnetic anomalies (LCMA) are sub-ion-gyroradius structures that have been shown to stand off the solar wind (SW) plasma from the Moon’s surface, forming shock-like discontinuities and reflecting incident SW protons. In this Letter, the results of high-resolution, 2D fully kinetic simulations show a bursty electron-only magnetic reconnection in the SW-LCMA interaction region, characterized by the quasiperiodic formation and ejection of magnetic islands and strong parallel electron flows along the X-point separator lines. The islands are observed to modify the magnetic pressure pileup and Hall electric field above the LCMA, leading to sharp increases in reflected protons that drive electromagnetic fluctuations downstream and short distances upstream in the SW.

79 ASTRONOMY AND ASTROPHYSICS↗

Realization of an Asymmetric Non‐Aqueous Redox Flow Battery through Molecular Design to Minimize Active Species Crossover and Decomposition

Abstract This communication presents a mechanism‐based approach to identify organic electrolytes for non‐aqueous redox flow batteries (RFBs). Symmetrical flow cell cycling of a pyridinium anolyte and a cyclopropenium catholyte resulted in extensive capacity fade due to competing decomposition of the pyridinium species. Characterization of this decomposition pathway enabled the rational design of next‐generation anolyte/catholyte pairs with dramatically enhanced cycling performance. Three factors were identified as critical for slowing capacity fade: (1) separating the anolyte–catholyte in an asymmetric flow cell using an anion exchange membrane (AEM); (2) moving from monomeric to oligomeric electrolytes to limit crossover through the AEM; and (3) removing the basic carbonyl moiety from the anolyte to slow the protonation‐induced decomposition pathway. Ultimately, these modifications led to a novel anolyte–catholyte pair that can be cycled in an AEM‐separated asymmetric RFB for 96 h with >95 % capacity retention at an open circuit voltage of 1.57 V.

Shrestha, Anuska↗

Design of a gas-solid-solid separator to remove ash from circulating fluidized bed reactors

Cyclones are one of the most common types of gas-solid separators used in circulating fluidized bed boilers. However, cyclones typically do not allow ash to leave the system through the cyclone exit, causing ash to build up in the fluidized bed and necessitating additional systems to remove ash that builds up in the bed. In this study, an alternative “disengager” gas-solid separator is proposed as a way of inherently separating small and large solids, resulting in a gas-solid-solid separation system where ash is allowed to leave the system along with gas while the desired fluidized bed material is retained. Unlike cyclones, which rely on centrifugal force to separate solids and gas, the disengager separates based on entrainment velocity of the particles. Upwards-flowing gas and particles strike a deflection plate and enter the disengaging chamber where particles with low terminal velocity such as ash fines flow with the gas, while larger particles such as sand fall to the bottom of the separator and are returned to the fluidized bed. In this study, several different proposed disengager configurations are simulated and compared to a typical cyclone using computational fluid dynamic (CFD) simulations. Furthermore, it was found that separation efficiency in the disengager is strongly influenced by the size of the deflection plate, rather than by the size of the unit itself. The predicted separation efficiency showed that compared to a cyclone, the disengager design allows significantly more ash to exit the system but retains a similar amount of desirable material. Additionally, the disengager was predicted to not suffer significantly more erosion that a cyclone.

42 ENGINEERING↗

Fiber Sorbents – A Versatile Platform for Sorption-Based Gas Separations

Increasing demand for high-purity fine chemicals and a drive for process intensification of large-scale separations have driven significant work on the development of highly engineered porous materials with promise for sorption-based separations. While sorptive separations in porous materials offer energy-efficient alternatives to longstanding thermal-based methods, the particulate nature of many of these sorbents has sometimes limited their large-scale deployment in high-throughput applications such as gas separations, for which the necessary high feed flow rates and gas velocities accrue prohibitive operational costs. These processability limitations have been historically addressed through powder shaping methods aimed at the fabrication of structured sorbent contactors based on pellets, beads or monoliths, commonly obtained as extrudates. These structures overcome limitations such as elevated pressure drops commonly recorded across powder adsorption beds but often accrue thermal limitations arising from elevated particle density and aggregation, which ultimately cap their maximum separation performance. Furthermore, the harsh mechanical strain to which powder particles are subjected during contactor fabrication, in the form of extrusion/compression forces, can result in partial pore occlusion and framework degradation, further limiting their performance. Here, we present the development of porous fiber sorbents as an alternative sorbent contactor design capable of addressing sorbent processability limitations while enabling an array of performance-maximizing heat integration capabilities. This new sorbent form factor leverages pre-existing know-how from hollow fiber spinning to produce fiber-shaped sorbent contactors through the phase inversion of known polymers in a process known as dry-jet/wet quenching. The process of phase inversion allows microporous sorbent particles to be latched onto a macroporous polymer matrix under mild processing conditions, thus making it compatible with soft porous materials prone to amorphization under traditional pelletization conditions. Sorbent fibers can be created with different geometries through control of the spinning apparatus and process, offering the possibility to produce monolithic and hollow fibers alike, the latter of which can be integrated with thermalization fluid flows. In this Account, we summarize our progress in the field of fiber sorbents from both design and application standpoints. We further guide the reader through the evolution of this field from the early inceptive work on zeolite hollow fibers to recent developments on MOF fibers. We highlight the versatile nature of fiber sorbents, both from the composition, fabrication and structure points of view, and further demonstrate how fiber sorbents offer alternative paths in tackling new and challenging chemical separation challenges like direct air capture (DAC), with a final perspective on the future of the field.

36 MATERIALS SCIENCE↗

Uranium Oxide Elutriation at Higher Flow Rates and Longer Times

This work supports the recovery of high-assay low-enriched uranium (i.e., U enriched in 235 U to >5% and <20%) from used nuclear fuel. One step includes elutriating triuranium octoxide (U 3 O 8 ) from a fluidized bed of alumina (Al 2 O 3 ) that fills the main reaction chamber. The experiments outlined in this report assessed the effects of the time and flow rate on the elutriation of U 3 O 8 powder from the reaction vessel. Testing at Oak Ridge National Laboratory (ORNL) was performed in a one-fourth scale experimental apparatus based on a pilot-scale system, which is the Material Recovery Pilot Plant (MRPP) at Idaho National Laboratory (INL). The experimental system comprises a supply of pressurized nitrogen (N 2 ), a gas heater, the reaction vessel, a knock-out pot, and a water trap where elutriated particles are collected for sampling. An Al 2 O 3 fluidized bed is contained within the reaction vessel. A bed of Al 2 O 3 with a nominal size of 180 µm and U 3 O 8 with a nominal size of 7 µm was fluidized at 2 different flow rates so the U 3 O 8 physically separated from the bed. The apparatus was simple and had no valves or sharp bends downstream from the reactor. The flow exited the reactor through a side port at the top. The flow passed through a knock-out vessel before the flow was passed through a bubbler to capture the U 3 O 8 . Seven tests were completed with U 3 O 8 using an elutriating gas of N 2 heated to 200°C. Tests were performed with bed loadings that comprise 25.0 g of U 3 O 8 , or 3.6 wt% of the bed mass. Except for one test performed at 12 SLPM, the gas flow rate through the bed for the other experiments was set at 10 SLPM as measured by the mass flow meter. Linear velocities through the bed were calculated as 28 and 34 cm/s, or 0.9 and 1.1 ft·s -1 . The duration of the tests ranged from 20 to 180 min, and one test had a flow cycled between 0 and 5 SLPM every 5 min for 40 min. The testing successfully demonstrated elutriation of U 3 O 8 from the fluidized bed as up to a 100% of the U 3 O 8 was transferred out of the bed. The removal of U 3 O 8 had a logarithmic dependence on flow rate. The design changes made the apparatus easier to clean, and the mass balances were better than those used in previous studies. Some elutriated Al 2 O 3 collected in the tubing that connected the reactor to the knock-out pot, but the amount of Al 2 O 3 was minimal for 10 SLPM. The flow at 12 SLPM showed considerable transport of Al 2 O 3 bed material into the reaction-vessel headspace and the transfer line to the knock-out pot.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Dynamic separation of gases using microsieves

Separation of light weight molecules, such as nitrogen, argon, and oxygen, from heavier compounds can have significant impacts on energy capture, environmental monitoring, or isotopic applications. Large-scale gas separation techniques, like gas centrifugation and membrane mitigation, can be problematic as they impart tremendous energy and induce high mechanical stress onto the instrumentation. Microsieves, also known as micronozzles or microfunnels, are developed to create physical barriers to separate specific isotopes and gases. Separation is achieved using a converging and diverging micronozzle to impose supersonic gas flow around a curved wall, and it has been used for the separation of heavy actinide isotopes in low weight gas as well as separation of low weight gas compositions of nitrogen and argon back in 1900s. However, systematic reviews of this unique technology are lacking. The application of the Laval style nozzle, which has a converging/diverging entrance fundamental to the micronozzle, is included in this review due to its importance in industrial applications in uranium (U) isotope refinement. Using advanced computational fluid dynamic (CFD) simulations, the extent of gas separation can be modelled. Herein, we first examine the literature and survey recent advances on fabrication techniques for creating curved micronozzles, methods and separation principles used to design devices. Furthermore, we then follow with highlights of CFD simulations applied to evaluate the separation effects using microsieves. Finally, identification of the gap and recommendation for future development and applications are suggested for using intrinsic molecular features and fluidic dynamics in formulating separation strategies.

30 Microfluidics↗

REBOUND: Reverse Engineering Bidirectional Outflow Under Non-Equilibrium Diffusion

Rare-earth elements (REEs) are essential for electronics, renewable energy, and defense technologies. However, the current supply of REEs relies on mining concentrated in a few countries and energy-intensive separations. DOE’s Basic Energy Sciences (BES) program has launched a grand challenge which aims to ensure a sustainable supply of critical REEs by developing innovative and environmentally friendly separation methods. As an alternative to costly and harmful traditional methods, the Non-Equilibrium Transport Driven Separations (NETS) initiative has created a microfluidic Y-channel co-flow method that applies external fields to exploit magneto- and electrohydrodynamic effects for separating dilute REE ions from complex feedstocks. Computational fluid dynamics (CFD) studies have identified a few operating conditions with promising ion selectivity and separation efficiency. However, challenges remain regarding Y-channel versatility across feedstocks and accurate incorporation of physical phenomena into CFD models. In this work, we develop a multi-fidelity modelling approach which integrates experimental results with CFD simulation to build a surrogate model for the dependence of separation efficiency to variation of design parameters. The surrogate model enables a reinforcement learning (RL) method to adaptively launch CFD and experimental runs, improving model fidelity around optimal Y-channel parameters.

36 MATERIALS SCIENCE↗

NREL's HydroGEN Photoreactor Platform Design and Characterization

This technical report details the design, fabrication, and construction of a flow-through photoreactor developed at The National Renewable Energy Laboratory (NREL). The photoreactor is designed to facilitate indoor or outdoor testing on photoelectrochemical or integrated solar water splitting devices. The photoreactor design described here is based on a previous demonstration. The photoreactor version described herein contains several advanced design features that include: 1) Chassis-chuck two-part design, 2) Separable counter electrode compartments (here, the counter electrodes are anodes), 3) Reduced electrode separation to reduce electrolyte resistance, 4) Improved flow pattern for bubble removal, 5) Mounting for spectroradiometer receptor, 6) Fresnel lens and collimating tube attachments. The photoreactor presented here was designed primarily to accommodate photocathode materials in an illuminated compartment and a dark anode as the oxygen evolution catalyst. The build materials were thus chosen to for acid electrolyte compatibility. Other embodiments could accommodate a photoanode with dark hydrogen evolution catalyst and other build materials for alkaline compatibility.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sample analysis systems and methods of use thereof

The invention generally relates to sample analysis systems and methods of use thereof. In certain aspects, the invention provides a system for analyzing a sample that includes an ion generator configured to generate ions from a sample. The system additionally includes an ion separator configured to separate at or above atmospheric pressure the ions received from the ion generator without use of laminar flowing gas, and a detector that receives and detects the separated ions.

42 ENGINEERING↗

Resource Recovery and Environmental Protection in Wyoming’s Greater Green River Basin Using Selective Nanostructured Membranes (Final Report)

Produced water (PW) represents a sizable waste stream that is co-generated with oil and natural gas production. In 2021 Wyoming ranked 8th and 9th, respectively in domestic oil and natural gas production. In 2017 Wyoming ranked as the 4th highest generator of PW in the U.S, accounting for 7% of the total volume generated. In the context of being the 3rd most arid state in the U.S., the value of water reuse becomes obvious. PW reuse, and resource recovery, in any form requires some level of treatment to remove particulates, residual (free, dispersed) hydrocarbons, organics, and salts. The level of treatment depends on the requirements of the reuse, or resource recovery, application. PW management systems in Wyoming employ in order of volume of PW managed the following management strategies: reinjection for enhanced oil recovery, surface discharge, deep well injection, evaporation ponds (impoundments), and commercial management/treatment. Complicating treatment efforts are the relatively high concentrations of organics (natural and synthetic), dispersed/free hydrocarbons, benzene-toluene-ethylbenzene, and xylenes (BTEX) compounds, biologicals, salts, and minerals. Hydrocarbons (dispersed/dissolved crude oils) and BTEX compounds, as well as synthetic organics, present economic and environmental concerns. The former represents lost revenue, while the latter results in negative environmental impacts like emissions from surface impoundments. The overall objective of this proposal was to synthesize superhydrophilic/oleophobic and superhydrophobic/oleophilic membranes for selectively concentrating and then separating BTEX compounds and oil and grease (O&G) from PW originating from the Greater Green River Basin (GGRB) in Wyoming. Three specific research aims were pursued to accomplish this overall objective. This final report details the development of the superhydrophobic and superhydrophilic membranes, as well as the design of the membrane module prototypes specifically. The technoeconomic assessment is separately reported in another document. 1. Aim #1 – Material optimization and performance evaluation of superhydrophilic/oleophobic and superhydrophobic/oleophilic membranes made by electrospinning/spraying. 2. Aim #2 – Design and construction of cross-flow membrane modules for selectively concentrating and then separating BTEX/oil from GGRB produced water. 3. Aim #3 – Techno-economic assessment of BTEX/oil recovery, and clean water production, using superhydrophilic/oleophobic and superhydrophobic/oleophilic membrane separation for GGRB PW. Superhydrophobic membranes were synthesized by electrospinning poly(vinylidene fluoride-co-hexafluoropropylene) (PVDF-HFP) nanofibers onto polyester (PET) substrates and electrospraying nano-carbon black/PVDF-HFP onto the nanofibrous layer. These membranes were characterized by high (>8000 liters per square meter per hour per bar (LMH/bar)) permeance values for pure hydrocarbon phases and a high hydrocarbon selectivity (>96%) when treating GGRB PW. All results were obtained when operating the membrane in a crossflow configuration representative of actual field operating conditions. Solvent/oil properties, specifically viscosity and total surface energy/tension, affected permeance across the membrane, which resulted in light mineral oil (394 LMH/bar) and o-xylene (1834 LMH/bar) being characterized by lower permeance values in the pure phase tests. Mixed phase fluxes between 40 to 80 LMH were obtained for the PW when operating the membrane at a feed pressure of 0.3 bar. Flux decreased as the mixed phase concentration in the feed decreased pointing to the importance of maximizing the collision efficiency between the emulsion and the membrane surface and maximizing the emulsion concentration in the feed and the turbulence within the feed channel. These tests demonstrated that the superhydrophobic membranes developed here are a viable hydrocarbon recovery method for GGRB PWs and should be pursued for testing in pilot-scale trials. Superhydrophilic membranes were successfully synthesized via electrospinning/spraying using polyacrylonitrile (PAN) nanofibers as a base nanofibrous matrix. Integration of polyaniline (PANI) into the nanofibrous matrix produced a superior membrane, for water filtration applications, relative to PAN alone and reduced graphene oxide (RGO) when integrated into the nanofibrous matrix. This conclusion was based on the PANI-PAN resistance to flux loss (fouling) when treating model solvent/oil solutions representative of GGRB PWs and field collected PW from the GGRB. The synthesized PAN membranes outperformed a commercially available PAN membrane designed for oil/water separation. This finding indicates that the surface chemical and physical characteristics of the electrospun membranes presents improved properties for filtration of challenging waters like GGRB PWs. The electrospun membranes therefore show promise overall as a substitute for conventionally polymerized membranes in PW treatment applications. The PANI-PAN membrane specifically presents superior performance characteristics for concentration O&G prior to treatment by the hydrocarbon recovery membrane and producing high-quality filtrate for reuse and/or additional treatment (desalination).

02 PETROLEUM↗

An extension of Thwaites’ method for turbulent boundary layers

Thwaites developed an approximate method for determining the evolution of laminar boundary layers. The approximation follows from an assumption that the growth of a laminar boundary layer in the presence of pressure gradients could be parameterized solely as a function of the Holstein–Bohlen flow parameter, thus reducing the von Kármán momentum integral to a first-order ordinary differential equation. This method is useful for the analysis of laminar flows, and in computational potential flow solvers to account for the viscous effects. In this work, an approximate method for determining the momentum thickness of a two-dimensional, turbulent boundary layer is proposed following Thwaites’ work. It is shown that the method provides good estimates of the momentum thickness for multiple boundary layers, including both favourable and adverse pressure gradient effects, up to the point of separation. In the limit of high Reynolds numbers, it is possible to derive a criterion for the onset of separation from the proposed model, which is shown to be in agreement with prior empirical observations. The sensitivity of the separation location with respect to upstream perturbations is also analysed through this model for the NASA/Boeing speed bump and the transonic Bachalo–Johnson bump.

17 WIND ENERGY↗

Electrochemical CO2 conversion to formic acid through the Andora Process

Electrochemical utilization or conversion of CO2 can be used to convert waste CO2 into targeted high-value and cost-effective commercial products. Conventional processes for CO2 conversion to fuel sources are challenging and expensive due to their requirement for high temperatures and pressures. Electrochemically, a catalyst can be used at a fixed applied potential to reduce CO2 to form a desired product with less competition at modest operating conditions. Formic acid (FA) has recently become a product of interest for CO2 conversion due to its potential for hydrogen storage and fuel cell applications. Electrochemically, FA is produced through a direct two-electron transfer process involving CO2 with a proton source, requiring less energy input and fewer reaction steps than traditional processes. However, there remain significant obstacles for the broader-scale implementation of CO2 on the market, including competition with unwanted CO and H, the presence of molecular oxygen, and stability of the electrocatalyst. In order to overcome some of these challenges, a new continuous flow Andora Process was developed consisting of two separate cells to decouple electrochemical reduction of the charge carrier with the FA production via an engineered catalyst. Recent experimental results show that with appropriate design changes and operating conditions to the flow system, FA production above 500 mM can be achieved along with efficiencies above 80%. TEA and LCA assessments of the Andora Process show a potential for the reduction in the GWP potential and cost saving compared to the current comparative formic acid production process.

60 APPLIED LIFE SCIENCES↗

Electrically Driven Ion Separations in Permeable Membranes

Ion separations are critical for producing materials such as Li and rare-earth element salts that are essential for applications including energy storage and production. This research aimed to develop and understand separations of singly and multiply charged ions as well as separations among different singly charged ions such as Li + and Na + . Thin layers of polyelectrolytes on membranes allow passage of singly charged ions at rates 1000 times faster than those for multiply charged ions. This allows purification of mixtures with equal concentrations of Li + and Mg 2+ to create a permeate solution with >99% Li + . Such membranes may also effectively isolate rare earth elements. To separate Li + and K + , this work developed opposing flow and current in membranes. The current overcomes the flow of K + to allow Li + to selectively pass through the membrane, leading to Li + /K + selectivities around 70. All of the research employs modelling to better understand the mechanism of these separations, which includes transport driven by gradients in electrical potential, concentration, and pressure. The interplay of these driving forces can create highly selective separations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Separating Sensible and Latent Cooling with Electrically Charged Rotating Vortexes and Vapor Capturing Air Handler Technology

The dehumidification of air is a highly energy-intensive process, typically requiring vapor compression systems to cool the air below its dew point temperature. This project investigated air dehumidification achieved by electrostatic and dielectrophoresis forces in electrosprays and centrifugal principles in cyclone separators. As a practical implementation of these processes, Electrospray Vortical Flow eXchanger (EVFX) prototypes were designed, fabricated, and tested. Electrosprays were employed to energize water droplets, which then captured water vapor via dielectrophoresis. A cyclone separator was positioned downstream of the electrosprays. This combination increased the dehumidification efficiency and controlled the flow of the water vapor condensate removed from the air.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Using Homogenized Models to Explore the Effect of Fracture Densities on Weathering

Despite its importance, only a few researchers have incorporated the effects of fracturing into models of reactive transport for rock weathering. Here we explore 2D simulations that describe weathering under conditions of diffusive and advective transport within heterogeneous media consisting of rocky blocks and fractures. In our simulations, the Darcy velocities vary in space and time and depend on weathering processes within the rock matrix. We explore simulations with saturated and unsaturated flow for weathering bedrock that consists of blocks separated by inert or weathered material. The simulations show that a simplified homogenized model can approximate exact solutions for some of the simulated columns and hills and can allow exploration of coupling between flow and reaction in fractured rock. These hillslope simulations document that, even in the presence of 2D water flow, i) an increase in fracture density results in faster weathering advance rates; and ii) the water table locates deeper for a rock system that is weathered and fractured rather than weathered and unfractured. Some of these patterns have also been observed for natural systems. But these simulations also highlight how simplified models that do not use appropriate averaging of heterogeneities can be inaccurate in predicting weathering rate for natural systems. For example, if water flows both vertically and laterally through the vadose zone of a hill, then a prediction of the depth of regolith that is based on modeling strictly unidirectional downward infiltration will be unrealistically large. Likewise, if the fracture density observed near the land surface is used in a model to predict depth of weathering for a system where the fracture density decreases downward, the model will overestimate regolith depth. Learning how to develop accurately homogenized models could thus enable better conceptual models and predictions of weathering advance in natural systems.

58 GEOSCIENCES↗