model.growth: Example: Growth model

In baccione-eawag/EawagSchoolTools: Miscellaneous Tools for Eawag Summer School in Environmental Systems Analysis

Description

Defines a simple growth model, which is used in the exercises for the Eawag Summer School in Environmental Systems Analysis.

Usage

model.growth(par,L)

Arguments

par	named vector of parameters
L	observation layout

Details

Calculates the growth of micro-organisms on a substrate in a batch reactor:

dC_M/dt = mu * C_S/(K+C_S) * C_M - b*C_M

dC_S/dt = -(mu/Y) * C_S/(K+C_S) * C_M

State variables are:

C_M	concentration of micro-organisms in the reactor
C_S	concentration of substrate in the reactor

Parameters are:

mu	maximum growth rate of micro-organisms
K	half-concentration of growth rate with respect to substrate
b	rate of death and respiration processes of micro-organisms
Y	yield of growth process
C_M_ini	initial concentration of micro-organisms
C_S_ini	initial concentration of substrate

The observation layout L contains the variable names and time points at which the model is to be evaluated. L should be a data frame with two columns. The first column contains the names of the variables to be evaluated, and the second column contains the time points (or other evaluation steps) for the model. See function checkLayout in package bayesBias for more details. For example, for the growth model:

var	t
C_M	1
C_M	2
C_M	3
...	...
C_S	1
C_S	2
C_S	3
...	...

Value

Returns a named vector of model output (concentrations). Names consist of the variable name and evaluation point (time), separated by an underscore (for example C_M_1).

Author(s)

Peter Reichert <[e-mail]>, Lukas M. Weber

See Also

model.monod, model.monod.external, and function checkLayout in package bayesBias

baccione-eawag/EawagSchoolTools documentation built on Dec. 19, 2021, 6:38 a.m.