Objective: Building a Bioreactor with intracellular reactions
|
|
- Robert Clark
- 6 years ago
- Views:
Transcription
1 BioCell Model Objective: Building a Bioreactor with intracellular reactions This example shows a more advanced bioreactor model where both intracellular and extracellular concentrations are considered. For an introduction to the basics of bioreactor modeling in REX, you may find this simpler BioReactor example more useful. We model the production of Xylitol from Xylose by the yeast Candida Mogii in a batch reactor. This model is fully based on the work of Tochampa et al. [Ref_1]. You may download the zipped rex file to view the model. Features Illustrated Building a Biocell model. Defining Cell or Cells and intracellular Compounds. Modelling compounds and reactions inside the cells. Fermentation Model Xylose available in the liquid solution is taken up by Candida Mogii yeast, where it is reduced to Xylitol. Part of the Xylitol is excreted by the cell, and the rest is used to generate cell mass. Glucose in the solution medium helps to increase the xylitol yield, as the glucose pathway to cell mass is a more direct way than xylitol to generate energy and biomass. The biomass growth rate is given by: [1] Where [BioMass] is the biomass concentration (g/l) and μ is the specific growth rate (1/t). The specific growth rate is modelled considering two substrates, Glucose in the medium and Xylitol inside the cell, as follows: [2] There, [Glucose] is the glucose concentration (g/l) and [Xylitol] in is the xylitol concentration within the cell.
2 Parameter μ Glc max is the maximum specific rate on glucose alone and μ Xit max is the maximum specific rate on xylitol; K s,glc and K s,xit are the Monod saturation constants for Glucose and Xylitol respectively and K r is the repression constant by glucose. The rates for glucose and xylose uptake are related by the biomass concentration as: [3] [4] where the specific uptakes of glucose (q Glc ) and xylose (q Xyl ) are expressed as follows: [5] [6] There q Glc max and q Xyl max are the maximum uptake rates of glucose and xylose respectively; k i,glc and k i,xyl are parameters related to the inhibition for the uptake of one sugar with respect to the presence of the other. Saturation constant K s,glc is assumed to be same as in [2]. From the two last equations, you can see that the presence of glucose decreases the uptake of xylose, and the presence of xylose decreases the uptake of glucose, because the transport of those species through the cell membrane is facilitated by the same diffusion system. It is assumed that all xylose taken up from the medium is converted inside the cell to xylitol. Thus formation of intracellular xylitol is same as the uptake of xylose, adjusted by their ratio of molecular weight M xylitol / M xylose = ( g/mol) / ( g.mol) = Given the specific uptake of Xylose from [6], the specific formation xylitol inside the cell is given by:
3 [7] Some of the intracellular xylitol is used to create Biomass. The biomass yield on xylitol Y BM,Xit thus tied the specific rates for biomass formation from xylitol (last term of Eq [2]) with respect to Xylitol consumption for cell formation: [8] Intracellular Xylitol is also transported outside across the cell membrane. The proposed formulation from Tochampa et al. [Ref_1] is estimated using a modification of Fick's law. The specific mass flux is proportional given by: [9] Where D is the diffusion coefficient of xylitol in the cell membrane, K is the partition coefficient, d is membrane thickness and a is the specific surface area of the cell. Thus, the Xylitol mass balance proposed by [Ref_1] is actually twofold: a mass balance for Xylitol inside the cells and another for Xylitol in the medium (that is, outside the cells): [10] [11] In the mass balance of xylitol inside the cell, ρ BM is the mass density of the cell. The last term of [10] represents the dilution effect of the species inside the cells due to the increase of biomass. A derivation of this is shown later in this document.
4 The above equations represent the model presented by Tochampa et al. Please note that all the equations here deal with mass concentrations. REX works directly with mass balances, but we will replicate the same model above. The equivalence of the REX model to the equations above is explained in more detail here. Basics of BioCell Modelling in REX REX allows to model reactions that could take place either inside or outside the cells. Also, more than one cell type can be considered. To enable this special model, you must go to the Chemistry node and right click on Enable Special Models Bio Cell Model: Once the Bio Cell model is enabled, the Cells node is shown in the Chemistry tree. In the first grid of that node, all compounds are listed and you must select which ones represent the cells. Suppose we mark compound C1 as a Cell: Then, in the second grid we can assign which compounds will be modelled inside the cell C1. You can have more than one cell defined, but a given compound can not be assigned to more than one cell.
5 Once the intracellular and extracellular compounds are assigned, the concentrations are calculated accordingly. Considering Mass Density as the unit for Concentrations in Units Configurations node, we describe the concentration calculation for the different types of species below: Cells: The concentration for any compound selected as Cell is calculated as the usual way M i / V rph, where M i is the mass of cells in the reactor and V rph is the volume of the liquid reaction phase. Compounds assigned to Cells: Those compounds are assumed to be located inside a given cell, thus their concentration are referred to the cell volume: M i / V c, where M i is the compound mass and V c is the volume of the cell that he compound belongs to. The Cell volume is calculated through its density value that must be provided in the Chemistry Compounds Properties node. Compounds that are not Cells nor assigned to Cells: Those compounds are in the liquid solution outside the cells, thus its concentration is calculated with respect to the liquid phase volume, excluding the volume occupied by the cells themselves: M i / [V rph - V c ]. There, M i is the mass of compound i in the liquid phase, V rph is the volume of the liquid phase and V c is the volume of the cell species. The mass balance equations used in REX are explained below. The equivalence between the REX model and the model in Tochampa et al is also described. In general, the mass balance for compound i for a pure batch reactor is given by: [12] Here, M i is the total mass of compound i in the reactor, V rph is the volume of the reaction phase, R r is the rate for reaction r, and v i,r is the stoichiometry coefficient (mass basis) for compound i in reaction r. For Biocell models, the REX equations are formulated differently to account for both the cells and the extracellular medium. For a batch reactor without inflows or outflows, it is given by: [13] Here we implicitly discriminate among three types of reactions, modelled in the corresponding terms on the Right Hand Side (RHS) of the above equation:
6 Biomass Reactions, indicated as by the bm subindex. These reactions are related to cell growth, and REX identifies them automatically: any reaction with a cell compound in its stoichiometry is assigned to this reaction type. The only exception is where a compound inside the cell it is also present in the reaction stoichiometry: for that case, the reaction is not considered in this reaction type but in the one described next. In the batch mass balance above, they are taken into account in the first RHS term, where the reaction rate is multiplied by the liquid phase volume V rph and by the stoichiometric coefficient v i,bm. Reactions taking place inside the Cell: Their rates are indicated by R ic,c in the second term of the above mass balance equation, where the subindex ic and c indicate that the reaction takes place inside the cell c. If a reaction stoichiometry has one or more compounds belonging inside a cell (identified in the Cells node), the REX internally assigns the reaction to this type. In the batch mass balance for compound i, the reaction rate is multiplied by the related stoichiometric coefficient of compound i and by the reaction phase volume V rph. An important issue is that in REX we use R ic,c that has units of [Mass] / [Time] / [Volume of liquid phase]. However, it is common to refer to specific rate r ic,c with units [Mass] / [Time] / [Mass of cells], as indicated in [1], [3] or [4]. The relationship between them is given by: where C c is the mass concentration of cell c in the reactor phase. Thus, when implementing a literature model based on specific rates r ic,c, we must add an order of one with respect to the cell to get R ic,c which is the rate type required in REX. Reactions taking place outside the Cells: Shown in the last term of the RHS of the mass balance above. There, R oc represents the reaction rate, that is multiplied by the stoichiometric coefficient and the liquid phase volume (V rph minus the volume V c taken by the cells c). REX identifies this reaction type when its stoichiometry does not have a cell or a compound inside the cell. In this other example, we showed the equivalence between the model based on concentration derivatives and the REX internal model where derivatives for total mass of compounds is used. We extend that equivalence for the models above where intracellular compounds are considered. For a species i marked as belonging inside a cell, all reactions that contain this species are considered to take place inside the cell. Thus [13] simplifies to: [14] [15] where V rph is the reaction phase volume.
7 Simplification: Even if there are multiple cells defined in the project, compound i can only be assigned to one of them. From here on, we refer only to the cell to which the compound belongs and discard any reference to the other cells, so the index c on the summation can be dropped. We now express the mass of compound i as a function of its concentration inside the cell C i and cell volume V BM. The latter is replaced by the ratio of M BM / ρ BM (Cell mass / Density): [16] Applying the product rule for the differential: [17] Solving for dc i / dt : [18] In the first term of the RHS, R ic,c is replaced by the specific rate r ic,c, according to [14]. In the second term, we replace M BM in the derivative by the product of the reactor volume V rph and the biomass concentration C c. In addition, we assume constant reaction phase volume and take it out of the derivative: [19] In the first RHS term, C c V rph = M BM, so it get eliminated by the M BM in the denominator. In the second term, we replace V rph / M BM, by C c -1 :
8 [20] Replacing dc c / dt by [1], considering that [Biomass] = C c, we have: [21] This last formula is the general expression, from which [10] is the particular equation for Xylitol_in in this example. The last term on the RHS takes into account the dilution of the compounds inside the biomass when it grows and is implicit in the REX model. The advantage of using Mass derivatives over Concentrations derivatives is that we need not assume constant reactor volume. Xylitol Model: Implementation in REX We now describe the model implementation in REX. You may download and import the zipped BioCell.rex file to follow this section. You will see that the Biocell Model has been enabled in the Chemistry node. Units have been selected in accordance with the work of [Ref_1] and Mass Density is chosen as the unit for Concentration. We enter the compounds. As we need to distinguish between the Xylitol inside and outside the cell, we load Xylitol_in and Xylitol_out as two compounds to track them individually: In the Cells node, BioMass is selected to represent the cell; Xylitol_in is the only compound we need to assign inside the cell:
9 Now we define the reactions in the system. For BioMass Growth, we need two reactions as the cell growth depends on Glucose and Xylitol reactants, according to equations [1] and [2]. As seen from the model for Biomass formation from Xylitol_in (second RHS term of equation [2]) and the equation for xylitol_in consumption (equation [8]), we can tie these equations by defining this reaction: [22] The value of the Y BM,Xit is unknown so it will estimated from the experimental data as in this other Example. In the Parameters node, we set the orders of that reaction which according to [1] and [2] is one with respect to both BioMass and Xylitol_in.
10 We must now define the LHHW sites that account for the inhibitions indicated in [2]. You can inspect them in the Kinetics Parameters node, where they are named InhBioMass-Xit and RepBioMass- Xit, and they are assigned to the BiomassGrowth-Xit reaction in node Kinetics Kinetics Sites. For glucose, the proposed model for conversion to biomass is not as direct. There is no simple proportionality constant between the equations for biomass growth from glucose ([1] and first RHS term of [2]) and the related model for glucose uptake (equations [3] and [5]). So a single reaction cannot model glucose to Biomass reaction as we could for the Xylitol to BioMass path. To address this, we define two reactions in REX: one where we enter BioMass as the only product and without any reactant, and another reaction with Glucose as reactant and no product. We will then use parameter relationships to tie the common parameters in these two reactions. The reaction definitions are shown below:
11 In the Kinetics Parameters node, we enter the orders for these reactions, that are one with respect to Glucose and BioMass in both. Then we enter the inhibition sites for each reaction, in accordance with [2] and [5]. Parameter K s,glc in both reactions is same, so we need to enforce a parameter relationship. As we can enforce parameter relationship for Preexponentials and Energies but not for the site constant in Parameter Relationships node, we rearrange equation [5] so that:
12 [23] Similar procedure is carried out for the inhibition of the first term in equation [2]. Thus, we can constrain the preexponentials of the glucose term of both sites in the Parameter Relationships node as shown below: As for the Xylose conversion to Xylitol inside the cell, based on [7] we can enter this reaction:
13 From equations [4] and [6], we enter order one with respect to Xylose and BioMass: An LHHW Site is defined in Kinetics Parameters node to account for the inhibition in [6] and assigned to this reaction in the Kinetics Kinetics Sites node. Finally, we have to model the diffusion of Xylitol through the cell membrane in accordance with [9]. We need to use an auxiliary variable to account for the concentration difference that is the driving force for diffusion. That is done by defining a PseudoCompound Xit_In-Out, as shown below: The diffusion equation [9] is now modeled as a reaction Xit_Diff. The orders of that reaction are one with respect to the above mentioned pseudocompound and with respect to the biomass. This reaction has no inhibition factor:
14 Parameter Estimation: Results The experimental data from [Ref_1] consists of concentration measurements for Glucose, Biomass, Xylitol and Xylose. The experiments combine sets with high and low initial Biomass concentration and 0, 7 and 14% of initial Glucose content. Measurements of Xylitol correspond to the concentration in the medium, so it is assigned to compound Xylitol_out. The experimental concentrations are entered in the phase tab of Measurements Sets node: In the Total Mass tab, we need to specify the initial amount of each compound and we calculate that from the provided initial concentrations. Thus, the initial amount of BioMass is set to initial
15 concentration multiplied by the volume. For the other species that are in the medium outside the cells, the calculation is done multiplying the initial concentration and the phase volume excluding the volume taken up by the Biomass at time=0. We assume that there is zero amount of xylitol inside the cell at time=0. Thus, the Total Mass tab for a given set is entered as below: The parameters that are estimated have open bounds in Estimation Parameters node. The energies are not estimated (fixed to zero) because all experimental data is at same temperature. Not all LHHW Sites preexponentials are estimated in accordance with the model. Stoichiometric coefficients are only estimated for the xylitol to biomass reaction. The results of the estimation can be seen in the Results tree. In particular, we show below the concentration predictions for a set with low biomass and 14% initial glucose:
16 You can see that Glucose is depleted early where it is converted to biomass. Interestingly enough, Xylose is also consumed but slowly during this initial stage because of glucose inhibition according to [6]. Xylose concentration actually increases a bit at the beginning due to low consumption while the volume of the medium is decreased by the formation of biomass. As Glucose concentration reduces, xylose consumption accelerates. Xylitol is formed inside the cell and its concentration outside evolves as shown in the chart. However, when the Xylose concentrations become low, so does Xylitol formation. Thus, a maximum Xylitol concentration is reached in the medium. Xylitol diffuses back to the cell and xylitol consumption to yield biomass becomes the predominant reaction. In essence, the product Xylitol is lost to biomass formation. The chart below compares Xylitol concentration outside and inside the cell: In the beginning, Xylitol inside the cell (pink line) increases fast as its consumption to produce biomass is inhibited by Glucose according to [8]. This increases the diffusion of Xylitol across the cell membrane and extracellular Xylitol concentration increases (yellow line) according to [9]. When Glucose is depleted, its inhibition on both the Xylose conversion to Xylitol ([6] and [7]) and the Xylitol conversion to BioMass [2] is gone. Rates of both Xylitol formation and consumption tends to accelerate as glucose inhibition is not present. However, as xylose reduces, conversion to xylitol becomes slower. External xylitol diffuses back to the cell and gets converted to biomass, thus explaining the xylitol maximum. Even though the REX model above and the model of [Ref_1] do a reasonable job of fitting the experimental data, the estimated kinetic parameters in these two models turned out to have
17 substantial differences. The reasons for this could be manifold and a few of them are highlighted below: Sets used for estimation: All experimental data was used for estimation in REX, whereas only the first three sets were used in [Ref_1] Weighting factors: The hybrid weighting strategy used in REX is quite different from the method used in [ Ref_1] Correlated Parameters and non-uniqueness: Given that both models are able to fit experimental data, parameter correlation could be a culprit. In LHHW systems, there is often strong correlation between the pre-exponentials for the numerator and denominator terms. Additional experiments or reformulation may be needed to get more unique estimates. Estimation solver: REX uses a gradient based (Generalized reduced gradient method) for the nonlinear estimation problem which is quite different from the genetic algorithm used in [Ref_1]. This example is a good reminder of the care needed in parameter estimation. Overparametrization can result in multiple parameter estimates giving an equally decent fit to the experimental data. Care should be taken in experiment design and the choice of model to get reliable estimates. References 1. Tochampa, W., Sitisansaneeyakul, S., Vanichsriratana, W., Srinophakun, P., Bakker, H., Chisti, Y., A Model of Xylitol Production by the Yeast Candida Moguii, Bioprocess Biosyst Eng, 2005 (28),
Objective: To Evaluate Predictions of Alternative Models
www.optience.com Methanol Synthesis Objective: To Evaluate Predictions of Alternative Models In this example, we propose and evaluate Mass Action and Langmuir Hinshelwood (LHHW) models for methanol synthesis
More informationObjective: To Estimate Kinetic Parameters for Toluene Reaction with Benzyl Alcohol
www.optience.com Toluene Benzylation Objective: To Estimate Kinetic Parameters for Toluene Reaction with Benzyl Alcohol In this example, an LHHW kinetic model for toluene reaction over solid catalyst in
More informationPolymerization Modeling
www.optience.com Polymerization Modeling Objective: Modeling Condensation Polymerization by using Functional Groups In this example, we develop a kinetic model for condensation polymerization that tracks
More informationModel for Propane Dehydrogenation
www.optience.com Model for Propane Dehydrogenation Objective: Develop a kinetic model for Propane Dehydrogenation In this example, we propose a simplified model for the catalytic dehydrogenation of propane
More informationPropylene Hydroformylation
www.optience.com Propylene Hydroformylation Objective: Building Fedbatch reactors with Pressure Control In this example, we simulate a multiphase batch reactor for Propylene Hydroformylation and also explain
More informationOxidative Coupling of Methane
www.optience.com Oxidative Coupling of Methane Objective: Parameter Estimation for Methane to Ethylene reaction In this example, we build a simplified reaction network for the Oxidative Coupling of Methane
More informationAGITATION AND AERATION
AGITATION AND AERATION Although in many aerobic cultures, gas sparging provides the method for both mixing and aeration - it is important that these two aspects of fermenter design be considered separately.
More informationADAPTIVE ALGORITHMS FOR ESTIMATION OF MULTIPLE BIOMASS GROWTH RATES AND BIOMASS CONCENTRATION IN A CLASS OF BIOPROCESSES
ADAPTIVE ALGORITHMS FOR ESTIMATION OF MULTIPLE BIOMASS GROWTH RATES AND BIOMASS ONENTRATION IN A LASS OF BIOPROESSES V. Lubenova, E.. Ferreira Bulgarian Academy of Sciences, Institute of ontrol and System
More informationChemical Reaction Engineering. Multiple Reactions. Dr.-Eng. Zayed Al-Hamamre
Chemical Reaction Engineering Multiple Reactions Dr.-Eng. Zayed Al-Hamamre 1 Content Types of Reactions Selectivity Reaction Yield Parallel Reactions Series Reactions Net Rates of Reaction Complex Reactions
More informationConvective Heat and Mass Transfer Prof. A.W. Date Department of Mechanical Engineering Indian Institute of Technology, Bombay
Convective Heat and Mass Transfer Prof. A.W. Date Department of Mechanical Engineering Indian Institute of Technology, Bombay Module No. # 01 Lecture No. # 32 Stefan Flow Model We are now familiar with
More informationDIFFERENTIAL EQUATION ANALYSIS IN BIOMEDICAL SCIENCE AND ENGINEERING
DIFFERENTIAL EQUATION ANALYSIS IN BIOMEDICAL SCIENCE AND ENGINEERING DIFFERENTIAL EQUATION ANALYSIS IN BIOMEDICAL SCIENCE AND ENGINEERING ORDINARY DIFFERENTIAL EQUATION APPLICATIONS WITH R William E.
More informationFigure 1: Ca2+ wave in a Xenopus oocyte following fertilization. Time goes from top left to bottom right. From Fall et al., 2002.
1 Traveling Fronts Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 2 When mature Xenopus oocytes (frog
More informationSurface Chemistry Tutorial
Surface Chemistry Tutorial Introduction Surface chemistry is often the most important and most overlooked aspect of reacting flow modeling. Surface rate expressions can be hard to find or not even exist
More information2. At quasi-steady state or equilibrium, the net in-flux of the carrier-substrate complex CS is balanced by the net out-flux of the free carrier C.
Facilitated Transport Instructor: Nam un Wang facilitmcd Process escription In facilitated transport, a carrier molecule C binds to the substrate to form a carrier-substrate complex C at the outer side
More informationCellular Metabolic Models
Cellular Metabolic Models. Cellular metabolism. Modeling cellular metabolism. Flux balance model of yeast glycolysis 4. Kinetic model of yeast glycolysis Cellular Metabolic Models Cellular Metabolism Basic
More informationBioinformatics: Network Analysis
Bioinformatics: Network Analysis Reaction Kinetics COMP 572 (BIOS 572 / BIOE 564) - Fall 2013 Luay Nakhleh, Rice University 1 Reaction kinetics is the study of how fast chemical reactions take place, what
More informationAdvanced Chemical Reaction Engineering Prof. H. S. Shankar Department of Chemical Engineering IIT Bombay. Lecture - 03 Design Equations-1
(Refer Slide Time: 00:19) Advanced Chemical Reaction Engineering Prof. H. S. Shankar Department of Chemical Engineering IIT Bombay Lecture - 03 Design Equations-1 We are looking at advanced reaction engineering;
More informationSystems II. Metabolic Cycles
Systems II Metabolic Cycles Overview Metabolism is central to cellular functioning Maintenance of mass and energy requirements for the cell Breakdown of toxic compounds Consists of a number of major pathways
More informationModeling Microbial Populations in the Chemostat
Modeling Microbial Populations in the Chemostat Hal Smith A R I Z O N A S T A T E U N I V E R S I T Y H.L. Smith (ASU) Modeling Microbial Populations in the Chemostat MBI, June 3, 204 / 34 Outline Why
More informationReactors. Reaction Classifications
Reactors Reactions are usually the heart of the chemical processes in which relatively cheap raw materials are converted to more economically favorable products. In other cases, reactions play essential
More informationSolutions to Problem Assignment 1 (Kinetics) Oxtoby Problems -- see Solutions Manual. Solutions to Lecture Problems I # 1-6
Solutions to Problem Assignment 1 (Kinetics) Oxtoby Problems -- see Solutions Manual Solutions to Lecture Problems I # 1-6 I-1. This problem illustrates how the differential rate law (DRL) can be simplified
More informationStationary phase. Time
An introduction to modeling of bioreactors Bengt Carlsson Dept of Systems and Control Information Technology Uppsala University August 19, 2002 Abstract This material is made for the course Wastewater
More informationChemical Reaction Engineering - Part 12 - multiple reactions Richard K. Herz,
Chemical Reaction Engineering - Part 12 - multiple reactions Richard K. Herz, rherz@ucsd.edu, www.reactorlab.net Multiple reactions are usually present So far we have considered reactors in which only
More informationComments Transient Energy Balances
Comments Transient Energy Balances General form of the stuff balance equation Rate of Rate Rate Rate of Rate of Accumulation In Out Generation Consumption F1 Q1 F6 F2 Q3 F3 Q2 F4 F5 2 Word form of the
More informationwhere a + b = 2 (this is the general case) These all come from the fact that this is an overall second order reaction.
Chapter 7 Problems Page of 6 //007 7. Hydrolysis of ethyl acetate is as follows: EtAc + OH - Ac - + EtOH. At 5 ºC, the disappearance of OH - is used to determine the extent of the reaction, leading to
More informationPhysics 142 Steady Currents Page 1. Steady Currents
Physics 142 Steady Currents Page 1 Steady Currents If at first you don t succeed, try, try again. Then quit. No sense being a damn fool about it. W.C. Fields Electric current: the slow average drift of
More informationTheoretical Models for Chemical Kinetics
Theoretical Models for Chemical Kinetics Thus far we have calculated rate laws, rate constants, reaction orders, etc. based on observations of macroscopic properties, but what is happening at the molecular
More information1/r plots: a brief reminder
L10-1 1/r plots: a brief reminder 1/r X target X L10-2 1/r plots: a brief reminder 1/r X target X L10-3 1/r plots: a brief reminder 1/r X target X Special Case: utocatalytic Reactions Let s assume a reaction
More informationStoichiometry. Please take out your notebooks
Stoichiometry Please take out your notebooks Stoichiometry stochio = Greek for element metry = measurement Stoichiometry is about measuring the amounts of elements and compounds involved in a reaction.
More informationWritten Exam 15 December Course name: Introduction to Systems Biology Course no
Technical University of Denmark Written Exam 15 December 2008 Course name: Introduction to Systems Biology Course no. 27041 Aids allowed: Open book exam Provide your answers and calculations on separate
More informationPlug flow Steady-state flow. Mixed flow
1 IDEAL REACTOR TYPES Batch Plug flow Steady-state flow Mixed flow Ideal Batch Reactor It has neither inflow nor outflow of reactants or products when the reaction is being carried out. Uniform composition
More informationPhoto-Biological Reactor for Organic Waste Consumption and Hydrogen Production
Photo-Biological Reactor for Organic Waste Consumption and Hydrogen Production L. F. de Souza 1 1 Universidade Federal do Paraná, Curitiba, Paraná, Brazil *Corresponding Author: Rua Cerro Azul 1498, Colombo
More informationA Gas Uniformly fills any container. Easily compressed. Mixes completely with any other gas. Exerts pressure on its surroundings.
Chapter 5 Gases Chapter 5 A Gas Uniformly fills any container. Easily compressed. Mixes completely with any other gas. Exerts pressure on its surroundings. Copyright Cengage Learning. All rights reserved
More informationBifurcation Analysis of Continuous Biochemical Reactor Models
Biotechnol. Prog. 2001, 17, 647 660 647 Bifurcation Analysis of Continuous Biochemical Reactor Models Yongchun Zhang and Michael A. Henson* Department of Chemical Engineering, Louisiana State University,
More informationDevelopment of Dynamic Models. Chapter 2. Illustrative Example: A Blending Process
Development of Dynamic Models Illustrative Example: A Blending Process An unsteady-state mass balance for the blending system: rate of accumulation rate of rate of = of mass in the tank mass in mass out
More informationChemical Processes. Transport Processes in the Environment. Transport Processes. Advection. Advection. t 1 t 2
hemical Processes by Stefan Trapp (taken from course 4) Transport Processes in the Environment a definition of our nomenklatura Transport Processes Only 3 types of transport processes exist ) by advection
More informationLecture 6 Asymptotic Structure for Four-Step Premixed Stoichiometric Methane Flames
Lecture 6 Asymptotic Structure for Four-Step Premixed Stoichiometric Methane Flames 6.-1 Previous lecture: Asymptotic description of premixed flames based on an assumed one-step reaction. basic understanding
More informationInflow Qin, Sin. S, X Outflow Qout, S, X. Volume V
UPPSALA UNIVERSITET AVDELNINGEN FÖR SYSTEMTEKNIK BC,PSA 9809, Last rev August 17, 2000 SIMULATION OF SIMPLE BIOREACTORS Computer laboratory work in Wastewater Treatment W4 1. Microbial growth in a "Monode"
More informationENZYME SCIENCE AND ENGINEERING PROF. SUBHASH CHAND DEPARTMENT OF BIOCHEMICAL ENGINEERING AND BIOTECHNOLOGY IIT DELHI
ENZYME SCIENCE AND ENGINEERING PROF. SUBHASH CHAND DEPARTMENT OF BIOCHEMICAL ENGINEERING AND BIOTECHNOLOGY IIT DELHI LECTURE 23 STEADY STATE ANALYSIS OF MASS TRANSFER & BIOCHEMICAL REACTION IN IME REACTORS
More informationCarbon labeling for Metabolic Flux Analysis. Nicholas Wayne Henderson Computational and Applied Mathematics. Abstract
Carbon labeling for Metabolic Flux Analysis Nicholas Wayne Henderson Computational and Applied Mathematics Abstract Metabolic engineering is rapidly growing as a discipline. Applications in pharmaceuticals,
More informationMITOCW enzyme_kinetics
MITOCW enzyme_kinetics In beer and wine production, enzymes in yeast aid the conversion of sugar into ethanol. Enzymes are used in cheese-making to degrade proteins in milk, changing their solubility,
More informationSteady-State Molecular Diffusion
Steady-State Molecular Diffusion This part is an application to the general differential equation of mass transfer. The objective is to solve the differential equation of mass transfer under steady state
More informationA. One-Substrate Reactions (1) Kinetic concepts
A. One-Substrate Reactions (1) Kinetic concepts (2) Kinetic analysis (a) Briggs-Haldane steady-state treatment (b) Michaelis constant (K m ) (c) Specificity constant (3) Graphical analysis (4) Practical
More informationTHE CHEMICAL REACTION EQUATION AND STOICHIOMETRY
9.1 Stoichiometry Stoichiometry provides a quantitative means of relating the amount of products produced by chemical reaction(s) to the amount of reactants. You should take the following steps in solving
More informationThermodynamic and Stochiometric Principles in Materials Balance
Thermodynamic and Stochiometric Principles in Materials Balance Typical metallurgical engineering problems based on materials and energy balance NiO is reduced in an open atmosphere furnace by excess carbon
More informationChapter 5 Ground Rules of Metabolism Sections 1-5
Chapter 5 Ground Rules of Metabolism Sections 1-5 5.1 A Toast to Alcohol Dehydrogenase In the liver, the enzyme alcohol dehydrogenase breaks down toxic ethanol to acetaldehyde, an organic molecule even
More informationReactions and Reactors
Reactions and Reactors ChE 400 - Reactive Process Engineering If we want to run a chemical process in order to convert some reactants (the reactor feed) to some product (the reactor effluent), we have
More informationProblems Points (Max.) Points Received
Chemical Engineering 142 Chemical Kinetics and Reaction Engineering Midterm 1 Tuesday, October 8, 2013 8:10 am-9:30 am The exam is 100 points total. Please read through the questions very carefully before
More informationMass Transfer in a Stirred Batch Reactor
Mass Transfer in a Stirred Batch Reactor In many processes, efficient reactor usage goes hand in hand with efficient mixing. The ability to accurately examine the effects of impeller placement, speed,
More informationSimulated niche partitioning by bacteria
Simulated niche partitioning by bacteria Steven S. Andrews and Adam P. Arkin Physical Biosciences Division Lawrence Berkeley National Laboratory 1 Cyclotron Road Berkeley, CA 94720 ssandrews@lbl.gov 1.
More informationOptimal Feeding Strategy for Bioreactors with Biomass Death
Optimal Feeding Strategy for Bioreactors with Biomass Death L. Bodizs, B. Srinivasan, D. Bonvin Laboratoire d Automatique, Ecole Polytechnique Féderale de Lausanne, CH-1015, Lausanne, Switzerland Abstract
More informationCHAPTER 1: ENZYME KINETICS AND APPLICATIONS
CHAPTER 1: ENZYME KINETICS AND APPLICATIONS EM 1 2012/13 ERT 317 BIOCHEMICAL ENGINEERING Course details Credit hours/units : 4 Contact hours : 3 hr (L), 3 hr (P) and 1 hr (T) per week Evaluations Final
More informationThe Rate Expression. The rate, velocity, or speed of a reaction
The Rate Expression The rate, velocity, or speed of a reaction Reaction rate is the change in the concentration of a reactant or a product with time. A B rate = - da rate = db da = decrease in concentration
More informationMATLAB Ordinary Differential Equation (ODE) solver for a simple example 1. Introduction
MATLAB Ordinary Differential Equation (ODE) solver for a simple example 1. Introduction Differential equations are a convenient way to express mathematically a change of a dependent variable (e.g. concentration
More informationImplementation of the King-Altman method in DynaFit
Implementation of the King-Altman method in DynaFit BioKin Technical Note TN-2015-03 Petr Kuzmič BioKin Ltd., Watertown, Massachusetts, USA http: // www. biokin. com Abstract The DynaFit software package
More informationUniversity of Groningen. Levulinic acid from lignocellulosic biomass Girisuta, Buana
University of Groningen Levulinic acid from lignocellulosic biomass Girisuta, Buana IMPRTANT NTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it. Please
More informationEngineering Model Reduction and Entropy-based Lyapunov Functions in Chemical Reaction Kinetics
Entropy 2010, 12, 772-797; doi:10.3390/e12040772 Article OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Engineering Model Reduction and Entropy-based Lyapunov Functions in Chemical Reaction
More informationBioethanol production sustainability: Outlook for improvement using computer-aided techniques
17 th European Symposium on Computer Aided Process Engineering ESCAPE17 V. Plesu and P.S. Agachi (Editors) 007 Elsevier B.V. All rights reserved. 1 Bioethanol production sustainability: Outlook for improvement
More informationII. The Relationship Between Photosynthesis and Respiration
Name Period Date Cellular Respiration Notes Outline I The Importance of Food A Food provides living things with the: B Food serves as a source of: C Food serves as a source of: II The Relationship Between
More informationMATHCAD SOLUTIONS TO THE CHEMICAL ENGINEERING PROBLEM SET
MATHCA SOLUTIONS TO THE CHEMICAL ENGINEERING PROBLEM SET Mathematical Software Session John J. Hwalek, epartment of Chemical Engineering University of Maine, Orono, Me 4469-5737 (hwalek@maine.maine.edu)
More informationFreezing Around a Pipe with Flowing Water
1 Introduction Freezing Around a Pipe with Flowing Water Groundwater flow can have a significant effect on ground freezing because heat flow via convection is often more effective at moving heat than conduction
More informationIntroduction to cardiac electrophysiology 1. Dr. Tóth András 2018
Introduction to cardiac electrophysiology 1. Dr. Tóth ndrás 2018 Topics Transmembran transport Donnan equilibrium Resting potential 1 Transmembran transport Major types of transmembran transport J: net
More information2. Review on Material Balances
2. Review on Material Balances Objectives After completing this chapter, students should be able to recall the principle of the conservation of mass recall the principle of the stoichiometry of chemical
More informationDr. Trent L. Silbaugh, Instructor Chemical Reaction Engineering Final Exam Study Guide
Chapter 1 Mole balances: Know the definitions of the rate of reaction, rate of disappearance and rate of appearance Know what a rate law is Be able to write a general mole balance and know what each term
More informationCollision Rates, Mean Free Path, and Diffusion. Chemistry B
Collision Rates, Mean Free Path, and Diffusion Chemistry 8-23B David Ronis McGill University In the previous section, we found the parameter b by computing the average force exerted on the walls of the
More informationRate laws, Reaction Orders. Reaction Order Molecularity. Determining Reaction Order
Rate laws, Reaction Orders The rate or velocity of a chemical reaction is loss of reactant or appearance of product in concentration units, per unit time d[p] = d[s] The rate law for a reaction is of the
More informationSlide 1 / 6. Free Response
Slide 1 / 6 Free Response Slide 2 / 6 1 To maximize efficiency, the rate of light harvesting at PSII has to equal the rates of electron transfer from water to the PSI reaction center. If electrons are
More informationStability Analysis of Nonlinear Hybrid System in Microbial Fed-batch Culture
The Fourth International Conference on Computational Systems Biology (ISB010) Suzhou, China, September 9 11, 010 Copyright 010 ORSC & APORC, pp. 363 37 Stability Analysis of Nonlinear Hybrid System in
More informationIntroduction to the Problem-Solving Methodology (PSM)
Problem Name: Problem Description: Date: P01 Simple Reactor Introduction to the Problem-Solving Methodology (PSM) Your Name: Problem Session Objectives To learn the five stages of the problem solving methodology.
More informationThermodynamic principles governing metabolic operation : inference, analysis, and prediction Niebel, Bastian
University of Groningen Thermodynamic principles governing metabolic operation : inference, analysis, and prediction Niebel, Bastian IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's
More informationPractical exercises INSA course: modelling integrated networks of metabolism and gene expression
Practical exercises INSA course: modelling integrated networks of metabolism and gene expression Hidde de Jong October 7, 2016 1 Carbon catabolite repression in bacteria All free-living bacteria have to
More informationPart I.
Part I bblee@unimp . Introduction to Mass Transfer and Diffusion 2. Molecular Diffusion in Gasses 3. Molecular Diffusion in Liquids Part I 4. Molecular Diffusion in Biological Solutions and Gels 5. Molecular
More informationFAQ about SIMATIC BATCH
FAQ about SIMATIC BATCH How do I Configure Scaled Parameters? FAQ Table of Contents Table of Contents... 2 1 General Information... 3 1.1 Task/motivation... 3 1.2 Definitions... 3 2 Configuring Scaled
More information/639 Final Solutions, Part a) Equating the electrochemical potentials of H + and X on outside and inside: = RT ln H in
580.439/639 Final Solutions, 2014 Question 1 Part a) Equating the electrochemical potentials of H + and X on outside and inside: RT ln H out + zf 0 + RT ln X out = RT ln H in F 60 + RT ln X in 60 mv =
More informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
More informationAE 205 Materials and Energy Balances Asst. Prof. Dr. Tippabust Eksangsri. Chapter 4 Stoichiometry and MB with Reactions
AE 205 Materials and Energy Balances Asst. Prof. Dr. Tippabust Eksangsri Chapter 4 Stoichiometry and MB with Reactions Stoichiometry Stoichiometry provides a quantitative means of relating the amount of
More informationSPA for quantitative analysis: Lecture 6 Modelling Biological Processes
1/ 223 SPA for quantitative analysis: Lecture 6 Modelling Biological Processes Jane Hillston LFCS, School of Informatics The University of Edinburgh Scotland 7th March 2013 Outline 2/ 223 1 Introduction
More informationPre-AP Biology Energy Unit Study Guide Part I
Pre-AP Biology Energy Unit Study Guide Part I The Law of conservation of matter/mass : matter can not be created or destroyed However, matter may be rearranged in space In chemical reactions, the mass
More informationPhysical Pharmacy. Diffusion
Physical Pharmacy Diffusion Diffusion Diffusion is defined as a process of mass transfer of individual molecules of a substance brought about by random molecular motion and associated with a driving force
More informationOutline. Definition and mechanism Theory of diffusion Molecular diffusion in gases Molecular diffusion in liquid Mass transfer
Diffusion 051333 Unit operation in gro-industry III Department of Biotechnology, Faculty of gro-industry Kasetsart University Lecturer: Kittipong Rattanaporn 1 Outline Definition and mechanism Theory of
More informationGrowth models for cells in the chemostat
Growth models for cells in the chemostat V. Lemesle, J-L. Gouzé COMORE Project, INRIA Sophia Antipolis BP93, 06902 Sophia Antipolis, FRANCE Valerie.Lemesle, Jean-Luc.Gouze@sophia.inria.fr Abstract A chemostat
More informationAA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow
AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit
More informationEnergy Transformation. Metabolism = total chemical reactions in cells.
Energy Transformation Metabolism = total chemical reactions in cells. metabole = change Metabolism is concerned with managing the material and energy resources of the cell -Catabolism -Anabolism -Catabolism
More informationA Producer-Consumer Model With Stoichiometry
A Producer-Consumer Model With Stoichiometry Plan B project toward the completion of the Master of Science degree in Mathematics at University of Minnesota Duluth Respectfully submitted by Laura Joan Zimmermann
More informationWorksheet 1.1. Chapter 1: Quantitative chemistry glossary
Worksheet 1.1 Chapter 1: Quantitative chemistry glossary Amount The number of moles of a substance present in a sample. Aqueous solution A solution with water as the solvent. Atmosphere The unit atmosphere
More informationSimple kinetics of enzyme action
Simple kinetics of enzyme action It is established that enzymes form a bound complex to their reactants (i.e. substrates) during the course of their catalysis and prior to the release of products. This
More informationChapter 7: Metabolic Networks
Chapter 7: Metabolic Networks 7.1 Introduction Prof. Yechiam Yemini (YY) Computer Science epartment Columbia University Introduction Metabolic flux analysis Applications Overview 2 1 Introduction 3 Metabolism:
More informationExam 1 Chemical Reaction Engineering 26 February 2001 Closed Book and Notes
Exam 1 Chemical Reaction Engineering 26 February 21 Closed Book and Notes (2%) 1. Derive the unsteady-state mole balance for a chemical species A for a packed bed reactor using the following steps: a)
More informationBIOLOGY 10/11/2014. An Introduction to Metabolism. Outline. Overview: The Energy of Life
8 An Introduction to Metabolism CAMPBELL BIOLOGY TENTH EDITION Reece Urry Cain Wasserman Minorsky Jackson Outline I. Forms of Energy II. Laws of Thermodynamics III. Energy and metabolism IV. ATP V. Enzymes
More informationBoat Reactor for Low Pressure Chemical Vapor Deposition
Boat Reactor for Low Pressure Chemical Vapor Deposition Introduction Chemical vapor deposition (CVD) is an important step in the process of manufacturing microchips. A common application is the deposition
More informationBimetallic ruthenium-copper nanoparticles embedded in. mesoporous carbon as an effective hydrogenation catalyst
Supporting Information Bimetallic ruthenium-copper nanoparticles embedded in mesoporous carbon as an effective hydrogenation catalyst Jiajia Liu, *a Li Li Zhang, b Jiatao Zhang, a Tao Liu, c and X. S.
More informationModule 1: Mole Balances, Conversion & Reactor Sizing (Chapters 1 and 2, Fogler)
CHE 309: Chemical Reaction Engineering Lecture-2 Module 1: Mole Balances, Conversion & Reactor Sizing (Chapters 1 and 2, Fogler) Module 1: Mole Balances, Conversion & Reactor Sizing Topics to be covered
More informationWUFI Workshop at NTNU /SINTEF Fundamentals
WUFI Workshop at NTNU /SINTEF 2008 Fundamentals Contents: From steady-state to transient Heat storage and -transport Moisture storage and -transport Calculation of coupled transport Model limitations 2
More informationTo increase the concentration of product formed in a PFR, what should we do?
To produce more moles of product per time in a flow reactor system, what can we do? a) Use less catalyst b) Make the reactor bigger c) Make the flow rate through the reactor smaller To increase the concentration
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) 1D Heat Equation: Derivation Current Semester 1 / 19 Introduction The derivation of the heat
More informationEXPERIMENT 1 DETERMINATION OF GAS DIFFUSION COEFFICIENT
EXPERIMENT 1 DETERMINATION OF GAS DIFFUSION COEFFICIENT Objective: The objective of this experiment is to calculate diffusion coefficient of a volatile organic compound in air by means of Chapman Enskog
More informationEnergy and the Cell. All living things need energy to survive and do work.
Energy and the Cell EQ: How do cells acquire energy? EQ: Why is the relationship between plants and animals essential to life? All living things need energy to survive and do work. Organisms who depend
More informationAP Biology Cellular Respiration
AP Biology Cellular Respiration The bonds between H and C represents a shared pair of electrons These are high-energy electrons This represents chemical potential energy Hydro-carbons posses a lot of chemical
More informationMultistability in the lactose utilization network of Escherichia coli
Multistability in the lactose utilization network of Escherichia coli Lauren Nakonechny, Katherine Smith, Michael Volk, Robert Wallace Mentor: J. Ruby Abrams Agenda Motivation Intro to multistability Purpose
More informationPrinciples of Synthetic Biology: Midterm Exam
Principles of Synthetic Biology: Midterm Exam October 28, 2010 1 Conceptual Simple Circuits 1.1 Consider the plots in figure 1. Identify all critical points with an x. Put a circle around the x for each
More information