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.

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.

A comparison of in situ constant and falling head permeameter tests to assess the distribution of clogging within horizontal subsurface flow constructed wetlands

Clogging is the main operational problem associated with horizontal subsurface flow constructed wetlands (HSSF CWs). The measurement of saturated hydraulic conductivity has proven to be a suitable technique to assess clogging within HSSF CWs. The vertical and horizontal distribution of hydraulic conductivity was assessed in two full-scale HSSF CWs by using two different in situ permeameter methods (falling head (FH) and constant head (CH) methods). Horizontal hydraulic conductivity profiles showed that both methods are correlated by a power function (FH= CH 0.7821, r 2=0.76) within the recorded range of hydraulic conductivities (0-70 m/day). However, the FH method provided lower values of hydraulic conductivity than the CH method (one to three times lower). Despite discrepancies between the magnitudes of reported readings, the relative distribution of clogging obtained via both methods was similar. Therefore, both methods are useful when exploring the general distribution of clogging and, specially, the assessment of clogged areas originated from preferential flow paths within full-scale HSSF CWs. Discrepancy between methods (either in magnitude and pattern) aroused from the vertical hydraulic conductivity profiles under highly clogged conditions. It is believed this can be attributed to procedural differences between the methods, such as the method of permeameter insertion (twisting versus hammering). Results from both methods suggest that clogging develops along the shortest distance between water input and output. Results also evidence that the design and maintenance of inlet distributors and outlet collectors appear to have a great influence on the pattern of clogging, and hence the asset lifetime of HSSF CWs. © Springer Science+Business Media B.V. 2011.

Clogging in subsurface-flow treatment wetlands:measurement, modeling and management

This paper reviews the state of the art in measuring, modeling, and managing clogging in subsurface-flow treatment wetlands. Methods for measuring in situ hydraulic conductivity in treatment wetlands are now available, which provide valuable insight into assessing and evaluating the extent of clogging. These results, paired with the information from more traditional approaches (e.g., tracer testing and composition of the clog matter) are being incorporated into the latest treatment wetland models. Recent finite element analysis models can now simulate clogging development in subsurface-flow treatment wetlands with reasonable accuracy. Various management strategies have been developed to extend the life of clogged treatment wetlands, including gravel excavation and/or washing, chemical treatment, and application of earthworms. These strategies are compared and available cost information is reported. © 2012 Elsevier Ltd.

Clogging in subsurface-flow treatment wetlands:occurrence and contributing factors

Clogging is a major operational and maintenance issue associated with the use of subsurface flow wetlands for wastewater treatment, and can ultimately limit the lifetime of the system. This review considers over two decades of accumulated knowledge regarding clogging in both vertical and horizontal subsurface flow treatment wetlands. The various physical, chemical and biological factors responsible for clogging are identified and discussed. The occurrence of clogging is placed into the context of various design and operational parameters such as wastewater characteristics, upstream treatment processes, intermittent or continuous operation, influent distribution, and media type. This information is then used to describe how clogging develops within, and subsequently impacts, common variants of subsurface flow treatment wetland typically used in the U.S., U.K., France and Germany. Comparison of these systems emphasized that both hydraulic loading rate and solids loading rate need to be considered when designing systems to operate robustly, i.e. hydraulic overloading makes horizontal-flow tertiary treatment systems in the U.K. more susceptible to clogging problems than vertical-flow primary treatment systems in France. Future research should focus on elucidating the underlying mechanisms of clogging as they relate to the design, operation, and maintenance of subsurface flow treatment wetlands. © 2010 Elsevier B.V.

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.

A three dimensional numerical investigation of hillslope flow processes

Physically based distributed models of catchment hydrology are likely to be made available as engineering tools in the near future. Although these models are based on theoretically acceptable equations of continuity, there are still limitations in the present modelling strategy. Of interest to this thesis are the current modelling assumptions made concerning the effects of soil spatial variability, including formations producing distinct zones of preferential flow. The thesis contains a review of current physically based modelling strategies and a field based assessment of soil spatial variability. In order to investigate the effects of soil nonuniformity a fully three dimensional model of variability saturated flow in porous media is developed. The model is based on a Galerkin finite element approximation to Richards equation. Accessibility to a vector processor permits numerical solutions on grids containing several thousand node points. The model is applied to a single hillslope segment under various degrees of soil spatial variability. Such variability is introduced by generating random fields of saturated hydraulic conductivity using the turning bands method. Similar experiments are performed under conditions of preferred soil moisture movement. The results show that the influence of soil variability on subsurface flow may be less significant than suggested in the literature, due to the integrating effects of three dimensional flow. Under conditions of widespread infiltration excess runoff, the results indicate a greater significance of soil nonuniformity. The recognition of zones of preferential flow is also shown to be an important factor in accurate rainfall-runoff modelling. Using the results of various fields of soil variability, experiments are carried out to assess the validity of the commonly used concept of `effective parameters'. The results of these experiments suggest that such a concept may be valid in modelling subsurface flow. However, the effective parameter is observed to be event dependent when the dominating mechanism is infiltration excess runoff.

The interaction of luminance and texture amplitude in surface depth perception

Previous studies have suggested separate channels for detection of first-order luminance modulations (LM) and second-order modulations of the local amplitude (AM) of a texture. Mixtures of LM and AM with different phase relationships appear very different: in-phase compounds (LM + AM) look like 3-D corrugated surfaces, while out-of-phase compounds (LM - AM) appear flat and/or transparent. This difference may arise because the in-phase compounds are consistent with multiplicative shading, while the out-of-phase compounds are not. We investigated the role of these modulation components in surface depth perception. We used a textured background with thin bars formed by local changes in luminance and/or texture amplitude. These stimuli appear as embossed surfaces with wide and narrow regions. Keeping the AM modulation depth fixed at a suprathreshold level, we determined the amount of luminance contrast required for observers to correctly indicate the width (narrow or wide) of 'raised' regions in the display. Performance (compared to the LM-only case) was facilitated by the presence of AM, but, unexpectedly, performance for LM - AM was as good as for LM + AM. Thus, these results suggest that there is an interaction between first-order and second-order mechanisms during depth perception based on shading cues, but the phase dependence is not yet understood.

Cell surface analysis of the basidiomycete yeast cryptococcus neoformans

Cell surface properties of the basidiomycete yeast Cryptococcus neoformans were investigated with a combination of novel and well proven approaches. Non-specific cell adhesion forces, as well as exposed carbohydrate and protein moieties potentially associated with specific cellular interaction, were analysed. Experimentation and analysis employed cryptococcal cells of different strains, capsular status and culture age. Investigation of cellular charge by particulate microelectrophoresis revealed encapsulated yeast forms of C. neoformans manifest a distinctive negative charge regardless of the age of cells involved; in turn, the neutral charge of acapsulate yeasts confirmed that the polysaccharide capsule, and not the cell wall, was responsible for this occurrence. Hydrophobicity was measured by MATH and HICH techniques, as well as by the attachment of polystyrene microspheres. All three techniques, where applicable, found C. neoformans yeast to be consistently hydrophilic; this state varied little regardless of strain and culture age. Cell surface carbohydrates and protein were investigated with novel fluorescent tagging protocols, flow cytometry and confocal microscopy. Cell surface carbohydrate was identified by controlled oxidation in association with biotin hydrazide and fluorescein-streptavidin tagging. Marked amounts of carbohydrate were measured and observed on the cell wall surface of cryptococcal yeasts. Furthermore, tagging of carbohydrates with selective fluorescent lectins supported the identification, measurement and observation of substantial amounts of mannose, glucose and N-acetyl-glucosamine. Cryptococcal cell surface protein was identified using sulfo-NHS-biotin with fluorescein-streptavidin, and then readily quantified by flow cytometry. Confocal imaging of surface exposed carbohydrate and protein revealed common localised areas of vivid fluorescence associated with buds, bud scars and nascent daughter cells. Carbohydrate and protein fluorescence often varied between strains, culture age and capsule status of cells examined. Finally, extension of protein tagging techniques resulted in the isolation and extraction of two biotinylated proteins from the yeast cell wall surface of an acapsulate strain of C.neoformans.

Effects of electric charge on osmotic flow across periodically arranged circular cylinders

An electrostatic model is developed for osmotic flow across a layer consisting of identical circular cylinders with a fixed surface charge, aligned parallel to each other so as to form an ordered hexagonal arrangement. The expression of the osmotic reflection coefficient is derived for spherical solutes with a fixed surface charge suspended in an electrolyte, based on low-Reynolds-number hydrodynamics and a continuum, point-charge description of the electric double layers. The repulsive electrostatic interaction between the surface charges with the same sign on the solute and the cylinders is shown to increase the exclusion region of solute from the cylinder surface, which enhances the osmotic flow. Applying the present model to the study of osmotic flow across the endothelial surface glycocalyx of capillary walls has revealed that this electrostatic model could account well for the reflection coefficients measured for charged macromolecules, such as albumin, in the physiological range of charge density and ion concentration.

Electrical charge effect on osmotic flow through pores

An electrostatic model for osmotic flow through circular cylindrical pores is developed to describe the reflection coefficient for the membrane transport in the presence of surface charges on the pore wall and the solute. For a spherical solute placed at an arbitrary radial position in the pore, the electrical potential was computed by a spectral element method applied to the Poisson-Boltzmann equation together with the condition of electrical neutrality. The interaction energy between the surface charges was used to estimate the osmotic reflection coefficient. The proposed model predicts that even for a small Debye length compared to the pore radius, the repulsive electrostatic interaction between the surface charges could significantly increase the osmotic flow through the pore.

Charge effects on the hindered transport of macromolecules across the endothelial surface glycocalyx layer

A fluid mechanical and electrostatic model for the transport of solute molecules across the vascular endothelial surface glycocalyx layer (EGL) was developed to study the charge effect on the diffusive and convective transport of the solutes. The solute was assumed to be a spherical particle with a constant surface charge density, and the EGL was represented as an array of periodically arranged circular cylinders of like charge, with a constant surface charge density. By combining the fluid mechanical analyses for the flow around a solute suspended in an electrolyte solution and the electrostatic analyses for the free energy of the interaction between the solute and cylinders based on a mean field theory, we estimated the transport coefficients of the solute across the EGL. Both of diffusive and convective transports are reduced compared to those for an uncharged system, due to the stronger exclusion of the solute that results from the repulsive electrostatic interaction. The model prediction for the reflection coefficient for serum albumin agreed well with experimental observations if the charge density in the EGL is ranged from approximately -10 to -30 mEq/l.

Flow across microvessel walls through the endothelial surface glycocalyx and the interendothelial cleft

A mathematical model is presented for steady fluid flow across microvessel walls through a serial pathway consisting of the endothelial surface glycocalyx and the intercellular cleft between adjacent endothelial cells, with junction strands and their discontinuous gaps. The three-dimensional flow through the pathway from the vessel lumen to the tissue space has been computed numerically based on a Brinkman equation with appropriate values of the Darcy permeability. The predicted values of the hydraulic conductivity Lp, defined as the ratio of the flow rate per unit surface area of the vessel wall to the pressure drop across it, are close to experimental measurements for rat mesentery microvessels. If the values of the Darcy permeability for the surface glycocalyx are determined based on the regular arrangements of fibres with 6nm radius and 8nm spacing proposed recently from the detailed structural measurements, then the present study suggests that the surface glycocalyx could be much less resistant to flow compared to previous estimates by the one-dimensional flow analyses, and the intercellular cleft could be a major determinant of the hydraulic conductivity of the microvessel wall.

Gravity spun polycaprolactone fibres for soft tissue engineering:interaction with fibroblasts and myoblasts in cell culture

Poly(ε-caprolactone) (PCL) fibres were produced by wet spinning from solutions in acetone under low shear (gravity flow) conditions. As-spun PCL fibres exhibited a mean strength and stiffness of 7.9 MPa and 0.1 GPa, respectively and a rough, porous surface morphology. Cold drawing to an extension of 500% resulted in increases in fibre strength (43 MPa) and stiffness (0.3 GPa) and development of an oriented, fibrillar surface texture. The proliferation rate of Swiss 3T3 mouse fibroblasts and C2C12 mouse myoblasts on as-spun, 500% cold-drawn and gelatin-modified PCL fibres was determined in cell culture to provide a basic measure of the biocompatibility of the fibres. Proliferation of both cell types was consistently higher on gelatin-coated fibres relative to as-spun fibres at time points below 7 days. Fibroblast growth rates on cold-drawn PCL fibres exceeded those on as-spun fibres but myoblast proliferation was similar on both substrates. After 1 day in culture, both cell types had spread and coalesced on the fibres to form a cell layer, which conformed closely to the underlying topography. The high fibre compliance combined with a potential for modifying the fibre surface chemistry with cell adhesion molecules and the surface architecture by cold drawing to enhance proliferation of fibroblasts and myoblasts, recommends further investigation of gravity-spun PCL fibres for 3-D scaffold production in soft tissue engineering. © 2005 Elsevier Ltd. All rights reserved.

Surface and subsurface Vickers indentation cracks in SiC, Si3N4, and sialon ceramics

The nature of subsurface cracks formed under and around Vickers hardness indentations is often assumed rather than identified. Subsurface cracks in four engineering ceramics are revealed using a penetrant technique, and flaw dimensions are recorded. The resulting data are used to investigate several aspects of indentation cracking, such as crack shape, functional relationships between indentation load and flaw dimensions, and the performance of indentation fracture toughness equations. An R curve is constructed for each of the materials. © 1995 The Institute of Materials.