#*-------------------------------------------------------------------
* EMSO Model Library (EML) Copyright (C) 2004 - 2007 ALSOC.
*
* This LIBRARY is free software; you can distribute it and/or modify
* it under the therms of the ALSOC FREE LICENSE as available at
* http://www.enq.ufrgs.br/alsoc.
*
* EMSO Copyright (C) 2004 - 2007 ALSOC, original code
* from http://www.rps.eng.br Copyright (C) 2002-2004.
* All rights reserved.
*
* EMSO is distributed under the therms of the ALSOC LICENSE as
* available at http://www.enq.ufrgs.br/alsoc.
*
*--------------------------------------------------------------------
* Sample file for model cstr
*--------------------------------------------------------------------
* 
* This sample file needs VRTherm (www.vrtech.com.br) to run.
*
*----------------------------------------------------------------------
* Author: Paula B. Staudt
* $Id: sample_cstr.mso 245 2007-04-19 18:37:24Z bicca $
*--------------------------------------------------------------------*#

using "reactors/cstr";

FlowSheet cstr_PropyleneGlycol

PARAMETERS
	PP   		as Plugin(Brief="Physical Properties",
		Type = "PP",
		Components  = ["1,2-propylene oxide", "methanol","water","1,2-propylene glycol" ],
		LiquidModel = "PR",
		VapourModel = "PR"
	);
	NComp 	as Integer;
	R 			as Real		(Default=8.3145, Unit='J/mol/K');

DEVICES
	s1   	as source;
	Reac 	as cstr;

SET
	NComp 			= PP.NumberOfComponents;
	Reac.NReac 	= 1;
	Reac.stoic 		= [-1, 0, -1, 1];
	
EQUATIONS

	Reac.Hr = -36400 * 'Btu/lbmol' - 7 * 'Btu/lbmol/degR' * (Reac.Outlet.T - 528 * 'degR');
	Reac.r = 16.96e12 / 'h' * exp(- 32400 * 'Btu/lbmol'/R/Reac.Outlet.T)
		* Reac.C(1);

SPECIFY
	s1.Outlet.F = 52126 * 'mol/h';
	s1.Outlet.P = 1     * 'atm';
	s1.Outlet.T = 534.67* 'degR';
	s1.Outlet.z = [0.0468, 0.0783, 0.8749, 0];

	Reac.Vr    		= 300   * 'gal';
	Reac.Outlet.F	= Reac.Inlet.F;

	"Adiabatic"
	Reac.q     		= 0     * 'J/s';
	Reac.Outlet.P	= 1     * 'atm';
	
CONNECTIONS

	s1.Outlet  to  Reac.Inlet;
	
INITIAL

	Reac.Outlet.T = 533 * 'degR';
	Reac.M = 100 * 'kmol';
	Reac.Outlet.z(1:3) = [0.0468, 0.0783, 0.874];

OPTIONS
	TimeStep 	= 0.1;
	TimeEnd 	= 5;

end