#*---------------------------------------------------------------------
* 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 of an equilibrium reactor
*----------------------------------------------------------------------
*
* This sample file needs VRTherm (www.vrtech.com.br) to run.
*
*----------------------------------------------------------------------
*
*   This is a sample based on a problem of the book "Introduction to
* Chemical Engineering Thermodynamics" by Smith, Van Ness and Abbott,
* 5th ed., 1996.
*
*----------------------------------------------------------------------
* Author: Rodolfo Rodrigues
* $Id$
*--------------------------------------------------------------------*#

using "reactors/equil";


#*---------------------------------------------------------------------
* Resolved example 15.5 - p.527 (Brazilian edition)
*
* System in gas fase contents CO, H2O, CO2, and H2.
*
*	Co + H2O -> CO2 + H2
*
* Initial conditions
*   1bar, 1100K, 1mol of CO, and 1mol of H2O
*
* Specify
*	Outlet.T
*	Outlet.P
*
* Return
*	Outlet.z
----------------------------------------------------------------------*#


FlowSheet sample_equil_vap
	PARAMETERS
	PP 		as Plugin(Brief="External physical properties", Type="PP",
			Components = ["carbon monoxide","water","carbon dioxide","hydrogen"],
			LiquidModel = "IdealLiquid",
			VapourModel = "Ideal");
	NComp	as Integer;
	NReac	as Integer;
	
	DEVICES
	Fin 	as source;
	R 		as equil_vap;

	CONNECTIONS
	Fin.Outlet 	to R.Inlet;
	
	SET
	NComp = PP.NumberOfComponents;
	
	R.NReac = 1;
	R.stoic(:,1) = [-1.0, -1.0, 1.0, 1.0];
	
	SPECIFY
	Fin.Outlet.F = 2*'mol/s';
	Fin.Outlet.z = [0.5, 0.5, 0.0, 0.0];
	Fin.Outlet.P = 1*'atm';
	Fin.Outlet.T = 1100*'K';
	
	R.Outlet.P = 1*'bar';
	R.Outlet.T = Fin.Outlet.T;
	R.Outlet.F = 2*'mol/s'*sqrt(R.Tank.Level/'m');

	R.Tank.Across = 5*'m^2';
	R.Tank.L = 1.5*'m';
	
	INITIAL
    R.Outletm.T = 1100*'K';
	R.M = [1.0, 1.0, 0.0, 0.0]*'mol';

	OPTIONS
	Dynamic = false;
end