using "types";

FlowSheet Pei_Ying_Lin_2020

	
	PARAMETERS
	
	qGmax as coefficient (Unit = 'mol/mol/h');
	ks1 as conc_mol; 
	Yxs as coefficient (Unit = 'mol/mol');
	fHAc as coefficient (Unit = 'mol/mol');
	fHBu as coefficient (Unit = 'mol/mol');
	fHEt as coefficient (Unit = 'mol/mol');
	fH2 as coefficient (Unit = 'mol/mol');
	
	VARIABLES

	SGli as conc_mol;
	SHAc as conc_mol;
	SHBu as conc_mol;
	SHEt as conc_mol;
	X1 as conc_mol;
	
	

	EQUATIONS
	
	"Biological Phase" 

	
	
	diff(SGli)=             -qGmax*(SGli/(ks1+SGli))*X1;
	diff(SHAc)= (1-Yxs)*fHAc*qGmax*(SGli/(ks1+SGli))*X1;
	diff(SHBu)= (1-Yxs)*fHBu*qGmax*(SGli/(ks1+SGli))*X1;
	diff(SHEt)= (1-Yxs)*fHEt*qGmax*(SGli/(ks1+SGli))*X1;
	diff(X1)  =          Yxs*qGmax*(SGli/(ks1+SGli))*X1;
	
	
	
	
	SET
	
	qGmax = 1.03*'mol/mol/h';
	ks1 = 0.0026*'mol/l'; 
	Yxs = 0.229*'mol/mol';
	fHAc = 0.464*'mol/mol';
	fHBu = 0.83*'mol/mol';
	fHEt = 0.0*'mol/mol';
	fH2 = 3.66*'mol/mol';
	
	INITIAL
	
	SGli = 0.0175*'mol/l';
	X1 = 0.0011*'mol/l';
	SHAc = 0.0005*'mol/l';
	SHBu = 0.0009*'mol/l';
	SHEt = 0.0006*'mol/l';
	
	
	
	OPTIONS

	TimeStart = 0;
	TimeStep = 0.5;
	TimeEnd = 30;
	TimeUnit = 'h';
end

Estimation Pei_Ying_Lin_2020_Est as Pei_Ying_Lin_2020
	ESTIMATE
	# PARAMETER  START     LOWER     UPPER     UNIT
	qGmax        1.03      0.0001    5        'mol/mol/h';
	ks1          0.0026    0.0001    1        'mol/l';
	Yxs          0.229     0.0001    1        'mol/mol';
	
	EXPERIMENTS
	# DATA FILE             WEIGTH
	"Pei_Ying_2020.dat"      1;
	
	OPTIONS
	
	NLPSolver(
		MaxIterations = 1000,
		File = "complex"
		#File = "ipopt_emso"
	);
end