source: mso/eml/reactors/cstr.mso @ 24

#*---------------------------------------------------------------------
* This file is property of the author and cannot be used, copyed
* or modified without permission.
*
* Copyright (C) 2004  the author
*----------------------------------------------------------------------
* Author: Paula B. Staudt
* $Id: cstr.mso 1 2006-06-20 17:33:53Z rafael $
*--------------------------------------------------------------------*#

using "streams";

Model cstr_basic
        PARAMETERS
ext     NComp as Integer;
        NReac as Integer;
        stoic(NComp, NReac) as Real (Brief = "Stoichiometric Matrix");

        VARIABLES
in      Inlet  as stream;
out     Outlet as stream_therm;
        q    as heat_rate;
        Vr   as volume (Brief="Reacting Volume");
        M    as mol (Brief = "Molar total amount");
        C(NComp)  as conc_mol(Brief="Components concentration");

        EQUATIONS
        "Molar fraction"
        Outlet.z * M = C * Vr;
        "Vapourisation Fraction"
        Outlet.v = Inlet.v;
end

Model cstr_yield as cstr_basic
        PARAMETERS
        yield(NReac) as fraction;
        
        EQUATIONS
        Outlet.z*Outlet.F = Inlet.z*Inlet.F * sumt(stoic*(1-yield));
        "Reactor Energy Balance"
        diff(M*Outlet.h) = Inlet.F*Inlet.h - Outlet.F*Outlet.h - q;
end


Model cstr as cstr_basic
        VARIABLES
        r(NReac) as reaction_mol;
        Hr(NReac) as heat_reaction;
        
        EQUATIONS
        "Component Molar Balance"
        diff(Outlet.z*M) = (Inlet.F*Inlet.z - Outlet.F*Outlet.z) + sumt(stoic*r)*Vr;

        "Reactor Energy Balance"
        diff(M*Outlet.h) = Inlet.F*Inlet.h - Outlet.F*Outlet.h +
                sum(Hr*sum(stoic*r))*Vr - q;
end