Catchment characteristics and soil and water chemistry of the San Francisco basin and subbasins, Ecuador

In this study, the Mean Transit Time and Mixing Model Analysis methods are combined to unravel the runoff generation process of the San Francisco River basin (73.5 km**2) situated on the Amazonian side of the Cordillera Real in the southernmost Andes of Ecuador. The montane basin is covered with cloud forest, sub-páramo, pasture and ferns. Nested sampling was applied for the collection of streamwater samples and discharge measurements in the main tributaries and outlet of the basin, and for the collection of soil and rock water samples. Weekly to biweekly water grab samples were taken at all stations in the period April 2007-November 2008. Hydrometric data, Mean Transit Time and Mixing Model Analysis allowed preliminary evaluation of the processes controlling the runoff in the San Francisco River basin. Results suggest that flow during dry conditions mainly consists of lateral flow through the C-horizon and cracks in the top weathered bedrock layer, and that all subcatchments have an important contribution of this deep water to runoff, no matter whether pristine or deforested. During normal to low precipitation intensities, when antecedent soil moisture conditions favour water infiltration, vertical flow paths to deeper soil horizons with subsequent lateral subsurface flow contribute most to streamflow. Under wet conditions in forested catchments, streamflow is controlled by near surface lateral flow through the organic horizon. Exceptionally, saturation excess overland flow occurs. By absence of the litter layer in pasture, streamflow under wet conditions originates from the A horizon, and overland flow.

Neutron-mapping polymer flow: Scattering, flow visualization, and molecular theory

Flows of complex fluids need to be understood at both macroscopic and molecular scales, because it is the macroscopic response that controls the fluid behavior, but the molecular scale that ultimately gives rise to rheological and solid-state properties. Here the flow field of an entangled polymer melt through an extended contraction, typical of many polymer processes, is imaged optically and by small-angle neutron scattering. The dual-probe technique samples both the macroscopic stress field in the flow and the microscopic configuration of the polymer molecules at selected points. The results are compared with a recent tube model molecular theory of entangled melt flow that is able to calculate both the stress and the single-chain structure factor from first principles. The combined action of the three fundamental entangled processes of reptation, contour length fluctuation, and convective constraint release is essential to account quantitatively for the rich rheological behavior. The multiscale approach unearths a new feature: Orientation at the length scale of the entire chain decays considerably more slowly than at the smaller entanglement length.

Investigating viscous fluid flow in an internal mixer using computational fluid dynamics

This thesis presents an effective methodology for the generation of a simulation which can be used to increase the understanding of viscous fluid processing equipment and aid in their development, design and optimisation. The Hampden RAPRA Torque Rheometer internal batch twin rotor mixer has been simulated with a view to establishing model accuracies, limitations, practicalities and uses. As this research progressed, via the analyses several 'snap-shot' analysis of several rotor configurations using the commercial code Polyflow, it was evident that the model was of some worth and its predictions are in good agreement with the validation experiments, however, several major restrictions were identified. These included poor element form, high man-hour requirements for the construction of each geometry and the absence of the transient term in these models. All, or at least some, of these limitations apply to the numerous attempts to model internal mixes by other researchers and it was clear that there was no generally accepted methodology to provide a practical three-dimensional model which has been adequately validated. This research, unlike others, presents a full complex three-dimensional, transient, non-isothermal, generalised non-Newtonian simulation with wall slip which overcomes these limitations using unmatched ridding and sliding mesh technology adapted from CFX codes. This method yields good element form and, since only one geometry has to be constructed to represent the entire rotor cycle, is extremely beneficial for detailed flow field analysis when used in conjunction with user defined programmes and automatic geometry parameterisation (AGP), and improves accuracy for investigating equipment design and operation conditions. Model validation has been identified as an area which has been neglected by other researchers in this field, especially for time dependent geometries, and has been rigorously pursued in terms of qualitative and quantitative velocity vector analysis of the isothermal, full fill mixing of generalised non-Newtonian fluids, as well as torque comparison, with a relatively high degree of success. This indicates that CFD models of this type can be accurate and perhaps have not been validated to this extent previously because of the inherent difficulties arising from most real processes.

Complementary methods to investigate the development of clogging within a horizontal sub-surface flow tertiary treatment wetland

A combination of experimental methods was applied at a clogged, horizontal subsurface flow (HSSF) municipal wastewater tertiary treatment wetland (TW) in the UK, to quantify the extent of surface and subsurface clogging which had resulted in undesirable surface flow. The three dimensional hydraulic conductivity profile was determined, using a purpose made device which recreates the constant head permeameter test in-situ. The hydrodynamic pathways were investigated by performing dye tracing tests with Rhodamine WT and a novel multi-channel, data-logging, flow through Fluorimeter which allows synchronous measurements to be taken from a matrix of sampling points. Hydraulic conductivity varied in all planes, with the lowest measurement of 0.1 md1 corresponding to the surface layer at the inlet, and the maximum measurement of 1550 md1 located at a 0.4m depth at the outlet. According to dye tracing results, the region where the overland flow ceased received five times the average flow, which then vertically short-circuited below the rhizosphere. The tracer break-through curve obtained from the outlet showed that this preferential flow-path accounted for approximately 80% of the flow overall and arrived 8 h before a distinctly separate secondary flow-path. The overall volumetric efficiencyof the clogged system was 71% and the hydrology was simulated using a dual-path, dead-zone storage model. It is concluded that uneven inlet distribution, continuous surface loading and high rhizosphere resistance is responsible for the clog formation observed in this system. The average inlet hydraulic conductivity was 2 md1, suggesting that current European design guidelines, which predict that the system will reach an equilibrium hydraulic conductivity of 86 md1, do not adequately describe the hydrology of mature systems.

Clogging in horizontal subsurface flow constructed wetlands

Horizontal Subsurface Flow Treatment Wetlands (HSSF TWs) are used by Severn Trent Water as a low-cost tertiary wastewater treatment for rural locations. Experience has shown that clogging is a major operational problem that reduces HSSF TW lifetime. Clogging is caused by an accumulation of secondary wastewater solids from upstream processes and decomposing leaf litter. Clogging occurs as a sludge layer where wastewater is loaded on the surface of the bed at the inlet. Severn Trent systems receive relatively high hydraulic loading rates, which causes overland flow and reduces the ability to mineralise surface sludge accumulations. A novel apparatus and method, the Aston Permeameter, was created to measure hydraulic conductivity in situ. Accuracy is ±30 %, which was considered adequate given that conductivity in clogged systems varies by several orders of magnitude. The Aston Permeameter was used to perform 20 separate tests on 13 different HSSF TWs in the UK and the US. The minimum conductivity measured was 0.03 m/d at Fenny Compton (compared with 5,000 m/d clean conductivity), which was caused by an accumulation of construction fines in one part of the bed. Most systems displayed a 2 to 3 order of magnitude variation in conductivity in each dimension. Statistically significant transverse variations in conductivity were found in 70% of the systems. Clogging at the inlet and outlet was generally highest where flow enters the influent distribution and exits the effluent collection system, respectively. Surface conductivity was lower in systems with dense vegetation because plant canopies reduce surface evapotranspiration and decelerate sludge mineralisation. An equation was derived to describe how the water table profile is influenced by overland flow, spatial variations in conductivity and clogging. The equation is calibrated using a single parameter, the Clog Factor (CF), which represents the equivalent loss of porosity that would reproduce measured conductivity according to the Kozeny-Carman Equation. The CF varies from 0 for ideal conditions to 1 for completely clogged conditions. Minimum CF was 0.54 for a system that had recently been refurbished, which represents the deviation from ideal conditions due to characteristics of non-ideal media such as particle size distribution and morphology. Maximum CF was 0.90 for a 15 year old system that exhibited sludge accumulation and overland flow across the majority of the bed. A Finite Element Model of a 15 m long HSSF TW was used to indicate how hydraulics and hydrodynamics vary as CF increases. It was found that as CF increases from 0.55 to 0.65 the subsurface wetted area increases, which causes mean hydraulic residence time to increase from 0.16 days to 0.18 days. As CF increases from 0.65 to 0.90, the extent of overland flow increases from 1.8 m to 13.1 m, which reduces hydraulic efficiency from 37 % to 12 % and reduces mean residence time to 0.08 days.

Linear stability of flow past two circular cylinders in a side-by-side arrangement

The linear stability of flow past two circular cylinders in a side-by-side arrangement is investigated theoretically, numerically and experimentally under the assumption of a two-dimensional flow field, in order to explore the origin of in-phase and antiphase oscillatory flows. Steady symmetric flow is realized at a small Reynolds number, but becomes unstable above a critical Reynolds number though the solution corresponding to the flow still satisfies the basic equations irrespective of the magnitude of the Reynolds number. We obtained the solution numerically and investigated its linear stability. We found that there are two kinds of unstable modes, i.e., antisymmetric and symmetric modes, which lead to in-phase and antiphase oscillatory flows, respectively. We determined the critical Reynolds numbers for the two modes and evaluated the critical distance at which the most unstable disturbance changes from the antisymmetric to the symmetric mode, or vice versa. ©2005 The Physical Society of Japan.

A finite element approach to modelling the hydrological regime in horizontal subsurface flow constructed wetlands for wastewater treatment

A Finite Element Analysis (FEA) model is used to explore the relationship between clogging and hydraulics that occurs in Horizontal Subsurface Flow Treatment Wetlands (HSSF TWs) in the United Kingdom (UK). Clogging is assumed to be caused by particle transport and an existing single collector efficiency model is implemented to describe this behaviour. The flow model was validated against HSSF TW survey results obtained from the literature. The model successfully simulated the influence of overland flow on hydrodynamics, and the interaction between vertical flow through the low permeability surface layer and the horizontal flow of the saturated water table. The clogging model described the development of clogging within the system but under-predicted the extent of clogging which occurred over 15 years. This is because important clogging mechanisms were not considered by the model, such as biomass growth and vegetation establishment. The model showed the usefulness of FEA for linking hydraulic and clogging phenomenon in HSSF TWs and could be extended to include treatment processes. © 2011 Springer Science+Business Media B.V.

An integrated flow and transport model to study the impact of mercury remediation strategies for East Fork Poplar Creek watershed, Oak Ridge, Tennessee

An integrated flow and transport model using MIKE SHE/MIKE 11 software was developed to predict the flow and transport of mercury, Hg(II), under varying environmental conditions. The model analyzed the impact of remediation scenarios within the East Fork Poplar Creek watershed of the Oak Ridge Reservation with respect to downstream concentration of mercury. The numerical simulations included the entire hydrological cycle: flow in rivers, overland flow, groundwater flow in the saturated and unsaturated zones, and evapotranspiration and precipitation time series. Stochastic parameters and hydrologic conditions over a five year period of historical hydrological data were used to analyze the hydrological cycle and to determine the prevailing mercury transport mechanism within the watershed. Simulations of remediation scenarios revealed that reduction of the highly contaminated point sources, rather than general remediation of the contaminant plume, has a more direct impact on downstream mercury concentrations.

Factors Controlling Surface Water Flow in a Low-gradient Subtropical Wetland

Surface water flow patterns in wetlands play a role in shaping substrates, biogeochemical cycling, and ecosystem characteristics. This paper focuses on the factors controlling flow across a large, shallow gradient subtropical wetland (Shark River Slough in Everglades National Park, USA), which displays vegetative patterning indicative of overland flow. Between July 2003 and December 2007, flow speeds at five sites were very low (s−1), and exhibited seasonal fluctuations that were correlated with seasonal changes in water depth but also showed distinctive deviations. Stepwise linear regression showed that upstream gate discharges, local stage gradients, and stage together explained 50 to 90% of the variance in flow speed at four of the five sites and only 10% at one site located close to a levee-canal combination. Two non-linear, semi-empirical expressions relating flow speeds to the local hydraulic gradient, water depths, and vegetative resistance accounted for 70% of the variance in our measured speed. The data suggest local-scale factors such as channel morphology, vegetation density, and groundwater exchanges must be considered along with landscape position and basin-scale geomorphology when examining the interactions between flow and community characteristics in low-gradient wetlands such as the Everglades.

Movement Effects on the Flow Physics and Nutrient Delivery in Engineered Valvular Tissues

Mechanical conditioning has been shown to promote tissue formation in a wide variety of tissue engineering efforts. However the underlying mechanisms by which external mechanical stimuli regulate cells and tissues are not known. This is particularly relevant in the area of heart valve tissue engineering (HVTE) owing to the intense hemodynamic environments that surround native valves. Some studies suggest that oscillatory shear stress (OSS) caused by steady flow and scaffold flexure play a critical role in engineered tissue formation derived from bone marrow derived stem cells (BMSCs). In addition, scaffold flexure may enhance nutrient (e.g. oxygen, glucose) transport. In this study, we computationally quantified the i) magnitude of fluid-induced shear stresses; ii) the extent of temporal fluid oscillations in the flow field using the oscillatory shear index (OSI) parameter, and iii) glucose and oxygen mass transport profiles. Noting that sample cyclic flexure induces a high degree of oscillatory shear stress (OSS), we incorporated moving boundary computational fluid dynamic simulations of samples housed within a bioreactor to consider the effects of: 1) no flow, no flexure (control group), 2) steady flow-alone, 3) cyclic flexure-alone and 4) combined steady flow and cyclic flexure environments. We also coupled a diffusion and convention mass transport equation to the simulated system. We found that the coexistence of both OSS and appreciable shear stress magnitudes, described by the newly introduced parameter OSI-t , explained the high levels of engineered collagen previously observed from combining cyclic flexure and steady flow states. On the other hand, each of these metrics on its own showed no association. This finding suggests that cyclic flexure and steady flow synergistically promote engineered heart valve tissue production via OSS, so long as the oscillations are accompanied by a critical magnitude of shear stress. In addition, our simulations showed that mass transport of glucose and oxygen is enhanced by sample movement at low sample porosities, but did not play a role in highly porous scaffolds. Preliminary in-house in vitro experiments showed that cell proliferation and phenotype is enhanced in OSI-t environments.

An Integrated Flow and Transport Model to Study the Impact of Mercury Remediation Strategies for East Fork Poplar Creek Watershed, Oak Ridge, Tennessee

An integrated flow and transport model using MIKE SHE/MIKE 11 software was developed to predict the flow and transport of mercury, Hg(II), under varying environmental conditions. The model analyzed the impact of remediation scenarios within the East Fork Poplar Creek watershed of the Oak Ridge Reservation with respect to downstream concentration of mercury. The numerical simulations included the entire hydrological cycle: flow in rivers, overland flow, groundwater flow in the saturated and unsaturated zones, and evapotranspiration and precipitation time series. Stochastic parameters and hydrologic conditions over a five year period of historical hydrological data were used to analyze the hydrological cycle and to determine the prevailing mercury transport mechanism within the watershed. Simulations of remediation scenarios revealed that reduction of the highly contaminated point sources, rather than general remediation of the contaminant plume, has a more direct impact on downstream mercury concentrations.

Do Changes in Current Flow as a Result of Arrays of Tidal Turbines Have an Effect on Benthic Communities?

Arrays of tidal energy converters have the potential to provide clean renewable energy for future generations. Benthic communities may, however, be affected by changes in current speeds resulting from arrays of tidal converters located in areas characterised by strong currents. Current speed, together with bottom type and depth, strongly influence benthic community distributions; however the interaction of these factors in controlling benthic dynamics in high energy environments is poorly understood. The Strangford Lough Narrows, the location of SeaGen, the world’s first single full-scale, grid-compliant tidal energy extractor, is characterised by spatially heterogenous high current flows. A hydrodynamic model was used to select a range of benthic community study sites that had median flow velocities between 1.5–2.4 m/s in a depth range of 25–30 m. 25 sites were sampled for macrobenthic community structure using drop down video survey to test the sensitivity of the distribution of benthic communities to changes in the flow field. A diverse range of species were recorded which were consistent with those for high current flow environments and corresponding to very tide-swept faunal communities in the EUNIS classification. However, over the velocity range investigated, no changes in benthic communities were observed. This suggested that the high physical disturbance associated with the high current flows in the Strangford Narrows reflected the opportunistic nature of the benthic species present with individuals being continuously and randomly affected by turbulent forces and physical damage. It is concluded that during operation, the removal of energy by marine tidal energy arrays in the far-field is unlikely to have a significant effect on benthic communities in high flow environments. The results are of major significance to developers and regulators in the tidal energy industry when considering the environmental impacts for site licences.

Design Optimization of Nozzle Shapes for Maximum Uniformity of Exit Flow

The objective of this study is to identify the optimal designs of converging-diverging supersonic and hypersonic nozzles that perform at maximum uniformity of thermodynamic and flow-field properties with respect to their average values at the nozzle exit. Since this is a multi-objective design optimization problem, the design variables used are parameters defining the shape of the nozzle. This work presents how variation of such parameters can influence the nozzle exit flow non-uniformities. A Computational Fluid Dynamics (CFD) software package, ANSYS FLUENT, was used to simulate the compressible, viscous gas flow-field in forty nozzle shapes, including the heat transfer analysis. The results of two turbulence models, k-e and k-ω, were computed and compared. With the analysis results obtained, the Response Surface Methodology (RSM) was applied for the purpose of performing a multi-objective optimization. The optimization was performed with ModeFrontier software package using Kriging and Radial Basis Functions (RBF) response surfaces. Final Pareto optimal nozzle shapes were then analyzed with ANSYS FLUENT to confirm the accuracy of the optimization process.

Movement effects on the flow physics and nutrient delivery in engineered valvular tissues

Mechanical conditioning has been shown to promote tissue formation in a wide variety of tissue engineering efforts. However the underlying mechanisms by which external mechanical stimuli regulate cells and tissues are not known. This is particularly relevant in the area of heart valve tissue engineering (HVTE) owing to the intense hemodynamic environments that surround native valves. Some studies suggest that oscillatory shear stress (OSS) caused by steady flow and scaffold flexure play a critical role in engineered tissue formation derived from bone marrow derived stem cells (BMSCs). In addition, scaffold flexure may enhance nutrient (e.g. oxygen, glucose) transport. In this study, we computationally quantified the i) magnitude of fluid-induced shear stresses; ii) the extent of temporal fluid oscillations in the flow field using the oscillatory shear index (OSI) parameter, and iii) glucose and oxygen mass transport profiles. Noting that sample cyclic flexure induces a high degree of oscillatory shear stress (OSS), we incorporated moving boundary computational fluid dynamic simulations of samples housed within a bioreactor to consider the effects of: 1) no flow, no flexure (control group), 2) steady flow-alone, 3) cyclic flexure-alone and 4) combined steady flow and cyclic flexure environments. We also coupled a diffusion and convention mass transport equation to the simulated system. We found that the coexistence of both OSS and appreciable shear stress magnitudes, described by the newly introduced parameter OSI-:τ: explained the high levels of engineered collagen previously observed from combining cyclic flexure and steady flow states. On the other hand, each of these metrics on its own showed no association. This finding suggests that cyclic flexure and steady flow synergistically promote engineered heart valve tissue production via OSS, so long as the oscillations are accompanied by a critical magnitude of shear stress. In addition, our simulations showed that mass transport of glucose and oxygen is enhanced by sample movement at low sample porosities, but did not play a role in highly porous scaffolds. Preliminary in-house in vitro experiments showed that cell proliferation and phenotype is enhanced in OSI-:τ: environments.^

Influence of land-use change on near-surface hydrological processes: Undisturbed forest to pasture

Soil compaction that follows the clearing of tropical forest for cattle pasture is associated with lower soil hydraulic conductivity and increased frequency and volume of overland flow. We investigated the frequency of perched water tables, overland flow and stormflow in an Amazon forest and in an adjacent 25-year-old pasture cleared from the same forest. We compared the results with the frequencies of these phenomena estimated from comparisons of rainfall intensity and soil hydraulic conductivity. The frequency of perched water tables based on rainfall intensity and soil hydraulic conductivity was expected to double in pasture compared with forest. This corresponded closely with an approximate doubling of the frequency of stormflow and overland flow in pasture. In contrast, the stormflow volume in pasture increased 17-fold. This disproportional increase of stormflow resulted from overland flow generation over large areas of pasture, while overland flow generation in the forest was spatially limited and was observed only very near the stream channel. In both catchments, stormflow was generated by saturation excess because of perched water tables and near-surface groundwater levels. Stormflow was occasionally generated in the forest by rapid return flow from macropores, while slow return flow from a continuous perched water table was more common in the pasture. These results suggest that deforestation for pasture alters fundamental mechanisms of stormflow generation and may increase runoff volumes over wide regions of Amazonia. (C) 2009 Elsevier B.V. All rights reserved.