1 | #*--------------------------------------------------------------------- |
---|

2 | * EMSO Model Library (EML) Copyright (C) 2004 - 2016 ALSOC. |
---|

3 | * |
---|

4 | * This LIBRARY is free software; you can distribute it and/or modify |
---|

5 | * it under the therms of the ALSOC FREE LICENSE as available at |
---|

6 | * http://www.enq.ufrgs.br/alsoc. |
---|

7 | * |
---|

8 | * EMSO Copyright (C) 2004 - 2016 ALSOC, original code |
---|

9 | * from http://www.rps.eng.br Copyright (C) 2002-2004. |
---|

10 | * All rights reserved. |
---|

11 | * |
---|

12 | * EMSO is distributed under the therms of the ALSOC LICENSE as |
---|

13 | * available at http://www.enq.ufrgs.br/alsoc. |
---|

14 | * |
---|

15 | *---------------------------------------------------------------------- |
---|

16 | * Author: Argimiro R. Secchi |
---|

17 | * COPPE/UFRJ - Group of Modeling, Simulation, Control, |
---|

18 | * and Optimization of Processes |
---|

19 | * $Id$ |
---|

20 | *--------------------------------------------------------------------*# |
---|

21 | |
---|

22 | FlowSheet OCM_din |
---|

23 | |
---|

24 | PARAMETERS |
---|

25 | MCO as Plugin(Type="OCFEM",Boundary="BOTH", |
---|

26 | InternalPoints=5, |
---|

27 | alfa=1, beta=1); |
---|

28 | |
---|

29 | np as Integer; |
---|

30 | nd as Integer; |
---|

31 | N as Integer; |
---|

32 | A(N) as Real(Brief="recebe valores de do vetor A"); |
---|

33 | B(N) as Real(Brief="recebe valores de do vetor B"); |
---|

34 | r(np) as Real(Brief="raizes de jacob"); |
---|

35 | |
---|

36 | Da as Real; |
---|

37 | Pe as Real; |
---|

38 | m as Real; |
---|

39 | |
---|

40 | SET |
---|

41 | np= MCO.NodalPoints; |
---|

42 | N= MCO.NodalPoints * MCO.NodalPoints; |
---|

43 | A = MCO.matrixA; |
---|

44 | B = MCO.matrixB; |
---|

45 | r = MCO.roots; |
---|

46 | |
---|

47 | Da = 5; |
---|

48 | m = 1; |
---|

49 | Pe = 2; |
---|

50 | |
---|

51 | VARIABLES |
---|

52 | y(np) as Real(Brief="variavel dependente"); |
---|

53 | |
---|

54 | mA(np,np) as Real(Brief="cria matriz A"); |
---|

55 | mB(np,np) as Real(Brief="cria matriz B"); |
---|

56 | |
---|

57 | dif1x(np) as Real(Brief="calcula 1derivada"); |
---|

58 | dif2x(np) as Real(Brief="calcula 2derivada"); |
---|

59 | |
---|

60 | EQUATIONS |
---|

61 | #comando para criar as matrizes |
---|

62 | for j in [1:np] do |
---|

63 | mA(:,j) = A((j-1)*np+[1:np]); |
---|

64 | mB(:,j) = B((j-1)*np+[1:np]); |
---|

65 | end |
---|

66 | |
---|

67 | #calculando as derivadas |
---|

68 | dif1x = sumt(mA*y); |
---|

69 | dif2x = sumt(mB*y); |
---|

70 | |
---|

71 | diff(y(2:np-1))*'s' + dif1x(2:np-1) = 1/Pe * dif2x(2:np-1) - Da*y(2:np-1)^m; |
---|

72 | |
---|

73 | -1/Pe * dif1x(1) = 1 - y(1); |
---|

74 | # y(1) = 1; |
---|

75 | |
---|

76 | dif1x(np) = 0; |
---|

77 | |
---|

78 | INITIAL |
---|

79 | y(2:np-1) = 0; |
---|

80 | |
---|

81 | OPTIONS |
---|

82 | TimeStep = 0.01; |
---|

83 | TimeEnd = 0.2; |
---|

84 | TimeUnit = 's'; |
---|

85 | Dynamic = true; |
---|

86 | end |
---|