Search NASA⌕ Search

SEARCH · Search NASA

Results for “Air simulations”

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 379 records · Page 21

Multiphysics Simulations of Nanofibrous High Power Targets for the Inference of Failure Modes

Fermilab s High Power Targetry Research and Development (HPT R&D) group have been developing and studying an electrospun nanofiber target concept to support the need for robust targets in future fixed target facilities. These nanofiber mats have demonstrated resistance to radiation damage, and the free motion of the individual fibers is expected to mitigate the cyclic stresses induced by a pulsed, high-power beam. To evaluate the efficacy of this concept, nanofiber mat samples have been sent by the HPT R&D group to the HiRadMat facility at CERN for prototypic thermal shock testing; the outcomes of these experiments show that the survivability of a nanofiber target depends on its construction parameters, in particular the packing density of the fibers. Samples with higher packing densities have consistently been destroyed by exposure to the HiRadMat beam, with a hole at the beam center visible, and layers of nanofibers peeled away from the center hole, whereas samples with lower densities have survived with limited damage. The exact reason for the failure of the higher density targets was unclear at the time of the original experiments, but the results of our recent multiphysics simulations which recreate the experiments support the hypothesis that the expansion and pressurization of the air inside the target after being heated by the pulsed beam is the cause; in a high density nanofiber mat, the motion of air through the pores of the mat is much more restricted, and induces a larger pressure on the fibers, blowing the mat apart. In this talk, we ll share the results of these simulations and discuss how they support this hypothesis.

Asztalos, Will↗

Aerodynamic effects on outer ash deposition rates in second generation atmospheric pressure oxy-coal combustion systems

Recent measurements in scenarios representative of second generation atmospheric pressure oxy-coal combustion systems have shown a significant increase in ash deposition rates in comparison to combustion in air. However, the causative mechanisms behind this increase have not been well understood. To fill this void, well-characterized experiments including fuel and deposit particle size distributions (PSDs) were coupled with highly resolved numerical simulations to isolate the aerodynamic effects impacting the deposition process such that mechanisms could be hypothesized. Three combustion scenarios (AIR, OXY27, OXY70) spanning a factor of three variation in flue gas volumetric flow rates were simulated and the deposition characteristics (impaction rates, deposit PSD, temperature, residence times, capture rates) tracked/predicted using a customized deposition module. The measured deposit PSD was significantly different from the PSD of the parent fuel indicating significant physio-chemical transformations (coalescence and particle growth in particular) at play. The use of simplistic modeling approaches (swelling parameter variations) to model particle growth led to inaccurate deposit ash PSD and rate predictions. A more satisfactory agreement between the measurements and simulations was obtained when the functional form of the parent fuel PSD (spread parameter) was modified to conform closely to the deposit PSD while still ensuring the fidelity of temperature and velocity predictions. Based on the measured ash compositions, identical capture criteria were employed across all three scenarios. Here this study further supports the theory that the ash deposition rates in these systems are dominated by aerodynamic effects with the ash PSD playing a dominant role.

42 ENGINEERING↗

CFD Analysis of RLUOB Zone 1 HEPA Filter Plenum and Testing Manifolds

The purpose of this computational fluid dynamics (CFD) analysis is to ensure that Camfil Farr’s upstream and downstream injection and sampling manifolds can meet or exceed the requirements outlined in the ASME AG-1 1997 a(2000) Code on Nuclear Air and Gas Treatment for the testing of HEPA and adsorbent filters. This paper will present a numerical simulation of airflow in the Radiological Laboratory Utility Office Building (RLUOB) zone 1 HEPA filter plenum and testing manifolds using the commercial CFD software ANSYS FLUENT 2020R1. The CFD analysis focuses on the investigation of the air flow distribution and air-aerosol mixing uniformity. The evaluation was done for all steps of the modeling process: grid generation, physics setup, simulation, and post-processing. The mass flow rate in each section of the zone 1 injection and sampling manifolds is also reported.

42 ENGINEERING↗

Air classification of forest residue for tissue and ash separation efficiency

The goal of this Case Study was to evaluate the performance of air classification of logging residues toward meeting conversion CMAs for carbon and ash contents, as compared to the static status quo Base Case system in which the residues are first dried and then ground in a hammer mill with a 6 mm screen and fines < 1.18 mm are removed. Also considered were moisture and ash impacts on throughput and Overall Operating Effectiveness (OOE), as well as delivered feedstock cost and minimum fuel selling price (MFSP). Laboratory data on the impacts of fan speed and moisture content on the separation efficiency of soil ash, needles and bark from white wood were received from FCIC Subtask 5.2. Average throughput and energy consumption data were obtained from the Bioenergy Feedstock National User Facility (BFNUF) for the same air classifier. These data were utilized to develop the necessary response surface equations to perform throughput analysis using discrete event simulation. Because the Base Case status quo system utilizes drying prior to grinding, we modeled the Case Study with drying prior to air classification and subsequent grinding of the separated white wood to isolate the individual quality and cost impacts of air classification relative to the Base Case system. Key takeaways from this Case Study are that air classification improves the quality of the final material, but increases the production cost, especially when lights are disposed; this becomes a tradeoff between increased conversion yield and the extra cost. Removing material should be done as early in the process as possible. Each operation that occurs prior to removing the material increases the cost of the disposed material and leads to wasted energy expenditures. As processes are included or modified, the impact on monetary and energy cost should be included in the decision process. Identifying alternative uses and the associated value for material that is separated from the feedstock stream will have a significant impact on the delivered cost of the material. Although there was an increase in production cost by adding air classification, there was a resulting benefit to the MFSP when meeting or exceeding all of the conversion CMAs; this was an important finding that should be explored further with FCIC Subtask 8.3 and potentially FCIC Task 6.

09 BIOMASS FUELS↗

Large Eddy Simulations of the Atmospheric Boundary Layer Over Satellite‐Sensed Sea Ice Maps

Surface heterogeneity in the marginal ice zone (MIZ) causes multiscale secondary atmospheric circulations that are challenging to model or observe. The absence or inadequate representation of these circulations in ocean‐atmosphere exchange schemes in climate models is partially responsible for the underestimation of Arctic sea ice loss. Observationally, such circulations obfuscate the interpretation of polar atmospheric chemistry measurements, among others. To address this open challenge, large‐eddy simulations are conducted over real‐world satellite‐sensed sea ice maps with an ice‐sea temperature contrast, as well as over idealizations of these maps that alter the ice pattern but conserve its fraction, showing that the ice fraction of a surface is not sufficient to predict the surface heat flux. In a second suite of simulations, three other heterogeneity metrics (representing the surface fragmentation, patch edge tortuosity, and patch size variability) are introduced to complement the ice fraction in describing the surface. Simulations varying these parameters suggest that they also significantly impact surface‐air interactions. A roughness contrast is then added to the surface temperature contrast, indicating that the contribution of roughness changes to the resulting atmospheric circulations is less pronounced than that of thermal heterogeneity. Based on these results, we illustrate, using a multi‐linear regression on these map features, that MIZ surface parameterizations in Earth Systems Models can be improved if they account for these various characteristics of the sea‐ice patterns.

54 ENVIRONMENTAL SCIENCES↗

2D Modeling of Plasma Streamer and Glow Phases at Ammonia-Air Flame Conditions

Streamer and glow plasma phases have been modeled at two thermochemical states of an ammonia-air flame: fresh reactants and burnt products. A new AMReX-based solver has been verified against benchmarks in the literature and has been used to perform these simulations. A Helmholtz-equation based photoionization model with parameters accounting for the presence of NH3 in air has been coupled with the solver to accurately model the streamer propagation phase. A detailed plasma kinetics mechanism has been compiled and used to predict the evolution of electrons, excited states, ions, and radicals during streamer propagation and glow formation. The propagation velocity of streamers was found increase by almost two-fold when the mixture was changed from the fresh reactants to the burnt products. Moreover, vibrational excitation was found to be limited to the streamer body, whereas ionization predominantly occurred at the streamer head, as is expected. Finally, the differences in the pathways of O and H radical production during the streamer propagation and glow phases have been briefly discussed.

ammonia-air flame↗

Cooling performance of an active-passive hybrid composite phase change material (HcPCM) finned heat sink: Constant operating mode

Here, the present study explores a hybrid thermal management technology based on air cooling and hybrid composite phase change material (HcPCM) filled finned heat sink for cooling performance of lower to medium heat flux dissipating electronic devices. Two-dimensional numerical simulations are conducted to study the conjugate heat transfer effects of three types of finned heat sink: air-cooled finned heat sink, HcPCM-cooled finned heat sink, and hybrid (air-HcPCM) cooled finned heat sink. A heat sink with a constant volume faction of plate-fins is designed in all cases and simultaneous effects of hybrid nanoparticles and air are investigated to keep the heat sink temperature at safe operating conditions between 40–60°C. The effect of air is incorporated into the heat sink by applying the convective heat transfer coefficient of h c = 10–100W /m 2 .K which tends to create the natural convection and forced convection heat transfer characteristics. The heat flux is varied from 25–40kW /m 2 in the current study. The hybrid nanoparticles of carbon additives (GO and MWCNTs) are dispersed into the RT-35HC, used as a PCM, with a volume fraction of 0% to 6%. Transient simulations are carried out using COMSOL Multiphysics to solve the governing equations for PCM based conjugate heat transfer model. The results showed that forced convection heat transfer improved the cooling performance of the hybrid heat sink compared to natural convection heat transfer. The addition of nanoparticles further enhanced thermal enhancement and uniform melting distribution of PCM inside the finned heat sink. The h c between 30 to 50W/m 2 .K shows optimized values for forced convection heat transfer operating conditions. The volume fraction of 2% of GO+MWCNTs nanoparticles in recommended or optimum concentration for uniform melting of PCM inside the finned heat sink.

36 MATERIALS SCIENCE↗

Energy Characterization of the UMD Electron Accelerator

Multiple different groups at Los Alamos National Laboratory have a cited interest in understanding the effects of ionizing radiation in materials. Of especial interest to the ISR division is understanding the effects in materials that are used in the space environment, where radiation is omnipresent and impossible to fully shield from. Development is under way to simulate the effects of radiation damage in detectors and electronics, primarily with particle transport codes such as Geant4 and RAMPART, a code that provides a user friendly interface to Geant4. However, these simulations need benchmarking to experimental tests, some of which have been performed at a variety of small electron accelerator facilities. One of these facilities is the Linear Electron Accelerator at the University of Maryland Radiation Facilities. This facility can provide electron beams on the order of 100 mA delivered in tight macropulse bunches, which allows for experiments where the behavior of electronics is monitored as a function of pulse number. While the beam rates per macropulse of this accelerator is monitored during an experiment via a current pick-off connected to a large graphite beamstop, the energy distribution of this electron beam is not measured during the experiment and has not been well characterized. The energy distribution of the beam is pertinent for these types of experiments, since dose is directly related to energy deposited in a medium. Further, these experiments are performed in open air and sample placement in the room can dictate what dose rate they receive per macropulse. We have created simulated dosemaps for the UMD vault room, but these simulations assume a given energy of the beam, both for calculating the dose and to simulate how much the beam diverges in the air. There have been two attempts now to characterize this accelerator beam energy with a plate type spectrometer, including the most recent experiment earlier this year to quantify the beam energy distribution. These results will also be compared to a dosimetry experiment performed by collaborators at the UMD facility, the details of which are described in Reference. Understanding the electron beam energy is important for future experiments, because the energy directly impacts the dose and radiation damage given to a sample and is important for properly simulating the electron beam when comparing simulation to experiment.

43 PARTICLE ACCELERATORS↗

Oxidation of porous stainless steel supports for metal-supported solid oxide fuel cells

Here, oxidation behavior of porous P434L ferritic stainless steel, used for the fabrication of metal-supported solid oxide fuel cells (MS-SOFC), is studied under anodic and cathodic atmospheres. Temperature- and atmosphere-dependence is determined for as-sintered and pre-oxidized stainless steel. Pre-oxidation reduced the long-term oxidation rate. For pre-oxidized samples, the oxidation rate in air exceeds that in humid hydrogen for temperatures above 700 °C. The influence of PrOx, LSCF-SDC, and Ni-SDC coatings is also examined. The coatings do not dramatically impact oxide scale growth. Oxidation in C-free and C-containing anodic atmospheres is similar. Addition of CO 2 , CH 4 , and CO to humidified hydrogen to simulate ethanol reformate does not significantly impact oxidation behavior. Cr transpiration in humid air is greatly reduced by the PrOx coating, and a PrCrO 3 reaction product is observed throughout the porous structure. A dense and protective chromia-based scale forms on steel samples during oxidation in all conditions. A thin silica enriched oxide layer also forms at the metal-scale interface. In general, the oxidation behavior at 700 °C is found to be acceptable.

25 ENERGY STORAGE↗

Validating simulated mountain wave impacts on hub-height wind speed using SoDAR observations

The ascent of stably stratified air over a mountain barrier can trigger the generation of mountain waves. Mountain waves occur frequently over the Columbia River Gorge in western North America and can impact wind power generation over the area. Therefore, predicting the details of mountain waves events (e.g., dominant wavelength, timing, and duration) can be very valuable for the wind energy community. In this study, the ability of the Weather Research and Forecasting (WRF) model to simulate mountain waves and their impact on hub-height wind speed is investigated. Our results suggest that the WRF model has moderate skill in simulating observed mountain wave. Further, given WRF predictions of wavelength range and wave period, the Fast Fourier Transform can calculate the simulated mountain wave impact on hub-height wind speed. The resulting wind speeds agree well with SoDAR observations in terms of both magnitude and pattern. Finally, for the simulated cases, WRF consistently predicts impacts of significant mountain wave events about an hour earlier than the actual observations. Overall, in this work, the sensitivities as well as uncertainties associated with our methodology are discussed in detail.

17 WIND ENERGY↗

Water balance response of permafrost-affected watersheds to changes in air temperatures

Abstract Observations show increases in river discharge to the Arctic Ocean especially in winter over the last decades but the physical mechanisms driving these changes are not yet fully understood. We hypothesize that even in the absence of a precipitation increase, permafrost degradation alone can lead to increased annual river runoff. To test this hypothesis we perform 12 millennium-long simulations over an idealized hypothetical watershed (1 km 2 ) using a distributed, physically based water balance model (Water flow and Balance Simulation Model, WaSiM). The model is forced by both a hypothetical warming defined by an air temperature increase of 7.5 ∘ C over 100 years, and a corresponding cooling scenario. To assess model sensitivity we vary soil saturated hydraulic conductivity and lateral subsurface flow configuration. Under the warming scenario, changes in subsurface water transport due to ground temperature changes result in a 7%–14% increase in annual runoff accompanied by a 6%–20% decrease in evapotranspiration. The increase in runoff is most pronounced in winter. Hence, the simulations demonstrate that changes in permafrost characteristics due to climate warming and associated changes in evapotranspiration provide a plausible mechanism for the observed runoff increases in Arctic watersheds. In addition, our experiments show that when lateral subsurface moisture transport is not included, as commonly done in global-scale Earth System Models, the equilibrium water balance in response to the warming or cooling is similar to the water balance in simulations where lateral subsurface transport is included. However, the transient changes in water balance components prior to reaching equilibrium differ greatly between the two. For example, for high saturated hydraulic conductivity only when lateral subsurface transport is considered, a period of decreased runoff occurs immediately after the warming. This period is characterized by a positive change in soil moisture storage caused by the soil moisture deficit developed during prior cooling.

Debolskiy, Matvey V. (ORCID:0000000296343627)↗

Summertime Near-Surface Temperature Biases Over the Central United States in Convection-Permitting Simulations

Convection-Permitting Model (CPM) simulations of the Central United States climate for the summer of 2011 are studied to understand the causes of warm biases in 2-m air temperature (T 2m ) and related underestimates of precipitation including that from mesoscale convective systems (MCSs). Based on 10 CPM simulations and 9 coarser-resolution model simulations, we quantify contributions from evaporative fraction (EF) and radiation to the T 2m bias with both types of models overestimating T 2m largely because they underestimate EF. The performance of CPMs in capturing MCS characteristics (frequency, rainfall, propagation) varies. The pre-summer precipitation bias has large correlation with mean summertime T 2m bias but the relationship between summertime MCS mean rainfall bias and T 2m bias is non-monotonic. Analysis of lifting condensation level deficit and convective available potential energy suggests that models with T 2m warm biases and low EF have too dry and stable boundary layers, inhibiting the formation of clouds, precipitation and MCSs. Among the CPMs with differing model formulations (e.g., transpiration, infiltration, cloud macrophysics and microphysics), evidence suggests that altering the land-surface model is more effective than altering the atmospheric model in reducing T 2m biases. In conclusion, these results demonstrate that land-atmosphere interactions play a very important role in determining the summertime climate of the Central United States.

2-m air temperature↗

A multiphase flow model of water droplets dielectrophoretic-induced air dehumidification phenomena

Air humidity in indoor spaces plays a critical role in human comfort and health. Dehumidification systems are used for building humidity controls, but they can take significant energy consumption, especially in geographic locations with high outdoor humidity and warm climates. Consequently, there is a growing demand for innovative dehumidification processes that consume minimal energy. Dielectrophoretic air dehumidification represents one such promising approach. However, it has not garnered significant attention due to the absence of engineering models and simulation tools capable of evaluating its performance and limitations at large-scale airflows. A new numerical multiphase CFD model, which is also experimentally validated, is developed in a customized Reacting Foam solver based on OpenFOAM® version 9. The newly developed model seeks to decrease substantial energy consumption and lower costs by leveraging the dielectrophoretic phenomenon to regulate moisture levels in the air. The solver integrates a hybrid Eulerian-Lagrangian framework to track the droplet's trajectory and growth rate while solving the continuum equations for the moist air. An electrospray produces electrically charged droplets, which grow during their in-flight trajectories as water vapor condenses onto their surfaces. The role of electrostatic forces in promoting vapor condensation within a high-gradient electrical field is investigated, and the dielectrophoretic vapor nucleation process on charged water droplets is discussed. The CFD model was validated against results from the literature and from proof-of-concept experiments conducted by the authors, which showed a 2 % air dehumidification with a single electrospray and airflow rate of 5 cubic feet per minute. The simulation results indicated that augmenting the number of electrically charged spray droplets increased the dehumidification of the air to 25 %. The initial mean droplet diameter, the orientation of the injector and relative humidity significantly influence the assessment of dehumidification. As a result, scaling up this approach to larger airflow volumes is identified as a potential future research direction.

42 ENGINEERING↗

Simulating Mixed‐Phase Open Cellular Clouds Observed During COMBLE: Evaluation of Parameterized Turbulence Closure

Marine cold-air outbreaks, or CAOs, are airmass transformations whereby relatively cold boundary layer (BL) air is transported over relatively warm water. To more deeply understand BL and mixed-phase cloud properties during CAO conditions, the Cold-Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) took place from late 2019 into early 2020. During COMBLE, the U.S. Department of Energy's first Atmospheric Radiation Measurement Mobile Facility (AMF1) was deployed to Andenes, Norway, far downstream (~1,000 km) from the Arctic pack ice. This study examines the two most intense CAOs sampled at the AMF1 site. The observed BL structures are open cellular with high (~3–5 km) and cold (–30 to –50 °C) cloud tops, and they often have pockets of high liquid water paths (LWPs; up to ~1,000 g m –2 ) associated with strong updrafts and enhanced turbulence. We use a high-resolution mesoscale model to explore how well four turbulence closure methods represent open cellular clouds. After applying a radar simulator to model outputs for direct evaluation, cloud top properties agree well with AMF1 observations (within ~10%), but radar reflectivity and LWP agreement is more variable. Results suggest that the turbulent Prandtl number may play an important role for the simulated BL and cloud properties. All simulations produce enhanced precipitation rates that are well-correlated with a cloud transition. Finally, the eddy-diffusivity/mass-flux approach produces the deepest cloud layer and therefore the largest and most coherent cellular structures. Furthermore, we recommend the use of a non-local turbulence closure approach to better capture turbulent processes in intense CAOs.

54 ENVIRONMENTAL SCIENCES↗

Measured influence of overhead HVAC on exposure to airborne contaminants from simulated speaking in a meeting and a classroom

Tracer gas experiments were conducted in a 158 m 3 room with overhead supply diffusers to study dispersion of contaminants from simulated speaking in physically distanced meeting and classroom configurations. The room was contained within a 237 m 3 cell with open plenum return to the HVAC system. Heated manikins at desks and a researcher operating the tracer release apparatus presented 8-9 thermal plumes. Experiments were conducted under conditions of no forced air and neutral, cooled, or heated air supplied at 980-1100 cmh, and with/out 20% outdoor air. CO 2 was released at the head of one manikin in each experiment to simulate small (< 5 µm diameter) respiratory aerosols. The metric of exposure relative to perfectly mixed (ERM) is introduced to quantify impacts, based on measurements at manikin heads and at three heights in the center and corners of the room. Chilled or neutral supply air provided good mixing with ERMs close to one. Thermal stratification during heating produced higher ERMs at most manikins: 25% were ≥ 2.5 and the highest were > 5× perfectly mixed conditions. Operation of two within-zone air cleaners together moving ≥400 cmh vertically in the room provided enough mixing to mitigate elevated exposure variations.

54 ENVIRONMENTAL SCIENCES↗

Enhancement of phase change material hysteresis model: A case study of modeling building envelope in EnergyPlus

Nowadays, buildings are expected to offer demand side services to the power grid to enhance the electrical load flexibility, which leads to the concepts of grid-interactive efficient buildings (GEBs). Phase change material (PCM)-based thermal energy storage has seen increasing attention in recent years for peak load shifting of grid-interactive efficient buildings (GEBs). Numerical models are critical tools for design and evaluation of PCM-integrated systems. Most industrial-grade PCMs are reported to melt/freeze over a temperature range instead of at a unique temperature. Such thermal hysteresis effect significantly affects the reliability of simulation results because not only the heat transfer process depends on melting and freezing temperatures, the PCM thermal properties change significantly during the phase change process as well. This study is aimed to develop a model for the PCMs used in the building envelope with the capability to accurately simulate hysteretic behaviors. Further, this model is based on a two-phase assumption and is implemented in a whole building energy performance simulation program (i.e., EnergyPlus). A comparison between numerical results and experimental data shows that during a complete phase transition, the two-phase model could achieve a good agreement with the experimental data. During a partial phase transition, the two-phase model could lead to significant improvements compared to other alternative PCM models, including the existing PCM model in EnergyPlus. Last, whole building simulations were performed to study this model's performance regarding heating/cooling loads and zone mean air temperature of a given building. The results show that the difference in hourly heating/cooling loads introduced by the models was less than 1% in design conditions, while significant changes were observed in both hourly heating/cooling loads and zone mean air temperature when the PCM envelope underwent partial phase transition processes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Impact of heating and cooling loads on battery energy storage system sizing in extreme cold climates

Efficient operation of battery energy storage systems requires that battery temperature remains within a specific range. Current techno-economic models neglect the parasitic loads heating and cooling operations have on these devices, assuming they operate at constant temperature. In this work, these effects are investigated considering the optimal sizing of battery energy storage systems when deployed in cold environments. Here, a peak shaving application is presented as a linear programming problem which is then formulated in the PYOMO optimization programming language. The building energy simulation software EnergyPlus is used to model the heating, ventilation, and air conditioning load of the battery energy storage system enclosure. Case studies are conducted for eight locations in the United States considering a nickel manganese cobalt oxide lithium ion battery type and whether the power conversion system is inside or outside the enclosure. The results show an increase of 42% to 300% in energy capacity size, 43% to 217% in power rating, and 43% to 296% increase in capital cost dependent on location. This analysis shows that the heating, ventilation, and air conditioning load can have a large impact on the optimal sizes and cost of a battery energy storage system and merit consideration in techno-economic studies.

25 ENERGY STORAGE↗

Hyporheic Exchange in Sand Dunes Under a Freely Deforming River Water Surface

Bedform-driven hyporheic exchange is conditioned by the head gradients at the sediment-water interface. Local exchange phenomena between the surface and the subsurface are often driven by the dynamic forces related to the velocity distributions around sediment waves. In open channels, the static forces are represented by the water depth, yet most computational fluid dynamics models use a rigid-lid approximation. Here, we investigated whether and when the deformation of the river’s free-surface influences bedform-driven hyporheic exchange. This was done through simulations of coupled open channel and hyporheic flows with the air-water interface modeled either as a free-surface or a rigid-lid across increasing subcritical Froude numbers from ~0.05 to ~0.95. The normalized hyporheic flux was higher when considering free-surface deformation across most of the range of Froude numbers considered. When the Froude number was larger than 0.6 and smaller than 0.85, both hydrostatic and nonhydrostatic-driven fluxes increased significantly compared to the rigid-lid approach and the total flux was about five times that predicted with the rigid-lid. These results indicate that bedform-driven hyporheic fluxes are underestimated in most studies that use a rigid-lid assumption in subcritical flows.

58 GEOSCIENCES↗