using "FDM";

FlowSheet FDM_din

	PARAMETERS
	L as Real;
	vs as Real;
	Pe as Real;
	qGmax as Real;
	ks1 as Real;
	Yxs as Real;
	fHAc as Real;
	fHBu as Real;
	fHPr as Real;
	fHEt as Real;
	fH2 as Real;
	R1 as Real;
	SHmd as Real;
	Kaac as Real;
	Kabu as Real;
	Kapr as Real;
	KaCO2 as Real;
	Scat as Real;
	San as Real;
	Kw as Real;
	pHLL as Real;
	pHUL as Real;
	n as Real;
	KpH as Real;
	Cy as Real;
	Cpr as Real;
	Cbu as Real;
	Cac as Real;
	Cbac as Real;
	kdecy as Real;
	Cxc as Real;
	KHCO2 as Real;
	KHH2 as Real;
	klaCO2 as Real;
	klaH2 as Real;
	T as Real;
	R as Real;
	Patm as Real;
	PH2O as Real;
	ro as Real;
	g as Real;
	eqCOD_SHAc as Real;
	eqCOD_SAcm as Real;
	eqCOD_SHBu as Real;
	eqCOD_SBum as Real;
	eqCOD_SHPr as Real;
	eqCOD_SPrm as Real;
	eqCOD_SH2 as Real;
	

	
	DEVICES
	#Carregando o modulo de diferencas finitas
	ex as MDF;
	
	SET
	#Especificando o valor de x0
	ex.xi = 0;
	#Especificando o valor de xf
	ex.xf = 1;
	#Especificando o valor de N
	ex.N = 100;


	L = 0.5;
	vs = 0.05;
	Pe = 150;
	qGmax = 1.03;
	ks1 = 0.0026;
	Yxs = 0.229;
	fHAc = 0.464;
	fHBu = 0.83;
	fHPr = 0;
	fHEt = 0;
	fH2 = 3.66;
	R1 = 0.0000000001;
	SHmd = 1;
	Kaac = 0.000017378;
	Kabu = 0.0000151;
	Kapr = 0.0000132;
	KaCO2 = 0.000000494;
	Scat = 0.02;
	San  = 0.02;
	Kw = 0.00000000000001;
	pHLL = 4;
	pHUL = 5.5;
	n = 2;
	KpH = 0.000018;
	Cy   = 0.0313;
	Cpr  = 0.268;
	Cbu  = 0.025;
	Cac  = 0.0313;
	Cbac = 0.0313;
	kdecy = 0.00083;
	Cxc = 0.02786;
	KHCO2 = 0.000000271;
	KHH2 =  0.00000000738;
	klaCO2 = 1.68;
	klaH2 = 0.4;
	T = 298;
	R = 8314.472;
	Patm = 106900;
	PH2O =   5570;
	ro = 1000;
	g = 9.81;
	eqCOD_SHAc=64;
	eqCOD_SAcm=64;
	eqCOD_SHBu=160;
	eqCOD_SBum=160;
	eqCOD_SHPr=112;
	eqCOD_SPrm=112;
	eqCOD_SH2=16;
	


	EQUATIONS
	#Equacao do modelo
	(L/vs)*diff(ex.y)*'s'   =    (1/Pe)*ex.dif2x   -    ex.dif1x   -(L/vs)*qGmax*(ex.y/(ks1+ex.y))*ex.y7*KpH^n/(KpH^n+(ex.y11)^n);
	(L/vs)*diff(ex.y2)*'s'  =    (1/Pe)*ex.dif22x  -    ex.dif12x  +(L/vs)*(1-Yxs)*fHAc*qGmax*(ex.y/(ks1+ex.y))*ex.y7*KpH^n/(KpH^n+(ex.y11)^n);
	(L/vs)*diff(ex.y3)*'s'  =    (1/Pe)*ex.dif23x  -    ex.dif13x  +(L/vs)*(1-Yxs)*fHBu*qGmax*(ex.y/(ks1+ex.y))*ex.y7*KpH^n/(KpH^n+(ex.y11)^n);
	(L/vs)*diff(ex.y4)*'s'  =    (1/Pe)*ex.dif24x  -    ex.dif14x  +(L/vs)*(1-Yxs)*fHPr*qGmax*(ex.y/(ks1+ex.y))*ex.y7*KpH^n/(KpH^n+(ex.y11)^n);
	(L/vs)*diff(ex.y5)*'s'  =    (1/Pe)*ex.dif25x  -    ex.dif15x  +(L/vs)*(1-Yxs)*fHEt*qGmax*(ex.y/(ks1+ex.y))*ex.y7*KpH^n/(KpH^n+(ex.y11)^n);
	(L/vs)*diff(ex.y6)*'s'  =    (1/Pe)*ex.dif26x  -    ex.dif16x  +(L/vs)*(1-Yxs)*fH2 *qGmax*(ex.y/(ks1+ex.y))*ex.y7*KpH^n/(KpH^n+(ex.y11)^n)-(L/vs)*klaH2*(ex.y6-KHH2*R*T*ex.y16);
	(L/vs)*diff(ex.y7)*'s'  = R1*(1/Pe)*ex.dif27x  - R1*ex.dif17x  +(L/vs)*         Yxs*qGmax*(ex.y/(ks1+ex.y))*ex.y7*KpH^n/(KpH^n+(ex.y11)^n);
	(L/vs)*diff(ex.y8)*'s'  =    (1/Pe)*ex.dif28x  -    ex.dif18x  +Kaac*(ex.y2-ex.y8)/(ex.y11);
	(L/vs)*diff(ex.y9)*'s'  =    (1/Pe)*ex.dif29x  -    ex.dif19x  +Kabu*(ex.y3-ex.y9)/(ex.y11);
	(L/vs)*diff(ex.y10)*'s' =    (1/Pe)*ex.dif210x -    ex.dif110x +Kapr*(ex.y4-ex.y10)/(ex.y11);
	(L/vs)*diff(ex.y11)*'s' =    (1/Pe)*ex.dif211x -    ex.dif111x -(Scat-San -ex.y8/64-ex.y9/160 -ex.y10/112+ex.y14)/2+0.5*sqrt((Scat-San -ex.y8/64-ex.y9/160 -ex.y10/112+ex.y14)^2+4*Kw);
	(L/vs)*diff(ex.y12)*'s' =    (1/Pe)*ex.dif212x -    ex.dif112x -(-Cy+(1-Yxs)*(fHPr*Cpr+fHBu*Cbu+fHAc*Cac)
																   +Yxs*Cbac)*qGmax*(L/vs)*(ex.y/(ks1+ex.y))*ex.y7*KpH^n/(KpH^n+(ex.y11)^n)
																   -(L/vs)*((-Cbac+Cxc)*kdecy*ex.y7)-(L/vs)*klaCO2*(ex.y13-KHCO2*R*T*ex.y15);
	(L/vs)*diff(ex.y13)*'s' =    (1/Pe)*ex.dif213x -    ex.dif113x +(ex.y12-ex.y14);
	(L/vs)*diff(ex.y14)*'s' =    (1/Pe)*ex.dif214x -    ex.dif114x +KaCO2*(ex.y12-ex.y14)/(ex.y11);
	(L/vs)*diff(ex.y15)*'s' =    (1/Pe)*ex.dif215x -    ex.dif115x -(L/vs)*(ex.y15* ex.y17/(0.005*ex.dx*0.005027))+(L/vs)*klaCO2*(ex.y13-KHCO2*R*T*ex.y15)*0.47/0.05;
	(L/vs)*diff(ex.y16)*'s' =    (1/Pe)*ex.dif216x -    ex.dif116x -(L/vs)*(ex.y16* ex.y17/(0.005*ex.dx*0.005027))+(L/vs)*klaH2*(ex.y6-KHH2*R*T*ex.y16)*0.47/0.05;
	(L/vs)*diff(ex.y17)*'s' = R1*(1/Pe)*ex.dif217x - R1*ex.dif117x +(L/vs)*(R*T/(Patm-PH2O+ro*g*(L-ex.x*L)))*ex.dx*0.005027*(klaCO2*(ex.y13-KHCO2*R*T*ex.y15)+klaH2*(ex.y6-KHH2*R*T*ex.y16)/16);
	(L/vs)*diff(ex.y18)*'s' =    (1/Pe)*ex.dif218x -    ex.dif118x;
	(L/vs)*diff(ex.y19)*'s' =    (1/Pe)*ex.dif219x -    ex.dif119x;
	(L/vs)*diff(ex.y20)*'s' =    (1/Pe)*ex.dif220x -    ex.dif120x;

		
	#Especificando as condicoes de contorno
	#Para x=xf
	ex.dif1x(ex.N+1)=0;
	ex.dif12x(ex.N+1)=0;
	ex.dif13x(ex.N+1)=0;
	ex.dif14x(ex.N+1)=0;
	ex.dif15x(ex.N+1)=0;
	ex.dif16x(ex.N+1)=0;
	ex.dif17x(ex.N+1)=0;
	ex.dif18x(ex.N+1)=0;
	ex.dif19x(ex.N+1)=0;
	ex.dif110x(ex.N+1)=0;
	ex.dif111x(ex.N+1)=0;
	ex.dif112x(ex.N+1)=0;
	ex.dif113x(ex.N+1)=0;
	ex.dif114x(ex.N+1)=0;
	ex.dif115x(ex.N+1)=0;
	ex.dif116x(ex.N+1)=0;
	ex.dif117x(ex.N+1)=0;
	ex.dif118x(ex.N+1)=0;
	ex.dif119x(ex.N+1)=0;
	ex.dif120x(ex.N+1)=0;


	#Para x=xi
#	-1/Pe * ex.dif1x(1) = 1 - ex.y(1);	
	ex.y(1)=10;
	ex.y2(1)=0;
	ex.y3(1)=0;
	ex.y4(1)=0;
	ex.y5(1)=0;
	ex.y6(1)=0;
	ex.y7(1)=0.01;
	ex.y8(1)=0;
	ex.y9(1)=0;
	ex.y10(1)=0;
	ex.y11(1)=0.0000001;
	ex.y12(1)=0.1275;
	ex.y13(1)=0;
	ex.y14(1)=0.1275;
	ex.y15(1)=0;
	ex.y16(1)=0;
	ex.y17(1)=0;
	ex.y18(1)=0;
	ex.y19(1)=0;
	ex.y20(1)=0;


	INITIAL
	#Condicao inicial
	ex.y(2:ex.N)=0;
	ex.y2(2:ex.N)=0;
	ex.y3(2:ex.N)=0;
	ex.y4(2:ex.N)=0;
	ex.y5(2:ex.N)=0;
	ex.y6(2:ex.N)=0;
	ex.y7(2:ex.N)=0.01;
	ex.y8(2:ex.N)=0;
	ex.y9(2:ex.N)=0;
	ex.y10(2:ex.N)=0;
	ex.y11(2:ex.N)=0.0000001;
	ex.y12(2:ex.N)=0;
	ex.y13(2:ex.N)=0;
	ex.y14(2:ex.N)=0;
	ex.y15(2:ex.N)=0;
	ex.y16(2:ex.N)=0;
	ex.y17(2:ex.N)=0;
	ex.y18(2:ex.N)=0;
	ex.y19(2:ex.N)=0;
	ex.y20(2:ex.N)=0;


	OPTIONS
	TimeStep = 1;
	TimeEnd = 150;
	TimeUnit = 's';
	Dynamic = true;
end