SEARCH · Search NASA

Results for “Numerical fluid dynamics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Multiphase computational fluid dynamics modeling of reacting flows in absorption columns for carbon capture

First-principles derived computational fluid dynamics (CFD) simulations have been proposed as a fundamental tool for investigating solvent-based CO 2 absorption in packed columns due to their ability to accurately represent the underlying nonlinear, multiscale dynamics. Numerous studies have previously utilized such CFD simulations to investigate hydrodynamics of columns with structured and random packings by assessing the key hydrodynamic metrics such as the interfacial and wetted areas. While mapping such metrics for different conditions is essential to the optimization of absorption columns, it is not sufficient, as the CO 2 capture rate depends also on the coupled, nonlinear dynamics from the underlying chemical reaction kinetics, thermodynamics, and heat-transfer rates. In this work, we present detailed CFD simulation results augmented by incorporating the effects of interfacial physical mass transfer of CO 2 , heat release from chemical reaction kinetics, and thermophysical property variations from resulting temperature gradients. We demonstrate the applicability of the proposed approach in numerically assessing the performance of packed columns by evaluating key hydrodynamic quantities, CO 2 absorption rates, and temperature rise in a reference column with packings that are structurally similar to the Sulzer Mellapak™ 250.Y packing, for different solvent inflow velocities and temperatures. Predictions from simulation results are found to be consistent with the trends in experimental observations from the literature, suggesting that the predictive capabilities of the simulation framework can be leveraged to guide the future development of absorber-column designs and optimized process flowsheets.

Absorption columns↗

DOE OSTI2022

The role of breakup and coalescence in fine-scale bubble-induced turbulence. I. Dynamics

We study the effect of bubble breakup and coalescence on fine-scale dynamics of bubbly turbulent flows using direct numerical simulations. We perform two different simulations of dilute bubbly flows of void fraction 0.5%: one with bubbles breaking up and coalescing and the other without these physical processes. The volume of the fluid method is used for simulating bubbles undergoing breakup and coalescence while the bubbles are treated as rigid spheres in the immersed boundary method simulation. The energy spectrum in both types of simulation, consistent with previous studies, exhibits a -3 slope. We follow a single infinitesimal fluid element as it evolves to understand velocity gradient dynamics using conditional mean trajectories. We note finite-time divergence when the fluid element evolves under the action of inertial and pressure dynamics. The inertial, pressure, and viscous velocity gradient dynamics, when considered individually, produce the same results for bubble-induced turbulence (BIT) as with the classical homogeneous isotropic turbulence (HIT). Yet when the overall velocity gradient dynamics is considered, BIT results in non-cyclic trajectories moving toward stable node and unstable saddle while classical HIT shows cyclic behavior in their trajectories that move toward the origin. Interestingly, both the volume of fluid and immersed boundary simulations produce similar results. Therefore, there are two main takeaways from this research. First, new velocity gradient models are needed for BIT as their velocity gradients behavior is entirely different from the HIT. Second, we can neglect the bubble topology, breakup, and coalescence while studying or modeling the fine-scale dynamics of BIT.

42 ENGINEERING↗

DOE OSTI2025

Multidimensional Modeling of Mixture Formation in a Hydrogen-Fueled Heavy-Duty Optical Engine With Direct Injection

Hydrogen (H 2 ), as a carbon-free fuel, is considered as one of the most promising solutions to reduce the carbon footprint of hard-to-decarbonize energy and transportation sectors. As such, hydrogen-fueled internal combustion engines (H 2 ICEs) have recently been receiving increasing attention, particularly in applications such as on-road/off-road heavy-duty transport and combined heat and power. The direct injection (DI) of gaseous hydrogen into the combustion chamber offers great potential for achieving high power density and high engine efficiency, while mitigating the risk of backfire and reducing pre-ignition. However, the numerical simulation of H 2 DI system remains a formidable challenge associated with the high computational cost of reproducing compressible supersonic flow and shocks in narrow injector passages and in near-nozzle regions. In general, there is a lack of well-established and validated practices for the modeling of high-pressure H 2 DI in large-bore engines. Here, to this end, this study focuses on computational fluid dynamics (CFD) modeling of the mixture formation process in a heavy-duty optical engine employing a medium-pressure H 2 DI system. Both large eddy simulations (LES) and Reynolds Averaged Navier–Stokes (RANS) simulations are performed and evaluated against optical data. Gaseous hydrogen is injected into the combustion chamber via a centrally located outward opening hollow-cone injector at a pressure of 40 bar. Simulations are carried out for two injection timings, namely, −120 and −60 °CA. The numerical predictions for H 2 distribution in different horizontal and vertical planes during the compression stroke are systematically compared against optical data obtained through planar laser-induced fluorescence (PLIF) measurements. Overall, the LES approach using the Dynamic Structure model is found to have good predictive capabilities for the early jet penetration in terms of length and shape, as well as the later H 2 distributions. However, the unsteady RANS approach with the renormalization group $k - ϵ$ model, which is widely used by industry to model heavy-duty ICEs, significantly underpredicts the H 2 mixing, even at similar mesh resolution to that used in LES. These results indicate that there is a need for the improvement of mixing submodels within the RANS approach when applied to H 2 DI simulations.

LES↗

DOE OSTI2023

Theory and simulations of linear and nonlinear two-dimensional Rayleigh–Taylor dynamics with variable acceleration

Interfacial Rayleigh–Taylor mixing is crucial to describing important natural and engineering processes, such as exploding supernovae, laser micromachining, hot spots in inertial confinement fusion, and optical telecommunications. These require the characterization of the time dependence of the driving acceleration. We compare our theoretical formulation based on group theory foundations with interface-capturing numerical simulations for linear and nonlinear two-dimensional Rayleigh–Taylor instabilities in a finite-sized domain with time-varying acceleration over broad ranges of Atwood numbers and acceleration exponents. Detailed corroboration between theory and simulations is provided for this foundational case. Both demonstrate the strong interfacial nature of Rayleigh–Taylor instabilities, which suggests that practical flow fields can be reconstructed from the derived fluid potential using the proposed theory. A robust agreement is also obtained for the early and late-time evolution of the amplitudes of the bubble and spike, which demonstrate that the Rayleigh–Taylor flow can transition to the mixing regime even for a single-mode initial perturbation. Corroboration with experiments of high energy density plasmas motivated by studies of supernovae is also achieved. In addition, a long-standing puzzle in Rayleigh–Taylor dynamics on the interplay between the acceleration, the shear, and the interface morphology in the theory and simulations is resolved by accounting for finite viscosity of the fluids. The characterization of Rayleigh–Taylor instabilities as a highly interfacial phenomenon provides valuable insight into its multiscale nature, which enhances the design and understanding of numerous processes of practical interest.

42 ENGINEERING↗

DOE OSTI2020

Modelling of nanometer scale dust grains in tokamak

Dust poses a serious threat to tokamak operation and safety. It is important to study the behaviour of dust grains under tokamak's discharge conditions, which depends heavily on their size and charge. Existing simulations mainly address issues on dust grains with radii larger than 1 μm, in which case, the drift effect due to electromagnetic fields can be safely ignored. For nanometer scale dust grains, however, the drift effect becomes significant and a new model based on guiding-centre system needs to be established. In this work, the NDS has been done under BOUT++ framework. The simulation contains two parts. Part one, NDS evaluates the charging and ablation processes of the dust grains. In the second part, the guiding-centre orbits of dust particles are tracked in tokamak plasmas, whose parameters are obtained from BOUT++, a highly desirable C++ code package for performing parallel plasma fluid simulations with an arbitrary number of equations in 3D curvilinear coordinates. The orbit of nanodust dynamics is described by guiding centre equations for simplicity, and these equations are numerically solved by conventional fourth-order Runge Kutta method. Simulations provide results such as trajectories and evolutions of dust particles with different sizes and velocities for different tokamak geometries. Results show tungsten dust grains with a radius of a few nanometers launched from outer midplane will oscillate before totally ablated in C-Mod. The oscillation in this case is driven by the ion drag force. Larger Nanodust with a radius of 100 nm, on the contrary, cannot be completely constrained by the electromagnetic field. In conclusion, the high plasma temperature and density in the seperatrix region causes severe dust ablation, resulting in total ablation within several ms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

On the Departure from Monin–Obukhov Surface Similarity and Transition to the Convective Mixed Layer

Large-eddy simulations are used to evaluate mean profile similarity in the convective boundary layer (CBL). Particular care is taken regarding the grid sensitivity of the profiles and the mitigation of inertial oscillations in the simulation spin-up. The nondimensional gradients Φ for wind speed and air temperature generally align with Monin–Obukhov similarity across cases but have a steeper slope than predicted within each profile. The same trend has been noted in several other recent studies. The Businger-Dyer relations are modified here with an exponential cutoff term to account for the decay in Φ to first-order approximation, yielding improved similarity from approximately 0.05z i to above 0.3z i , where z i is the CBL depth. The necessity for the exponential correction is attributed to an extended transition from surface scaling to zero gradient in the mixed layer, where the departure from Monin–Obukhov similarity may be negligible at the surface but becomes substantial well below the conventional surface layer height of 0.1 z i .

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Programmable topotaxis of magnetic rollers in time-varying fields

We discuss how spatially uniform, time-periodic magnetic fields can be designed to power and direct the migration of ferromagnetic spheres up (or down) local gradients in the topography of a solid substrate. Our results are based on a dynamical model that considers the time-varying magnetic torques on the particle and its motion through the fluid at low Reynolds number. We use both analytical theory and numerical simulation to design magnetic fields that maximize the migration velocity up (or down) an inclined plane. We show how “topotaxis” of spherical particles relies on differences in the hydrodynamic resistance to rotation about axes parallel and perpendicular to the plane. Importantly, the designed fields can drive multiple independent particles to move simultaneously in different directions as determined by gradients in their respective environments. Experiments on ferromagnetic spheres provide evidence for topotactic motions up inclined substrates. The ability to program the autonomous navigation of driven particles within anisotropic environments is relevant to the design of colloidal robots.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Physics-Reinforced Machine Learning Algorithms for Multiscale Closure Model Discovery

The central objective of this project was to address the challenge of modeling and simulating complex multiscale turbulence phenomena by leveraging physics-guided machine learning (PGML) and hybrid modeling approaches. By integrating physics-based methods with data-driven models, the research focused on achieving robust and scalable solutions for geophysical turbulence, enhancing numerical weather prediction and climate research tools. The project resulted in significant advancements in computational modeling paradigms, predictive tools for reduced-order modeling, and innovative algorithms for fluid dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Glovebox Fires Phase 1: Simulations of Open Burner Experiments

Understanding the duration of a fire required to breach a glovebox via glove ports and/or windows within a facility is important in developing mitigation strategies for safety. To facilitate this understanding, both simulation and experiments are utilized. The approach is to first validate the Computation Fluid Dynamics (CFD) code, Fire Dynamics Simulator (FDS), with experimental data collected by New Mexico Tech (NMT) and then perform simulations of full-scale rooms containing gloveboxes to assess numerous scenarios which would otherwise be cost-prohibitive experimentally. This report provides comparison to first-phase experiments involving a fire without a glovebox. A subsequent report will provide comparison to experiments involving a glovebox. The results indicate good agreement with FDS tending to over predict the pre-mixed and diffusion-mode tests by 6% and 10%, respectively.

42 ENGINEERING↗

DOE OSTI2021

Lagrangian particle model for 3D simulation of pellets and SPI fragments in tokamaks

A 3D numerical model for the ablation of pellets and shattered pellet injection (SPI) fragments in tokamaks in the plasma disruption mitigation and fueling parameter space has been developed based on the Lagrangian particle code [R. Samulyak, X. Wang, H.-S. Chen, Lagrangian Particle Method for Compressible Fluid Dynamics, J. Comput. Phys., 362 (2018), 1-19]. The pellet code implements the low magnetic Reynolds number MHD equations, kinetic models for the electronic heating, a pellet surface ablation model, an equation of state that supports multiple ionization states, radiation, and a model for grad-B drift of the ablated material across the magnetic field. The Lagrangian particle algorithm is highly adaptive, capable of simulating a large number of fragments in 3D while eliminating numerical difficulties of dealing with the tokamak background plasma. The code has achieved good agreement with theory for spherically symmetric ablation flows. Axisymmetric simulations of neon and deuterium pellets in magnetic fields ranging from 1 to 6 Tesla have been compared with previous simulations using the FronTier code, and very good agreement has also been obtained. Furthermore, the main physics contribution of the paper is a detailed study of the influence of 3D effects, in particular grad-B drift, on pellet ablation rates and properties of ablation clouds. Smaller reductions of ablation rates in magnetic fields compared to axially symmetric simulations have been demonstrated because the ablated material is not confined to narrowing channels in the presence of grad-B drift. Contribution of various factors in the grad-B drift model has also been quantified.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

NeuroSEM: A hybrid framework for simulating multiphysics problems by coupling PINNs and spectral elements

Multiphysics problems that are characterized by complex interactions among fluid dynamics, heat transfer, structural mechanics, and electromagnetics, are inherently challenging due to their coupled nature. While experimental data on certain state variables may be available, integrating these data with numerical solvers remains a significant challenge. Physics-informed neural networks (PINNs) have shown promising results in various engineering disciplines, particularly in handling noisy data and solving inverse problems in partial differential equations (PDEs). However, their effectiveness in forecasting nonlinear phenomena in multiphysics regimes, particularly involving turbulence, is yet to be fully established. Here, this study introduces NeuroSEM, a hybrid framework integrating PINNs with the highfidelity Spectral Element Method (SEM) solver, Nektar++. NeuroSEM leverages the strengths of both PINNs and SEM, providing robust solutions for multiphysics problems. PINNs are trained to assimilate data and model physical phenomena in specific subdomains, which are then integrated into the Nektar++ solver. We demonstrate the efficiency and accuracy of NeuroSEM for thermal convection in cavity flow and flow past a cylinder. The framework effectively handles data assimilation by addressing those subdomains and state variables where the data is available. We applied NeuroSEM to the Rayleigh-B´enard convection system, including cases with missing thermal boundary conditions and noisy datasets. Finally, we applied the proposed NeuroSEM framework to real particle image velocimetry (PIV) data to capture flow patterns characterized by horseshoe vortical structures. Our results indicate that NeuroSEM accurately models the physical phenomena and assimilates the data within the specified subdomains. The framework’s plug-and-play nature facilitates its extension to other multiphysics or multiscale problems. Furthermore, NeuroSEM is optimized for efficient execution on emerging integrated GPU-CPU architectures. This hybrid approach enhances the accuracy and efficiency of simulations, making it a powerful tool for tackling complex engineering challenges in various scientific domains.

42 ENGINEERING↗

DOE OSTI2022

On the Response of a Herschel–Bulkley Fluid Due to a Moving Plate

In this paper, we study the boundary-layer flow of a Herschel–Bulkley fluid due to a moving plate; this problem has been experimentally investigated by others, where the fluid was assumed to be Carbopol, which has similar properties to cement. The computational fluid dynamics finite volume method from the open-source toolbox/library OpenFOAM is used on structured quad grids to solve the mass and the linear momentum conservation equations using the solver “overInterDyMFoam” customized with non-Newtonian viscosity libraries. The governing equations are solved numerically by using regularization methods in the context of the overset meshing technique. The results indicate that there is a good comparison between the experimental data and the simulations. The boundary layer thicknesses are predicted within the uncertainties of the measurements. The simulations indicate strong sensitivities to the rheological properties of the fluid.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Midfidelity model verification for a point-absorbing wave energy converter with linear power take-off

In the preliminary design stage of a wave energy converter (WEC), researchers need fast and reliable simulation tools. High-fidelity numerical models are usually employed to study the wave-structure interaction, but the computational cost is demanding. As an alternative, midfidelity models can provide simulations in the order of real time. In this study, we operate Uppsala University’s WEC in a relatively mild sea state and model it using WEC-Sim. The model is verified based on OpenFOAM simulations. To analyze the ability of the midfidelity model to capture WEC dynamics, we investigate the system separately with 1, 2, and 3 degrees of freedom. We examine the contribution of viscous phenomena, and study both linear and weakly nonlinear solutions provided by WEC-Sim. Our results indicate that the viscous effects can be neglected in heave and surge motion, but not for pitch. We also find that the weakly nonlinear WEC-Sim solution successfully agrees with the computational fluid dynamics, whereas the linear solution could suggest misleading results.

16 TIDAL AND WAVE POWER↗

DOE OSTI2020

Numerical Investigation of a Central Fuel Property Hypothesis Under Boosted Spark-Ignition Conditions

In the present work, a central fuel property hypothesis (CFPH), which states that fuel properties are sufficient to provide an indication of a fuel’s performance irrespective of its chemical composition, was numerically investigated. In particular, the objective of the study was to determine whether Research Octane Number (RON) and Motor Octane Number (MON), as fuel properties, are sufficient to describe a fuel’s knock-limited performance under boosted spark-ignition (SI) conditions within the framework of CFPH. To this end, four TPRF-bioblendstock surrogates having different compositions but matched RON (=98) and MON (=90), were first generated using a non-linear regression model based on artificial neural network (ANN). Additionally, three unconventional bioblendstocks were included in the analysis: di-isobutylene (DIB), isobutanol, and Anisole. Skeletal reaction mechanisms were generated for the TPRF-DIB, TPRF-isobutanol, and TPRF-anisole blends from a detailed kinetic mechanism. Thereafter, numerical simulations were performed for the fuel surrogates using the skeletal mechanisms and a virtual cooperative fuel research (CFR) engine model, under a representative boosted operating condition. In the computational fluid dynamics (CFD) model, the G-equation approach was employed to track the turbulent flame front and the well-stirred reactor model combined with the multi-zone binning strategy was used to capture auto-ignition in the end-gas. In addition, laminar flame speed (LFS) was tabulated for each blend as a function of pressure, temperature, and equivalence ratio a priori, and the lookup tables were used to prescribe laminar flame speed as an input to the G-equation model. Parametric spark timing sweeps were performed for each fuel blend to determine the corresponding knock-limited spark advance (KLSA) and 50% burn point (CA50) at the respective KLSA timing. It was observed that despite same RON, MON, and engine operating conditions, the TPRF-anisole blend exhibited markedly different knock-limited performance from the other three blends. This deviation from the octane index (OI) expectation was shown to be caused by differences in laminar flame speed. However, it was found that relatively large fuel-specific differences in LFS (>20%) would have to be present to cause any appreciable deviation from the OI framework. Otherwise, RON and MON would still be robust enough to predict a fuel’s knock-limited performance.

42 ENGINEERING↗

DOE OSTI2020

Numerical Investigation of a Central Fuel Property Hypothesis under Boosted Spark-ignition Conditions

In the present work, a central fuel property hypothesis (CFPH), which states that fuel properties are sufficient to provide an indication of a fuel's performance irrespective of its chemical composition, was numerically investigated In particular, the objective of the study was to determine whether Research Octane Number (RON) and Motor Octane Number (MON), as fuel properties, are sufficient to describe a fuel's knock-limited performance under boosted spark-ignition (SI) conditions within the framework of CFPH. To this end, four TPRF-bioblendstock surrogates having different compositions but matched RON (= 98) and MON (= 90), were first generated using a non-linear regression model based on artificial neural network (ANN). Three unconventional bioblendstocks were included in the analysis: Di-isobutylene (DIB), Isobutanol and Anisole. Skeletal reaction mechanisms were generated for the TPRF-DIB, TPRF-isobutanol and TPRF-anisole blends from a detailed kinetic mechanism. Thereafter, numerical simulations were performed for the fuel surrogates using the skeletal mechanisms and a virtual cooperative fuel research (CFR) engine model, under a representative boosted operating condition. In the computational fluid dynamics (CFD) model, the G-equation approach was employed to track the turbulent flame front and the well-stirred reactor model combined with multi-zone binning strategy was used to capture auto-ignition in the end-gas. In addition, laminar flame speed was tabulated for each blend as a function of pressure, temperature and equivalence ratio a priori, and the lookup tables were used to prescribe laminar flame speed as an input to the G-equation model. Parametric spark timing sweeps were performed for each fuel blend to determine the corresponding knock-limited spark advance (KLSA) and 50% burn point (CA50) at the respective KLSA timing. It was observed that despite same RON, MON and engine operating conditions, the TPRF-Anisole blend exhibited markedly different knock-limited performance from the other three blends. This deviation from the octane index (DI) expectation was shown to be caused by differences in laminar flame speed (LFS). However; it was found that relatively large fuel-specific differences in LFS (> 20%) would have to be present to cause any appreciable deviation from the 01 framework. Otherwise, RON and MON would still be robust enough to predict a fuel's knock-limited performance.

Pal, Pinaki↗

DOE OSTI2020

Numerical Modeling and Parametric Study of a Dual Purpose Underground Thermal Battery

Buildings’ thermal energy systems for providing space heating/cooling consume a considerable amount of electricity nationwide. Integrating thermal energy storage with buildings’ thermal systems has the potential to shave the peak electric demand and to overcome the mismatch between the intermittent renewable power and the fluctuating daily demand for electricity. A novel dual-purpose underground thermal battery (DPUTB) has been developed to provide diurnal thermal energy storage in conjunction with a dual-source heat pump (DSHP)system. A two-dimensional numerical model is developed to simulate the transient performance of the DPUTB. This model has similar accuracy as a detailed three-dimensional model developed with commercial heat transfer and fluid dynamics simulation program but with 1000 times faster computational speed. A parametric study is performed to identify a design of the DPUTB that can provide 14 kWh thermal storage capacity for each charging and discharging cycle. A system simulation has been developed to assess the performance of the designed DPUTB when it is integrated with a DSHP for conditioning a typical residential building in summer and winter. Simulation results indicate that the DPUTB can meet all the design requirements in both summer and winter with the simple rule-based control strategy and the overall integrated DPUTB and DSHP system can shift the electric load for meeting the thermal demand from the peak hours to off-peak hours during a day, which would result in considerable energy cost savings if a Time of Use electricity tariff is applicable.

25 ENERGY STORAGE↗

DOE OSTI2025

Development of an Advanced Multiphysics Simulation Capability for Radiant's Microreactor Design

Argonne National Laboratory and Idaho National Laboratory, through a Department of Energy Gateway for Accelerated Innovation in Nuclear Voucher, supported key analysis needs of Radiant related to (i) air jacket thermal fluid performance, (ii) evaluation of decay heat source terms defining air jacket technical requirements, and (iii) assessment of modeling methodologies employed for core analysis. All of these activities center on numerical simulation of various aspects of Kaleidos using the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, the Cardinal multiphysics application, the OpenMC Monte Carlo code, and the Nek5000 computational fluid dynamics (CFD) code. This project builds upon an earlier Nuclear Energy Advanced Modeling and Simulation (NEAMS) Thermal-Hydraulic (T/H) Center of Excellence (CoE) project focused on initial demonstration of Cardinal multiphysics simulation of High Temperature Gas Reactors (HTGRs) and now focuses on Radiant’s Kaleidos concept.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Hyporheic Exchange Driven by Submerged Rigid Vegetation: A Modeling Study

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

54 ENVIRONMENTAL SCIENCES↗