Last change
on this file since 390 was
390,
checked in by Argimiro Resende Secchi, 15 years ago

Adding AUTO DAE examples.

File size:
1.1 KB

Rev  Line  

[390]  1  using "types"; 

 2  

 3  FlowSheet cstr_dae 

 4  

 5  PARAMETERS 

 6  fi as Real; 

 7  beta as Real; 

 8  delta as Real; 

 9  q as Real; 

 10  alfa as Real; 

 11  S as Real; 

 12  csi as Real; 

 13  delta1 as Real; 

 14  delta2 as Real; 

 15  gama as Real; 

 16  x1f as Real; 

 17  x2f as Real; 

 18  x3f as Real; 

 19  x4f as Real; 

 20  qc as Real; 

 21  

 22  VARIABLES 

 23  u1 as Real; 

 24  u2 as Real; 

 25  u3 as Real; 

 26  u4 as Real; 

 27  u5 as Real; 

 28  u6 as Real; 

 29  

 30  SET 

 31  fi = 0.04; 

 32  beta = 8; 

 33  delta = 0.78; 

 34  q = 1; 

 35  alfa = 0.19; 

 36  S = 1.015; 

 37  csi = 0.32; 

 38  delta1 = 10; 

 39  delta2 = 0.952; 

 40  gama = 27.85; 

 41  x1f = 1; 

 42  x2f = 0; 

 43  x3f = 0; 

 44  x4f = 1; 

 45  qc = 0; 

 46  

 47  EQUATIONS 

 48  

 49  q*(x1fu1)u1*u5*fi  diff(u1)*'s' = 0; 

 50  q*(x2fu2)u2*fi*S*u6+u1*fi*u5  diff(u2)*'s' = 0; 

 51  q*(x3fu3)+delta*(u4u3)+beta*fi*(u1*u5+alfa*u2*u6*S)  diff(u3)*'s' = 0; 

 52  delta1*(qc*(x4fu4)+delta*delta2*(u3u4))  diff(u4)*'s' = 0; 

 53  u5 + exp(u3/(1+(u3/gama))) = 0; 

 54  u6 + exp(csi*u3/(1+(u3/gama))) = 0; 

 55  

 56  GUESS 

 57  u1 = 0.6419; 

 58  u2 = 0.3270; 

 59  u3 = 2.9101; 

 60  u4 = 2.9101; 

 61  u5 = 13.9403; 

 62  u6 = 2.3236; 

 63  

 64  INITIAL 

 65  u1 = 0.6; 

 66  u2 = 0.3; 

 67  u3 = 3; 

 68  u4 = 3; 

 69  

 70  OPTIONS 

 71  TimeEnd = 10; 

 72  TimeStep = 0.2; 

 73  Dynamic = false; 

 74  end 

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