#*---------------------------------------------------------------------
* This file is property of the author and cannot be used, copyed
* or modified without permission.
*
* Copyright (C) 2002-2004  the author
*----------------------------------------------------------------------
* Author: Rafael de Pelegrini Soares
* $Id: sample_bio.mso 1 2006-06-20 17:33:53Z rafael $
*----------------------------------------------------------------------
*
* EMSO sample file.
*
* FlowSheet with a bio-reactor based on Monod model.
*
*-------------------------------------------------------------------*#

using "types";

FlowSheet bio
 	PARAMETERS
	miMax   as frequency(Brief="Monod parameter", Unit="1/h", Default=0.53);
	D       as frequency(Brief="Dilution rate F/V", Unit="1/h");
	Km      as conc_mass(Brief="Monod parameter", Default=0.12);
	K1      as inv_conc_mass(Brief="Monod parameter", Default=0.4545);
	Y	as coefficient(Brief="Yield biomass/substrate", Default=0.4);
	x2f     as conc_mass(Brief="Substrate feed concentration", Default=4);
 
	VARIABLES
	biomass   as conc_mass (Brief="Biomass concentration");
	substrate as conc_mass (Brief="Substrate concentration");
	mi        as frequency (Brief="Specific growth rate", Unit="1/h");
	
 	EQUATIONS
	"Specific growth rate"
	mi = miMax*substrate/(Km+substrate+(K1*substrate^2));
	
	"Biomass production"
	diff(biomass)/biomass = mi-D;
	
	"Substrate consumption"
	diff(substrate) = D*(x2f-substrate)-(biomass*mi/Y);
	
 	SET
 	D      = 0.3  *  "1/h";
	
  	INITIAL
	biomass = 1  *  "kg/m^3";
	substrate = 0.5  *  "kg/m^3";
 
	OPTIONS
	time = [0:0.5:10]*"h";
	outputLevel = "high";
	LASolver = "dense"; # dense or sparse
end