/[ascend]/trunk/models/dyn_flash.a4l
ViewVC logotype

Diff of /trunk/models/dyn_flash.a4l

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1 by aw0a, Fri Oct 29 20:54:12 2004 UTC revision 576 by johnpye, Tue May 9 03:42:08 2006 UTC
# Line 249  END NOTES; Line 249  END NOTES;
249    
250      METHOD seqmod;      METHOD seqmod;
251          RUN state.specify;          RUN state.specify;
252          state.y[cd.components].fixed := FALSE;          FREE state.y[cd.components];
253          IF equilibrated THEN          IF equilibrated THEN
254          state.T.fixed :=FALSE;          FREE state.T;
255          ELSE          ELSE
256          state.T.fixed :=TRUE;          FIX state.T;
257          END IF;          END IF;
258          split_fraction['vapor'].fixed :=TRUE;          FIX split_fraction['vapor'];
259          state.phase_fraction['vapor'].fixed :=FALSE;          FREE state.phase_fraction['vapor'];
260          dMi_dt[cd.other_components].fixed   :=TRUE;          FIX dMi_dt[cd.other_components];
261          Mi[cd.components].fixed         :=FALSE;          FREE Mi[cd.components];
262          dM_dt.fixed     :=TRUE;          FIX dM_dt;
263          M.fixed         :=FALSE;          FREE M;
264          dU_dt.fixed     :=TRUE;          FIX dU_dt;
265          U.fixed         :=FALSE;          FREE U;
266          Vol.fixed       :=TRUE;          FIX Vol;
267          vol_liq.fixed   :=TRUE;          FIX vol_liq;
268          IF dynamic THEN          IF dynamic THEN
269          dMi_dt[cd.other_components].fixed   :=FALSE;          FREE dMi_dt[cd.other_components];
270          Mi[cd.other_components].fixed       :=TRUE;          FIX Mi[cd.other_components];
271          dM_dt.fixed     :=FALSE;          FREE dM_dt;
272          M.fixed     :=TRUE;          FIX M;
273          dU_dt.fixed     :=FALSE;          FREE dU_dt;
274          U.fixed     :=TRUE;          FIX U;
275          Qin.fixed   :=TRUE;          FIX Qin;
276          IF equilibrated THEN          IF equilibrated THEN
277              state.P.fixed       :=FALSE;              FREE state.P;
278          ELSE          ELSE
279              state.T.fixed       :=FALSE;              FREE state.T;
280          END IF;          END IF;
281          vol_liq.fixed   :=FALSE;          FREE vol_liq;
282          liqout.flow.fixed :=TRUE;          FIX liqout.flow;
283          vapout.flow.fixed :=TRUE;          FIX vapout.flow;
284          split_fraction['vapor'].fixed :=FALSE;          FREE split_fraction['vapor'];
285          ELSE          ELSE
286          dMi_dt[cd.other_components] :=0 {mole/sec};          dMi_dt[cd.other_components] :=0 {mole/sec};
287          dM_dt :=0 {mole/sec};          dM_dt :=0 {mole/sec};
# Line 509  END NOTES; Line 509  END NOTES;
509    
510      METHOD seqmod;      METHOD seqmod;
511          RUN state.specify;          RUN state.specify;
512          state.y[cd.components].fixed := FALSE;          FREE state.y[cd.components];
513          IF equilibrated THEN          IF equilibrated THEN
514          state.T.fixed :=FALSE;          FREE state.T;
515          ELSE          ELSE
516          state.T.fixed :=TRUE;          FIX state.T;
517          END IF;          END IF;
518          state.phase_fraction['vapor'].fixed :=FALSE;          FREE state.phase_fraction['vapor'];
519          cmo_ratio.fixed :=TRUE;          FIX cmo_ratio;
520          dMi_dt[cd.other_components].fixed :=TRUE;          FIX dMi_dt[cd.other_components];
521          Mi[cd.other_components].fixed :=FALSE;          FREE Mi[cd.other_components];
522          dM_dt.fixed :=TRUE;          FIX dM_dt;
523          M.fixed :=FALSE;          FREE M;
524          dU_dt.fixed :=TRUE;          FIX dU_dt;
525          U.fixed :=FALSE;          FREE U;
526          Vol.fixed :=TRUE;          FIX Vol;
527          vol_liq.fixed :=TRUE;          FIX vol_liq;
528          tauv.fixed :=FALSE;          FREE tauv;
529          taul.fixed :=FALSE;          FREE taul;
530              liqout.flow.fixed :=FALSE;              FREE liqout.flow;
531          vapout.flow.fixed :=FALSE;          FREE vapout.flow;
532          IF dynamic THEN          IF dynamic THEN
533          dMi_dt[cd.other_components].fixed :=FALSE;          FREE dMi_dt[cd.other_components];
534          Mi[cd.other_components].fixed :=TRUE;          FIX Mi[cd.other_components];
535          dM_dt.fixed :=FALSE;          FREE dM_dt;
536          M.fixed :=TRUE;          FIX M;
537          dU_dt.fixed :=FALSE;          FREE dU_dt;
538          U.fixed :=TRUE;          FIX U;
539          Qin.fixed :=TRUE;          FIX Qin;
540          IF equilibrated THEN          IF equilibrated THEN
541              state.P.fixed       :=FALSE;              FREE state.P;
542          ELSE          ELSE
543              state.T.fixed       :=FALSE;              FREE state.T;
544          END IF;          END IF;
545          vol_liq.fixed :=FALSE;          FREE vol_liq;
546          tauv.fixed :=TRUE;          FIX tauv;
547          taul.fixed :=TRUE;          FIX taul;
548          cmo_ratio.fixed :=FALSE;          FREE cmo_ratio;
549          ELSE          ELSE
550          dMi_dt[cd.other_components] :=0 {mole/sec};          dMi_dt[cd.other_components] :=0 {mole/sec};
551          dM_dt :=0 {mole/sec};          dM_dt :=0 {mole/sec};
# Line 560  END NOTES; Line 560  END NOTES;
560      END specify;      END specify;
561    
562      METHOD reset_to_adiabatic;      METHOD reset_to_adiabatic;
563          cmo_ratio.fixed := FALSE;          FREE cmo_ratio;
564          Qin.fixed := TRUE;          FIX Qin;
565          Qin := 0{W};          Qin := 0{W};
566      END reset_to_adiabatic;      END reset_to_adiabatic;
567    
# Line 951  END NOTES; Line 951  END NOTES;
951    
952      METHOD seqmod;      METHOD seqmod;
953          RUN state.specify;          RUN state.specify;
954          state.y[cd.components].fixed := FALSE;          FREE state.y[cd.components];
955          IF equilibrated THEN          IF equilibrated THEN
956          state.T.fixed :=FALSE;          FREE state.T;
957          ELSE          ELSE
958          state.T.fixed :=TRUE;          FIX state.T;
959          END IF;          END IF;
960          state.phase_fraction['vapor'].fixed :=FALSE;          FREE state.phase_fraction['vapor'];
961          q.fixed :=TRUE;          FIX q;
962          dMi_dt[cd.other_components].fixed :=TRUE;          FIX dMi_dt[cd.other_components];
963          Mi[cd.other_components].fixed :=FALSE;          FREE Mi[cd.other_components];
964          dM_dt.fixed :=TRUE;          FIX dM_dt;
965          M.fixed :=FALSE;          FREE M;
966          dU_dt.fixed :=TRUE;          FIX dU_dt;
967          U.fixed :=FALSE;          FREE U;
968          Vol.fixed :=TRUE;          FIX Vol;
969          vol_liq.fixed :=TRUE;          FIX vol_liq;
970          tauv.fixed :=FALSE;          FREE tauv;
971          taul.fixed :=FALSE;          FREE taul;
972          IF dynamic THEN          IF dynamic THEN
973          dMi_dt[cd.other_components].fixed :=FALSE;          FREE dMi_dt[cd.other_components];
974          Mi[cd.other_components].fixed :=TRUE;          FIX Mi[cd.other_components];
975          dM_dt.fixed :=FALSE;          FREE dM_dt;
976          M.fixed :=TRUE;          FIX M;
977          dU_dt.fixed :=FALSE;          FREE dU_dt;
978          U.fixed :=TRUE;          FIX U;
979          Qin.fixed :=TRUE;          FIX Qin;
980          IF equilibrated THEN          IF equilibrated THEN
981              state.P.fixed       :=FALSE;              FREE state.P;
982          ELSE          ELSE
983              state.T.fixed       :=FALSE;              FREE state.T;
984          END IF;          END IF;
985          vol_liq.fixed :=FALSE;          FREE vol_liq;
986          tauv.fixed :=TRUE;          FIX tauv;
987          taul.fixed :=TRUE;          FIX taul;
988          q.fixed :=FALSE;          FREE q;
989          ELSE          ELSE
990          dMi_dt[cd.other_components] :=0 {mole/sec};          dMi_dt[cd.other_components] :=0 {mole/sec};
991          dM_dt :=0 {mole/sec};          dM_dt :=0 {mole/sec};
# Line 1002  END NOTES; Line 1002  END NOTES;
1002    
1003      METHOD reset_to_adiabatic;      METHOD reset_to_adiabatic;
1004          RUN seqmod;          RUN seqmod;
1005          q.fixed := FALSE;          FREE q;
1006          Qin.fixed := TRUE;          FIX Qin;
1007          Qin := 0{W};          Qin := 0{W};
1008      END reset_to_adiabatic;      END reset_to_adiabatic;
1009    
# Line 1245  END NOTES; Line 1245  END NOTES;
1245    
1246      METHOD seqmod;      METHOD seqmod;
1247          RUN state.specify;          RUN state.specify;
1248          state.y[cd.components].fixed := FALSE;          FREE state.y[cd.components];
1249          state.phase_fraction['vapor'].fixed :=FALSE;          FREE state.phase_fraction['vapor'];
1250          IF equilibrated THEN          IF equilibrated THEN
1251          state.T.fixed :=FALSE;          FREE state.T;
1252          ELSE          ELSE
1253          state.T.fixed :=TRUE;          FIX state.T;
1254          END IF;          END IF;
1255          RUN split.seqmod;          RUN split.seqmod;
1256          split.split[1].fixed :=FALSE;          FREE split.split[1];
1257          liqout.flow.fixed :=FALSE;          FREE liqout.flow;
1258          vap_distillate.flow.fixed :=FALSE;          FREE vap_distillate.flow;
1259          liq_distillate.flow.fixed :=FALSE;          FREE liq_distillate.flow;
1260          IF partial THEN          IF partial THEN
1261          liq_distillate.flow.fixed :=TRUE;          FIX liq_distillate.flow;
1262          ELSE          ELSE
1263          vap_distillate.flow.fixed :=TRUE;          FIX vap_distillate.flow;
1264          END IF;          END IF;
1265          reflux_ratio.fixed :=TRUE;          FIX reflux_ratio;
1266          dMi_dt[cd.other_components].fixed :=TRUE;          FIX dMi_dt[cd.other_components];
1267          Mi[cd.other_components].fixed :=FALSE;          FREE Mi[cd.other_components];
1268          dM_dt.fixed :=TRUE;          FIX dM_dt;
1269          M.fixed :=FALSE;          FREE M;
1270          dU_dt.fixed :=TRUE;          FIX dU_dt;
1271          U.fixed :=FALSE;          FREE U;
1272          Qin.fixed :=FALSE;          FREE Qin;
1273          Vol.fixed :=TRUE;          FIX Vol;
1274          vol_liq.fixed :=TRUE;          FIX vol_liq;
1275          IF dynamic THEN          IF dynamic THEN
1276          dMi_dt[cd.other_components].fixed :=FALSE;          FREE dMi_dt[cd.other_components];
1277          Mi[cd.other_components].fixed :=TRUE;          FIX Mi[cd.other_components];
1278          dM_dt.fixed :=FALSE;          FREE dM_dt;
1279          M.fixed :=TRUE;          FIX M;
1280          dU_dt.fixed :=FALSE;          FREE dU_dt;
1281          U.fixed :=TRUE;          FIX U;
1282          Qin.fixed :=TRUE;          FIX Qin;
1283          IF equilibrated THEN          IF equilibrated THEN
1284              state.P.fixed       :=FALSE;              FREE state.P;
1285          ELSE          ELSE
1286              state.T.fixed       :=FALSE;              FREE state.T;
1287          END IF;          END IF;
1288          vol_liq.fixed   :=FALSE;          FREE vol_liq;
1289          reflux_ratio.fixed :=FALSE;          FREE reflux_ratio;
1290          state.phase_fraction['vapor'].fixed :=FALSE;          FREE state.phase_fraction['vapor'];
1291          liqout.flow.fixed :=TRUE;          FIX liqout.flow;
1292  (*      liqout.Details.state.slack_PhaseDisappearance['liquid1'].fixed :=TRUE;*)  (*      FIX liqout.Details.state.slack_PhaseDisappearance['liquid1'];*)
1293          vap_distillate.flow.fixed :=TRUE;          FIX vap_distillate.flow;
1294          liq_distillate.flow.fixed :=TRUE;          FIX liq_distillate.flow;
1295          ELSE          ELSE
1296          dMi_dt[cd.other_components] :=0 {mole/sec};          dMi_dt[cd.other_components] :=0 {mole/sec};
1297          dM_dt :=0 {mole/sec};          dM_dt :=0 {mole/sec};
# Line 1542  END NOTES; Line 1542  END NOTES;
1542    
1543      METHOD seqmod;      METHOD seqmod;
1544          RUN state.specify;          RUN state.specify;
1545          state.y[cd.components].fixed := FALSE;          FREE state.y[cd.components];
1546          state.phase_fraction['vapor'].fixed :=FALSE;          FREE state.phase_fraction['vapor'];
1547          IF equilibrated THEN          IF equilibrated THEN
1548          state.T.fixed :=FALSE;          FREE state.T;
1549          ELSE          ELSE
1550          state.T.fixed :=TRUE;          FIX state.T;
1551          END IF;          END IF;
1552          RUN split.seqmod;          RUN split.seqmod;
1553          split.split[1].fixed :=FALSE;          FREE split.split[1];
1554          vapout.flow.fixed :=FALSE;          FREE vapout.flow;
1555          vap_bottom.flow.fixed :=FALSE;          FREE vap_bottom.flow;
1556          liq_bottom.flow.fixed :=FALSE;          FREE liq_bottom.flow;
1557          IF partial THEN          IF partial THEN
1558          vap_bottom.flow.fixed :=TRUE;          FIX vap_bottom.flow;
1559          ELSE          ELSE
1560          liq_bottom.flow.fixed :=TRUE;          FIX liq_bottom.flow;
1561          END IF;          END IF;
1562          reboil_ratio.fixed :=TRUE;          FIX reboil_ratio;
1563          dMi_dt[cd.other_components].fixed :=TRUE;          FIX dMi_dt[cd.other_components];
1564          Mi[cd.other_components].fixed :=FALSE;          FREE Mi[cd.other_components];
1565          dM_dt.fixed :=TRUE;          FIX dM_dt;
1566          M.fixed :=FALSE;          FREE M;
1567          dU_dt.fixed :=TRUE;          FIX dU_dt;
1568          U.fixed :=FALSE;          FREE U;
1569          Qin.fixed :=FALSE;          FREE Qin;
1570          Vol.fixed :=TRUE;          FIX Vol;
1571          vol_liq.fixed :=TRUE;          FIX vol_liq;
1572          IF dynamic THEN          IF dynamic THEN
1573          dMi_dt[cd.other_components].fixed :=FALSE;          FREE dMi_dt[cd.other_components];
1574          Mi[cd.other_components].fixed :=TRUE;          FIX Mi[cd.other_components];
1575          dM_dt.fixed :=FALSE;          FREE dM_dt;
1576          M.fixed :=TRUE;          FIX M;
1577          dU_dt.fixed :=FALSE;          FREE dU_dt;
1578          U.fixed :=TRUE;          FIX U;
1579          Qin.fixed :=TRUE;          FIX Qin;
1580          IF equilibrated THEN          IF equilibrated THEN
1581              state.P.fixed       :=FALSE;              FREE state.P;
1582          ELSE          ELSE
1583              state.T.fixed       :=FALSE;              FREE state.T;
1584          END IF;          END IF;
1585          vol_liq.fixed   :=FALSE;          FREE vol_liq;
1586          reboil_ratio.fixed :=FALSE;          FREE reboil_ratio;
1587          state.phase_fraction['vapor'].fixed :=FALSE;          FREE state.phase_fraction['vapor'];
1588          vapout.flow.fixed :=TRUE;          FIX vapout.flow;
1589          vap_bottom.flow.fixed :=TRUE;          FIX vap_bottom.flow;
1590          liq_bottom.flow.fixed :=TRUE;          FIX liq_bottom.flow;
1591          ELSE          ELSE
1592          dMi_dt[cd.other_components] :=0 {mole/sec};          dMi_dt[cd.other_components] :=0 {mole/sec};
1593          dM_dt :=0 {mole/sec};          dM_dt :=0 {mole/sec};

Legend:
Removed from v.1  
changed lines
  Added in v.576

john.pye@anu.edu.au
ViewVC Help
Powered by ViewVC 1.1.22