
MBW:Dynamic Biodegradation model with FlocculationFrom MathBioMBW: Flocculated Growth of Microbes Article Review by Robert Griffin Hale Brandt, Bernd W., and Sebastiaan A. L. M. Kooijman. "Two Parameters Account for the Flocculated Growth of Microbes in Biodegradation Assays." Biotechnology and Bioengineering 70.6 (2000): 67784.
ContentsSummaryBrandt and Kooijan focus on modeling the dynamics of biodegradation of organic compounds by microbes. Previous models and data assume individual cell suspension and neglect the formation of flocs or aggregated microbes. Flocs can greatly reduce biodegradation because cells near the center of the floc can either be considered dead biomass or considered to have a reduced consumption. Brandt and Kooijan use a simplistic model to include the effects of flocs. It is assumed that flocs have two phases and are modeled as spheres. Initially small flocs increase biomass at an exponentially rate. At a certain volume, the floc transition to a linear growth rate when a dead ‘kernel’ forms in the center of the biomass. It is also assumed that the substrate is held at a constant concentration exterior to the floc. The transition is a result of diffusionlimited substrate, which dictates the radius of living microbes. Beyond the maximum radius, cells no longer receive adequate nutrients and die. Eventually the flocs will grow to a critical radius and rupture. The dead biomass is released and the surface of critical floc is then assumed to separate into smaller flocs and the cycle repeats. Context/Biological phenomenon under considerationBiodegradation is the main method for the processing of household chemicals in wastewater treatment plants. Commonly, the chemicals are added to a large tank of ‘sludge’ that contains live microbial organisms. These microbes feed on the chemicals, metabolizing them to a stable or less harmful state. Many studies of the exact processes within the plants’ tanks have been done. The researchers try to model exact growth rates of the microbes and subsequent breakdown of the substrate being controlled. Often however, the models used assume that the microbes are dispersed throughout the tank evenly; as in an idealized suspension. By assuming true suspension, these models are ignoring the fact that microbes often form coagulates or floc. The process of these clusters forming and growing is known as flocculation. Flocculation has obvious, and sofar unaccounted, effects on the dynamics of a wastewater treatment facility.
Once the kernel has formed, the increase in volume is determined solely by the growth of the living exterior layer. Thus, growth now slows to a roughly linear rate, similar to the tapering of a logistic growth model. The spherical floc will continue to grow, albeit at a reduced speed, until a second critical size is reached. At this point, the floc is no longer physically stable and will fragment. The dead core is released into the environment, and the thin living layer breaks up into multiple small pieces, each of which form its own new floc. At this point the process begins again for the small, kernelless flocs. The critical size for selfdestruction is determined by many factors: the adhesiveness of neighboring cells and their density, and most importantly the sheer forces felt from perturbations in the water, due to stirring and fluid convection. HistoryIn today's hyperengineered and critically designed world, the applications for fluid dynamic models are huge. Applications range from sewage water treatment and large scale water distribution to raw water collection and largescale reservoirs. Each of these areas contain large lists of variables that affect the flow and circulation of the water. Often, the systems are so large and complex that any type of controlled experiment would be far too expensive to construct. Thus, to save construction costs and valuable resources, computer simulated models are used instead, to best investigate the often hardtoobserve phenomena. Using the computer models, engineers can work together with the mathematicians to best increase efficiency in cost management and also control risks to individuals in the plants and those in the general public. Before any largescale system is installed, rigorous tests can be applied via simulations to determine whether the design is worth constructing. In addition, an engineer can also simulation the effect of minor alterations to an already existing infrastructure. This will prove invaluable in the coming years as we struggle to increase efficiency to save valuable resources and energy. Using computer models, we will be able to systematically identify the weak points of a system and update them in the most effective ways. Current models take into account many factors on all scales; ranging from the micro to the macro. For example, effective water treatment requires a strong grasp of the biological systems developing, the chemical reactions and mass diffusion, as well as the hydraulic conditions caused by the treatment system being used. A true insight into how the hydrodynamics and chemical equations effect eachother can lead to novel treatment standards and practices that can help to save time, money, and lives.Today, the chemical, mechanical and manufacturing industries use computational fluid dynamic models to test their products and systems. The water treatment sector has lagged behind slightly, but has recently seen a large movement towards improved models and applications. In the United States and abroad, fluid dynamics models have begun to lead to increased understanding, and innovation. However, most of these models glaze over or ignore completely the effects of flocculation on the results they are looking for. Often flocculation is removed from the equations because of it is difficult to accurately describe mathematically and measure in experiments. The inherent fractal nature of floc formation makes accurate spatial models of the phenomena computationally taxing. To adapt these models and make them temporally dynamic would require more computing power than is available to the standard engineer, or even larger corporations. The size of the data sets, and dimensions of the variables would be staggering to consider using current models. Such variables are floc size distribution, shape, composition, density, and structural stability. In addition, there are various behaviors and processes related to flocculation including aggregation, sedimentation, fragmentation, deposition and suspension. For this reason, a novel approach to describing flocculation is required. Parameter Definitionsdiffusion coefficient energy density maximum energy density scaled functional response scaled functional response average scaled response
specific energy conductance: maintenance rate coefficient:
massspecific maintenance rate scaled radius at division scaled maximum thickness of living layer: scaled total radius of floc:
radius at birth radius at division maximum thickness of living layer of floc total radius of floc
amount of reserves in living biomass
amount of living biomass in floc amount of dead biomass in floc number of daughters from a floc volumespecific maximum assimilation rate volumespecific maintenance rate specific growth rate of cells specific growth rate of flocks time time at birth interdivision time volume at division volume of living biomass total volume of biomass volume of dead biomass scaled substrate concentration: substrate concentration at saturation coefficient of scaled functional response amount of biomass / volume floc structural biomass / volume dead biomass / volume minimum substrate concentration for support substrate needed per biomass formed substrate needed per reserve formed ModelThe purpose of the model is to simulate the dynamics of Flocs, and their effects on reactor dynamics. It is therefore first necessary to approximate the growth of these flocs. Floc GrowthNormal flocs are irregular in shape but for the purpose of this model they will be estimated as spherical. It will also be assumed that the growth of the floc is slow compared to diffusion such that steady states can be achieved. The volumespecfic biomass will additionally be held constant. There are two types of flocs: Small flocs were all microbes receive adequate nutrition via diffusion and large flocs where the core consists of dead biomass. If the floc grows exponentially as seen in the following equation. This represents small flocs where all cells receive adequate nutrition.
However if the floc has a living mantel and a core consisting of a dead biomass the volume changes as follows
DivisionNow that the growth rate of a floc has been modeled it is necessary to describe the life cycle of a floc. A floc will continue to grow until it reaches a critical volume and then the floc will rupture. This rupture is due to mechanical instability. The critical size is a function of the shear forces and floc porosity, age of floc. For this model a single variable accounts for these parameters denoting the maximum volume. It is assumed that the living mantle dissociates into daughter flocs with out changing thickness. The number of daughter flocs is a result of the total volume and the living mantle. The living mantle is dependent on the concentration of the substrate. A figure of this cycle can be seen below..
A plot showing the ratio of to as a function of critical size can be seen below. .Substrate Concentration and Reactor DynamicsNow that the growth and life of a floc are defined the model explores the steady states of the floc. The steady states will help determine the maximum thickness for the living mantle, and thus the amount a metabolic activity occurring in a floc. Cells require substrate that is transported via diffusion. The substrate concentration profile in a sphere as a result of diffusion is as follows.
It is important to assume that the concentration of the substrate is constant and the growth is slow compared the substrate profile.
This final equation uses the assumptions previously stated and the MarrPitt cellular growth model to simulate the dynamics of flocs in biodegradation. Cellular GrowthThus far, the growth of a floc has been defined as a scaled function of a known cellular growth rate: the MarrPitt model However, the MarrPitt is overly simple. In a more accurate model for growth rate of the flocs, the specific growth rate is found using the Dynamic Energy Budget (DEB) Theory (Kooijman,2000). According to DEB, cells can only grow while reserves are sufficient. Both the energy required of cellular maintenance and growth detract from the reserves. Also, a parameter known as the 'investment ratio' is introduced to account for the varying resources devoted to either growth or maintenance. When taking this into account, the growth rate of the living cells now becomes:
Using this new rate definition, the change in floc radius can be recomputed as follows:
Which is obviously more complex than the original time rate of the radius of the floc. This is due to the variation of parameters previously held constant. To allow easier computation, though, we can use the same simplifications as in the previous sections, and the above equation becomes: . The only difference being that is now the DEB specific growth rate as defined above. AnalysisTaking into account the effects of flocculation on reactor dynamics greatly improves the accuracy of predictions being made. This accuracy can deepen our understanding of the requirements and restrictions made on manufacturing and safety standards. The limitations placed on the growth of actively reacting biomass are based on observations made by many different laboratories. This slowed growth rate helps to improve the predictions made of reactors. Instead of assuming continued exponential growth of cells proportional to the nutrients available, the model smartly rectifies the growth of flocs as their size grows beyond sustainability. Other improvements over current models are made in the interaction between cells and the substrate being consumed. Older models, assuming cell suspension, have used various functional responses based on the concentration of substrate and other parameters. However, few accounted for the need for substrates to diffuse through layers of the flocs to reach internal cells. In this way, the rate of reaction is reduced significantly, and the new model accounts for this (fig. 5). The model does have some shortcomings, and variables that are ignored or held to be constant. Just like any model, steadystate assumptions lead to vastly simpler answers, but at a cost to accuracy. The shape of the flocs is assumed to be spherical; the irregularities ignored. This assumption affects all further aspects of the model, including growth rate, size at division and consumption rate. Another assumption made is that the volume at division is held constant, based on hydrodynamic conditions and other characteristics of the flocs. In reality, division volume will vary with respect to many characteristics of the reactor. Also, it will vary in space, as shear forces will be much greater in some locations. Overall, the assumptions made are not implausible. By considering the effects of flocculation, the model greatly improves upon its contemporaries. Brandt and Kooijman did not run any experiments to test their model, but used existing data sets to verify their model. Further StudyFurther research could be done on flocculation. The effects of the shape and size of the containing vessel should be studied, as hydrodynamic forces vary greatly in differing vehicles. Figuring out the metabolic activity of the kernel could also prove helpful. The paper assumed a steady state distribution of substrate such that the kernel had absolutely no metabolic activity. In reality, this may not be the truth, as the interior cells may realize their reduced supply of nutrients and slow growth, thus reducing the expenses of survival. This would allow the kernel to have a continued effect on chemical decomposition. Currently, it would prove very difficult, if not impossible, to explicitly measure some of the states of a floc, as they are often on the microscale. But, as nanotechnology improves, perhaps soon we will be able to answer some of the unknowns left in the understanding of reactor dynamics. External LinksComputational fluid dynamics with flocculation 