Changes in / [71:63]


Ignore:
Location:
/mso
Files:
4 deleted
23 edited

Legend:

Unmodified
Added
Removed

  • /mso/eml/heat_exchangers

    • Property svn:ignore deleted

  • /mso/eml/heat_exchangers/HEX_Engine.mso

    r71 r63  
    3838*               Basic_Pdrop     : Hot and Cold Side Pressure Drop
    3939*----------------------------------------------------------------------
    40 *       - Main_DoublePipe : Double Pipe Heat Exchanger Block
    41 *               DoublePipe_HeatTransfer
    42 *               DoublePipe_PressureDrop
    43 *----------------------------------------------------------------------
    4440* Author: Gerson Balbueno Bicca
    4541* $Id$
    … …  
    281277Cold    as Basic_Pdrop;# Cold Stream
    282278end     
    283 
    284 Model DoublePipe_HeatTransfer
    285 #=====================================================================
    286 #       Heat Transfer
    287 #===================================================================== 
    288 PARAMETERS
    289 As              as area                         (Brief="Cross Sectional Area for Flow",Default=0.05,Lower=1e-8);
    290 Dh      as length                       (Brief="Hydraulic Diameter of Pipe for Heat Transfer",Lower=1e-8);
    291 
    292 VARIABLES
    293 Re              as positive                     (Brief="Reynolds Number",Default=100,Lower=1);
    294 hcoeff  as heat_trans_coeff (Brief="Film Coefficient",Default=1,Lower=1e-12, Upper=1e6);
    295 PR              as positive                     (Brief="Prandtl Number",Default=0.5,Lower=1e-8);
    296 Phi     as positive                     (Brief="Phi Correction",Default=1,Lower=1e-3);
    297 Vmean   as velocity                     (Brief="Tube Velocity",Lower=1e-8);
    298 end
    299 
    300 Model DoublePipe_PressureDrop
    301 #=====================================================================
    302 #       Pressure Drop
    303 #===================================================================== 
    304 PARAMETERS
    305 Dh      as length               (Brief="Hydraulic Diameter of Pipe for Pressure Drop",Lower=1e-6);
    306 
    307 VARIABLES
    308 Pdrop   as press_delta  (Brief="Pressure Drop",Default=0.01, Lower=1e-10);
    309 fi      as fricfactor   (Brief="Friction Factor", Default=0.05, Lower=1e-10, Upper=2000);
    310 Re              as positive             (Brief="Reynolds Number",Default=100,Lower=1);
    311 end     
    312 
    313 Model Main_DoublePipe
    314 VARIABLES
    315 HeatTransfer as DoublePipe_HeatTransfer;
    316 PressureDrop as DoublePipe_PressureDrop;
    317 end

  • /mso/eml/heat_exchangers/HeatExchangerDetailed.mso

    r71 r63  
    165165        then
    166166
    167 "Hot Stream Average Heat Capacity"
     167"Heat Capacity Hot Stream"
    168168        Properties.Hot.Average.Cp       =               PP.LiquidCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
    169 
    170 "Hot Stream Inlet Heat Capacity"
    171169        Properties.Hot.Inlet.Cp         =               PP.LiquidCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
    172 
    173 "Hot Stream Outlet Heat Capacity"
    174170        Properties.Hot.Outlet.Cp        =               PP.LiquidCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    175171
    176 "Hot Stream Average Mass Density"
     172"Mass Density Hot Stream"
    177173        Properties.Hot.Average.rho      =               PP.LiquidDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
    178 
    179 "Hot Stream Inlet Mass Density"
    180174        Properties.Hot.Inlet.rho        =               PP.LiquidDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
    181 
    182 "Hot Stream Outlet Mass Density"       
    183175        Properties.Hot.Outlet.rho       =               PP.LiquidDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    184176
    185 "Hot Stream Average Viscosity"
     177"Viscosity Hot Stream"
    186178        Properties.Hot.Average.Mu       =               PP.LiquidViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);     
    187 
    188 "Hot Stream Inlet Viscosity"
    189179        Properties.Hot.Inlet.Mu         =               PP.LiquidViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);       
    190 
    191 "Hot Stream Outlet Viscosity"
    192180        Properties.Hot.Outlet.Mu        =               PP.LiquidViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);     
    193181
    194 "Hot Stream Average Conductivity"
    195         Properties.Hot.Average.K        =               PP.LiquidThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);   
    196 
    197 "Hot Stream Inlet Conductivity"
    198         Properties.Hot.Inlet.K          =               PP.LiquidThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);     
    199 
    200 "Hot Stream Outlet Conductivity"
    201         Properties.Hot.Outlet.K         =               PP.LiquidThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);   
    202 
    203 "Hot Stream Heat Capacity at Wall Temperature"
    204         Properties.Hot.Wall.Cp          =               PP.LiquidCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
    205 
    206 "Hot Stream Viscosity  at Wall Temperature"
    207         Properties.Hot.Wall.Mu          =               PP.LiquidViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);     
    208 
    209 "Hot Stream Conductivity at Wall Temperature"
    210         Properties.Hot.Wall.K           =               PP.LiquidThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);   
    211        
    212 
    213         else
    214 
    215 "Hot Stream Average Heat Capacity"
     182"Conductivity Hot Stream"
     183        Properties.Hot.Average.K        =       PP.LiquidThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);   
     184        Properties.Hot.Inlet.K  =               PP.LiquidThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);     
     185        Properties.Hot.Outlet.K         =       PP.LiquidThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);   
     186
     187"Heat Capacity Hot Stream"
     188        Properties.Hot.Wall.Cp          =       PP.LiquidCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
     189
     190"Viscosity Hot Stream"
     191        Properties.Hot.Wall.Mu          =       PP.LiquidViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);     
     192
     193"Conductivity Hot Stream"
     194        Properties.Hot.Wall.K           =       PP.LiquidThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);   
     195       
     196
     197        else
     198
     199"Heat Capacity Hot Stream"
    216200        Properties.Hot.Average.Cp       =               PP.VapourCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
    217 
    218 "Hot Stream Inlet Heat Capacity"
    219201        Properties.Hot.Inlet.Cp         =               PP.VapourCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
    220 
    221 "Hot Stream Outlet Heat Capacity"
    222202        Properties.Hot.Outlet.Cp        =               PP.VapourCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    223203
    224 "Hot Stream Average Mass Density"
     204"Mass Density Hot Stream"
    225205        Properties.Hot.Average.rho      =               PP.VapourDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
    226 
    227 "Hot Stream Inlet Mass Density"
    228206        Properties.Hot.Inlet.rho        =               PP.VapourDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
    229 
    230 "Hot Stream Outlet Mass Density"
    231207        Properties.Hot.Outlet.rho       =               PP.VapourDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    232208
    233 "Hot Stream Average Viscosity"
     209"Viscosity Hot Stream"
    234210        Properties.Hot.Average.Mu       =               PP.VapourViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
    235 
    236 "Hot Stream Inlet Viscosity"
    237211        Properties.Hot.Inlet.Mu         =               PP.VapourViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
    238 
    239 "Hot Stream Outlet Viscosity"
    240212        Properties.Hot.Outlet.Mu        =               PP.VapourViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    241213
    242 "Hot Stream Average Conductivity"
    243         Properties.Hot.Average.K        =               PP.VapourThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);   
    244 
    245 "Hot Stream Inlet Conductivity"
    246         Properties.Hot.Inlet.K          =               PP.VapourThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);     
    247        
    248 "Hot Stream Outlet Conductivity"
    249         Properties.Hot.Outlet.K         =               PP.VapourThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);   
    250 
    251 "Hot Stream Heat Capacity at Wall Temperature"
    252         Properties.Hot.Wall.Cp          =               PP.VapourCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
    253 
    254 "Hot Stream Viscosity at Wall Temperature"
    255         Properties.Hot.Wall.Mu          =               PP.VapourViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
    256 
    257 "Hot Stream Conductivity at Wall Temperature"
    258         Properties.Hot.Wall.K           =               PP.VapourThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);   
     214"Conductivity Hot Stream"
     215        Properties.Hot.Average.K        =       PP.VapourThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);   
     216        Properties.Hot.Inlet.K  =               PP.VapourThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);     
     217        Properties.Hot.Outlet.K         =       PP.VapourThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);   
     218
     219"Heat Capacity Hot Stream"
     220        Properties.Hot.Wall.Cp  =       PP.VapourCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
     221
     222"Viscosity Hot Stream"
     223        Properties.Hot.Wall.Mu  =       PP.VapourViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
     224
     225"Conductivity Hot Stream"
     226        Properties.Hot.Wall.K   =       PP.VapourThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);   
    259227
    260228
    … …  
    455423
    456424"Shell Side Phi correction"
    457         Shell.HeatTransfer.Phi  = HE.PhiCorrection(Properties.Hot.Average.Mu,Properties.Hot.Wall.Mu);
     425        Shell.HeatTransfer.Phi = HE.PhiCorrection(Properties.Hot.Average.Mu,Properties.Hot.Wall.Mu);
    458426       
    459427"Tube Side Phi correction"
    … …  
    463431       
    464432"Shell Side Phi correction"
    465         Shell.HeatTransfer.Phi  = HE.PhiCorrection(Properties.Cold.Average.Mu,Properties.Cold.Wall.Mu);
     433        Shell.HeatTransfer.Phi = HE.PhiCorrection(Properties.Cold.Average.Mu,Properties.Cold.Wall.Mu);
    466434       
    467435"Tube Side Phi correction"
    … …  
    481449       
    482450"Hot Wall Temperature"
    483         Properties.Hot.Wall.Twall       = (Properties.Hot.Average.T+Properties.Cold.Average.T)/2;
     451        Properties.Hot.Wall.Twall  = (Properties.Hot.Average.T+Properties.Cold.Average.T)/2;
    484452
    485453"ColdWall Temperature"
    … …  
    487455
    488456"Tube Side Velocity"
    489         Tubes.HeatTransfer.Vtube        = Properties.Cold.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Properties.Cold.Average.rho*Ntt);
     457        Tubes.HeatTransfer.Vtube  = Properties.Cold.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Properties.Cold.Average.rho*Ntt);
    490458
    491459"Tube Side Reynolds Number"
    492         Tubes.HeatTransfer.Re           = (Properties.Cold.Average.rho*Tubes.HeatTransfer.Vtube*Ditube)/Properties.Cold.Average.Mu;
     460        Tubes.HeatTransfer.Re     =     (Properties.Cold.Average.rho*Tubes.HeatTransfer.Vtube*Ditube)/Properties.Cold.Average.Mu;
    493461       
    494462"Tube Side Prandtl Number"
    495         Tubes.HeatTransfer.PR           = ((Properties.Cold.Average.Cp/Properties.Cold.Average.Mw)*Properties.Cold.Average.Mu)/Properties.Cold.Average.K;
     463        Tubes.HeatTransfer.PR = ((Properties.Cold.Average.Cp/Properties.Cold.Average.Mw)*Properties.Cold.Average.Mu)/Properties.Cold.Average.K;
    496464
    497465"Tube Side Prandtl Number at Wall Temperature"
    498         Tubes.HeatTransfer.PRw          = ((Properties.Cold.Wall.Cp/Properties.Cold.Average.Mw)*Properties.Cold.Wall.Mu)/Properties.Cold.Wall.K;
     466        Tubes.HeatTransfer.PRw = ((Properties.Cold.Wall.Cp/Properties.Cold.Average.Mw)*Properties.Cold.Wall.Mu)/Properties.Cold.Wall.K;
    499467
    500468"Tube Side Film Coefficient"
    501         Tubes.HeatTransfer.htube        = HE.TubeFilmCoeff(Tubes.HeatTransfer.Re,Tubes.HeatTransfer.PR,Properties.Cold.Average.K)*Tubes.HeatTransfer.Phi;
     469        Tubes.HeatTransfer.htube= HE.TubeFilmCoeff(Tubes.HeatTransfer.Re,Tubes.HeatTransfer.PR,Properties.Cold.Average.K)*Tubes.HeatTransfer.Phi;
    502470
    503471"Shell Side Prandtl Number"
    504         Shell.HeatTransfer.PR           = ((Properties.Hot.Average.Cp/Properties.Hot.Average.Mw)*Properties.Hot.Average.Mu)/Properties.Hot.Average.K;
     472        Shell.HeatTransfer.PR = ((Properties.Hot.Average.Cp/Properties.Hot.Average.Mw)*Properties.Hot.Average.Mu)/Properties.Hot.Average.K;
    505473
    506474"Shell Side Prandtl Number at Wall Temperature"
    507         Shell.HeatTransfer.PRw          = ((Properties.Hot.Wall.Cp/Properties.Hot.Average.Mw)*Properties.Hot.Wall.Mu)/Properties.Hot.Wall.K;
     475        Shell.HeatTransfer.PRw = ((Properties.Hot.Wall.Cp/Properties.Hot.Average.Mw)*Properties.Hot.Wall.Mu)/Properties.Hot.Wall.K;
    508476
    509477"Tube Side Pressure Drop"
    … …  
    529497
    530498"Shell Pressure Drop Outlet Nozzle"
    531         Shell.PressureDrop.Pdnozzle_out = HE.DeltaPshellNozzleOut(Properties.Hot.Outlet.rho,Properties.Hot.Inlet.Fw);
     499        Shell.PressureDrop.Pdnozzle_out =HE.DeltaPshellNozzleOut(Properties.Hot.Outlet.rho,Properties.Hot.Inlet.Fw);
    532500
    533501"Velocity Shell Side Outlet Nozzle"
    … …  
    653621       
    654622VARIABLES
    655 
    656623DT0     as temp_delta   (Brief="Temperature Difference at Inlet",Lower=1);
    657624DTL             as temp_delta   (Brief="Temperature Difference at Outlet",Lower=1);
    … …  
    13281295VARIABLES
    13291296
    1330 LMTD(Nshell)            as temp_delta   (Brief="Logarithmic Mean Temperature Difference",Lower=5);
    1331 Fc(Nshell)                      as positive             (Brief="LMTD Correction Factor",Lower=0.5);
    1332 MTD(Nshell)                     as temp_delta   (Brief="Mean Temperature Difference",Lower=5);
     1297LMTD(Nshell)            as temp_delta   (Brief="Logarithmic Mean Temperature Difference",Lower=10);
     1298Fc(Nshell)                      as positive             (Brief="LMTD Correction Factor",Lower=0.75);
     1299MTD(Nshell)                     as temp_delta   (Brief="Mean Temperature Difference");
    13331300Unity(Nshell)           as HeatExchangerDetailed_Basic;
    13341301
    … …  
    15481515       
    15491516"Overall Heat Transfer Coefficient"
    1550         Unity(i).Details.U*(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)))=1;
     1517        Unity(i).Details.U=1/(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)));
    15511518
    15521519"Exchange Surface Area"
    … …  
    16061573
    16071574"Exchange Surface Area"
    1608 #       Unity(i).Details.Q   = Unity(i).Details.U*Unity(i).Details.A*Fc(i)*LMTD(i);
    1609         Unity(i).Details.Q   = Unity(i).Details.U*Pi*Dotube*Ntt*Ltube*Fc(i)*LMTD(i);
     1575        Unity(i).Details.Q   = Unity(i).Details.U*Unity(i).Details.A*Fc(i)*LMTD(i);
    16101576       
    16111577"Mean Temperature Difference"   

  • /mso/eml/heat_exchangers/HeatExchangerDiscretized.mso

    r71 r63  
    6060        then   
    6161       
    62 "Cold Stream Average Heat Capacity"
     62"Heat Capacity Cold Stream"
    6363        Properties.Cold.Average.Cp      =       PP.LiquidCp(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
    64 
    65 "Cold Stream Inlet Heat Capacity"
    6664        Properties.Cold.Inlet.Cp        =       PP.LiquidCp(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
    67 
    68 "Cold Stream Outlet Heat Capacity"
    6965        Properties.Cold.Outlet.Cp       =       PP.LiquidCp(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    7066
    71 "Cold Stream Average Mass Density"
     67"Mass Density Cold Stream"
    7268        Properties.Cold.Average.rho =   PP.LiquidDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
    73 
    74 "Cold Stream Inlet Mass Density"
    7569        Properties.Cold.Inlet.rho       =       PP.LiquidDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
    76 
    77 "Cold Stream Outlet Mass Density"
    7870        Properties.Cold.Outlet.rho      =       PP.LiquidDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    7971
    80 "Cold Stream Average Viscosity"
    81         Properties.Cold.Average.Mu      =       PP.LiquidViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
    82 
    83 "Cold Stream inlet Viscosity"
    84         Properties.Cold.Inlet.Mu        =       PP.LiquidViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
    85        
    86 "Cold Stream Outlet Viscosity"
    87         Properties.Cold.Outlet.Mu       =       PP.LiquidViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    88 
    89 "Cold Stream Average Conductivity"
    90         Properties.Cold.Average.K       =       PP.LiquidThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
    91 
    92 "Cold Stream Inlet Conductivity"       
    93         Properties.Cold.Inlet.K         =       PP.LiquidThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
    94 
    95 "Cold Stream Outlet Conductivity"
    96         Properties.Cold.Outlet.K        =       PP.LiquidThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    97 
    98 "Cold Stream Heat Capacity at Wall Temperature"
     72"Viscosity Cold Stream"
     73        Properties.Cold.Average.Mu =    PP.LiquidViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
     74        Properties.Cold.Inlet.Mu =              PP.LiquidViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
     75        Properties.Cold.Outlet.Mu =     PP.LiquidViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
     76
     77"Conductivity Cold Stream"
     78        Properties.Cold.Average.K =     PP.LiquidThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
     79        Properties.Cold.Inlet.K =               PP.LiquidThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
     80        Properties.Cold.Outlet.K =              PP.LiquidThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
     81
     82"Heat Capacity Cold Stream"
    9983        Properties.Cold.Wall.Cp         =       PP.LiquidCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
    10084       
    101 "Cold Stream Viscosity at Wall Temperature"
    102         Properties.Cold.Wall.Mu         =       PP.LiquidViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
    103 
    104 "Cold Stream Conductivity at Wall Temperature"
    105         Properties.Cold.Wall.K          =       PP.LiquidThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
     85"Viscosity Cold Stream"
     86        Properties.Cold.Wall.Mu =       PP.LiquidViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
     87
     88"Conductivity Cold Stream"
     89        Properties.Cold.Wall.K =        PP.LiquidThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
    10690
    10791
    10892        else
    10993
    110 "Cold Stream Average Heat Capacity"
     94"Heat Capacity Cold Stream"
    11195        Properties.Cold.Average.Cp      =       PP.VapourCp(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
    112 
    113 "Cold Stream Inlet Heat Capacity"       
    11496        Properties.Cold.Inlet.Cp        =       PP.VapourCp(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
    115 
    116 "Cold Stream Outlet Heat Capacity"     
    11797        Properties.Cold.Outlet.Cp       =       PP.VapourCp(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    11898
    119 "Cold Stream Average Mass Density"
    120         Properties.Cold.Average.rho =   PP.VapourDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
    121 
    122 "Cold Stream Inlet Mass Density"
     99"Mass Density Cold Stream"
     100        Properties.Cold.Average.rho     =       PP.VapourDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
    123101        Properties.Cold.Inlet.rho       =       PP.VapourDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
    124 
    125 "Cold Stream Outlet Mass Density"       
    126102        Properties.Cold.Outlet.rho      =       PP.VapourDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    127103
    128 "Cold Stream Average Viscosity "
    129         Properties.Cold.Average.Mu      =       PP.VapourViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
    130 
    131 "Cold Stream Inlet Viscosity " 
    132         Properties.Cold.Inlet.Mu        =       PP.VapourViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
    133 
    134 "Cold Stream Outlet Viscosity "
    135         Properties.Cold.Outlet.Mu       =       PP.VapourViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    136 
    137 "Cold Stream Average Conductivity "
    138         Properties.Cold.Average.K       =       PP.VapourThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
    139 
    140 "Cold Stream Inlet Conductivity "
    141         Properties.Cold.Inlet.K         =       PP.VapourThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
    142 
    143 "Cold Stream Outlet Conductivity "
    144         Properties.Cold.Outlet.K        =       PP.VapourThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    145        
    146 "Cold Stream Heat Capacity at Wall Temperature"
     104"Viscosity Cold Stream"
     105        Properties.Cold.Average.Mu =    PP.VapourViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
     106        Properties.Cold.Inlet.Mu =              PP.VapourViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
     107        Properties.Cold.Outlet.Mu =             PP.VapourViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
     108
     109"Conductivity Cold Stream"
     110        Properties.Cold.Average.K =     PP.VapourThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
     111        Properties.Cold.Inlet.K =               PP.VapourThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
     112        Properties.Cold.Outlet.K =              PP.VapourThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
     113       
     114"Heat Capacity Cold Stream"
    147115        Properties.Cold.Wall.Cp         =       PP.VapourCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
    148116
    149117
    150 "Cold Stream Viscosity at Wall Temperature"
    151         Properties.Cold.Wall.Mu         =       PP.VapourViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
    152 
    153 "Cold Stream Conductivity at Wall Temperature"
    154         Properties.Cold.Wall.K          =       PP.VapourThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
    155        
    156        
    157        
    158 end
    159 
     118"Viscosity Cold Stream"
     119        Properties.Cold.Wall.Mu =       PP.VapourViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
     120
     121"Conductivity Cold Stream"
     122        Properties.Cold.Wall.K =        PP.VapourThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
     123       
     124       
     125       
     126end
    160127
    161128if Inlet.Hot.v equal 0
    … …  
    163130        then
    164131
    165 "Hot Stream Average Heat Capacity"
     132"Heat Capacity Hot Stream"
    166133        Properties.Hot.Average.Cp       =               PP.LiquidCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
    167 
    168 "Hot Stream Inlet Heat Capacity"
    169134        Properties.Hot.Inlet.Cp         =               PP.LiquidCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
    170 
    171 "Hot Stream Outlet Heat Capacity"
    172135        Properties.Hot.Outlet.Cp        =               PP.LiquidCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    173136
    174 "Hot Stream Average Mass Density"
     137"Mass Density Hot Stream"
    175138        Properties.Hot.Average.rho      =               PP.LiquidDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
    176 
    177 "Hot Stream Inlet Mass Density"
    178139        Properties.Hot.Inlet.rho        =               PP.LiquidDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
    179 
    180 "Hot Stream Outlet Mass Density"       
    181140        Properties.Hot.Outlet.rho       =               PP.LiquidDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    182141
    183 "Hot Stream Average Viscosity"
     142"Viscosity Hot Stream"
    184143        Properties.Hot.Average.Mu       =               PP.LiquidViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);     
    185 
    186 "Hot Stream Inlet Viscosity"
    187144        Properties.Hot.Inlet.Mu         =               PP.LiquidViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);       
    188 
    189 "Hot Stream Outlet Viscosity"
    190145        Properties.Hot.Outlet.Mu        =               PP.LiquidViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);     
    191146
    192 "Hot Stream Average Conductivity"
     147"Conductivity Hot Stream"
    193148        Properties.Hot.Average.K        =               PP.LiquidThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);   
    194 
    195 "Hot Stream Inlet Conductivity"
    196         Properties.Hot.Inlet.K          =               PP.LiquidThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);     
    197 
    198 "Hot Stream Outlet Conductivity"
     149        Properties.Hot.Inlet.K  =               PP.LiquidThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);     
    199150        Properties.Hot.Outlet.K         =               PP.LiquidThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);   
    200151
    201 "Hot Stream Heat Capacity at Wall Temperature"
    202         Properties.Hot.Wall.Cp          =               PP.LiquidCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
    203 
    204 "Hot Stream Viscosity  at Wall Temperature"
    205         Properties.Hot.Wall.Mu          =               PP.LiquidViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);     
    206 
    207 "Hot Stream Conductivity at Wall Temperature"
    208         Properties.Hot.Wall.K           =               PP.LiquidThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);   
     152"Heat Capacity Hot Stream"
     153        Properties.Hot.Wall.Cp  =               PP.LiquidCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
     154
     155"Viscosity Hot Stream"
     156        Properties.Hot.Wall.Mu  =               PP.LiquidViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);     
     157
     158"Conductivity Hot Stream"
     159        Properties.Hot.Wall.K   =               PP.LiquidThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);   
    209160       
    210161
    211162        else
    212163
    213 "Hot Stream Average Heat Capacity"
     164"Heat Capacity Hot Stream"
    214165        Properties.Hot.Average.Cp       =               PP.VapourCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
    215 
    216 "Hot Stream Inlet Heat Capacity"
    217166        Properties.Hot.Inlet.Cp         =               PP.VapourCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
    218 
    219 "Hot Stream Outlet Heat Capacity"
    220167        Properties.Hot.Outlet.Cp        =               PP.VapourCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    221168
    222 "Hot Stream Average Mass Density"
     169"Mass Density Hot Stream"
    223170        Properties.Hot.Average.rho      =               PP.VapourDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
    224 
    225 "Hot Stream Inlet Mass Density"
    226171        Properties.Hot.Inlet.rho        =               PP.VapourDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
    227 
    228 "Hot Stream Outlet Mass Density"
    229172        Properties.Hot.Outlet.rho       =               PP.VapourDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    230173
    231 "Hot Stream Average Viscosity"
     174"Viscosity Hot Stream"
    232175        Properties.Hot.Average.Mu       =               PP.VapourViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
    233 
    234 "Hot Stream Inlet Viscosity"
    235176        Properties.Hot.Inlet.Mu         =               PP.VapourViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
    236 
    237 "Hot Stream Outlet Viscosity"
    238177        Properties.Hot.Outlet.Mu        =               PP.VapourViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    239178
    240 "Hot Stream Average Conductivity"
     179"Conductivity Hot Stream"
    241180        Properties.Hot.Average.K        =               PP.VapourThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);   
    242 
    243 "Hot Stream Inlet Conductivity"
    244         Properties.Hot.Inlet.K          =               PP.VapourThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);     
    245        
    246 "Hot Stream Outlet Conductivity"
     181        Properties.Hot.Inlet.K  =               PP.VapourThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);     
    247182        Properties.Hot.Outlet.K         =               PP.VapourThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);   
    248183
    249 "Hot Stream Heat Capacity at Wall Temperature"
    250         Properties.Hot.Wall.Cp          =               PP.VapourCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
    251 
    252 "Hot Stream Viscosity at Wall Temperature"
    253         Properties.Hot.Wall.Mu          =               PP.VapourViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
    254 
    255 "Hot Stream Conductivity at Wall Temperature"
    256         Properties.Hot.Wall.K           =               PP.VapourThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);   
     184"Heat Capacity Hot Stream"
     185        Properties.Hot.Wall.Cp  =               PP.VapourCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
     186
     187"Viscosity Hot Stream"
     188        Properties.Hot.Wall.Mu  =               PP.VapourViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
     189
     190"Conductivity Hot Stream"
     191        Properties.Hot.Wall.K   =               PP.VapourThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);   
    257192
    258193
    … …  
    338273"TEMA E Shell Effectiveness"
    339274        Eft = HE.EshellEffectiveness(Details.Cr,Details.NTU);
    340 
     275       
    341276end
    342277
    … …  
    357292        MTD   = Fc*LMTD;
    358293       
    359 "LMTD Correction Factor"
    360         Fc = HE.EshellCorrectionFactor(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T);
     294#"LMTD Correction Factor"
     295#       Fc = HE.EshellCorrectionFactor(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T);
    361296
    362297"Temperature Difference at Inlet"
    … …  
    810745       
    811746"Overall Heat Transfer Coefficient"
    812         Unity(i).Details.U*(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)))=1;
     747        Unity(i).Details.U      =       1/(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)));
    813748
    814749"Exchange Surface Area"
    … …  
    1020955        Unity(Nb+1).Shell.HeatTransfer.Sm = HE.CrossFlowArea(Unity(Nb+1).Baffles.Lso);
    1021956
    1022 #"LMTD Correction Factor"
    1023 #       Unity.Fc = HE.EshellCorrectionFactor(Unity(1).Inlet.Hot.T,Unity(Nb+1).Outlet.Hot.T,Unity(Nb+1).Inlet.Cold.T,Unity(1).Outlet.Cold.T);
     957"LMTD Correction Factor"
     958        Unity.Fc = HE.EshellCorrectionFactor(Unity(1).Inlet.Hot.T,Unity(Nb+1).Outlet.Hot.T,Unity(Nb+1).Inlet.Cold.T,Unity(1).Outlet.Cold.T);
    1024959
    1025960if side equal 1
    … …  
    12991234       
    13001235"Overall Heat Transfer Coefficient"
    1301         Unity(i).Details.U*(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)))=1;
     1236        Unity(i).Details.U      =       1/(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)));
    13021237
    13031238"Exchange Surface Area"

  • /mso/eml/heat_exchangers/HeatExchangerSimplified.mso

    r71 r63  
    5353
    5454
     55
    5556if Inlet.Cold.v equal 0
    56        
    5757        then   
    58        
    59 "Cold Stream Average Heat Capacity"
    60         Properties.Cold.Average.Cp      =       PP.LiquidCp(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
    61 
    62 "Cold Stream Inlet Heat Capacity"
    63         Properties.Cold.Inlet.Cp        =       PP.LiquidCp(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
    64 
    65 "Cold Stream Outlet Heat Capacity"
    66         Properties.Cold.Outlet.Cp       =       PP.LiquidCp(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    67 
    68 "Cold Stream Average Mass Density"
    69         Properties.Cold.Average.rho =   PP.LiquidDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
    70 
    71 "Cold Stream Inlet Mass Density"
    72         Properties.Cold.Inlet.rho       =       PP.LiquidDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
    73 
    74 "Cold Stream Outlet Mass Density"
    75         Properties.Cold.Outlet.rho      =       PP.LiquidDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    76 
    77 "Cold Stream Average Viscosity"
    78         Properties.Cold.Average.Mu      =       PP.LiquidViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
    79 
    80 "Cold Stream inlet Viscosity"
    81         Properties.Cold.Inlet.Mu        =       PP.LiquidViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
    82        
    83 "Cold Stream Outlet Viscosity"
    84         Properties.Cold.Outlet.Mu       =       PP.LiquidViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    85 
    86 "Cold Stream Average Conductivity"
    87         Properties.Cold.Average.K       =       PP.LiquidThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
    88 
    89 "Cold Stream Inlet Conductivity"       
    90         Properties.Cold.Inlet.K         =       PP.LiquidThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
    91 
    92 "Cold Stream Outlet Conductivity"
    93         Properties.Cold.Outlet.K        =       PP.LiquidThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    94 
    95 "Cold Stream Heat Capacity at Wall Temperature"
    96         Properties.Cold.Wall.Cp         =       PP.LiquidCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
    97        
    98 "Cold Stream Viscosity at Wall Temperature"
    99         Properties.Cold.Wall.Mu         =       PP.LiquidViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
    100 
    101 "Cold Stream Conductivity at Wall Temperature"
    102         Properties.Cold.Wall.K          =       PP.LiquidThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
     58"Heat Capacity Cold Stream"
     59        Properties.Cold.Average.Cp              =       PP.LiquidCp(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
     60        Properties.Cold.Inlet.Cp                =       PP.LiquidCp(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
     61        Properties.Cold.Outlet.Cp               =       PP.LiquidCp(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
     62
     63"Mass Density Cold Stream"
     64        Properties.Cold.Average.rho     =       PP.LiquidDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
     65        Properties.Cold.Inlet.rho               =       PP.LiquidDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
     66        Properties.Cold.Outlet.rho              =       PP.LiquidDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
     67
     68"Viscosity Cold Stream"
     69        Properties.Cold.Average.Mu              =       PP.LiquidViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
     70        Properties.Cold.Inlet.Mu                =       PP.LiquidViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
     71        Properties.Cold.Outlet.Mu               =       PP.LiquidViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
     72
     73"Conductivity Cold Stream"
     74        Properties.Cold.Average.K               =       PP.LiquidThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
     75        Properties.Cold.Inlet.K                 =       PP.LiquidThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
     76        Properties.Cold.Outlet.K                =       PP.LiquidThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
     77
     78"Heat Capacity Cold Stream"
     79        Properties.Cold.Wall.Cp                 =       PP.LiquidCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
     80       
     81"Viscosity Cold Stream"
     82        Properties.Cold.Wall.Mu                 =       PP.LiquidViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
     83
     84"Conductivity Cold Stream"
     85        Properties.Cold.Wall.K                  =       PP.LiquidThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
    10386
    10487
    10588        else
    10689
    107 "Cold Stream Average Heat Capacity"
     90"Heat Capacity Cold Stream"
    10891        Properties.Cold.Average.Cp      =       PP.VapourCp(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
    109 
    110 "Cold Stream Inlet Heat Capacity"       
    11192        Properties.Cold.Inlet.Cp        =       PP.VapourCp(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
    112 
    113 "Cold Stream Outlet Heat Capacity"     
    11493        Properties.Cold.Outlet.Cp       =       PP.VapourCp(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    11594
    116 "Cold Stream Average Mass Density"
    117         Properties.Cold.Average.rho =   PP.VapourDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
    118 
    119 "Cold Stream Inlet Mass Density"
    120         Properties.Cold.Inlet.rho       =       PP.VapourDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
    121 
    122 "Cold Stream Outlet Mass Density"       
    123         Properties.Cold.Outlet.rho      =       PP.VapourDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    124 
    125 "Cold Stream Average Viscosity "
    126         Properties.Cold.Average.Mu      =       PP.VapourViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
    127 
    128 "Cold Stream Inlet Viscosity " 
    129         Properties.Cold.Inlet.Mu        =       PP.VapourViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
    130 
    131 "Cold Stream Outlet Viscosity "
    132         Properties.Cold.Outlet.Mu       =       PP.VapourViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    133 
    134 "Cold Stream Average Conductivity "
    135         Properties.Cold.Average.K       =       PP.VapourThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
    136 
    137 "Cold Stream Inlet Conductivity "
    138         Properties.Cold.Inlet.K         =       PP.VapourThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
    139 
    140 "Cold Stream Outlet Conductivity "
    141         Properties.Cold.Outlet.K        =       PP.VapourThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    142        
    143 "Cold Stream Heat Capacity at Wall Temperature"
    144         Properties.Cold.Wall.Cp         =       PP.VapourCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
    145 
    146 
    147 "Cold Stream Viscosity at Wall Temperature"
    148         Properties.Cold.Wall.Mu         =       PP.VapourViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
    149 
    150 "Cold Stream Conductivity at Wall Temperature"
    151         Properties.Cold.Wall.K          =       PP.VapourThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
     95"Mass Density Cold Stream"
     96        Properties.Cold.Average.rho     =       PP.VapourDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
     97        Properties.Cold.Inlet.rho               =       PP.VapourDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
     98        Properties.Cold.Outlet.rho              =       PP.VapourDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
     99
     100"Viscosity Cold Stream"
     101        Properties.Cold.Average.Mu              =       PP.VapourViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
     102        Properties.Cold.Inlet.Mu                =       PP.VapourViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
     103        Properties.Cold.Outlet.Mu               =       PP.VapourViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
     104
     105"Conductivity Cold Stream"
     106        Properties.Cold.Average.K               =       PP.VapourThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
     107        Properties.Cold.Inlet.K                 =       PP.VapourThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
     108        Properties.Cold.Outlet.K                =       PP.VapourThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
     109       
     110"Heat Capacity Cold Stream"
     111        Properties.Cold.Wall.Cp                 =       PP.VapourCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
     112
     113
     114"Viscosity Cold Stream"
     115        Properties.Cold.Wall.Mu                 =       PP.VapourViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
     116
     117"Conductivity Cold Stream"
     118        Properties.Cold.Wall.K                  =       PP.VapourThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
    152119       
    153120       
    … …  
    159126        then
    160127
    161 "Hot Stream Average Heat Capacity"
     128"Heat Capacity Hot Stream"
    162129        Properties.Hot.Average.Cp       =               PP.LiquidCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
    163 
    164 "Hot Stream Inlet Heat Capacity"
    165130        Properties.Hot.Inlet.Cp         =               PP.LiquidCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
    166 
    167 "Hot Stream Outlet Heat Capacity"
    168131        Properties.Hot.Outlet.Cp        =               PP.LiquidCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    169132
    170 "Hot Stream Average Mass Density"
     133"Mass Density Hot Stream"
    171134        Properties.Hot.Average.rho      =               PP.LiquidDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
    172 
    173 "Hot Stream Inlet Mass Density"
    174135        Properties.Hot.Inlet.rho        =               PP.LiquidDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
    175 
    176 "Hot Stream Outlet Mass Density"       
    177136        Properties.Hot.Outlet.rho       =               PP.LiquidDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    178137
    179 "Hot Stream Average Viscosity"
     138"Viscosity Hot Stream"
    180139        Properties.Hot.Average.Mu       =               PP.LiquidViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);     
    181 
    182 "Hot Stream Inlet Viscosity"
    183140        Properties.Hot.Inlet.Mu         =               PP.LiquidViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);       
    184 
    185 "Hot Stream Outlet Viscosity"
    186141        Properties.Hot.Outlet.Mu        =               PP.LiquidViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);     
    187142
    188 "Hot Stream Average Conductivity"
     143"Conductivity Hot Stream"
    189144        Properties.Hot.Average.K        =               PP.LiquidThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);   
    190 
    191 "Hot Stream Inlet Conductivity"
    192145        Properties.Hot.Inlet.K          =               PP.LiquidThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);     
    193 
    194 "Hot Stream Outlet Conductivity"
    195146        Properties.Hot.Outlet.K         =               PP.LiquidThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);   
    196147
    197 "Hot Stream Heat Capacity at Wall Temperature"
     148"Heat Capacity Hot Stream"
    198149        Properties.Hot.Wall.Cp          =               PP.LiquidCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
    199150
    200 "Hot Stream Viscosity  at Wall Temperature"
     151"Viscosity Hot Stream"
    201152        Properties.Hot.Wall.Mu          =               PP.LiquidViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);     
    202153
    203 "Hot Stream Conductivity at Wall Temperature"
     154"Conductivity Hot Stream"
    204155        Properties.Hot.Wall.K           =               PP.LiquidThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);   
    205156       
    … …  
    207158        else
    208159
    209 "Hot Stream Average Heat Capacity"
     160"Heat Capacity Hot Stream"
    210161        Properties.Hot.Average.Cp       =               PP.VapourCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
    211 
    212 "Hot Stream Inlet Heat Capacity"
    213162        Properties.Hot.Inlet.Cp         =               PP.VapourCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
    214 
    215 "Hot Stream Outlet Heat Capacity"
    216163        Properties.Hot.Outlet.Cp        =               PP.VapourCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    217164
    218 "Hot Stream Average Mass Density"
     165"Mass Density Hot Stream"
    219166        Properties.Hot.Average.rho      =               PP.VapourDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
    220 
    221 "Hot Stream Inlet Mass Density"
    222167        Properties.Hot.Inlet.rho        =               PP.VapourDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
    223 
    224 "Hot Stream Outlet Mass Density"
    225168        Properties.Hot.Outlet.rho       =               PP.VapourDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    226169
    227 "Hot Stream Average Viscosity"
     170"Viscosity Hot Stream"
    228171        Properties.Hot.Average.Mu       =               PP.VapourViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
    229 
    230 "Hot Stream Inlet Viscosity"
    231172        Properties.Hot.Inlet.Mu         =               PP.VapourViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
    232 
    233 "Hot Stream Outlet Viscosity"
    234173        Properties.Hot.Outlet.Mu        =               PP.VapourViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    235174
    236 "Hot Stream Average Conductivity"
     175"Conductivity Hot Stream"
    237176        Properties.Hot.Average.K        =               PP.VapourThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);   
    238 
    239 "Hot Stream Inlet Conductivity"
    240177        Properties.Hot.Inlet.K          =               PP.VapourThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);     
    241        
    242 "Hot Stream Outlet Conductivity"
    243178        Properties.Hot.Outlet.K         =               PP.VapourThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);   
    244179
    245 "Hot Stream Heat Capacity at Wall Temperature"
     180"Heat Capacity Hot Stream"
    246181        Properties.Hot.Wall.Cp          =               PP.VapourCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
    247182
    248 "Hot Stream Viscosity at Wall Temperature"
     183"Viscosity Hot Stream"
    249184        Properties.Hot.Wall.Mu          =               PP.VapourViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
    250185
    251 "Hot Stream Conductivity at Wall Temperature"
     186"Conductivity Hot Stream"
    252187        Properties.Hot.Wall.K           =               PP.VapourThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);   
    253188
    254189
    255190end
    256 
    257191
    258192#=====================================================================
    … …  
    483417
    484418Model HeatExchanger_NTU         as Heatex_Basic_NTU
    485        
    486 PARAMETERS
    487 
    488 Side as Integer (Brief="Flow Direction",Lower=0,Upper=1);
    489 
    490 SET
    491 
    492 Side = HE.FlowDir(); # Return Flow Direction
    493 
    494 EQUATIONS
    495 
    496 if Details.Cr equal 0
    497        
    498         then   
     419
     420EQUATIONS
    499421"Effectiveness"
    500         Eft = 1-exp(-Details.NTU);
    501        
    502         else
    503 
    504 if Side equal 0
    505 
    506         then
    507 "Effectiveness in Cocurrent Flow"
    508         Eft = (1-exp(-Details.NTU*(1+Details.Cr)))/(1+Details.Cr);
    509        
    510         else
    511 
    512 if Details.Cr equal 1
    513        
    514         then
    515 "Effectiveness in Counter Flow"
    516         Eft = Details.NTU/(1+Details.NTU);
    517        
    518         else
    519 "Effectiveness in Counter Flow"
    520         Eft*(1-Details.Cr*exp(-Details.NTU*(1-Details.Cr))) = (1-exp(-Details.NTU*(1-Details.Cr)));
    521        
    522 end
    523 
    524 end
    525 
    526 
    527 end
    528 
     422        Eft=HE.Effectiveness(Details.Cr,Details.NTU);
    529423       
    530424end

  • /mso/eml/reactors/pfr.mso

    r71 r63  
    4545        q(NDisc)    as heat_rate;
    4646        M(NComp, NDisc)    as mol (Brief = "Molar holdup");
    47         Mt(NDisc)    as mol (Brief = "Molar holdup");
    48         C(NComp, NDisc)  as conc_mol(Brief="Components concentration", Lower=-1e-6);
     47        C(NComp, NDisc)  as conc_mol(Brief="Components concentration");
    4948        E(NDisc) as energy (Brief="Total Energy Holdup on element");
    5049       
    … …  
    7170       
    7271        for z in [1:NDisc]
    73                 for c in [1:NComp-1]
     72                for c in [1:NComp]
    7473                        "Component Molar Balance"
    7574                        diff(M(c,z)) = (str(z).F*str(z).z(c) - str(z+1).F*str(z+1).z(c))
    … …  
    8281
    8382                "Energy Holdup"
    84                 E(z) = Mt(z)*str(z+1).h - str(z+1).P*Across*L/NDisc;
     83                E(z) = sum(M(:,z))*str(z+1).h - str(z+1).P*Across*L/NDisc;
    8584
    8685                "mass flow is considered constant"
    8786                str(z+1).F*vol(z+1) = str(z).F*vol(z); # FIXME: is this correct? No (constant velocity: only for equimolar)
    88                 #rho(z+1)*vel(z+1) = rho(z)*vel(z);  # FIXME: this is correct! But does not converge.
     87#               rho(z+1)*vel(z+1) = rho(z)*vel(z);  # FIXME: this is correct! But does not converge.
    8988       
    9089                "Molar concentration"
    9190                C(:,z) * Across*L/NDisc = M(:,z);
    92                
    93                 "Sum of M"
    94                 Mt(z) = sum(M(:,z));
    95                
    96                 "Geometrical constraint"
    97                 Across*L/NDisc = Mt(z) * vol(z);
    9891
    9992                "Molar fraction"
    100                 str(z+1).z * Mt(z) = M(:,z);
    101         end
     93                str(z+1).z = C(:,z) * vol(z+1);
     94        end     
    10295
    10396        for z in [1:NDisc+1]
    … …  
    10699               
    107100                "Specific Mass"
    108                 rho(z) = PP.VapourDensity(str(z).T, str(z).P, str(z).z);
     101                rho(z) = sum(str(z).z*Mw)/vol(z);
    109102
    110103                "Velocity"

  • /mso/eml/stage_separators/condenser.mso

    r71 r63  
    101101ext PP as CalcObject;
    102102ext NComp as Integer;
     103        DP as press_delta (Brief="Pressure Drop in the condenser");
    103104
    104105        VARIABLES
    … …  
    106107out     OutletL as stream_therm; #(Brief="Liquid outlet stream");
    107108in      Q as heat_rate (Brief="Heat supplied");
    108         DP as press_delta (Brief="Pressure Drop in the condenser");
    109109
    110110        EQUATIONS

  • /mso/eml/stage_separators/tank.mso

    r71 r63  
    131131        Fout = k*sqrt(h);               
    132132end
    133 
    134 Model tank_feed
    135 
    136         PARAMETERS
    137 ext PP as CalcObject;
    138 ext NComp as Integer;
    139         Across as area (Brief="Tank cross section area", Default=2);
    140        
    141         VARIABLES
    142 in      Feed as stream;
    143 in      Inlet as stream;
    144 out     Outlet as stream_therm;
    145 
    146 in      Q as heat_rate (Brief="Rate of heat supply");
    147         Level    as length(Brief="Tank level");
    148         M(NComp) as mol (Brief="Molar Holdup in the tank");
    149         E as energy (Brief="Total Energy Holdup on tank");
    150         vL as volume_mol (Brief="Liquid Molar Volume");
    151 
    152         EQUATIONS
    153         "Mass balance"
    154         diff(M) = Feed.F*Feed.z + Inlet.F*Inlet.z - Outlet.F*Outlet.z;
    155        
    156         "Energy balance"
    157         diff(E) = Feed.F*Feed.h + Inlet.F*Inlet.h - Outlet.F*Outlet.h + Q;
    158 
    159         "Energy Holdup"
    160         E = sum(M)*Outlet.h;
    161 
    162         "Mechanical Equilibrium"
    163         Inlet.P = Outlet.P;
    164        
    165         "Liquid Volume"
    166         vL = PP.LiquidVolume(Outlet.T, Outlet.P, Outlet.z);
    167        
    168         "Composition"
    169         M = Outlet.z*sum(M);
    170        
    171         "Level of liquid phase"
    172         Level = sum(M)*vL/Across;
    173        
    174         "Vapourisation Fraction"
    175         Outlet.v = Inlet.v;
    176        
    177 end

  • /mso/eml/types.mso

    r71 r63  
    2222positive as Real (Brief = "Positive General Constant", Default=1.0, Lower=-1e-6);
    2323negative as Real (Brief = "Negative General Constant", Default=-1.0, Upper=0.0);
    24 fraction as positive (Brief = "Fraction" , Default=0.5, Upper=1.00001);
     24fraction as positive (Brief = "Fraction" , Default=0.5, Upper=1);
    2525percent as Real (Brief = "Percent", Default=50, Lower=0, Upper=100);
    2626control_signal as Real (Brief = "Control Signal", Default=1, Lower=-1e9, Upper=1e9);

  • /mso/sample/heat_exchangers/Eshell_Detailed_LMTD.mso

    r71 r63  
    22#===============================================================
    33#   Heat Exchanger TEMA E Shell - LMTD Method
     4#  converge após várias tentativas
    45#===============================================================
    56

  • /mso/sample/heat_exchangers/Eshell_Detailed_NTU.mso

    r71 r63  
    3535#       Shell Geometrical Parameters
    3636#=====================================================================   
    37 exchanger.Tpass                         = 4;
     37exchanger.Tpass                         = 2;
    3838exchanger.Dishell               = 0.75  *"m";
    3939exchanger.Lcf                   = 0.043 *"m";
    … …  
    5353exchanger.Ditube                = 0.013395  *"m";
    5454exchanger.Dotube                = 0.015875  *"m";
    55 exchanger.Kwall                         = 0.057         *"kW/m/K";
     55exchanger.Kwall                         = 0.57          *"kW/m/K";
    5656exchanger.Donozzle_Tube     = 0.203     *"m";
    5757exchanger.Dinozzle_Tube         = 0.203         *"m";

  • /mso/sample/heat_exchangers/Eshell_Discretized_LMTD.mso

    r71 r63  
    100100#   Simulation Options
    101101#============================================
    102  mode           = "steady";
    103 guessFile   = "E_Shell_Discretized_NTU";
     102 mode                   = "steady";
     103
    104104end

  • /mso/sample/heat_exchangers/Eshell_Discretized_NTU.mso

    r71 r63  
    9797#============================================
    9898 mode = "steady";
    99 guessFile = "E_Shell_Discretized_NTU";
     99
    100100end

  • /mso/sample/heat_exchangers/Multipass_Detailed.mso

    r71 r63  
    9999#============================================
    100100 mode                   = "steady";
    101 guessFile  = "Multipass_LMTD";
     101
    102102end

  • /mso/sample/heat_exchangers/NTU_Method.mso

    r71 r63  
    1111        streamcold      as streamTP;
    1212
     13
     14       
    1315CONNECTIONS
    1416
    … …  
    1820PARAMETERS
    1921
    20         PP                      as CalcObject   (File="vrpp");
     22        PP                      as CalcObject   (File="vrpp.dll");
     23        HE                      as CalcObject   (File="heatex.dll");
    2124        NComp           as Integer;
    2225       
    … …  
    3033        PP.Components           = ["water"];
    3134        NComp                           = PP.NumberOfComponents;
    32         exchanger.HE.HotSide          = "Shell";
     35        HE.HotSide          = "Shell";
    3336       
    34         exchanger.HE.TurbulentFlow     = "Petukhov";
    35 #       exchanger.HE.TurbulentFlow     = "SiederTate";
     37#       HE.TurbulentFlow     = "Petukhov";
     38        HE.TurbulentFlow     = "SiederTate";
    3639       
    3740#=====================================================================

  • /mso/sample/heat_exchangers/sampleEshell_LMTD.mso

    r71 r63  
    2323
    2424        exchanger.HE.LMTDcorrection     = "Bowmann";
    25 #       exchanger.HE.LMTDcorrection     = "Fakeri";
     25#       HE.LMTDcorrection       = "Fakeri";
    2626       
    2727

  • /mso/sample/heat_exchangers/sampleNTU.mso

    r71 r63  
    99NComp as Integer;
    1010       
     11
     12       
    1113DEVICES
    12 
    1314exchanger       as HeatExchanger_NTU;
    1415streamhot       as streamTP;
    … …  
    1920PP.LiquidModel          = "PR";
    2021PP.VapourModel          = "PR";
    21 PP.Components           = ["water","n-butane", "benzene", "n-octane" ];
    22 NComp                           = PP.NumberOfComponents;
    23 
    24 exchanger.HE.FlowDirection = "Cocurrent";
     22PP.Components   = ["water","n-butane", "benzene", "n-octane" ];
     23NComp                                           = PP.NumberOfComponents;
     24exchanger.HE.FlowDirection    = "Cocurrent";
    2525
    2626CONNECTIONS
    … …  
    3232
    3333exchanger.Details.U                             = 210*"W/(m^2*K)";
    34 exchanger.PressureDrop.Hot.Pdrop        = 0.2*"kPa";
     34exchanger.PressureDrop.Hot.Pdrop        = 0.1*"kPa";
    3535exchanger.PressureDrop.Cold.Pdrop   = 0.2*"kPa";
    36 
    3736streamhot.F             = 20    * "kmol/h";
    3837streamhot.T             = 450     * "K";
    3938streamhot.P             = 120    * "kPa";
    4039streamhot.z             = [1,0,0,0];
    41 
    4240streamcold.F            = 10 * "kmol/h";
    4341streamcold.P            = 120 * "kPa";
    4442streamcold.T            = 300 * "K";   
    4543streamcold.z            = [0,0.5, 0.1, 0.4];
    46 
    47 exchanger.Details.A = 20*"m^2";
    48 
     44exchanger.Outlet.Cold.T = 340*"K";
     45       
    4946OPTIONS
    50 
    51  mode = "steady";
    52 
    53 end
     47 mode   = "steady";
     48 end

  • /mso/sample/heat_exchangers/samples1.mso

    r71 r63  
    2525SET
    2626
    27         PP.LiquidModel          = "PR";
    28         PP.VapourModel          = "PR";
     27        PP.LiquidModel          = "RK";
     28        PP.VapourModel          = "RK";
    2929        PP.Components           = ["water"];
    3030        NComp                           = PP.NumberOfComponents;

  • /mso/sample/miscellaneous/sample_streams.mso

    r71 r63  
    6666        sv.v = 1.0; #vapor stream
    6767
    68         EQUATIONS
    6968        "Liquid-Vapour equilibrium"
    7069        PP.VapourFugacityCoefficient(sv.T, sv.P, sv.z) * sv.z =

  • /mso/sample/optimization/flash_opt.mso

    r71 r63  
    5454        FREE
    5555        fl.OutletL.T;
    56         #fl.OutletL.P;
    5756       
    5857        EQUATIONS
    59         fl.OutletL.T < 320 * "K";
     58        fl.OutletL.T < 315 * "K";
    6059        fl.OutletL.T > 300 * "K";
    6160

  • /mso/sample/reactors/sample_pfr.mso

    r71 r63  
    3535        NComp = PP.NumberOfComponents;
    3636       
    37         Reac.NDisc = 10;
     37        Reac.NDisc = 20;
    3838        Reac.Across = 0.7 * "in^2";
    3939        Reac.L = 2.28 * "m";
    … …  
    7171        for z in [2:Reac.NDisc+1]
    7272                Reac.str(z).T = Reac.Inlet.T;
    73                 Reac.str(z).z(1:NComp-1) = Reac.Inlet.z(1:NComp-1);
     73                Reac.str(z).z = Reac.Inlet.z;
    7474        end
    7575       

  • /mso/sample/stage_separators/sample_condenser.mso

    r71 r63  
    8282        s1.v = 1.0;
    8383        Q = 3.71e6 * "kJ/h";
     84
     85        SET
    8486        c1.DP = 100 * "kPa";
    8587
Note: See TracChangeset for help on using the changeset viewer.