source: trunk/sample/miscellaneous/ocfem/OCFEM_ss.mso @ 976

#*---------------------------------------------------------------------
* EMSO Model Library (EML) Copyright (C) 2004 - 2016 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 - 2016 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.
*
*----------------------------------------------------------------------
* Author: Argimiro R. Secchi
* COPPE/UFRJ - Group of Modeling, Simulation, Control,
*                               and Optimization of Processes
* $Id$
*--------------------------------------------------------------------*#

using "MODEL_FE";

FlowSheet OCFEM_ss
        
        DEVICES
        #Carregando o modulo do MCO em elementos finitos        
        d as MCO_EF;
        
        PARAMETERS 

        Da as Real;
        Pe as Real; 
        m as Real;      

        SET

        Da = 5;
        m = 1;
        Pe = 2;
        
        SPECIFY
        #Define os valores das fronteiras dos elementos finitos(vetor h) e valor da variável
        #independente para cálculo da interpolação

        d.h=[0,1/4,2/4,3/4,1];
        
        EQUATIONS
        
        #Equacao do problema
        for k in [1:d.ne] do
                for i in [2:d.np-1] do
                        d.dif1x(i,k)  = (1/Pe) * d.dif2x(i,k) - Da*(d.y(i,k))^m;
                end
        end
        
        #Especificando as condicoes de contorno
        #Para x=x0
#               (-1/Pe) * d.dif1x(1,1) = 1 - d.y(1,1);
                d.y(1,1) = 1;
        
        #Para x=xf
                d.dif1x(d.np,d.ne) = 0;

        OPTIONS
        Dynamic = false;
        
end