Last change
on this file since 22 was
22,
checked in by Rafael de Pelegrini Soares, 16 years ago

Added optimization samples


Property svn:eolstyle set to
native

Property svn:keywords set to
Id

File size:
1.2 KB

Line  

1  using "stage_separators/flash"; 

2  

3  # exemplo baseado nos dados do artigo do Gani 1985. 

4  FlowSheet FlashSteadyTest 

5  PARAMETERS 

6  PP as CalcObject(Brief="Physical Properties",File="vrpp"); 

7  NComp as Integer; 

8  

9  VARIABLES 

10  Q as heat_rate (Brief="Heat supplied"); 

11  leves as fraction; 

12  

13  SET 

14  PP.Components = ["1,3butadiene", "isobutene", "npentane", "1pentene", "1hexene", "benzene"]; 

15  PP.LiquidModel = "PR"; 

16  PP.VapourModel = "PR"; 

17  NComp = PP.NumberOfComponents; 

18  

19  DEVICES 

20  fl as flash_Steady; 

21  s1 as stream_therm; 

22  

23  CONNECTIONS 

24  s1 to fl.Inlet; 

25  Q to fl.Q; 

26  

27  EQUATIONS 

28  leves = fl.OutletV.z(1)+fl.OutletV.z(2); 

29  

30  SPECIFY 

31  s1.F = 496.3 * "kmol/h"; 

32  s1.T = 338 * "K"; 

33  s1.P = 507.1 * "kPa"; 

34  s1.v = 0.1380; 

35  s1.z = [0.2379,0.3082,0.09959,0.1373,0.08872,0.1283]; 

36  

37  fl.OutletL.P = 2.5 * "atm"; 

38  

39  #Q = 0 * "kJ/h"; 

40  fl.OutletL.T = 315.06 * "K"; 

41  

42  OPTIONS 

43  mode = "steady"; 

44  outputLevel = "all"; 

45  end 

46  

47  #* 

48  * Maximização da pureza dos leves no topo do flash. 

49  * Achou : 0.83 frente a 0.79 da simulação do FlowSheet 

50  *# 

51  Optimization FlashOpt1 as FlashSteadyTest 

52  MAXIMIZE 

53  leves; 

54  

55  FREE 

56  fl.OutletL.T; 

57  

58  EQUATIONS 

59  fl.OutletL.T < 315 * "K"; 

60  fl.OutletL.T > 300 * "K"; 

61  

62  OPTIONS 

63  mode = "steady"; 

64  outputLevel = "high"; 

65  end 

66  

Note: See
TracBrowser
for help on using the repository browser.