Periodic venting of MABR lumen allows high removal rates and high gas-transfer efficiencies.

The membrane-aerated biofilm reactor (MABR) is a novel treatment technology that employs gas-supplying membranes to deliver oxygen directly to a biofilm growing on the membrane surface. When operated with closed-end membranes, the MABR provides 100-percent oxygen transfer efficiencies (OTE), resulting in significant energy savings. However, closed-end MABRs are more sensitive to back-diffusion of inert gases, such as nitrogen. Back-diffusion reduces the average oxygen transfer rates (OTR), consequently decreasing the average contaminant removal fluxes (J). We hypothesized that venting the membrane lumen periodically would increase the OTR and J. Using an experimental flow cell and mathematical modeling, we showed that back-diffusion gas profiles developed over relatively long timescales. Thus, very short ventings could re-establish uniform gas profiles for relatively long time periods. Using modeling, we systematically explored the effect of the venting interval (time between ventings). At moderate venting intervals, opening the membrane for 20 s every 30 min, the venting significantly increased the average OTR and J without substantially impacting the OTEs. When the interval was short enough, in this case shorter than 20 min, the OTR was actually higher than for continuous open-end operation. Our results show that periodic venting is a promising strategy to combine the advantages of open-end and closed end operation, maximizing both the OTR and OTE.

Combined biofouling and scaling in membrane feed channels: a new modeling approach.

A mathematical model was developed for combined fouling due to biofilms and mineral precipitates in membrane feed channels with spacers. Finite element simulation of flow and solute transport in two-dimensional geometries was coupled with a particle-based approach for the development of a composite (cells and crystals) foulant layer. Three fouling scenarios were compared: biofouling only, scaling only and combined fouling. Combined fouling causes a quicker flux decline than the summed flux deterioration when scaling and biofouling act independently. The model results indicate that the presence of biofilms leads to more mineral formation due to: (1) an enhanced degree of saturation for salts next to the membrane and within the biofilm; and (2) more available surface for nucleation to occur. The impact of biofilm in accelerating gypsum precipitation depends on the composition of the feed water (eg the presence of NaCl) and the kinetics of crystal nucleation and growth. Interactions between flow, solute transport and biofilm-induced mineralization are discussed.

Modeling the nutrient removal process in aerobic granular sludge system by coupling the reactor- and granule-scale models.

We developed a model for nutrient removal in an aerobic granular sludge system. This model can quantitatively describe the start-up of the system by coupling a model for studying the population dynamics of the granules in the reactor (reactor-scale model) and a model for studying the microbial community structure in the granules (granule-scale model). The reactor-scale model is used for simulation for 10 days from the start, during which the granule size is relatively small; the granule-scale model is used after Day 10. The present approach proposes the output data of the reactor-scale model after 10 days as initial conditions for the granule-scale model. The constructed model satisfactorily describes experimental data in various spatial and temporal scales, which were obtained in this study by performing the anaerobic-aerobic-anoxic cycles using a sequencing batch reactor. Simulations using this model quantitatively predicted that the stability of nutrient removal process depended largely on the dissolved oxygen (DO) concentration, and the DO setpoint adaptation could improve the nutrient removal performance.

Modeling phototrophic biofilms in a plug-flow reactor.

The use of phototrophic biofilms in wastewater treatment has been recognized as a potential option for development of new reactor configurations. For better understanding of these systems, a numerical model was developed including relevant microbial processes. As a novelty, this model was implemented in COMSOL Multiphysics, a modern computational environment for complex dynamic models. A two-dimensional biofilm model was used to study the spatial distribution of microbial species within the biofilm and along the length of the reactor. The biofilm model was coupled with a one-dimensional plug-flow bulk liquid model. The impact of different operational conditions on the chemical oxygen demand (COD) and ammonia conversions was assessed. The model was tuned by varying two parameters: the half-saturation coefficient for light use by phototrophs and the oxygen mass transfer coefficient. The mass transfer coefficient was found to be determining for the substrate conversion rate. Simulations indicate that heterotrophs would overgrow nitrifiers and phototrophs within the biofilm until a low biodegradable COD value in the wastewater is reached (organic loading rate <2.32 gCOD/(m(2) d)). This limits the proposed positive effect of treating wastewater with a combination of algae and heterotrophs/autotrophs. Mechanistic models like this one are made for understanding the microbial interactions and their influence on the reactor performance.

Spacer geometry and particle deposition in spiral wound membrane feed channels.

Deposition of microspheres mimicking bacterial cells was studied experimentally and with a numerical model in feed spacer membrane channels, as used in spiral wound nanofiltration (NF) and reverse osmosis (RO) membrane systems. In-situ microscopic observations in membrane fouling simulators revealed formation of specific particle deposition patterns for different diamond and ladder feed spacer orientations. A three-dimensional numerical model combining fluid flow with a Lagrangian approach for particle trajectory calculations could describe very well the in-situ observations on particle deposition in flow cells. Feed spacer geometry, positioning and cross-flow velocity sensitively influenced the particle transport and deposition patterns. The deposition patterns were not influenced by permeate production. This combined experimental-modeling approach could be used for feed spacer geometry optimization studies for reduced (bio)fouling.

A mass-spring model unveils the morphogenesis of phototrophic Diatoma biofilms.

Diatoms often dominate planktonic communities in the ocean and phototrophic biofilms in streams and rivers, greatly contributing to global biogeochemical fluxes. In pelagic ecosystems, these microscopic algae can form chain-like microcolonies, which seem advantageous for nutrient uptake and protect against grazing, and at the same time reduce sinking. Despite the capability of many diatoms to form chains, their contribution to the architecture of phototrophic biofilms remains elusive. Here we propose a computational model to simulate the growth and behaviour of Diatoma chains in contrasting flow environments. This mass-spring mechanical model captures the natural behaviour of Diatoma chains well, emphasising the relevance of chain growth and entanglement for biofilm morphogenesis. The model qualitatively describes formation of intricate dome-shaped structures and of dreadlock-type streamers as observed in nature in multidirectional and unidirectional flow, respectively. The proposed model is a useful tool to study the effect of fluid dynamics on biofilm morphogenesis.

Numerical modelling of tooth enamel subsurface lesion formation induced by dental plaque.

Using a one-dimensional mathematical model that couples tooth demineralisation and remineralisation with metabolic processes occurring in the dental plaque, two mechanisms for subsurface lesion formation were evaluated. It was found that a subsurface lesion can develop only as the result of alternating periods of demineralisation (acid attack during sugar consumption) and remineralisation (resting period) in tooth enamel with uniform mineral composition. It was also shown that a minimum plaque thickness that can induce an enamel lesion exists. The subsurface lesion formation can also be explained by assuming the existence of a fluoride-containing layer at the tooth surface that decreases enamel solubility. A nearly constant thickness of the surface layer was obtained with both proposed mechanisms. Sensitivity analysis showed that surface layer formation is strongly dependent on the length of remineralisation and demineralisation cycles. The restoration period is very important and the numerical simulations support the observation that often consumption of sugars is a key factor in caries formation. The calculated profiles of mineral content in enamel are similar to those observed experimentally. Most probably, both studied mechanisms interact in vivo in the process of caries development, but the simplest explanation for subsurface lesion formation remains the alternation between demineralisation and remineralisation cycles without any pre-imposed gradients.

The granule size distribution in an anammox-based granular sludge reactor affects the conversion--implications for modeling.

Mathematical models are useful tools to optimize the performance of granular sludge reactors. In these models, typically a uniform granule size is assumed for the whole reactor, even though in reality the granules follow a size distribution and the granule size as such affects the process performance. This study assesses the effect of the granule size distribution on the performance of a granular sludge reactor in which autotrophic nitrogen removal is realized through one-stage partial nitritation-anammox. A comparison is made between different approaches to deal with particle size distributions in one-dimensional biofilm models, from the use of a single characteristic diameter to applying a multiple compartment model. The results show a clear impact on the conversion efficiency of the way in which particle size distribution is modeled, resulting from the effect of the granule size on the competition between nitrite oxidizing and anammox bacteria and from the interaction between granules of different sizes in terms of the exchange of solutes. Whereas the use of a uniform granule size is sufficient in case only the overall reactor behavior needs to be assessed, taking into account the detailed granule size distribution is required to study the solute exchange between particles of different sizes. For the latter purpose, the application of the widespread software package Aquasim is limited and the development of dedicated software applications is required.

Effect of granule size on autotrophic nitrogen removal in a granular sludge reactor.

Autotrophic nitrogen removal through sequential partial nitritation and anammox reactions can be achieved in biofilm reactors by controlling the oxygen concentration in the bulk liquid in such a way that nitrite oxidizers are outcompeted by anammox bacteria. In the case of granular sludge reactors, the granule size may influence the optimal range of oxygen concentration, as has been confirmed in the present study by means of numerical simulations. The range of oxygen concentrations for which combined partial nitritation and anammox conversion is established becomes broader for larger particles and with increasing influent ammonium concentrations. At the same time the likelihood of nitrite accumulation in the reactor effluent also increases.

Mathematical model for microbial fuel cells with anodic biofilms and anaerobic digestion.

This study describes the integration of IWA's anaerobic digestion model (ADM1) within a computational model of microbial fuel cells (MFCs). Several populations of methanogenic and electroactive microorganisms coexist suspended in the anolyte and in the biofilm attached to the anode. A number of biological, chemical and electrochemical reactions occur in the bulk liquid, in the biofilm and at the electrode surface, involving glucose, organic acids, H2 and redox mediators. Model output includes the evolution in time of important measurable MFC parameters (current production, consumption of substrates, suspended and attached biomass growth). Two- and three-dimensional model simulations reveal the importance of current and biomass heterogeneous distribution over the planar anode surface. Voltage- and power-current characteristics can be calculated at different moments in time to evaluate the limiting regime in which the MFC operates. Finally, model simulations are compared with experimental results showing that, in a batch MFC, smaller electrical resistance of the circuit leads to selection of electroactive bacteria. Higher coulombic yields are so obtained because electrons from substrate are transferred to anode rather than following the methanogenesis pathway. In addition to higher currents, faster COD consumption rates are so achieved. The potential of this general modelling framework is in the understanding and design of more complex cases of wastewater-fed microbial fuel cells.

Experimental and simulation analysis of community structure of nitrifying bacteria in a membrane-aerated biofilm.

Until now, only few attempts have been made to assess biofilm models simulating microenvironments in a biofilm. As a first step, we compare the microenvironment observed in a membrane aerated biofilm (MAB) to that derived from a two-dimensional computational model with individual ammonia oxidizing bacteria (AOB) and nitrite oxidizing bacteria (NOB) embedded in a continuum EPS matrix. Gradients of oxygen were determined by means of microelectrodes. The change in nitrifying bacterial populations with the biofilm depth was quantified using fluorescence in situ hybridization (FISH) in combination with a confocal laser scanning microscopy (CLSM). Microelectrode measurements revealed that oxic and anoxic or anaerobic regions exist within the MAB. The oxygen profile predicted by the model showed good agreement with that obtained by microelectrode measurements. The oxic part of the biofilm was dominated by NSO190 probe-hybridized AOB, which formed relatively large clusters of cells directly on the membrane surface, and by the NOB belonging to genus Nitrobacter sp. On the other hand, NOB belonging to genus Nitrospira sp. were abundant at the oxic-anoxic interface. The model prediction regarding AOB and Nitrobacter sp. distribution was consistent with the experimental counterpart. Measurements of AOB cluster size distribution showed that colonies are slightly larger adjacent to the membrane than at the inner part of the biofilm. The sizes predicted by the current model are larger than those obtained in the experiment, leading to the arguments that some factors not contained in the model would affect the cluster size.

Microbial motility involvement in biofilm structure formation--a 3D modelling study.

A computational model explaining formation of mushroom-like biofilm colonies is proposed in this study. The biofilm model combines for the first time cell growth with twitching motility in a three-dimensional individual-based approach. Model simulations describe the tendency of motile cells to form flat biofilms spreading out on the substratum, in contrast with the immotile variants that form only round colonies. These computational results are in good qualitative agreement with the experimental data obtained from Pseudomonas aeruginosa biofilms grown in flowcells. Simulations reveal that motile cells can possess a serious ecological advantage by becoming less affected by mass transfer limitations. Twitching motility alone appears to be insufficient to generate mushroom-like biofilm structures with caps on stalks. Rather, a substrate limitation-induced detachment of motile cells followed by reattachment could explain this intriguing effect leading to higher-level biofilm structure.

Kinetic model of a granular sludge SBR: influences on nutrient removal.

A mathematical model was developed that can be used to describe an aerobic granular sludge reactor, fed with a defined influent, capable of simultaneously removing COD, nitrogen and phosphate in one sequencing batch reactor (SBR). The model described the experimental data from this complex system sufficiently. The effect of process parameters on the nutrient removal rates could therefore be reliably evaluated. The influence of oxygen concentration, temperature, granule diameter, sludge loading rate, and cycle configuration were analyzed. Oxygen penetration depth in combination with the position of the autotrophic biomass played a crucial role in the conversion rates of the different components and thus on overall nutrient removal efficiencies. The ratio between aerobic and anoxic volume in the granule strongly determines the N-removal efficiency as it was shown by model simulations with varying oxygen concentration, temperature, and granule size. The optimum granule diameter for maximum N- and P-removal in the standard case operating conditions (DO 2 mg L(-1), 20 degrees C) was found between 1.2 and 1.4 mm and the optimum COD loading rate was 1.9 kg COD m(-3) day(-1). When all ammonia is oxidized, oxygen diffuses to the core of the granule inhibiting the denitrification process. In order to optimize the process, anoxic phases can be implemented in the SBR-cycle configuration, leading to a more efficient overall N-removal. Phosphate removal efficiency mainly depends on the sludge age; if the SRT exceeds 30 days not enough biomass is removed from the system to keep effluent phosphate concentrations low.

Multidimensional modelling to investigate interspecies hydrogen transfer in anaerobic biofilms.

Anaerobic digestion is a multistep process, mediated by a functionally and phylogenetically diverse microbial population. One of the crucial steps is oxidation of organic acids, with electron transfer via hydrogen or formate from acetogenic bacteria to methanogens. This syntrophic microbiological process is strongly restricted by a thermodynamic limitation on the allowable hydrogen or formate concentration. In order to study this process in more detail, we developed an individual-based biofilm model which enables to describe the processes at a microbial resolution. The biochemical model is the ADM1, implemented in a multidimensional domain. With this model, we evaluated three important issues for the syntrophic relationship: (i) Is there a fundamental difference in using hydrogen or formate as electron carrier? (ii) Does a thermodynamic-based inhibition function produced substantially different results from an empirical function? and; (iii) Does the physical co-location of acetogens and methanogens follow directly from a general model. Hydrogen or formate as electron carrier had no substantial impact on model results. Standard inhibition functions or thermodynamic inhibition function gave similar results at larger substrate field grid sizes (> 10 microm), but at smaller grid sizes, the thermodynamic-based function reduced the number of cells with long interspecies distances (> 2.5 microm). Therefore, a very fine grid resolution is needed to reflect differences between the thermodynamic function, and a more generic inhibition form. The co-location of syntrophic bacteria was well predicted without a need to assume a microbiological based mechanism (e.g., through chemotaxis) of biofilm formation.

Global impact and application of the anaerobic ammonium-oxidizing (anammox) bacteria.

In the anaerobic ammonium oxidation (anammox) process, ammonia is oxidized with nitrite as primary electron acceptor under strictly anoxic conditions. The reaction is catalysed by a specialized group of planctomycete-like bacteria. These anammox bacteria use a complex reaction mechanism involving hydrazine as an intermediate. The reactions are assumed to be carried out in a unique prokaryotic organelle, the anammoxosome. This organelle is surrounded by ladderane lipids, which make the organelle nearly impermeable to hydrazine and protons. The localization of the major anammox protein, hydrazine oxidoreductase, was determined via immunogold labelling to be inside the anammoxosome. The anammox bacteria have been detected in many marine and freshwater ecosystems and were estimated to contribute up to 50% of oceanic nitrogen loss. Furthermore, the anammox process is currently implemented in water treatment for the low-cost removal of ammonia from high-strength waste streams. Recent findings suggested that the anammox bacteria may also use organic acids to convert nitrate and nitrite into dinitrogen gas when ammonia is in short supply.

Role of nitrogen oxides in the metabolism of ammonia-oxidizing bacteria.

Ammonia-oxidizing bacteria (AOB) can use oxygen and nitrite as electron acceptors. Nitrite reduction by Nitrosomonas is observed under three conditions: (i) hydrogen-dependent denitrification, (ii) anoxic ammonia oxidation with nitrogen dioxide (NO(2)) and (iii) NO(x)-induced aerobic ammonia oxidation. NO(x) molecules play an important role in the conversion of ammonia and nitrite by AOB. Absence of nitric oxide (NO), which is generally detectable during ammonia oxidation, severely impairs ammonia oxidation by AOB. The lag phase of recovery of aerobic ammonia oxidation was significantly reduced by NO(2) addition. Acetylene inhibition tests showed that NO(2)-dependent and oxygen-dependent ammonia oxidation can be distinguished. Addition of NO(x) increased specific activity of ammonia oxidation, growth rate and denitrification capacity. Together, these findings resulted in a hypothetical model on the role of NO(x) in ammonia oxidation: the NO(x) cycle.

Multidimensional modelling of anaerobic granules.

A multispecies, two- and three-dimensional model was developed, based on a previously published planar biofilm model, and the biochemical structure of the ADM1. Several soluble substrates diffuse and react in the granule. Local pH is calculated from acid-base equilibria and charge balance. The model uses individual-based representation of biomass particles within the granule (biofilm), and describes spreading by an iterative pushing technique. The overall computational domain consists of one granule, and is divided into a grid with Cartesian coordinates. The number of grid elements does not limit the number of biomass particles, and it is not necessary to use grid-spreading techniques, such as cellular automata, which result in Cartesian artefacts. The model represents both microscopic and macroscopic features in granule structure, previously observed using in-situ molecular techniques, and can be effectively used to interpret these results.

Comparing biofilm models for a single species biofilm system.

A benchmark problem was defined to evaluate the performance of different mathematical biofilm models. The biofilm consisted of heterotrophic bacteria degrading organic substrate and oxygen. Mathematical models tested ranged from simple analytical to multidimensional numerical models. For simple and more or less flat biofilms it was shown that analytical biofilm models provide very similar results compared to more complex numerical solutions. When considering a heterogeneous biofilm morphology it was shown that the effect of an increased external mass transfer resistance was much more significant compared to the effect of an increased surface area inside the biofilm.

Modelling a spatially heterogeneous biofilm and the bulk fluid: selected results from benchmark problem 2 (BM2).

The numerical simulation of mass transfer and conversion in spatially heterogeneous biofilms on the meso-scale requires an accurate description of the hydrodynamics in the biofilm systems and of spatial effects. This leads to systems of three-dimensional nonlinear partial differential equations that are numerically very expensive to solve and to data requirements that are not easy to meet. In this paper several modeling approaches to reduce the physical complexity and, hence, accelerate the computation are compared. They range from a mere reduction of dimensionality by lumping the problem along a secondary flow direction to global mass balances or empirical correlations, at the core of which a one-dimensional boundary value problem must be solved. It is found that even strongly simplified models can describe the qualitative behaviour of the model with regard to variations in the geometrical and hydrodynamic model parameters quite well. In order to obtain also quantitatively reliable results the hydrodynamics must be considered in an appropriate manner.

Results from the multi-species benchmark problem 3 (BM3) using two-dimensional models.

In addition to the one-dimensional solutions of a multi-species benchmark problem (BM3) presented earlier (Rittmann et al., 2004), we offer solutions using two-dimensional (2-D) models. Both 2-D models (called here DN and CP) used numerical solutions to BM3 based on a similar mathematical framework of the one-dimensional AQUASIM-built models submitted by Wanner (model W) and Morgenroth (model M1), described in detail elsewhere (Rittmann et al., 2004). The CP model used differential equations to simulate substrate gradients and biomass growth and a particle-based approach to describe biomass division and biofilm growth. The DN model simulated substrate and biomass using a cellular automaton approach. For several conditions stipulated in BM3, the multidimensional models provided very similar results to the 1-D models in terms of bulk substrate concentrations and fluxes into the biofilm. The similarity can be attributed to the definition of BM3, which restricted the problem to a flat biofilm in contact with a completely mixed liquid phase, and therefore, without any salient characteristics to be captured in a multidimensional domain. On the other hand, the models predicted significantly different accumulations of the different types of biomass, likely reflecting differences in the way biomass spread within the biofilm is simulated.