#*---------------------------------------------------------------------
* 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_din
	
	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
			diff(d.y(i,k))*'s' + 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;

	INITIAL
	d.y(2:d.np-1,1:d.ne) = 0;

	OPTIONS
	TimeStep = 0.02;
	TimeEnd = 0.5;
	TimeUnit = 's';
	Dynamic = true;
	
end