| 1 | using "stage_separators/batch_dist"; |
|---|
| 2 | using "controllers/PIDs"; |
|---|
| 3 | |
|---|
| 4 | FlowSheet TesteBatch_PID |
|---|
| 5 | |
|---|
| 6 | PARAMETERS |
|---|
| 7 | PP as CalcObject(Brief="Physical Properties",File="vrpp"); |
|---|
| 8 | NComp as Integer; |
|---|
| 9 | P_max as pressure(Brief="Maximum Pressure on distillator"); |
|---|
| 10 | P_min as pressure(Brief="Minimum Pressure on distillator"); |
|---|
| 11 | |
|---|
| 12 | VARIABLES |
|---|
| 13 | P_ad as fraction(Brief="Dimensionless pressure"); |
|---|
| 14 | P_set as pressure(Brief="Minimum Pressure on distillator"); |
|---|
| 15 | |
|---|
| 16 | SET |
|---|
| 17 | PP.Components = [ "isobutane", "n-pentane" ]; |
|---|
| 18 | PP.LiquidModel = "PR"; |
|---|
| 19 | PP.VapourModel = "PR"; |
|---|
| 20 | NComp = PP.NumberOfComponents; |
|---|
| 21 | |
|---|
| 22 | DEVICES |
|---|
| 23 | batch as Diff_Dist; |
|---|
| 24 | reflux as stream_therm; |
|---|
| 25 | feed as stream_therm; |
|---|
| 26 | pidP as PID_Ideal_AWBT; |
|---|
| 27 | |
|---|
| 28 | CONNECTIONS |
|---|
| 29 | reflux to batch.InletL; |
|---|
| 30 | feed to batch.Inlet; |
|---|
| 31 | |
|---|
| 32 | EQUATIONS |
|---|
| 33 | if batch.Level > 1E-3 then |
|---|
| 34 | batch.Q = 3.7743e6 * "kJ/h"; |
|---|
| 35 | else |
|---|
| 36 | batch.Q = 0 * "kJ/h"; |
|---|
| 37 | end |
|---|
| 38 | |
|---|
| 39 | "Pressure Controller" |
|---|
| 40 | pidP.Parameters.tau = 0*"s"; |
|---|
| 41 | pidP.Parameters.tauSet = 0*"s"; |
|---|
| 42 | pidP.Parameters.alpha = 0.3; |
|---|
| 43 | pidP.Parameters.bias = 0; |
|---|
| 44 | pidP.Parameters.gamma = 1; |
|---|
| 45 | pidP.Parameters.beta = 1; |
|---|
| 46 | pidP.Options.action = -1; |
|---|
| 47 | pidP.Options.clip = 1; |
|---|
| 48 | pidP.Options.autoMan = 0; |
|---|
| 49 | pidP.Parameters.intTime = 0.1*"min"; |
|---|
| 50 | pidP.Parameters.gain = 1; |
|---|
| 51 | pidP.Parameters.derivTime = 0.02*"min"; |
|---|
| 52 | pidP.Ports.setPoint = (P_set-P_min)/(P_max-P_min); |
|---|
| 53 | pidP.Ports.input = P_ad; |
|---|
| 54 | P_ad = (batch.P-P_min)/(P_max-P_min); |
|---|
| 55 | batch.OutletV.F = 500 * "kmol/h" * pidP.Ports.output; |
|---|
| 56 | |
|---|
| 57 | SPECIFY |
|---|
| 58 | reflux.F = 0 * "kmol/h"; |
|---|
| 59 | reflux.T = 328 * "K"; |
|---|
| 60 | reflux.P = 180 * "kPa"; |
|---|
| 61 | reflux.z = [0.5, 0.5]; |
|---|
| 62 | reflux.v = 0; |
|---|
| 63 | |
|---|
| 64 | feed.F = 0 * "kmol/h"; |
|---|
| 65 | feed.T = 328 * "K"; |
|---|
| 66 | feed.P = 180 * "kPa"; |
|---|
| 67 | feed.z = [0.5, 0.5]; |
|---|
| 68 | feed.v = 0; |
|---|
| 69 | |
|---|
| 70 | P_set = 300 * "kPa"; |
|---|
| 71 | |
|---|
| 72 | SET |
|---|
| 73 | batch.V = 3 * "m^3"; |
|---|
| 74 | batch.Across = 1 * "m^2"; |
|---|
| 75 | P_max = 500 * "kPa"; |
|---|
| 76 | P_min = 100 * "kPa"; |
|---|
| 77 | |
|---|
| 78 | INITIAL |
|---|
| 79 | batch.T = 298 *"K"; |
|---|
| 80 | batch.Level = 2.5 * "m"; |
|---|
| 81 | batch.x(2) = 0.2; |
|---|
| 82 | |
|---|
| 83 | OPTIONS |
|---|
| 84 | relativeAccuracy = 1e-5; |
|---|
| 85 | |
|---|
| 86 | time = [0:0.01:6]*"min"; |
|---|
| 87 | end |
|---|
