Changeset 125 for branches/newlanguage


Ignore:
Timestamp:
Jan 22, 2007, 5:57:16 PM (15 years ago)
Author:
Rafael de Pelegrini Soares
Message:

Updated some models for the new language proposal

Location:
branches/newlanguage
Files:
9 edited

Legend:

Unmodified
Added
Removed

  • branches/newlanguage/eml/pressure_changers/pump.mso

    r75 r125  
    3333       
    3434        PARAMETERS
    35 ext NComp       as Integer                      (Brief = "Number of chemical components", Lower = 1);
    36 ext PP                  as CalcObject           (Brief = "External Physical Properties");
     35        outer PP        as Plugin                       (Brief = "External Physical Properties");
     36        outer NComp     as Integer                      (Brief = "Number of chemical components", Lower = 1);
    3737        Mw(NComp)       as molweight            (Brief = "Molar Weight");
    3838        Eff             as positive             (Default = 0.72, Brief = "Pump Efficiency");
    … …  
    147147Model pump
    148148        PARAMETERS
    149 ext PP as CalcObject;
    150 ext NComp as Integer;
     149        outer PP as Plugin;
     150        outer NComp as Integer;
    151151       
    152152        VARIABLES
    153153in      Inlet as stream;
    154 out     Outlet as stream_therm;
     154out     Outlet as streamPH;
    155155       
    156156        dP as press_delta (Brief="Pump head");
    … …  
    164164        Outlet.P = Inlet.P + dP;
    165165       
    166         "Temperature"
    167         Outlet.T = Inlet.T;
    168 
    169         "Vapourisation Fraction"
    170         Outlet.v = Inlet.v;
     166        "FIXME: pump potency"
     167        Outlet.h = Inlet.h;
    171168end

  • branches/newlanguage/eml/stage_separators/column.mso

    r102 r125  
    4646Model Section_Column
    4747        PARAMETERS
    48 ext PP as CalcObject;
    49 ext NComp as Integer;
     48        outer PP as Plugin;
     49        outer NComp as Integer;
    5050        NTrays as Integer(Brief="Number of trays", Default=2);
    5151        topdown as Integer(Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1);
    … …  
    8686Model Distillation_kettle_cond
    8787        PARAMETERS
    88 ext PP as CalcObject;
    89 ext NComp as Integer;
     88        outer PP as Plugin;
     89        outer NComp as Integer;
    9090        NTrays as Integer(Brief="Number of trays", Default=2);
    9191        topdown as Integer(Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1);
    … …  
    157157Model Distillation_thermosyphon_subcooling
    158158        PARAMETERS
    159 ext PP as CalcObject;
    160 ext NComp as Integer;
     159        outer PP as Plugin;
     160        outer NComp as Integer;
    161161        NTrays as Integer(Brief="Number of trays", Default=2);
    162162        topdown as Integer(Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1);
    … …  
    233233Model Distillation_thermosyphon_cond
    234234        PARAMETERS
    235 ext PP as CalcObject;
    236 ext NComp as Integer;
     235        outer PP as Plugin;
     236        outer NComp as Integer;
    237237        NTrays as Integer(Brief="Number of trays", Default=2);
    238238        topdown as Integer(Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1);
    … …  
    305305Model Distillation_kettle_subcooling
    306306        PARAMETERS
    307 ext PP as CalcObject;
    308 ext NComp as Integer;
     307        outer PP as Plugin;
     308        outer NComp as Integer;
    309309        NTrays as Integer(Brief="Number of trays", Default=2);
    310310        topdown as Integer(Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1);
    … …  
    372372Model Rectifier
    373373        PARAMETERS
    374 ext PP as CalcObject;
    375 ext NComp as Integer;
     374        outer PP as Plugin;
     375        outer NComp as Integer;
    376376        NTrays as Integer(Brief="Number of trays", Default=2);
    377377        topdown as Integer(Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1);
    … …  
    428428Model Rectifier_subcooling
    429429        PARAMETERS
    430 ext PP as CalcObject;
    431 ext NComp as Integer;
     430        outer PP as Plugin;
     431        outer NComp as Integer;
    432432        NTrays as Integer(Brief="Number of trays", Default=2);
    433433        topdown as Integer(Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1);
    … …  
    483483Model Refluxed_Stripping
    484484        PARAMETERS
    485 ext PP as CalcObject;
    486 ext NComp as Integer;
     485        outer PP as Plugin;
     486        outer NComp as Integer;
    487487        NTrays as Integer(Brief="Number of trays", Default=2);
    488488        topdown as Integer(Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1);
    … …  
    539539Model Refluxed_Stripping_subcooling
    540540        PARAMETERS
    541 ext PP as CalcObject;
    542 ext NComp as Integer;
     541        outer PP as Plugin;
     542        outer NComp as Integer;
    543543        NTrays as Integer(Brief="Number of trays", Default=2);
    544544        topdown as Integer(Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1);
    … …  
    594594Model Refluxed_Absorption
    595595        PARAMETERS
    596 ext PP as CalcObject;
    597 ext NComp as Integer;
     596        outer PP as Plugin;
     597        outer NComp as Integer;
    598598        NTrays as Integer(Brief="Number of trays", Default=2);
    599599        topdown as Integer(Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1);
    … …  
    651651Model Refluxed_Absorption_subcooling
    652652        PARAMETERS
    653 ext PP as CalcObject;
    654 ext NComp as Integer;
     653        outer PP as Plugin;
     654        outer NComp as Integer;
    655655        NTrays as Integer(Brief="Number of trays", Default=2);
    656656        topdown as Integer(Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1);
    … …  
    706706Model Reboiled_Stripping_kettle
    707707        PARAMETERS
    708 ext PP as CalcObject;
    709 ext NComp as Integer;
     708        outer PP as Plugin;
     709        outer NComp as Integer;
    710710        NTrays as Integer(Brief="Number of trays", Default=2);
    711711        topdown as Integer(Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1);
    … …  
    762762Model Reboiled_Stripping_thermosyphon
    763763        PARAMETERS
    764 ext PP as CalcObject;
    765 ext NComp as Integer;
     764        outer PP as Plugin;
     765        outer NComp as Integer;
    766766        NTrays as Integer(Brief="Number of trays", Default=2);
    767767        topdown as Integer(Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1);
    … …  
    820820Model Reboiled_Absorption_kettle
    821821        PARAMETERS
    822 ext PP as CalcObject;
    823 ext NComp as Integer;
     822        outer PP as Plugin;
     823        outer NComp as Integer;
    824824        NTrays as Integer(Brief="Number of trays", Default=2);
    825825        topdown as Integer(Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1);
    … …  
    876876Model Reboiled_Absorption_thermosyphon
    877877        PARAMETERS
    878 ext PP as CalcObject;
    879 ext NComp as Integer;
     878        outer PP as Plugin;
     879        outer NComp as Integer;
    880880        NTrays as Integer(Brief="Number of trays", Default=2);
    881881        topdown as Integer(Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1);
    … …  
    930930Model ReactiveDistillation
    931931        PARAMETERS
    932 ext PP as CalcObject;
    933 ext NComp as Integer;
     932        outer PP as Plugin;
     933        outer NComp as Integer;
    934934        NTrays as Integer(Brief="Number of trays", Default=2);
    935935       

  • branches/newlanguage/eml/stage_separators/condenser.mso

    r72 r125  
    4343Model condenser
    4444        PARAMETERS
    45 ext PP as CalcObject;
    46 ext NComp as Integer;
     45        outer PP as Plugin;
     46        outer NComp as Integer;
    4747        V as volume (Brief="Condenser total volume");
    4848        Across as area (Brief="Cross Section Area of reboiler");
    4949
    5050        VARIABLES
    51 in      InletV as stream; #(Brief="Vapour inlet stream");
    52 out     OutletL as stream_therm; #(Brief="Liquid outlet stream");
    53 out     OutletV as stream_therm; #(Brief="Vapour outlet stream");
     51in      InletV as stream(Brief="Vapour inlet stream");
     52out     OutletL as liquid_stream(Brief="Liquid outlet stream");
     53out     OutletV as vapour_stream(Brief="Vapour outlet stream");
    5454in      Q as heat_rate (Brief="Heat supplied");
    5555
    … …  
    101101        "Level of liquid phase"
    102102        Level = ML*vL/Across;
    103        
    104         "Vapourisation Fraction"
    105         OutletL.v = 0.0;
    106         OutletV.v = 1.0;
    107103end
    108104
    … …  
    113109Model condenserSteady
    114110        PARAMETERS
    115 ext PP as CalcObject;
    116 ext NComp as Integer;
     111        outer PP as Plugin;
     112        outer NComp as Integer;
    117113
    118114        VARIABLES
    119 in      InletV as stream; #(Brief="Vapour inlet stream");
    120 out     OutletL as stream_therm; #(Brief="Liquid outlet stream");
     115in      InletV as stream(Brief="Vapour inlet stream");
     116out     OutletL as liquid_stream(Brief="Liquid outlet stream");
    121117in      Q as heat_rate (Brief="Heat supplied");
    122118        DP as press_delta (Brief="Pressure Drop in the condenser");
    … …  
    142138Model condenserReact
    143139        PARAMETERS
    144 ext PP as CalcObject;
    145 ext NComp as Integer;
     140        outer PP as Plugin;
     141        outer NComp as Integer;
    146142        V as volume (Brief="Condenser total volume");
    147143        Across as area (Brief="Cross Section Area of reboiler");
    … …  
    152148
    153149        VARIABLES
    154 in      InletV as stream;                       #(Brief="Vapour inlet stream");
    155 out     OutletL as stream_therm;        #(Brief="Liquid outlet stream");
    156 out     OutletV as stream_therm;        #(Brief="Vapour outlet stream");
     150in      InletV as stream(Brief="Vapour inlet stream");
     151out     OutletL as liquid_stream(Brief="Liquid outlet stream");
     152out     OutletV as vapour_stream(Brief="Vapour outlet stream");
    157153
    158154        M(NComp) as mol (Brief="Molar Holdup in the tray");
    … …  
    208204        Level = ML*vL/Across;
    209205       
    210         "Vapourisation Fraction"
    211         OutletL.v = 0.0;
    212         OutletV.v = 1.0;
    213        
    214206        "Chemical Equilibrium"
    215207        PP.LiquidFugacityCoefficient(OutletL.T, OutletL.P, OutletL.z)*OutletL.z =
    … …  
    217209
    218210        sum(OutletL.z)=sum(OutletV.z);
    219 
    220211end

  • branches/newlanguage/eml/stage_separators/reboiler.mso

    r72 r125  
    4848Model reboiler
    4949        PARAMETERS
    50 ext PP as CalcObject;
    51 ext NComp as Integer;
     50        outer PP as Plugin;
     51        outer NComp as Integer;
    5252        Across as area (Brief="Cross Section Area of reboiler");
    5353        V as volume (Brief="Total volume of reboiler");
    5454
    5555        VARIABLES
    56 in      Inlet as stream; # (Brief="Feed Stream");
    57 in      InletL as stream; # (Brief="Liquid inlet stream");
    58 out     OutletL as stream_therm; # (Brief="Liquid outlet stream");
    59 out     OutletV as stream_therm; # (Brief="Vapour outlet stream");
     56in      Inlet as stream(Brief="Feed Stream");
     57in      InletL as stream(Brief="Liquid inlet stream");
     58out     OutletL as liquid_stream(Brief="Liquid outlet stream");
     59out     OutletV as vapour_stream(Brief="Vapour outlet stream");
    6060in      Q as heat_rate (Brief="Heat supplied");
    6161
    … …  
    112112        "Level of liquid phase"
    113113        Level = ML*vL/Across;
    114                
    115         "vaporization fraction"
    116         OutletV.v = 1.0;
    117         OutletL.v = 0.0;
    118114end
    119115
    … …  
    123119Model reboilerSteady
    124120        PARAMETERS
    125 ext PP as CalcObject;
    126 ext NComp as Integer;
     121        outer PP as Plugin;
     122        outer NComp as Integer;
    127123        DP as press_delta (Brief="Pressure Drop in the reboiler");
    128124
    129125        VARIABLES
    130 in      InletL as stream; #(Brief="Liquid inlet stream");
    131 out     OutletV as stream_therm; #(Brief="Vapour outlet stream");
     126in      InletL as stream(Brief="Liquid inlet stream");
     127out     OutletV as vapour_stream(Brief="Vapour outlet stream");
    132128in      Q as heat_rate (Brief="Heat supplied");
    133129        vV as volume_mol (Brief="Vapour Molar volume");
    … …  
    150146        "Pressure"
    151147        DP = InletL.P - OutletV.P;
    152 
    153         "Vapourisation Fraction"
    154         OutletV.v = 1.0;
    155148end
    156149
    … …  
    162155Model reboilerSteady_fakeH
    163156        PARAMETERS
    164 ext PP as CalcObject;
    165 ext NComp as Integer;
     157        outer PP as Plugin;
     158        outer NComp as Integer;
    166159        DP as press_delta (Brief="Pressure Drop in the reboiler");
    167160        k as Real (Brief = "Flow Constant", Unit="mol/J");
    168161       
    169162        VARIABLES
    170 in      InletL as stream; #(Brief="Liquid inlet stream");
    171 out     OutletV as stream; #(Brief="Vapour outlet stream");
     163in      InletL as stream(Brief="Liquid inlet stream");
     164out     OutletV as stream(Brief="Vapour outlet stream");
    172165in      Q as heat_rate (Brief="Heat supplied");
    173166
    … …  
    198191Model reboilerReact
    199192        PARAMETERS
    200 ext PP as CalcObject;
    201 ext NComp as Integer;
     193        outer PP as Plugin;
     194        outer NComp as Integer;
    202195        Across as area (Brief="Cross Section Area of reboiler");
    203196        V as volume (Brief="Total volume of reboiler");
    … …  
    208201
    209202        VARIABLES
    210 in      Inlet as stream;                        #(Brief="Feed Stream");
    211 in      InletL as stream;                       #(Brief="Liquid inlet stream");
    212 out     OutletL as stream_therm;        #(Brief="Liquid outlet stream");
    213 out     OutletV as stream_therm;        #(Brief="Vapour outlet stream");
     203in      Inlet as stream(Brief="Feed Stream");
     204in      InletL as stream(Brief="Liquid inlet stream");
     205out     OutletL as liquid_stream(Brief="Liquid outlet stream");
     206out     OutletV as vapour_stream(Brief="Vapour outlet stream");
    214207
    215208        Q as heat_rate (Brief="Heat supplied");
    … …  
    260253        Vol = ML*vL;
    261254       
    262         "vaporization fraction "
    263         OutletV.v = 1.0;
    264         OutletL.v = 0.0;
    265        
    266255        "Mechanical Equilibrium"
    267256        OutletL.P = OutletV.P;
    … …  
    271260       
    272261        "Geometry Constraint"
    273         V = ML*vL + MV*vV;             
     262        V = ML*vL + MV*vV;
    274263
    275264        "Chemical Equilibrium"

  • branches/newlanguage/eml/stage_separators/tank.mso

    r72 r125  
    3939
    4040        PARAMETERS
    41 ext PP as CalcObject;
    42 ext NComp as Integer;
     41        outer PP as Plugin;
     42        outer NComp as Integer;
    4343        Across as area (Brief="Tank cross section area", Default=2);
    4444       
    4545        VARIABLES
    4646in      Inlet as stream;
    47 out     Outlet as stream_therm;
     47out     Outlet as liquid_stream;
    4848
    4949in      Q as heat_rate (Brief="Rate of heat supply");
    … …  
    7474        "Level of liquid phase"
    7575        Level = sum(M)*vL/Across;
    76        
    77         "Vapourisation Fraction"
    78         Outlet.v = Inlet.v;
    7976end
    8077
    … …  
    8784
    8885        PARAMETERS
    89 ext PP as CalcObject;
    90 ext NComp as Integer;
     86        outer PP as Plugin;
     87        outer NComp as Integer;
    9188        radius as length(Brief="Tank radius");
    9289        L as length(Brief="Tank length");
    … …  
    9491        VARIABLES
    9592in      Inlet as stream;
    96 out     Outlet as stream_therm;
     93out     Outlet as liquid_stream;
    9794
    9895in      Q as heat_rate (Brief="Rate of heat supply");
    … …  
    128125        "Level of liquid phase"
    129126        Level = sum(M)*vL/Across;
    130        
    131         "Vapourisation Fraction"
    132         Outlet.v = Inlet.v;
    133127end
    134128
    … …  
    154148
    155149        PARAMETERS
    156 ext PP as CalcObject;
    157 ext NComp as Integer;
     150        outer PP as Plugin;
     151        outer NComp as Integer;
    158152        Across as area (Brief="Tank cross section area", Default=2);
    159153       
    … …  
    161155in      Feed as stream;
    162156in      Inlet as stream;
    163 out     Outlet as stream_therm;
     157out     Outlet as liquid_stream;
    164158
    165159in      Q as heat_rate (Brief="Rate of heat supply");
    … …  
    190184        "Level of liquid phase"
    191185        Level = sum(M)*vL/Across;
    192        
    193         "Vapourisation Fraction"
    194         Outlet.v = Inlet.v;
    195        
    196186end

  • branches/newlanguage/eml/streams.mso

    r124 r125  
    8181end
    8282
     83Model streamPH as stream
     84        PARAMETERS
     85        outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
     86        outer NComp as Integer (Brief = "Number of chemical components", Lower = 1);
     87       
     88        VARIABLES
     89        x(NComp) as fraction(Brief = "Liquid Molar Fraction");
     90        y(NComp) as fraction(Brief = "Vapour Molar Fraction");
     91       
     92        EQUATIONS
     93        "Flash Calculation"
     94        [v, x, y] = PP.FlashPH(P, h, z);
     95        "Enthalpy"
     96        h = (1-v)*PP.LiquidEnthalpy(T, P, x) +
     97                v*PP.VapourEnthalpy(T, P, y);
     98end
     99
    83100Model source
    84101        ATTRIBUTES

  • branches/newlanguage/sample/miscellaneous/sample_arrays.mso

    r83 r125  
    4040
    4141FlowSheet SampleArrays
    42 
    4342        DEVICES
    4443        dev1 as Arrays;
    … …  
    4847
    4948        OPTIONS
    50         time = [0:0.1:2*3.1415];
     49        TimeStep = 0.1;
     50        TimeEnd = 3.1415;
    5151end

  • branches/newlanguage/sample/miscellaneous/sample_if.mso

    r124 r125  
    4545        y = 0;
    4646       
     47        OPTIONS
     48        TimeStep = 0.1;
     49        TimeEnd = 10.0;
     50        TimeUnit = "s";
    4751end
    4852

  • branches/newlanguage/sample/stage_separators/sample_column.mso

    r112 r125  
    3232FlowSheet SectionColumn_Test_with2tray
    3333        PARAMETERS
    34         PP      as CalcObject(Brief="Physical Properties",File="vrpp");
     34        PP      as Plugin(Brief="Physical Properties",File="vrpp");
    3535        NComp   as Integer;
    3636
    … …  
    4343        DEVICES
    4444        sec as Section_Column;
    45         feed as stream_therm;
    46         reb as stream_therm;
    47         cond as stream_therm;
     45        feed as liquid_stream;
     46        reb as vapour_stream;
     47        cond as liquid_stream;
    4848        zero as stream;
    4949       
    … …  
    5959        feed.P = 168.3 * "kPa";
    6060        feed.z = [0.5, 0.5];
    61         feed.v = 0;
    6261
    6362        zero.F = 0 * "kmol/h";
    … …  
    7170        cond.P = 150 * "kPa";
    7271        cond.T = 281.75 * "K";
    73         cond.v = 0.0;
    7472        cond.z = [0.6664, 0.3336];
    7573
    … …  
    7876        reb.T = 328.12 * "K";
    7977        reb.z = [0.001848, 0.9982];
    80         reb.v = 1.0;
    8178       
    8279        sec.trays.Emv = 1;
    … …  
    10299        #relativeAccuracy = 1e-3;
    103100        NLASolver = "sundials";
    104         time = [0:10:1000];
     101        TimeStep = 10;
     102        TimeEnd = 1000;
    105103       
    106104# After running few seconds of transient the steady-state
    … …  
    114112FlowSheet SectionColumn_Test_with8tray
    115113        PARAMETERS
    116         PP      as CalcObject(Brief="Physical Properties",File="vrpp");
     114        PP      as Plugin(Brief="Physical Properties",File="vrpp");
    117115        NComp   as Integer;
    118116
    … …  
    126124        DEVICES
    127125        sec as Section_Column;
    128         feed as stream_therm;
    129         reb as stream_therm;
    130         cond as stream_therm;
     126        feed as liquid_stream;
     127        reb as vapour_stream;
     128        cond as liquid_stream;
    131129        zero as stream;
    132130       
    … …  
    145143        feed.P = 168.3 * "kPa";
    146144        feed.z = [0.5, 0.5];
    147         feed.v = 0;
    148145
    149146        zero.F = 0 * "kmol/h";
    … …  
    157154        cond.P = 150 * "kPa";
    158155        cond.T = 281.75 * "K";
    159         cond.v = 0.0;
    160156        cond.z = [0.6664, 0.3336];
    161157
    … …  
    164160        reb.T = 328.12 * "K";
    165161        reb.z = [0.001848, 0.9982];
    166         reb.v = 1.0;
    167162       
    168163        sec.trays.Emv = 1;
    … …  
    186181       
    187182        OPTIONS
    188         relativeAccuracy = 1e-5;
    189         time = [0:1:100];
     183        RelativeAccuracy = 1e-5;
     184        TimeStep = 1;
     185        TimeEnd = 100;
    190186        #guessFile="SectionColumn_Test_with8tray.rlt";
    191187        #mode = "steady";
    … …  
    195191FlowSheet Distillation_kettle_cond_Test
    196192        PARAMETERS
    197         PP      as CalcObject(Brief="Physical Properties",File="vrpp");
     193        PP      as Plugin(Brief="Physical Properties",File="vrpp");
    198194        NComp   as Integer;
    199195       
    … …  
    212208        DEVICES
    213209        col as Distillation_kettle_cond;
    214         feed as streamTP;
     210        feed as source;
    215211        zero as stream;
    216212       
    … …  
    224220       
    225221        SPECIFY
    226         feed.F = 113.4 * "kmol/h";
    227         feed.T = 291 * "K";
    228         feed.P = 168.3 * "kPa";
    229         feed.z = 1/NComp;
     222        feed.Outlet.F = 113.4 * "kmol/h";
     223        feed.Outlet.T = 291 * "K";
     224        feed.Outlet.P = 168.3 * "kPa";
     225        feed.Outlet.z = 1/NComp;
    230226       
    231227        zero.F = 0 * "kmol/h";
    … …  
    277273
    278274        OPTIONS
    279         relativeAccuracy = 1e-3;
    280         time = [0:0.01:1, 2:50];
     275        RelativeAccuracy = 1e-3;
     276        TimeStep = 0.01;
     277        TimeEnd = 50;
     278        #time = [0:0.01:1, 2:50];
    281279        #guessFile="Distillation_kettle_cond_Test.rlt";
    282280        #mode = "steady";       
    … …  
    285283FlowSheet Column_ctrl
    286284        PARAMETERS
    287         PP      as CalcObject(Brief="Physical Properties",File="vrpp");
     285        PP      as Plugin(Brief="Physical Properties",File="vrpp");
    288286        NComp   as Integer;
    289287       
    … …  
    324322        DEVICES
    325323        col as Distillation_kettle_cond;
    326         feed as streamTP;
     324        feed as source;
    327325        zero as stream;
    328326        TCcond as PIDIncr_Ideal_AW;
    … …  
    421419
    422420        SPECIFY
    423         feed.F = 113.4 * "kmol/h";
    424         feed.T = 291 * "K";
    425         feed.P = 168.3 * "kPa";
    426         feed.z = 1/NComp;
     421        feed.Outlet.F = 113.4 * "kmol/h";
     422        feed.Outlet.T = 291 * "K";
     423        feed.Outlet.P = 168.3 * "kPa";
     424        feed.Outlet.z = 1/NComp;
    427425       
    428426        zero.F = 0 * "kmol/h";
    … …  
    486484
    487485        OPTIONS
    488         relativeAccuracy = 1e-3;
    489         time = [0:0.01:0.1, 0.11:0.01:2]*"h";
     486        RelativeAccuracy = 1e-3;
     487        TimeStep = 0.01;
     488        TimeEnd = 2;
     489        TimeUnit = "h";
     490        #time = [0:0.01:0.1, 0.11:0.01:2]*"h";
    490491        #initialFile = "Column_ctrl.rlt";
    491492        #guessFile = "Column_ctrl.rlt";
Note: See TracChangeset for help on using the changeset viewer.