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.

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.

Reconstruction of an unsteady flow from incomplete boundary data:Part I. Theory

Purpose – To propose and investigate a stable numerical procedure for the reconstruction of the velocity of a viscous incompressible fluid flow in linear hydrodynamics from knowledge of the velocity and fluid stress force given on a part of the boundary of a bounded domain. Design/methodology/approach – Earlier works have involved the similar problem but for stationary case (time-independent fluid flow). Extending these ideas a procedure is proposed and investigated also for the time-dependent case. Findings – The paper finds a novel variation method for the Cauchy problem. It proves convergence and also proposes a new boundary element method. Research limitations/implications – The fluid flow domain is limited to annular domains; this restriction can be removed undertaking analyses in appropriate weighted spaces to incorporate singularities that can occur on general bounded domains. Future work involves numerical investigations and also to consider Oseen type flow. A challenging problem is to consider non-linear Navier-Stokes equation. Practical implications – Fluid flow problems where data are known only on a part of the boundary occur in a range of engineering situations such as colloidal suspension and swimming of microorganisms. For example, the solution domain can be the region between to spheres where only the outer sphere is accessible for measurements. Originality/value – A novel variational method for the Cauchy problem is proposed which preserves the unsteady Stokes operator, convergence is proved and using recent for the fundamental solution for unsteady Stokes system, a new boundary element method for this system is also proposed.

Modelling and simulation of biomass gasification in a circulating fluidized bed reactor

Computational Fluid Dynamics (CFD) has found great acceptance among the engineering community as a tool for research and design of processes that are practically difficult or expensive to study experimentally. One of these processes is the biomass gasification in a Circulating Fluidized Bed (CFB). Biomass gasification is the thermo-chemical conversion of biomass at a high temperature and a controlled oxygen amount into fuel gas, also sometime referred to as syngas. Circulating fluidized bed is a type of reactor in which it is possible to maintain a stable and continuous circulation of solids in a gas-solid system. The main objectives of this thesis are four folds: (i) Develop a three-dimensional predictive model of biomass gasification in a CFB riser using advanced Computational Fluid Dynamic (CFD) (ii) Experimentally validate the developed hydrodynamic model using conventional and advanced measuring techniques (iii) Study the complex hydrodynamics, heat transfer and reaction kinetics through modelling and simulation (iv) Study the CFB gasifier performance through parametric analysis and identify the optimum operating condition to maximize the product gas quality. Two different and complimentary experimental techniques were used to validate the hydrodynamic model, namely pressure measurement and particle tracking. The pressure measurement is a very common and widely used technique in fluidized bed studies, while, particle tracking using PEPT, which was originally developed for medical imaging, is a relatively new technique in the engineering field. It is relatively expensive and only available at few research centres around the world. This study started with a simple poly-dispersed single solid phase then moved to binary solid phases. The single solid phase was used for primary validations and eliminating unnecessary options and steps in building the hydrodynamic model. Then the outcomes from the primary validations were applied to the secondary validations of the binary mixture to avoid time consuming computations. Studies on binary solid mixture hydrodynamics is rarely reported in the literature. In this study the binary solid mixture was modelled and validated using experimental data from the both techniques mentioned above. Good agreement was achieved with the both techniques. According to the general gasification steps the developed model has been separated into three main gasification stages; drying, devolatilization and tar cracking, and partial combustion and gasification. The drying was modelled as a mass transfer from the solid phase to the gas phase. The devolatilization and tar cracking model consist of two steps; the devolatilization of the biomass which is used as a single reaction to generate the biomass gases from the volatile materials and tar cracking. The latter is also modelled as one reaction to generate gases with fixed mass fractions. The first reaction was classified as a heterogeneous reaction while the second reaction was classified as homogenous reaction. The partial combustion and gasification model consisted of carbon combustion reactions and carbon and gas phase reactions. The partial combustion considered was for C, CO, H2 and CH4. The carbon gasification reactions used in this study is the Boudouard reaction with CO2, the reaction with H2O and Methanation (Methane forming reaction) reaction to generate methane. The other gas phase reactions considered in this study are the water gas shift reaction, which is modelled as a reversible reaction and the methane steam reforming reaction. The developed gasification model was validated using different experimental data from the literature and for a wide range of operating conditions. Good agreement was observed, thus confirming the capability of the model in predicting biomass gasification in a CFB to a great accuracy. The developed model has been successfully used to carry out sensitivity and parametric analysis. The sensitivity analysis included: study of the effect of inclusion of various combustion reaction; and the effect of radiation in the gasification reaction. The developed model was also used to carry out parametric analysis by changing the following gasifier operating conditions: fuel/air ratio; biomass flow rates; sand (heat carrier) temperatures; sand flow rates; sand and biomass particle sizes; gasifying agent (pure air or pure steam); pyrolysis models used; steam/biomass ratio. Finally, based on these parametric and sensitivity analysis a final model was recommended for the simulation of biomass gasification in a CFB riser.

Stochastically driven instability in rotating shear flows

Origin of hydrodynamic turbulence in rotating shear flows is investigated. The particular emphasis is on flows whose angular velocities decrease but specific angular momenta increase with increasing radial coordinate. Such flows are Rayleigh stable, but must be turbulent in order to explain observed data. Such a mismatch between the linear theory and observations/experiments is more severe when any hydromagnetic/magnetohydrodynamic instability and the corresponding turbulence therein is ruled out. The present work explores the effect of stochastic noise on such hydrodynamic flows. We focus on a small section of such a flow which is essentially a plane shear flow supplemented by the Coriolis effect. This also mimics a small section of an astrophysical accretion disk. It is found that such stochastically driven flows exhibit large temporal and spatial correlations of perturbation velocities, and hence large energy dissipations, that presumably generate instability. A range of angular velocity profiles (for the steady flow), starting with the constant angular momentum to that of the constant circular velocity are explored. It is shown that the growth and roughness exponents calculated from the contour (envelope) of the perturbed flows are all identical, revealing a unique universality class for the stochastically forced hydrodynamics of rotating shear flows. This work, to the best of our knowledge, is the first attempt to understand origin of instability and turbulence in the three-dimensional Rayleigh stable rotating shear flows by introducing additive stochastic noise to the underlying linearized governing equations. This has important implications in resolving the turbulence problem in astrophysical hydrodynamic flows such as accretion disks.

Modelling a biochemical reaction with computational fluid dynamics

Earlier investigations (Cartland Glover et al., 2004) into the use of computational fluid dynamics (CFD) for the modelling of gas-liquid and gas-liquid-solid flow allowed a simple biochemical reaction model to be implemented. A single plane mesh was used to represent the transport and reaction of molasses, the mould Aspergillus niger and citric acid in a bubble column with a height to diameter aspect ratio of 20:1. Two specific growth rates were used to examine the impact that biomass growth had on the local solids concentration and the effect this had on the local hydrodynamics of the bubble column.

Concurrent multiscale modelling of atomistic and hydrodynamic processes in liquids

Fluctuations of liquids at the scales where the hydrodynamic and atomistic descriptions overlap are considered. The importance of these fluctuations for atomistic motions is discussed and examples of their accurate modelling with a multi-space-time-scale fluctuating hydrodynamics scheme are provided. To resolve microscopic details of liquid systems, including biomolecular solutions, together with macroscopic fluctuations in space-time, a novel hybrid atomistic-fluctuating hydrodynamics approach is introduced. For a smooth transition between the atomistic and continuum representations, an analogy with two-phase hydrodynamics is used that leads to a strict preservation of macroscopic mass and momentum conservation laws. Examples of numerical implementation of the new hybrid approach for the multiscale simulation of liquid argon in equilibrium conditions are provided. © 2014 The Author(s) Published by the Royal Society.

A CFD study of biomass pyrolysis in a downer reactor equipped with a novel gas-solid separator - I:hydrodynamic performance

This study presents the first part of a CFD study on the performance of a downer reactor for biomass pyrolysis. The reactor was equipped with a novel gas-solid separation method, developed by the co-authors from the ICFAR (Canada). The separator, which was designed to allow for fast separation of clean pyrolysis gas, consisted of a cone deflector and a gas exit pipe installed inside the downer reactor. A multi-fluid model (Eulerian-Eulerian) with constitutive relations adopted from the kinetic theory of granular flow was used to simulate the multiphase flow. The effects of the various parameters including operation conditions, separator geometry and particle properties on the overall hydrodynamics and separation efficiency were investigated. The model prediction of the separator efficiency was compared with experimental measurements. The results revealed distinct hydrodynamic features around the cone separator, allowing for up to 100% separation efficiency. The developed model provided a platform for the second part of the study, where the biomass pyrolysis is simulated and the product quality as a function of operating conditions is analyzed. Crown Copyright © 2014 Published by Elsevier B.V. All rights reserved.

Predicting the effect of mixing on oxygen transfer and nutrient removal in activated sludge basins

Activated sludge basins (ASBs) are a key-step in wastewater treatment processes that are used to eliminate biodegradable pollution from the water discharged to the natural environment. Bacteria found in the activated sludge consume and assimilate nutrients such as carbon, nitrogen and phosphorous under specific environmental conditions. However, applying the appropriate agitation and aeration regimes to supply the environmental conditions to promote the growth of the bacteria is not easy. The agitation and aeration regimes that are applied to activated sludge basins have a strong influence on the efficacy of wastewater treatment processes. The major aims of agitation by submersible mixers are to improve the contact between biomass and wastewater and the prevention of biomass settling. They induce a horizontal flow in the oxidation ditch, which can be quantified by the mean horizontal velocity. Mean values of 0.3-0.35 m s-1 are recommended as a design criteria to ensure best conditions for mixing and aeration (Da Silva, 1994). To give circulation velocities of this order of magnitude, the positioning and types of mixers are chosen from the plant constructors' experience and the suppliers' data for the impellers. Some case studies of existing plants have shown that measured velocities were not in the range that was specified in the plant design. This illustrates that there is still a need for design and diagnosis approach to improve process reliability by eliminating or reducing the number of short circuits, dead zones, zones of inefficient mixing and poor aeration. The objective of the aeration is to facilitate the quick degradation of pollutants by bacterial growth. To achieve these objectives a wastewater treatment plant must be adequately aerated; thus resulting in 60-80% of all energetic consummation being dedicated to the aeration alone (Juspin and Vasel, 2000). An earlier study (Gillot et al., 1997) has illustrated the influence that hydrodynamics have on the aeration performance as measure by the oxygen transfer coefficient. Therefore, optimising the agitation and aeration systems can enhance the oxygen transfer coefficient and consequently reduce the operating costs of the wastewater treatment plant. It is critically important to correctly estimate the mass transfer coefficient as any errors could result in the simulations of biological activity not being physically representative. Therefore, the transfer process was rigorously examined in several different types of process equipment to determine the impact that different hydrodynamic regimes and liquid-side film transfer coefficients have on the gas phase and the mass transfer of oxygen. To model the biological activity occurring in ASBs, several generic biochemical reaction models have been developed to characterise different biochemical reaction processes that are known as Activated Sludge Models, ASM (Henze et al., 2000). The ASM1 protocol was selected to characterise the impact of aeration on the bacteria consuming and assimilating ammonia and nitrate in the wastewater. However, one drawback of ASM protocols is that the hydrodynamics are assumed to be uniform by the use of perfectly mixed, plug flow reactors or as a number of perfectly mixed reactors in series. This makes it very difficult to identify the influence of mixing and aeration on oxygen mass transfer and biological activity. Therefore, to account for the impact of local gas-liquid mixing regime on the biochemical activity Computational Fluid Dynamics (CFD) was used by applying the individual ASM1 reaction equations as the source terms to a number of scalar equations. Thus, the application of ASM1 to CFD (FLUENT) enabled the investigation of the oxygen transfer efficiency and the carbon & nitrogen biological removal in pilot (7.5 cubic metres) and plant scale (6000 cubic metres) ASBs. Both studies have been used to validate the effect that the hydrodynamic regime has on oxygen mass transfer (the circulation velocity and mass transfer coefficient) and the effect that this had on the biological activity on pollutants such as ammonia and nitrate (Cartland Glover et al., 2005). The work presented here is one part to of an overall approach for improving the understanding of ASBs and the impact that they have in terms of the hydraulic and biological performance on the overall wastewater treatment process. References CARTLAND GLOVER G., PRINTEMPS C., ESSEMIANI K., MEINHOLD J., (2005) Modelling of wastewater treatment plants ? How far shall we go with sophisticated modelling tools? 3rd IWA Leading-Edge Conference & Exhibition on Water and Wastewater Treatment Technologies, 6-8 June 2005, Sapporo, Japan DA SILVA G. (1994). Eléments d'optimisation du transfert d'oxygène par fines bulles et agitateur séparé en chenal d'oxydation. PhD Thesis. CEMAGREF Antony ? France. GILLOT S., DERONZIER G., HEDUIT A. (1997). Oxygen transfer under process conditions in an oxidation ditch equipped with fine bubble diffusers and slow speed mixers. WEFTEC, Chicago, USA. HENZE M., GUJER W., MINO T., van LOOSDRECHT M., (2000). Activated Sludge Models ASM1, ASM2, ASM2D and ASM3, Scientific and Technical Report No. 9. IWA Publishing, London, UK. JUSPIN H., VASEL J.-L. (2000). Influence of hydrodynamics on oxygen transfer in the activated sludge process. IWA, Paris - France.

Hydrogen production by steam reforming of DME in a large scale CFB reactor. Part I:computational model and predictions

This study presents a computational fluid dynamic (CFD) study of Dimethyl Ether steam reforming (DME-SR) in a large scale Circulating Fluidized Bed (CFB) reactor. The CFD model is based on Eulerian-Eulerian dispersed flow and solved using commercial software (ANSYS FLUENT). The DME-SR reactions scheme and kinetics in the presence of a bifunctional catalyst of CuO/ZnO/Al2O3+ZSM-5 were incorporated in the model using in-house developed user-defined function. The model was validated by comparing the predictions with experimental data from the literature. The results revealed for the first time detailed CFB reactor hydrodynamics, gas residence time, temperature distribution and product gas composition at a selected operating condition of 300 °C and steam to DME mass ratio of 3 (molar ratio of 7.62). The spatial variation in the gas species concentrations suggests the existence of three distinct reaction zones but limited temperature variations. The DME conversion and hydrogen yield were found to be 87% and 59% respectively, resulting in a product gas consisting of 72 mol% hydrogen. In part II of this study, the model presented here will be used to optimize the reactor design and study the effect of operating conditions on the reactor performance and products.

Multiscale modelling:approaches and challenges

Multiscale systems that are characterized by a great range of spatial–temporal scales arise widely in many scientific domains. These range from the study of protein conformational dynamics to multiphase processes in, for example, granular media or haemodynamics, and from nuclear reactor physics to astrophysics. Despite the diversity in subject areas and terminology, there are many common challenges in multiscale modelling, including validation and design of tools for programming and executing multiscale simulations. This Theme Issue seeks to establish common frameworks for theoretical modelling, computing and validation, and to help practical applications to benefit from the modelling results. This Theme Issue has been inspired by discussions held during two recent workshops in 2013: ‘Multiscale modelling and simulation’ at the Lorentz Center, Leiden (http://www.lorentzcenter.nl/lc/web/2013/569/info.php3?wsid=569&venue=Snellius), and ‘Multiscale systems: linking quantum chemistry, molecular dynamics and microfluidic hydrodynamics’ at the Royal Society Kavli Centre. The objective of both meetings was to identify common approaches for dealing with multiscale problems across different applications in fluid and soft matter systems. This was achieved by bringing together experts from several diverse communities.

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.

Two-phase flow analogy as an effective boundary condition for modelling liquids at atomistic resolution

A hybrid Molecular Dynamics/Fluctuating Hydrodynamics framework based on the analogy with two-phase hydrodynamics has been extended to dynamically tracking the feature of interest at all-atom resolution. In the model, the hydrodynamics description is used as an effective boundary condition to close the molecular dynamics solution without resorting to standard periodic boundary conditions. The approach is implemented in a popular Molecular Dynamics package GROMACS and results for two biomolecular systems are reported. A small peptide dialanine and a complete capsid of a virus porcine circovirus 2 in water are considered and shown to reproduce the structural and dynamic properties compared to those obtained in theory, purely atomistic simulations, and experiment.

Geometrical optimization of a fast pyrolysis bubbling fluidized bed reactor using comutational fluid dynamics

This paper analyzes the physical phenomena that take place inside an 1 kg/h bubbling fluidized bed reactor located at Aston University and presents a geometrically modified version of it, in order to improve certain hydrodynamic and gas flow characteristics. The bed uses, in its current operation, 40 L/min of N2 at 520 °C fed through a distributor plate and 15 L/min purge gas stream, i.e., N2 at 20 °C, via the feeding tube. The Eulerian model of FLUENT 6.3 is used for the simulation of the bed hydrodynamics, while the k - ε model accounts for the effect of the turbulence field of one phase on the other. The three-dimensional simulation of the current operation of the reactor showed that a stationary bubble was formed next to the feeding tube. The size of the permanent bubble reaches up to the splash zone of the reactor, without any fluidizaton taking place underneath the feeder. The gas flow dynamics in the freeboard of the reactor is also analyzed. A modified version of the reactor is presented, simulated, and analyzed, together with a discussion on the impact of the flow dynamics on the fast pyrolysis of biomass. © 2010 American Chemical Society.