#*--------------------------------------------------------------------- * 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