The biological performance assessment (BioPA) toolset developed by Pacific Northwest National Laboratory (PNNL) estimates the relative biological performance of fish passage at a hydroelectric power turbine unit. The tool is based on the use of computational fluid dynamics (CFD) and fish biological response relationships. The recent release, BioPA-v3, is based on directly computed trajectory and collision of material Lagrangian particles using CFD simulation codes rather than the prior version that relies on Tecplot to compute streamtrace trajectories. Before modifying the toolset, a series of validation tests were performed at the various steps of modification in the toolset. Validation is a critical step of any numerical investigation that reflects the accuracy and reliability of the predicted results. It raises the confidence level of the user to use the modified version of the BioPA toolset. Several test cases were simulated and compared, where available, to observed data. The trajectory and collision of the small spherical and cylindrical particles in a water flume were compared to in-house experiments. The CFD predicted collision rate and flow field compared well with experimental observation for vane array and large cylinder as target bodies. Next, the CFD-predicted flow field and hydraulic performance of a laboratory-scale model of a Francis turbine was also successfully validated. Note that the trajectory of the particles is significantly affected by the flow field in such extreme conditions. In addition to the particle trajectories and flow field, the collision detection method employed in the CFD simulations was also successfully validated. The CFD predicted impact velocity, collision time, velocity, and trajectory of a sphere excellently matched with analytical value for a bouncing ball in the elastic collision. A similar approach was also tested and successfully validated for a collision of sphere with a 45? inclined plane. After successfully validating different cases, the BioPA toolset was modified to use direct output of the CFD prediction and the new version can be used in evaluating biological performance at hydroelectric turbines.

The systematic response of coastal ecosystems to inundation and salinity exposure is fundamental to their ecology and biogeochemical function. Here we observe and model freshwater‐seawater interactions in a first‐order stream—floodplain system where tidal access was recently restored. Subsurface flow and transport modeling were used to quantify and better understand the interplay of processes, properties, and conditions that control water level and salinity in the floodplain to the tidal stream. Water levels in the stream were highly correlated with tidal forcing, which resulted in episodic inundation of the floodplain at quasi‐monthly frequency. The tidal stream is the only source of salinity to the floodplain, yet shallow groundwater salinity was considerably higher than average stream salinity. The low‐permeability clay floodplain soils limit lateral groundwater flow and transport, resulting in floodplain groundwater and salinity dynamics driven almost exclusively by infiltration during inundation events. As inundation occurs during high tide, estuarine waters reach the floodplain with minor attenuation in salinity from the stream's freshwater discharge. Infiltration and salinity exposure are topography controlled and regulated by ponding depth and duration, seasonal ground saturation, and depth to water table. The model suggests that floodplain salinity is currently in an early stage of transition from pre‐restoration freshwater conditions and will not reach equilibrium for ~20 years. These findings have broad relevance for understanding how and over what time scales coastal ecosystems will respond to increasing seawater exposure from sea level rise, ocean‐originating storms, and changes in natural and man‐made barriers.

This report documents a numerical modeling study conducted to assess the pressure environment and relative biological performance of the existing turbine runner design at the U.S. Army Corps of Engineers (USACE) McNary Dam to establish a baseline for comparison to planned replacement runner designs. Computational fluid dynamics (CFD) models of the turbine unit constructed for this study were analyzed using a streamtrace-based technique called the Biological Performance Assessment (BioPA), which was developed at the Pacific Northwest National Laboratory (PNNL).

For the past two decades a slowly growing consensus has been that rising temperature is being caused by increases in carbon dioxide (CO2) in the Earth's atmosphere. What is less consensual is what to do about it. Solutions range from reducing emissions to various ways to capture CO2 (from effluents of power plants, as plants, etc.) and store it away from the atmosphere. The problem is of such magnitude that it is likely that we will need many methods. Also some methods would cause draconian harm to global economies. This work is about the storage part of the reducing of CO2 in the atmosphere. There have been many methods proposed for storing CO2: in vegetation, at the bottom of oceans, in land farms, and in subsurface formations. The latter is particularly prominent given the large estimated volume in porous formation in the Earth, and the 50-year old experience in injecting CO2 for enhanced oil recovery. This report is about the science behind subsurface storage, called here geological carbon storage or GCS. Effective GCS requires that billions of metric tons of CO2 be injected into storage reservoirs annually. GCS feasibility requires balancing the opposing impacts of large emplacement (injection) rates without compromising caprock seal integrity during injection. Once injected, the challenge is to assure the injected CO2 remains permanently stored. Under most conditions, CO2 is in a supercritical state, or scCO2. scCO2 is immiscible with the subsurface brine and has a factor of ten smaller viscosity, leading to the development of two-phase fluid flow and scCO2 bypassing. After the injection of large volumes of scCO2, the pore pressure, state of stress, chemical, thermal, and biological steady-state or equilibrium conditions in the subsurface are disrupted, which may lead to unpredictable behavior because the responses to these changes are often non-linear. Current challenges for efficient and reliable scCO2 storage are: (1) sustaining large injection rates; (2) using available pore space efficiently, and (3) controlling undesired or unexpected behavior, such as scCO2 leakage. Various aspects of these challenges are addressed by collaborative projects within Center for Frontiers in Subsurface Energy Security (CFSES). Challenge 1 Sustaining large storage rates To offset the annual CO2 emissions of in the United States, ~7 billion metric tons of CO2 must be injected into the storage reservoirs annually. Injected CO2 forms a plume that starts spreading, displacing originally resident brine, though inefficiently so. Because of its low viscosity, and buoyant forces, a CO2 plume rises to the reservoir rock-caprock interface. Injecting large volumes of scCO2 at sufficiently high rates into subsurface storage formations leads to increases in pore pressures, potential expansion of the reservoir rock, and may result in fracturing of the reservoir rock and/or caprock. CFSES research uses existing and new experimental and modeling approaches to identify chemical-physical controls on the permeability and pore pressure dynamics, anticipating geomechanical (fracturing) events, predicting multi-phase flow patterns and trapping of CO2. Challenge 2 – Using pore space with unprecedented efficiency Since large volumes of scCO2 are to be stored, the volume available in storage reservoirs must be used efficiently. Storage efficiency is defined as a fraction of pore space occupied by carbon dioxide. Current estimates are that less than 5 percent of available pore volume is available for storage because scCO2, having a smaller density and viscosity than reservoir brine, tends to form viscous fingers, and bypass available pore space. This flow regime results in a CO2 plume spreading over large areas in subsurface. The research goal in Challenge 2 is to advance the fundamental understanding of two-phase flow through porous media, including molecular-to meter- scale experimental and modeling approaches. This work leads to novel strategies to ensure that injected CO2 occupies more than 5% of the pore space, ideally up to 50%. Challenge 3 – Controlling undesired or unexpected behavior Following the injection of scCO2, the subsurface reservoir re-equilibrates to different physical and chemical states. This re-equilibration in space and time could lead to unexpected and/or undesired behavior: fracturing of reservoir rock and/or caprock, development of preferential flow paths, CO2 leakage (to the surface or overlaying fresh water reservoirs), induced seismicity, and activation of faults. One of the main research goals of CFSES was to advance the fundamental understanding of how individual chemical and physical processes are coupled. We particularly focused on the chemical-mechanical coupling, and flow-mechanical coupling. We collected experimental and modeling data to determine the temporal and spatial scales of the coupled processes, which may lead to pore collapse caused by creep, activation of fractures and faults, chemically-induced fracture propagation, and self-focusing of scCO2 flow.

Columbia River water elevations and flows in the Hanford Reach affect the environment and facilities along the shoreline, including movement of contaminants in groundwater, fish habitat, and infrastructure subject to flooding. This report describes the hydraulic simulation of hypothetical flood flows using the best available topographic and bathymetric data for the Hanford Reach and the Modular Aquatic Simulation System in 1 Dimension (MASS1) hydrodynamic model. The MASS1 model of the Hanford Reach was previously calibrated to field measurements of water surface elevations. The current model setup can be used for other studies of flow, water levels, and temperature in the Reach. The existing MASS1 channel geometry and roughness and other model configuration inputs for the Hanford Reach were used for this study, and previous calibration and validation results for the model are reprinted here for reference. The flood flows for this study were simulated by setting constant flow rates obtained from the U.S. Army Corps of Engineers (USACE) for the Columbia, Snake, and Yakima Rivers, and a constant water level at McNary Dam, and then running the model to steady state. The discharge levels simulated were all low-probability events; for example, a 100-year flood is one that would occur on average every 100 years, or put another way, in any given year there is a 1% chance that a discharge of that level or higher will occur. The simulated floods and their corresponding Columbia River discharges were 100-year (445,000 cfs), 500-year (520,000 cfs), and the USACE-defined Standard Project Flood (960,000 cfs). The resulting water levels from the steady-state floods can be viewed as “worst case” outcomes for the respective discharge levels. The MASS1 output for water surface elevations was converted to the North American Vertical Datum of 1988 and projected across the channel and land surface to enable mapping of the floodplain for each scenario. Floodplain maps show that for the 100-year and 500-year discharge levels, flooding is mainly confined to the topographic trench that is the river channel. The flooded area for the Standard Project Flood extends out of the channel area in some places, particularly in the 100-F Area. All of the output from the simulations have been archived and are available for future investigations in the Hanford Reach.

This report contains summarized progress on the construction, validation, and calibration of a hydrodynamic and water temperature model, and individual-based sturgeon model, for the Transboundary Reach of the Columbia River and Lake Roosevelt. A 1-dimensional hydrodynamic and water temperature model, the Modular Aquatic Simulation System in 1-dimension (MASS1), was constructed using existing bathymetry data from the U.S.-Canada international border to Grand Coulee Dam, and validated using velocity, temperature, and water elevation data collected by staff of the CTCR during three distinct hydrodynamic periods occurring in early spring, late spring, and summer in the study area. Parameters of MASS1 were then calibrated so that model simulations matched empirical data of water surface elevations collected in 2016. Several years of MASS1 data were archived and are available for use with the sturgeon IBM: 1975–2001, 2004, 2005, 2007–2015, and April–August 2016. An individual-based simulation model was concomitantly constructed to simulate four of the early life stages of white sturgeon: spawning, embryo incubation, free embryos, and early larvae. Submodels for each of the four life stages contain mathematical algorithms that primarily describe the growth, development, and movement of individuals based on outputs from MASS1. Sturgeon simulations were then run for two historical years where limited recruitment may have occurred (1997 and 2011) and two years in which recruitment likely did not occur (2004 and 2005).

Accurate modeling of the velocity field in the forebay of a hydroelectric power station is important for both power generation and fish passage, and is able to be increasingly well represented by computational fluid dynamics (CFD) simulations. Acoustic Doppler Current Profiler (ADCP) are investigated herein as a method of validating the numerical flow solutions, particularly in observed and calculated regions of non -homogeneous flow velocity. By using a numerical model of an ADCP operating in a velocity field calculated using CFD, the errors due to the spatial variation of the flow velocity are quantified. The numerical model of the ADCP is referred to herein as a Virtual ADCP (VADCP).Two applications of the VADCP are modeled in the numerical analyses presented. Firstly the virtual measurement error of the VADCP is calculated for a single instrument adjacent to the short converging intake of a powerhouse. Secondly, the flow discharge through the forebay is estimated from a transect of VADCP instruments at different distances from the powerhouse. The influence of instrument location and orientation are investigated for both cases.A velocity error of up to 94% of the reference velocity is calculated for a VADCP modeled adjacent to an operating intake and is shown to decrease with distance from the powerhouse. Qualitative agreement is observed between the calculated VADCP velocities and reference velocities by a horizontal offset distance of 18 m upstream of the powerhouse. (C) 2015 Elsevier Ltd. All rights reserved.

Four sets of nonreactive solute transport experiments were conducted with micromodels. Each set consisted of three experiments with one variable, i.e., flow velocity, grain diameter, pore-aspect ratio, and flow-focusing heterogeneity. The data sets were offered to pore-scale modeling groups to test their numerical simulators. Each set consisted of two learning experiments, for which all results were made available, and one challenge experiment, for which only the experimental description and base input parameters were provided. The experimental results showed a nonlinear dependence of the transverse dispersion coefficient on the Peclet number, a negligible effect of the pore-aspect ratio on transverse mixing, and considerably enhanced mixing due to flow focusing. Five pore-scale models and one continuum-scale model were used to simulate the experiments. Of the pore-scale models, two used a pore-network (PN) method, two others are based on a lattice Boltzmann (LB) approach, and one used a computational fluid dynamics (CFD) technique. The learning experiments were used by the PN models to modify the standard perfect mixing approach in pore bodies into approaches to simulate the observed incomplete mixing. The LB and CFD models used the learning experiments to appropriately discretize the spatial grid representations. For the continuum modeling, the required dispersivity input values were estimated based on published nonlinear relations between transverse dispersion coefficients and Peclet number. Comparisons between experimental and numerical results for the four challenge experiments show that all pore-scale models were all able to satisfactorily simulate the experiments. The continuum model underestimated the required dispersivity values, resulting in reduced dispersion. The PN models were able to complete the simulations in a few minutes, whereas the direct models, which account for the micromodel geometry and underlying flow and transport physics, needed up to several days on supercomputers to resolve the more complex problems.

Pore-scale models are useful for studying relationships between fundamental processes and phenomena at larger (i.e., Darcy) scales. However, the size of domains that can be simulated with explicit pore-scale resolution is limited by computational and observational constraints. Direct numerical simulation of pore-scale flow and transport is typically performed on millimeter-scale volumes at which X-ray computed tomography (XCT), often used to characterize pore geometry, can achieve micrometer resolution. In contrast, laboratory experiments that measure continuum properties are typically performed on decimeter-scale columns. At this scale, XCT resolution is coarse (tens to hundreds of micrometers) and prohibits characterization of small pores and grains. We performed simulations of pore-scale processes over a decimeter-scale volume of natural porous media with a wide range of grain sizes, and compared to results of column experiments using the same sample. Simulations were conducted using high-performance codes executed on a supercomputer. Two approaches to XCT image segmentation were evaluated, a binary (pores and solids) segmentation and a ternary segmentation that resolved a third category (porous solids with pores smaller than the imaged resolution). We used a multiscale Stokes-Darcy simulation method to simulate the combination of Stokes flow in large open pores and Darcy-like flow in porous solid regions. Flow and transport simulations based on the binary segmentation were inconsistent with experimental observations because of overestimation of large connected pores. Simulations based on the ternary segmentation provided results that were consistent with experimental observations, demonstrating our ability to successfully model pore-scale flow over a column-scale domain.

Rapid pressure changes in hydroelectric turbine flows can cause barotrauma that can be hazardous to the passage of fish, in particular migratory juvenile salmonids. Although numerous laboratory tests have evaluated the effect of rapid decompression in fish species of relevance, numerical modeling studies offer the advantage of predicting, for new turbine designs, the potential risks of mortality and injury from rapid pressure change during turbine passage. However, rapid pressure change is only one of several hydraulic risks encountered by fish during turbine passage in addition to blade strike, shear, and turbulence. To better understand the role of rapid pressure changes, the present work focuses on the application of a computational fluid dynamics based method for evaluating the risk of pressure-related mortality to fish passing through an early 1960s era original hydroelectric Kaplan turbine at Wanapum Dam (Columbia River, Washington), and a modern advanced Kaplan turbine installed in 2005. The results show that the modeling approach acceptably reproduced the nadir pressure distributions compared to field data previously collected at the site using an autonomous sensor. Our findings show that the new advanced-design unit performs better, in terms of reduced barotrauma risk to fish from exposure to low pressures, than the original turbine unit. The outcomes allow for comparative analyses of turbine designs and operations prior to installation, an advantage that can potentially be integrated in the process of designing new turbine units to achieve superior environmental performance. Overall, the results show that modern turbine designs can achieve the multiple objectives of increasing power generation, lowering cavitation potential, and reducing barotrauma risks to passing fish.

Public Utility District No. 2 of Grant County (GCPUD) operates the Priest Rapids Dam (PRD), a hydroelectric facility on the Columbia River in Washington State. The dam contains 10 Kaplan-type turbine units that are now more than 50 years old. Plans are underway to refit these aging turbines with new runners. The Columbia River at PRD is a migratory pathway for several species of juvenile and adult salmonids, so passage of fish through the dam is a major consideration when upgrading the turbines. In this paper, a method for turbine biological performance assessment (BioPA) is demonstrated. Using this method, a suite of biological performance indicators is computed based on simulated data from a CFD model of a proposed turbine design. Each performance indicator is a measure of the probability of exposure to a certain dose of an injury mechanism. Using known relationships between the dose of an injury mechanism and frequency of injury (dose–response) from laboratory or field studies, the likelihood of fish injury for a turbine design can be computed from the performance indicator. By comparing the values of the indicators from proposed designs, the engineer can identify the more-promising alternatives. We present an application of the BioPA methodmore » for baseline risk assessment calculations for the existing Kaplan turbines at PRD that will be used as the minimum biological performance that a proposed new design must achieve.« less

We introduce a method for hydro turbine biological performance assessment (BioPA) to bridge the gap between field and laboratory studies on fish injury and turbine engineering design. Using this method, a suite of biological performance indicators is computed based on simulated data from a computational fluid dynamics (CFD) model of a proposed hydro turbine design. Each performance indicator is a measure of the probability of exposure to a certain dose of an injury mechanism. If the relationship between the dose of an injury mechanism (stressor) and frequency of injury (dose–response) is known from laboratory or field studies, the likelihood of fish injury for a turbine design can be computed from the performance indicator. By comparing the values of the indicators from various turbine designs, engineers and biologists can identify the more-promising designs and operating conditions to minimize hydraulic conditions hazardous to passing fish. In this paper, the BioPA method is applied to estimate barotrauma induced mortal injury rates for Chinook salmon exposed to rapid pressure changes in Kaplan-type hydro turbines. Following the description of the general method, application of the BioPA to estimate the probability of mortal injury from exposure to rapid decompression is illustrated using a Kaplan hydro turbine at the John Day Dam on the Columbia River in the Pacific Northwest region of the USA. The estimated rates of mortal injury increased from 0.3% to 1.7% as discharge through the turbine increased from 334 to 564m3/s for fish assumed to be acclimated to a depth of 5m. The majority of pressure nadirs occurred immediately below the runner blades, with the lowest values in the gap at the blade tips and just below the leading edge of the blades. Such information can help engineers focus on problem areas when designing new turbine runners to be more fish-friendly than existing units.

Over the past two decades, there have been many studies describing injury mechanisms associated with turbine passage, the response of various fish species to these mechanisms, and the probability of survival through dams. Although developing tools to design turbines that improve passage survival has been difficult and slow, a more robust quantification of the turbine environment has emerged through integrating physical model data, fish survival data, and computational fluid dynamics (CFD) studies. Grant County Public Utility District (GCPUD) operates the Priest Rapids Dam (PRD), a hydroelectric facility on the Columbia River in Washington State. The dam contains 10 Kaplan-type turbine units that are now almost 50 years old. The Utility District plans to refit all of these aging turbines with new turbines. The Columbia River at PRD is a migratory pathway for several species of juvenile and adult salmonids, so passage of fish through the dam is a major consideration when replacing the turbines. In this presentation, a method for turbine biological performance assessment (BioPA) is introduced. Using this method, a suite of biological performance indicators is computed based on simulated data from a CFD model of a proposed turbine design. Each performance indicator is a measure of the probabilitymore » of exposure to a certain dose of an injury mechanism. Using known relationships between the dose of an injury mechanism and frequency of injury (dose–response) from laboratory or field studies, the likelihood of fish injury for a turbine design can be computed from the performance indicator. By comparing the values of the indicators from proposed designs, the engineer can identify the more-promising alternatives. We will present application of the BioPA method for baseline risk assessment calculations for the existing Kaplan turbines at PRD that will be used as the minimum biological performance that a proposed new design must achieve.« less

This Handbook of Data Reduction Procedures, Workbooks, and Exchange Templates is designed to support the Oncor geodatabase for the Columbia Estuary Ecosystem Restoration Program (CEERP). The following data categories are covered: water-surface elevation and temperature, sediment accretion rate, photo points, herbaceous wetland vegetation cover, tree plots and site summaries, fish catch and density, fish size, fish diet, fish prey, and Chinook salmon genetic stock identification. The handbook is intended for use by scientists collecting monitoring and research data for the CEERP. The ultimate goal of Oncor is to provide quality, easily accessible, geospatial data for synthesis and evaluation of the collective performance of CEERP ecosystem restoration actions at a program scale.

Pacific Northwest National Laboratory (PNNL) conducted this project for the U.S. Army Corps of Engineers, Portland District (Corps). The purpose of the project is to develop a geospatial, web-accessible database (called “Oncor”) for action effectiveness and related data from monitoring and research efforts for the Columbia Estuary Ecosystem Restoration Program (CEERP). The intent is for the Oncor database to enable synthesis and evaluation, the results of which can then be applied in subsequent CEERP decision-making. This is the first annual report in what is expected to be a 3- to 4-year project, which commenced on February 14, 2012.