Modeling Ion Transport across Thin-Film Composite Membranes During Saltwater Electrolysis.

Affordable thin-film composite (TFC) membranes are a potential alternative to more expensive ion exchange membranes in saltwater electrolyzers used for hydrogen gas production. We used a solution-friction transport model to study how the induced potential gradient controls ion transport across the polyamide (PA) active layer and support layers of TFC membranes during electrolysis. The set of parameters was simplified by assigning the same size-related partition and friction coefficients for all salt ions through the membrane active layer. The model was fit to experimental ion transport data from saltwater electrolysis with 600 mM electrolytes at a current density of 10 mA cm-2. When the electrolyte concentration and current density were increased, the transport of major charge carriers was successfully predicted by the model. Ion transport calculated using the model only minimally changed when the negative active layer charge density was varied from 0 to 600 mM, indicating active layer charge was not largely responsible for controlling ion crossover during electrolysis. Based on model simulations, a sharp pH gradient was predicted to occur within the supporting layer of the membrane. These results can help guide membrane design and operation conditions in water electrolyzers using TFC membranes.

Relative Insignificance of Polyamide Layer Selectivity for Seawater Electrolysis Applications.

Low-cost polyamide thin-film composite (TFC) membranes are being explored as alternatives to cation exchange membranes for seawater electrolysis. An optimal membrane should have a low electrical resistance to minimize applied potentials needed for water electrolysis and be able to block chloride ions present in a seawater catholyte from reaching the anode. The largest energy loss associated with a TFC membrane was the Nernstian overpotential of 0.74 V (equivalent to 37 Ω cm2 at 20 mA cm-2), derived from the pH difference between the anolyte and catholyte and not the membrane ohmic overpotential. Based on analysis using electrochemical impedance spectroscopy, the pristine TFC membrane contributed only 5.00 Ω cm2 to the ohmic resistance. Removing the polyester support layer reduced the resistance by 79% to only 1.04 Ω cm2, without altering the salt ion transport between the electrolytes. Enlarging the pore size (∼5 times) in the polyamide active layer minimally impacted counterion transport across the membrane during electrolysis, but it increased the total concentration of chloride transported by 60%. Overall, this study suggests that TFC membranes with thinner but mechanically strong supporting layers and size-selective active layers should reduce energy consumption and the potential for chlorine generation for seawater electrolyzers.

Quantitative chemical biosensing by bacterial chemotaxis in microfluidic chips.

Whole-cell bacterial bioreporters are proposed as alternatives to chemical analysis of, for example, pollutants in environmental compartments. Commonly based on reporter gene induction, bioreporters produce a detectable signal within 30 min to a few hours after exposure to the chemical target, which is impractical for applications aiming at a fast response. In an attempt to attain faster readout but maintain flexibility of chemical targeting, we explored the concept for quantitative chemical sensing by bacterial chemotaxis. Chemotaxis was quantified from enrichment of cells across a 600 µm-wide chemical gradient stabilized by parallel flow in a microfluidic chip, further supported by transport and chemotaxis steady state and kinetic modelling. As proof-of-concept, we quantified Escherichia coli chemotaxis towards serine, aspartate and methylaspartate as a function of attractant concentration and exposure time. E. coli chemotaxis enrichment increased sharply between 0 and 10 µM serine, before saturating at 100 µM. The chemotaxis accumulation rate was maximal at 10 µM serine, leading to observable cell enrichment within 5 min. The potential application for biosensing of environmental toxicants was investigated by quantifying chemotaxis of Cupriavidus pinatubonensis JMP134 towards the herbicide 2,4-dichlorophenoxyacetate. Our results show that bacterial chemotaxis can be quantified on a scale of minutes and may be used for developing faster bioreporter assays.

Mechanisms of nitrous oxide (N2 O) formation and reduction in denitrifying biofilms.

Nitrous oxide (N2 O) is a potent greenhouse gas that can be formed in wastewater treatment processes by ammonium oxidizing and denitrifying microorganisms. While N2 O emissions from suspended growth systems have been extensively studied, and some recent studies have addressed emissions from nitrifying biofilms, much less is known about N2 O emissions from denitrifying biofilm processes. This research used modeling to evaluate the mechanisms of N2 O formation and reduction in denitrifying biofilms. The kinetic model included formation and consumption of key denitrification species, including nitrate (NO3-), nitrite (NO2-), nitric oxide (NO), and N2 O. The model showed that, in presence of excess of electron donor, denitrifying biofilms have two distinct layers of activity: an outer layer where there is net production of N2 O and an inner layer where there is net consumption. The presence of oxygen (O2 ) had an important effect on N2 O emission from suspended growth systems, but a smaller effect on biofilm systems. The effects of NO3- and O2 differed significantly based on the biofilm thickness. Overall, the effects of biofilm thickness and bulk substrate concentrations on N2 O emissions are complex and not always intuitive. A key mechanism for denitrifying biofilms is the diffusion of N2 O and other intermediates from one zone of the biofilm to another. This leads to zones of N2 O formation or consumption transformations that would not exist in suspended growth systems.

Predicting N2O emissions from nitrifying and denitrifying biofilms: a modeling study.

Wastewater treatment plants can be significant sources of nitrous oxide (N2O), a potent greenhouse gas. While our understanding of N2O emissions from suspended-growth processes has advanced significantly, less is known about emissions from biofilm processes. Biofilms may behave differently due to their substrate gradients and microbial stratification. In this study, we used mathematical modeling to explore the mechanisms of N2O emissions from nitrifying and denitrifying biofilms. Our ammonia-oxidizing bacteria biofilm model suggests that N2O emissions from biofilm can be significantly greater than from suspended-growth systems. The driving factor is the diffusion of hydroxylamine, a nitrification intermediate, from the aerobic to the anoxic regions of the biofilm. The presence of nitrite-oxidizing bacteria further increased emissions. For denitrifying biofilms, our results suggest that emissions are generally greater than for suspended-growth systems. However, the magnitude of the difference depends on the bulk dissolved oxygen, chemical oxygen demand, and nitrate concentrations, as well as the biofilm thickness. Overall, the accumulation and diffusion of key intermediates, i.e. hydroxylamine and nitrite, distinguish biofilms from suspended-growth systems. Our research suggests that the mechanisms of N2O emissions from biofilms are much more complex than suspended-growth systems, and that emissions may be higher in many cases.

Short-range guiding can result in the formation of circular aggregates in myxobacteria populations.

Myxobacteria are social bacteria that upon starvation form multicellular fruiting bodies whose shape in different species can range from simple mounds to elaborate tree-like structures. The formation of fruiting bodies is a result of collective cell movement on a solid surface. In the course of development, groups of flexible rod-shaped cells form streams and move in circular or spiral patterns to form aggregation centers that can become sites of fruiting body formation. The mechanisms of such cell movement patterns are not well understood. It has been suggested that myxobacterial development depends on short-range contact-mediated interactions between individual cells, i.e. cell aggregation does not require long-range signaling in the population. In this study, by means of a computational mass-spring model, we investigate what types of short-range interactions between cells can result in the formation of streams and circular aggregates during myxobacterial development. We consider short-range head-to-tail guiding between individual cells, whereby movement direction of the head of one cell is affected by the nearby presence of the tail of another cell. We demonstrate that stable streams and circular aggregates can arise only when the trailing cell, in addition to being steered by the tail of the leading cell, is able to speed up to catch up with it. It is suggested that necessary head-to-tail interactions between cells can arise from physical adhesion, response to a diffusible substance or slime extruded by cells, or pulling by motility engine pili. Finally, we consider a case of long-range guiding between cells and show that circular aggregates are able to form without cells increasing speed. These findings present a possibility to discriminate between short-range and long-range guiding mechanisms in myxobacteria by experimentally measuring distribution of cell speeds in circular aggregates.

Assessing microbial competition in a hydrogen-based membrane biofilm reactor (MBfR) using multidimensional modeling.

The membrane biofilm reactor (MBfR) is a novel technology that safely delivers hydrogen to the base of a denitrifying biofilm via gas-supplying membranes. While hydrogen is an effective electron donor for denitrifying bacteria (DNB), it also supports sulfate-reducing bacteria (SRB) and methanogens (MET), which consume hydrogen and create undesirable by-products. SRB and MET are only competitive for hydrogen when local nitrate concentrations are low, therefore SRB and MET primarily grow near the base of the biofilm. In an MBfR, hydrogen concentrations are greatest at the base of the biofilm, making SRB and MET more likely to proliferate in an MBfR system than a conventional biofilm reactor. Modeling results showed that because of this, control of the hydrogen concentration via the intramembrane pressure was a key tool for limiting SRB and MET development. Another means is biofilm management, which supported both sloughing and erosive detachment. For the conditions simulated, maintaining thinner biofilms promoted higher denitrification fluxes and limited the presence of SRB and MET. The 2-d modeling showed that periodic biofilm sloughing helped control slow-growing SRB and MET. Moreover, the rough (non-flat) membrane assembly in the 2-d model provided a special niche for SRB and MET that was not represented in the 1-d model. This study compared 1-d and 2-d biofilm model applicability for simulating competition in counter-diffusional biofilms. Although more computationally expensive, the 2-d model captured important mechanisms unseen in the 1-d model.

Hydroxylamine diffusion can enhance N2O emissions in nitrifying biofilms: a modeling study.

Wastewater treatment plants can be significant sources of nitrous oxide (N2O), a potent greenhouse gas. However, little is known about N2O emissions from biofilm processes. We adapted an existing suspended-growth mathematical model to explore N2O emissions from nitrifying biofilms. The model included N2O formation by ammonia-oxidizing bacteria (AOB) via the hydroxylamine and the nitrifier denitrification pathways. Our model suggested that N2O emissions from nitrifying biofilms could be significantly greater than from suspended growth systems under similar conditions. The main cause was the formation and diffusion of hydroxylamine, an AOB nitrification intermediate, from the aerobic to the anoxic regions of the biofilm. In the anoxic regions, hydroxylamine oxidation by AOB provided reducing equivalents used solely for nitrite reduction to N2O, since there was no competition with oxygen. For a continuous system, very high and very low dissolved oxygen (DO) concentrations resulted in lower emissions, while intermediate values led to higher emissions. Higher bulk ammonia concentrations and greater biofilm thicknesses increased emissions. The model effectively predicted N2O emissions from an actual pilot-scale granular sludge reactor for sidestream nitritation, but significantly underestimated the emissions when the NH2OH diffusion coefficient was assumed to be minimal. This numerical study suggests an unexpected and important role of hydroxylamine in N2O emission in biofilms.

Variable cell morphology approach for individual-based modeling of microbial communities.

An individual-based, mass-spring modeling framework has been developed to investigate the effect of cell properties on the structure of biofilms and microbial aggregates through Lagrangian modeling. Key features that distinguish this model are variable cell morphology described by a collection of particles connected by springs and a mechanical representation of deformable intracellular, intercellular, and cell-substratum links. A first case study describes the colony formation of a rod-shaped species on a planar substratum. This case shows the importance of mechanical interactions in a community of growing and dividing rod-shaped cells (i.e., bacilli). Cell-substratum links promote formation of mounds as opposed to single-layer biofilms, whereas filial links affect the roundness of the biofilm. A second case study describes the formation of flocs and development of external filaments in a mixed-culture activated sludge community. It is shown by modeling that distinct cell-cell links, microbial morphology, and growth kinetics can lead to excessive filamentous proliferation and interfloc bridging, possible causes for detrimental sludge bulking. This methodology has been extended to more advanced microbial morphologies such as filament branching and proves to be a very powerful tool in determining how fundamental controlling mechanisms determine diverse microbial colony architectures.

Effect of heterotrophic growth on autotrophic nitrogen removal in a granular sludge reactor.

This study deals with the influence of heterotrophic growth on autotrophic nitrogen removal from wastewater in a granular sludge reactor. A mathematical model was set-up including autotrophic and heterotrophic growth and decay in the granules from a partial nitritation-anammox process. A distinction between heterotrophic bacteria was made based on the electron acceptor (dissolved oxygen, nitrite or nitrate) on which they grow, while the nitrogen gas produced was 'labelled' to retrieve its origin, from anammox or heterotrophic bacteria. Taking into account heterotrophic growth resulted in a lower initial nitrogen removal, but in a higher steady state nitrogen removal compared with a model in which heterotrophic growth was neglected. The anammox activity is related with the fact that heterotrophs initially use nitrite as electron acceptor, but when they switch to nitrate the produced nitrite can be used by anammox bacteria. Increased anammox activity in the presence of heterotrophs, therefore, resulted in a marginally increased N2 production at steady state. Heterotrophic denitrification of nitrate to nitrite also explains why small amounts of organic substrate present in the influent positively affect the maximum nitrogen removal capacity. However, the process efficiency deteriorates once the amount of organic substrate in the influent exceeds a certain threshold. The bulk oxygen concentration and the granule size have a dual effect on the autotrophic nitrogen removal efficiency. Besides, the maximum nitrogen removal efficiency decreases and the corresponding optimal bulk oxygen concentration increases with increasing granule size.

Biofilm development and the dynamics of preferential flow paths in porous media.

A two-dimensional pore-scale numerical model was developed to evaluate the dynamics of preferential flow paths in porous media caused by bioclogging. The liquid flow and solute transport through the pore network were coupled with a biofilm model including biomass attachment, growth, decay, lysis, and detachment. Blocking of all but one flow path was obtained under constant liquid inlet flow rate and biomass detachment caused by shear forces only. The stable flow path formed when biofilm detachment balances growth, even with biomass weakened by decay. However, shear forces combined with biomass lysis upon starvation could produce an intermittently shifting location of flow channels. Dynamic flow pathways may also occur when combined liquid shear and pressure forces act on the biofilm. In spite of repeated clogging and unclogging of interconnected pore spaces, the average permeability reached a quasi-constant value. Oscillations in the medium permeability were more pronounced for weaker biofilms.

Unique stratification of biofilm density in heterotrophic membrane-aerated biofilms: An experimental and modeling study.

We consistently find a band of high cell density develop within heterotrophic membrane-aerated biofilms. This study reports and attempts to explain this unique behavior. Biofilm density affects volumetric reaction rates, biofilm growth rates, substrate diffusion, and mechanical behavior. Yet the mechanisms and dynamics of biofilm density development are poorly understood. In this study, a membrane-aerated biofilm, where O2 was supplied from the base of the biofilm and acetate from the bulk liquid, was used to explore spatial and temporal patterns of density development. Biofilm density was assessed by optical coherence tomography. After inoculation, the biofilm quickly increased in thickness, with a low density throughout. However, as the biofilm reached a stable thickness of around 1000 µm, a high-density layer developed in the biofilm interior. The layer slowly expanded over time. Oxygen microprofiles in the biofilm showed this layer coincided with the most metabolically active zone, resulting from counter-diffusing O2 and acetate. The formation of this dense layer appeared to be related to changes in growth rates. Initially, high growth rates throughout the biofilm presumably led to fast-growing, low-density biofilms. As the biofilm became thicker, and as substrates became limiting in the biofilm interior, growth rates decreased, resulting in new growth at a higher density. A 1-D mathematical model with variable biofilm density was developed by linking the rates of extracellular polymeric substances (EPS) production to the growth rate. The model captured the initial fast growth at a low density, followed by a slower, substrate-limited growth in the biofilm interior, producing a dense band within the biofilm. Together, these results suggest that low growth rates can lead to high-density zones within the interior of counter-diffusional biofilms. These findings should also be relevant to conventional, co-diffusional biofilms, although differences in density may be less obvious.

Downscaling Industrial-Scale Syngas Fermentation to Simulate Frequent and Irregular Dissolved Gas Concentration Shocks.

In large-scale syngas fermentation, strong gradients in dissolved gas (CO, H2) concentrations are very likely to occur due to locally varying mass transfer and convection rates. Using Euler-Lagrangian CFD simulations, we analyzed these gradients in an industrial-scale external-loop gas-lift reactor (EL-GLR) for a wide range of biomass concentrations, considering CO inhibition for both CO and H2 uptake. Lifeline analyses showed that micro-organisms are likely to experience frequent (5 to 30 s) oscillations in dissolved gas concentrations with one order of magnitude. From the lifeline analyses, we developed a conceptual scale-down simulator (stirred-tank reactor with varying stirrer speed) to replicate industrial-scale environmental fluctuations at bench scale. The configuration of the scale-down simulator can be adjusted to match a broad range of environmental fluctuations. Our results suggest a preference for industrial operation at high biomass concentrations, as this would strongly reduce inhibitory effects, provide operational flexibility and enhance the product yield. The peaks in dissolved gas concentration were hypothesized to increase the syngas-to-ethanol yield due to the fast uptake mechanisms in C. autoethanogenum. The proposed scale-down simulator can be used to validate such results and to obtain data for parametrizing lumped kinetic metabolic models that describe such short-term responses.

Mass transfer enhancement in moving biofilm structures.

Biofilms are layers of microbial cells growing on an interface and they can form highly complex structures adapted to a wide variety of environmental conditions. Biofilm streamers have a small immobile base attached to the support and a flexible tail elongated in the flow direction, which can vibrate in fast flows. Herein we report numerical results for the role of the periodical movement of biofilm streamers on the nutrient uptake and in general on the solute mass transfer enhancement due to flow-induced oscillations. We developed what to our knowledge is a novel two-dimensional fluid-structure interaction model coupled to unsteady solute mass transport and solved the model using the finite element method with a moving mesh. Results demonstrate that the oscillatory movement of the biofilm tail significantly increases the substrate uptake. The mass transfer coefficient is the highest in regions close to the streamer tip. The reason for substrate transfer enhancement is the increase in speed of tip movement relative to the surrounding liquid, thereby reducing the thickness of the mass transfer boundary layer. In addition, we show that the relative mass transfer enhancement in unsteady conditions compared with the rigid static structure is larger at higher flow velocities, and this relative increase favors a more flexible structure.

The effect of biofilm permeability on bio-clogging of porous media.

A 3D Biofilm model, appropriate for complex porous media support structures, is successfully modified such that non-zero permeability of biofilms structures is enabled. A systematic study is then conducted into the influence of biofilm permeability on overall biomass growth rate. This reveals a significant influence at large biofilm concentrations; even when the permeability of the biomass is 1.25% of that of the free pore space, biomass accumulation increased by a factor of ∼3 over 40 h. The effect is shown to be retained when allowing for biomass detachment or erosion as a consequence of adjacent velocity shear. We conclude that biofilm permeability should be included in biofilm models and that further experimental work is required to better describe the link between biofilm permeability and local microstructure.

Mathematical modelling of tooth demineralisation and pH profiles in dental plaque.

A mathematical model of dental plaque has been developed in order to investigate the processes leading to dental caries. The one-dimensional time-dependent model integrates existing knowledge on biofilm processes (mass transfer, microbial composition, microbial conversions and substrate availability) with tooth demineralisation kinetics. This work is based on the pioneering studies of Dibdin who, nearly two decades ago, build a mathematical model roughly describing the metabolic processes taking place in dental plaque. We extended Dibdin's model with: multiple microbial species (aciduric and non-aciduric Streptococci, Actinomyces and Veillonella), more metabolic processes (i.e., aerobic and anaerobic glucose conversion, low and high glucose uptake affinity pathways, formation and consumption of storage compounds), ion transport by Nernst-Planck equations, and we coupled the obtained pH and chemical component gradients inside the plaque with tooth demineralisation. The new model implementation was complemented with faster and more rigorous numerical methods for the model solution. Model results confirm the protective effect of Veillonella due to lactate consumption. Interestingly, on short term, the storage compounds may not necessarily have a negative effect on demineralisation. Individual feeding patterns can also be easily studied with this model. For example, slow ("social") consumption of sugar-containing drinks proves to be more harmful than drinking the same amount over a short period of time.

Structured morphological modeling as a framework for rational strain design of Streptomyces species.

Successful application of a computational model for rational design of industrial Streptomyces exploitation requires a better understanding of the relationship between morphology-dictated by microbial growth, branching, fragmentation and adhesion-and product formation. Here we review the state-of-the-art in modeling of growth and product formation by filamentous microorganisms and expand on existing models by combining a morphological and structural approach to realistically model and visualize a three-dimensional pellet. The objective is to provide a framework to study the effect of morphology and structure on natural product and enzyme formation and yield. Growth and development of the pellet occur via the processes of apical extension, branching and cross-wall formation. Oxygen is taken to be the limiting component, with the oxygen concentration at the tips regulating growth kinetics and the oxygen profile within the pellet affecting the probability of branching. Biological information regarding the processes of differentiation and branching in liquid cultures of the model organism Streptomyces coelicolor has been implemented. The model can be extended based on information gained in fermentation trials for different production strains, with the aim to provide a test drive for the fermentation process and to pre-assess the effect of different variables on productivity. This should aid in improving Streptomyces as a production platform in industrial biotechnology.

New approaches to characterizing and understanding biofouling of spiral wound membrane systems.

Historically, biofouling research on spiral wound membrane systems is typically problem solving oriented. Membrane modules are studied as black box systems, investigated by autopsies. Biofouling is not a simple process. Many factors influence each other in a non-linear fashion. These features make biofouling a subject which is not easy to study using a fundamental scientific approach. Nevertheless to solve or minimize the negative impacts of biofouling, a clear understanding of the interacting basic principles is needed. Recent research into microbiological characterizing of biofouling, small scale test units, application of in situ visualization methods, and model approaches allow such an integrated study of biofouling.

No model is perfect, but some are useful.

Agent-based model should inform the action plan to curb algal blooms in Lake Erie.

Controlling the hydraulic resistance of membrane biofilms by engineering biofilm physical structure.

The application of membrane technology for water treatment and reuse is hampered by the development of a microbial biofilm. Biofilm growth in micro-and ultrafiltration (MF/UF) membrane modules, on both the membrane surface and feed spacer, can form a secondary membrane and exert resistance to permeation and crossflow, increasing energy demand and decreasing permeate quantity and quality. In recent years, exhaustive efforts were made to understand the chemical, structural and hydraulic characteristics of membrane biofilms. In this review, we critically assess which specific structural features of membrane biofilms exert resistance to forced water passage in MF/UF membranes systems applied to water and wastewater treatment, and how biofilm physical structure can be engineered by process operation to impose less hydraulic resistance ("below-the-pain threshold"). Counter-intuitively, biofilms with greater thickness do not always cause a higher hydraulic resistance than thinner biofilms. Dense biofilms, however, had consistently higher hydraulic resistances compared to less dense biofilms. The mechanism by which density exerts hydraulic resistance is reported in the literature to be dependant on the biofilms' internal packing structure and EPS chemical composition (e.g., porosity, polymer concentration). Current reports of internal porosity in membrane biofilms are not supported by adequate experimental evidence or by a reliable methodology, limiting a unified understanding of biofilm internal structure. Identifying the dependency of hydraulic resistance on biofilm density invites efforts to control the hydraulic resistance of membrane biofilms by engineering internal biofilm structure. Regulation of biofilm internal structure is possible by alteration of key determinants such as feed water nutrient composition/concentration, hydraulic shear stress and resistance and can engineer biofilm structural development to decrease density and therein hydraulic resistance. Future efforts should seek to determine the extent to which the concept of "biofilm engineering" can be extended to other biofilm parameters such as mechanical stability and the implication for biofilm control/removal in engineered water systems (e.g., pipelines and/or, cooling towers) susceptible to biofouling.