Changeset 68 for mso


Ignore:
Timestamp:
Nov 27, 2006, 6:37:59 PM (16 years ago)
Author:
gerson bicca
Message:

added double pipe heat exchanger model and updated eml/heat_exchangers

Location:
mso
Files:
4 added
15 edited

Legend:

Unmodified
Added
Removed
  • mso/eml/heat_exchangers

    • Property svn:ignore set to
      *.dll
  • mso/eml/heat_exchangers/HEX_Engine.mso

    r45 r68  
    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*----------------------------------------------------------------------
    4044* Author: Gerson Balbueno Bicca
    4145* $Id$
     
    277281Cold    as Basic_Pdrop;# Cold Stream
    278282end     
     283
     284Model DoublePipe_HeatTransfer
     285#=====================================================================
     286#       Heat Transfer
     287#===================================================================== 
     288PARAMETERS
     289As              as area                         (Brief="Cross Sectional Area for Flow",Default=0.05,Lower=1e-8);
     290Dh      as length                       (Brief="Hydraulic Diameter of Pipe for Heat Transfer",Lower=1e-8);
     291
     292VARIABLES
     293Re              as positive                     (Brief="Reynolds Number",Default=100,Lower=1);
     294hcoeff  as heat_trans_coeff (Brief="Film Coefficient",Default=1,Lower=1e-12, Upper=1e6);
     295PR              as positive                     (Brief="Prandtl Number",Default=0.5,Lower=1e-8);
     296Phi     as positive                     (Brief="Phi Correction",Default=1,Lower=1e-3);
     297Vmean   as velocity                     (Brief="Tube Velocity",Lower=1e-8);
     298end
     299
     300Model DoublePipe_PressureDrop
     301#=====================================================================
     302#       Pressure Drop
     303#===================================================================== 
     304PARAMETERS
     305Dh      as length               (Brief="Hydraulic Diameter of Pipe for Pressure Drop",Lower=1e-6);
     306
     307VARIABLES
     308Pdrop   as press_delta  (Brief="Pressure Drop",Default=0.01, Lower=1e-10);
     309fi      as fricfactor   (Brief="Friction Factor", Default=0.05, Lower=1e-10, Upper=2000);
     310Re              as positive             (Brief="Reynolds Number",Default=100,Lower=1);
     311end     
     312
     313Model Main_DoublePipe
     314VARIABLES
     315HeatTransfer as DoublePipe_HeatTransfer;
     316PressureDrop as DoublePipe_PressureDrop;
     317end
  • mso/eml/heat_exchangers/HeatExchangerDetailed.mso

    r45 r68  
    165165        then
    166166
    167 "Heat Capacity Hot Stream"
     167"Hot Stream Average Heat Capacity"
    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"
    169171        Properties.Hot.Inlet.Cp         =               PP.LiquidCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
     172
     173"Hot Stream Outlet Heat Capacity"
    170174        Properties.Hot.Outlet.Cp        =               PP.LiquidCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    171175
    172 "Mass Density Hot Stream"
     176"Hot Stream Average Mass Density"
    173177        Properties.Hot.Average.rho      =               PP.LiquidDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
     178
     179"Hot Stream Inlet Mass Density"
    174180        Properties.Hot.Inlet.rho        =               PP.LiquidDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
     181
     182"Hot Stream Outlet Mass Density"       
    175183        Properties.Hot.Outlet.rho       =               PP.LiquidDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    176184
    177 "Viscosity Hot Stream"
     185"Hot Stream Average Viscosity"
    178186        Properties.Hot.Average.Mu       =               PP.LiquidViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);     
     187
     188"Hot Stream Inlet Viscosity"
    179189        Properties.Hot.Inlet.Mu         =               PP.LiquidViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);       
     190
     191"Hot Stream Outlet Viscosity"
    180192        Properties.Hot.Outlet.Mu        =               PP.LiquidViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);     
    181193
    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"
     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"
    200216        Properties.Hot.Average.Cp       =               PP.VapourCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
     217
     218"Hot Stream Inlet Heat Capacity"
    201219        Properties.Hot.Inlet.Cp         =               PP.VapourCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
     220
     221"Hot Stream Outlet Heat Capacity"
    202222        Properties.Hot.Outlet.Cp        =               PP.VapourCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    203223
    204 "Mass Density Hot Stream"
     224"Hot Stream Average Mass Density"
    205225        Properties.Hot.Average.rho      =               PP.VapourDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
     226
     227"Hot Stream Inlet Mass Density"
    206228        Properties.Hot.Inlet.rho        =               PP.VapourDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
     229
     230"Hot Stream Outlet Mass Density"
    207231        Properties.Hot.Outlet.rho       =               PP.VapourDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    208232
    209 "Viscosity Hot Stream"
     233"Hot Stream Average Viscosity"
    210234        Properties.Hot.Average.Mu       =               PP.VapourViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
     235
     236"Hot Stream Inlet Viscosity"
    211237        Properties.Hot.Inlet.Mu         =               PP.VapourViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
     238
     239"Hot Stream Outlet Viscosity"
    212240        Properties.Hot.Outlet.Mu        =               PP.VapourViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    213241
    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);   
     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);   
    227259
    228260
     
    423455
    424456"Shell Side Phi correction"
    425         Shell.HeatTransfer.Phi = HE.PhiCorrection(Properties.Hot.Average.Mu,Properties.Hot.Wall.Mu);
     457        Shell.HeatTransfer.Phi  = HE.PhiCorrection(Properties.Hot.Average.Mu,Properties.Hot.Wall.Mu);
    426458       
    427459"Tube Side Phi correction"
     
    431463       
    432464"Shell Side Phi correction"
    433         Shell.HeatTransfer.Phi = HE.PhiCorrection(Properties.Cold.Average.Mu,Properties.Cold.Wall.Mu);
     465        Shell.HeatTransfer.Phi  = HE.PhiCorrection(Properties.Cold.Average.Mu,Properties.Cold.Wall.Mu);
    434466       
    435467"Tube Side Phi correction"
     
    449481       
    450482"Hot Wall Temperature"
    451         Properties.Hot.Wall.Twall  = (Properties.Hot.Average.T+Properties.Cold.Average.T)/2;
     483        Properties.Hot.Wall.Twall       = (Properties.Hot.Average.T+Properties.Cold.Average.T)/2;
    452484
    453485"ColdWall Temperature"
     
    455487
    456488"Tube Side Velocity"
    457         Tubes.HeatTransfer.Vtube  = Properties.Cold.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Properties.Cold.Average.rho*Ntt);
     489        Tubes.HeatTransfer.Vtube        = Properties.Cold.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Properties.Cold.Average.rho*Ntt);
    458490
    459491"Tube Side Reynolds Number"
    460         Tubes.HeatTransfer.Re     =     (Properties.Cold.Average.rho*Tubes.HeatTransfer.Vtube*Ditube)/Properties.Cold.Average.Mu;
     492        Tubes.HeatTransfer.Re           = (Properties.Cold.Average.rho*Tubes.HeatTransfer.Vtube*Ditube)/Properties.Cold.Average.Mu;
    461493       
    462494"Tube Side Prandtl Number"
    463         Tubes.HeatTransfer.PR = ((Properties.Cold.Average.Cp/Properties.Cold.Average.Mw)*Properties.Cold.Average.Mu)/Properties.Cold.Average.K;
     495        Tubes.HeatTransfer.PR           = ((Properties.Cold.Average.Cp/Properties.Cold.Average.Mw)*Properties.Cold.Average.Mu)/Properties.Cold.Average.K;
    464496
    465497"Tube Side Prandtl Number at Wall Temperature"
    466         Tubes.HeatTransfer.PRw = ((Properties.Cold.Wall.Cp/Properties.Cold.Average.Mw)*Properties.Cold.Wall.Mu)/Properties.Cold.Wall.K;
     498        Tubes.HeatTransfer.PRw          = ((Properties.Cold.Wall.Cp/Properties.Cold.Average.Mw)*Properties.Cold.Wall.Mu)/Properties.Cold.Wall.K;
    467499
    468500"Tube Side Film Coefficient"
    469         Tubes.HeatTransfer.htube= HE.TubeFilmCoeff(Tubes.HeatTransfer.Re,Tubes.HeatTransfer.PR,Properties.Cold.Average.K)*Tubes.HeatTransfer.Phi;
     501        Tubes.HeatTransfer.htube        = HE.TubeFilmCoeff(Tubes.HeatTransfer.Re,Tubes.HeatTransfer.PR,Properties.Cold.Average.K)*Tubes.HeatTransfer.Phi;
    470502
    471503"Shell Side Prandtl Number"
    472         Shell.HeatTransfer.PR = ((Properties.Hot.Average.Cp/Properties.Hot.Average.Mw)*Properties.Hot.Average.Mu)/Properties.Hot.Average.K;
     504        Shell.HeatTransfer.PR           = ((Properties.Hot.Average.Cp/Properties.Hot.Average.Mw)*Properties.Hot.Average.Mu)/Properties.Hot.Average.K;
    473505
    474506"Shell Side Prandtl Number at Wall Temperature"
    475         Shell.HeatTransfer.PRw = ((Properties.Hot.Wall.Cp/Properties.Hot.Average.Mw)*Properties.Hot.Wall.Mu)/Properties.Hot.Wall.K;
     507        Shell.HeatTransfer.PRw          = ((Properties.Hot.Wall.Cp/Properties.Hot.Average.Mw)*Properties.Hot.Wall.Mu)/Properties.Hot.Wall.K;
    476508
    477509"Tube Side Pressure Drop"
     
    497529
    498530"Shell Pressure Drop Outlet Nozzle"
    499         Shell.PressureDrop.Pdnozzle_out =HE.DeltaPshellNozzleOut(Properties.Hot.Outlet.rho,Properties.Hot.Inlet.Fw);
     531        Shell.PressureDrop.Pdnozzle_out = HE.DeltaPshellNozzleOut(Properties.Hot.Outlet.rho,Properties.Hot.Inlet.Fw);
    500532
    501533"Velocity Shell Side Outlet Nozzle"
     
    621653       
    622654VARIABLES
     655
    623656DT0     as temp_delta   (Brief="Temperature Difference at Inlet",Lower=1);
    624657DTL             as temp_delta   (Brief="Temperature Difference at Outlet",Lower=1);
     
    12951328VARIABLES
    12961329
    1297 LMTD(Nshell)            as temp_delta   (Brief="Logarithmic Mean Temperature Difference",Lower=10);
    1298 Fc(Nshell)                      as positive             (Brief="LMTD Correction Factor",Lower=0.75);
    1299 MTD(Nshell)                     as temp_delta   (Brief="Mean Temperature Difference");
     1330LMTD(Nshell)            as temp_delta   (Brief="Logarithmic Mean Temperature Difference",Lower=5);
     1331Fc(Nshell)                      as positive             (Brief="LMTD Correction Factor",Lower=0.5);
     1332MTD(Nshell)                     as temp_delta   (Brief="Mean Temperature Difference",Lower=5);
    13001333Unity(Nshell)           as HeatExchangerDetailed_Basic;
    13011334
     
    15151548       
    15161549"Overall Heat Transfer Coefficient"
    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)));
     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;
    15181551
    15191552"Exchange Surface Area"
     
    15731606
    15741607"Exchange Surface Area"
    1575         Unity(i).Details.Q   = Unity(i).Details.U*Unity(i).Details.A*Fc(i)*LMTD(i);
     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);
    15761610       
    15771611"Mean Temperature Difference"   
  • mso/eml/heat_exchangers/HeatExchangerDiscretized.mso

    r45 r68  
    6060        then   
    6161       
    62 "Heat Capacity Cold Stream"
     62"Cold Stream Average Heat Capacity"
    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"
    6466        Properties.Cold.Inlet.Cp        =       PP.LiquidCp(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
     67
     68"Cold Stream Outlet Heat Capacity"
    6569        Properties.Cold.Outlet.Cp       =       PP.LiquidCp(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    6670
    67 "Mass Density Cold Stream"
     71"Cold Stream Average Mass Density"
    6872        Properties.Cold.Average.rho =   PP.LiquidDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
     73
     74"Cold Stream Inlet Mass Density"
    6975        Properties.Cold.Inlet.rho       =       PP.LiquidDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
     76
     77"Cold Stream Outlet Mass Density"
    7078        Properties.Cold.Outlet.rho      =       PP.LiquidDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    7179
    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"
     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"
    8399        Properties.Cold.Wall.Cp         =       PP.LiquidCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
    84100       
    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);
     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);
    90106
    91107
    92108        else
    93109
    94 "Heat Capacity Cold Stream"
     110"Cold Stream Average Heat Capacity"
    95111        Properties.Cold.Average.Cp      =       PP.VapourCp(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
     112
     113"Cold Stream Inlet Heat Capacity"       
    96114        Properties.Cold.Inlet.Cp        =       PP.VapourCp(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
     115
     116"Cold Stream Outlet Heat Capacity"     
    97117        Properties.Cold.Outlet.Cp       =       PP.VapourCp(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    98118
    99 "Mass Density Cold Stream"
    100         Properties.Cold.Average.rho     =       PP.VapourDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
     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"
    101123        Properties.Cold.Inlet.rho       =       PP.VapourDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
     124
     125"Cold Stream Outlet Mass Density"       
    102126        Properties.Cold.Outlet.rho      =       PP.VapourDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    103127
    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"
     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"
    115147        Properties.Cold.Wall.Cp         =       PP.VapourCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
    116148
    117149
    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        
    126 end
     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       
     158end
     159
    127160
    128161if Inlet.Hot.v equal 0
     
    130163        then
    131164
    132 "Heat Capacity Hot Stream"
     165"Hot Stream Average Heat Capacity"
    133166        Properties.Hot.Average.Cp       =               PP.LiquidCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
     167
     168"Hot Stream Inlet Heat Capacity"
    134169        Properties.Hot.Inlet.Cp         =               PP.LiquidCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
     170
     171"Hot Stream Outlet Heat Capacity"
    135172        Properties.Hot.Outlet.Cp        =               PP.LiquidCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    136173
    137 "Mass Density Hot Stream"
     174"Hot Stream Average Mass Density"
    138175        Properties.Hot.Average.rho      =               PP.LiquidDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
     176
     177"Hot Stream Inlet Mass Density"
    139178        Properties.Hot.Inlet.rho        =               PP.LiquidDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
     179
     180"Hot Stream Outlet Mass Density"       
    140181        Properties.Hot.Outlet.rho       =               PP.LiquidDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    141182
    142 "Viscosity Hot Stream"
     183"Hot Stream Average Viscosity"
    143184        Properties.Hot.Average.Mu       =               PP.LiquidViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);     
     185
     186"Hot Stream Inlet Viscosity"
    144187        Properties.Hot.Inlet.Mu         =               PP.LiquidViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);       
     188
     189"Hot Stream Outlet Viscosity"
    145190        Properties.Hot.Outlet.Mu        =               PP.LiquidViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);     
    146191
    147 "Conductivity Hot Stream"
     192"Hot Stream Average Conductivity"
    148193        Properties.Hot.Average.K        =               PP.LiquidThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);   
    149         Properties.Hot.Inlet.K  =               PP.LiquidThermalConductivity(Inlet.Hot.T,Inlet.Hot.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"
    150199        Properties.Hot.Outlet.K         =               PP.LiquidThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);   
    151200
    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);   
     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);   
    160209       
    161210
    162211        else
    163212
    164 "Heat Capacity Hot Stream"
     213"Hot Stream Average Heat Capacity"
    165214        Properties.Hot.Average.Cp       =               PP.VapourCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
     215
     216"Hot Stream Inlet Heat Capacity"
    166217        Properties.Hot.Inlet.Cp         =               PP.VapourCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
     218
     219"Hot Stream Outlet Heat Capacity"
    167220        Properties.Hot.Outlet.Cp        =               PP.VapourCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    168221
    169 "Mass Density Hot Stream"
     222"Hot Stream Average Mass Density"
    170223        Properties.Hot.Average.rho      =               PP.VapourDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
     224
     225"Hot Stream Inlet Mass Density"
    171226        Properties.Hot.Inlet.rho        =               PP.VapourDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
     227
     228"Hot Stream Outlet Mass Density"
    172229        Properties.Hot.Outlet.rho       =               PP.VapourDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    173230
    174 "Viscosity Hot Stream"
     231"Hot Stream Average Viscosity"
    175232        Properties.Hot.Average.Mu       =               PP.VapourViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
     233
     234"Hot Stream Inlet Viscosity"
    176235        Properties.Hot.Inlet.Mu         =               PP.VapourViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
     236
     237"Hot Stream Outlet Viscosity"
    177238        Properties.Hot.Outlet.Mu        =               PP.VapourViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    178239
    179 "Conductivity Hot Stream"
     240"Hot Stream Average Conductivity"
    180241        Properties.Hot.Average.K        =               PP.VapourThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);   
    181         Properties.Hot.Inlet.K  =               PP.VapourThermalConductivity(Inlet.Hot.T,Inlet.Hot.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"
    182247        Properties.Hot.Outlet.K         =               PP.VapourThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);   
    183248
    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);   
     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);   
    192257
    193258
     
    273338"TEMA E Shell Effectiveness"
    274339        Eft = HE.EshellEffectiveness(Details.Cr,Details.NTU);
    275        
     340
    276341end
    277342
     
    292357        MTD   = Fc*LMTD;
    293358       
    294 #"LMTD Correction Factor"
    295 #       Fc = HE.EshellCorrectionFactor(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T);
     359"LMTD Correction Factor"
     360        Fc = HE.EshellCorrectionFactor(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T);
    296361
    297362"Temperature Difference at Inlet"
     
    745810       
    746811"Overall Heat Transfer Coefficient"
    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)));
     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;
    748813
    749814"Exchange Surface Area"
     
    9551020        Unity(Nb+1).Shell.HeatTransfer.Sm = HE.CrossFlowArea(Unity(Nb+1).Baffles.Lso);
    9561021
    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);
     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);
    9591024
    9601025if side equal 1
     
    12341299       
    12351300"Overall Heat Transfer Coefficient"
    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)));
     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;
    12371302
    12381303"Exchange Surface Area"
  • mso/eml/heat_exchangers/HeatExchangerSimplified.mso

    r45 r68  
    5353
    5454
    55 
    5655if Inlet.Cold.v equal 0
     56       
    5757        then   
    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);
    86 
    87 
    88         else
    89 
    90 "Heat Capacity Cold Stream"
     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);
     103
     104
     105        else
     106
     107"Cold Stream Average Heat Capacity"
    91108        Properties.Cold.Average.Cp      =       PP.VapourCp(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
     109
     110"Cold Stream Inlet Heat Capacity"       
    92111        Properties.Cold.Inlet.Cp        =       PP.VapourCp(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
     112
     113"Cold Stream Outlet Heat Capacity"     
    93114        Properties.Cold.Outlet.Cp       =       PP.VapourCp(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
    94115
    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);
     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);
    119152       
    120153       
     
    126159        then
    127160
    128 "Heat Capacity Hot Stream"
     161"Hot Stream Average Heat Capacity"
    129162        Properties.Hot.Average.Cp       =               PP.LiquidCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
     163
     164"Hot Stream Inlet Heat Capacity"
    130165        Properties.Hot.Inlet.Cp         =               PP.LiquidCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
     166
     167"Hot Stream Outlet Heat Capacity"
    131168        Properties.Hot.Outlet.Cp        =               PP.LiquidCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    132169
    133 "Mass Density Hot Stream"
     170"Hot Stream Average Mass Density"
    134171        Properties.Hot.Average.rho      =               PP.LiquidDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
     172
     173"Hot Stream Inlet Mass Density"
    135174        Properties.Hot.Inlet.rho        =               PP.LiquidDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
     175
     176"Hot Stream Outlet Mass Density"       
    136177        Properties.Hot.Outlet.rho       =               PP.LiquidDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    137178
    138 "Viscosity Hot Stream"
     179"Hot Stream Average Viscosity"
    139180        Properties.Hot.Average.Mu       =               PP.LiquidViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);     
     181
     182"Hot Stream Inlet Viscosity"
    140183        Properties.Hot.Inlet.Mu         =               PP.LiquidViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);       
     184
     185"Hot Stream Outlet Viscosity"
    141186        Properties.Hot.Outlet.Mu        =               PP.LiquidViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);     
    142187
    143 "Conductivity Hot Stream"
     188"Hot Stream Average Conductivity"
    144189        Properties.Hot.Average.K        =               PP.LiquidThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);   
     190
     191"Hot Stream Inlet Conductivity"
    145192        Properties.Hot.Inlet.K          =               PP.LiquidThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);     
     193
     194"Hot Stream Outlet Conductivity"
    146195        Properties.Hot.Outlet.K         =               PP.LiquidThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);   
    147196
    148 "Heat Capacity Hot Stream"
     197"Hot Stream Heat Capacity at Wall Temperature"
    149198        Properties.Hot.Wall.Cp          =               PP.LiquidCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
    150199
    151 "Viscosity Hot Stream"
     200"Hot Stream Viscosity  at Wall Temperature"
    152201        Properties.Hot.Wall.Mu          =               PP.LiquidViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);     
    153202
    154 "Conductivity Hot Stream"
     203"Hot Stream Conductivity at Wall Temperature"
    155204        Properties.Hot.Wall.K           =               PP.LiquidThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);   
    156205       
     
    158207        else
    159208
    160 "Heat Capacity Hot Stream"
     209"Hot Stream Average Heat Capacity"
    161210        Properties.Hot.Average.Cp       =               PP.VapourCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
     211
     212"Hot Stream Inlet Heat Capacity"
    162213        Properties.Hot.Inlet.Cp         =               PP.VapourCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
     214
     215"Hot Stream Outlet Heat Capacity"
    163216        Properties.Hot.Outlet.Cp        =               PP.VapourCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    164217
    165 "Mass Density Hot Stream"
     218"Hot Stream Average Mass Density"
    166219        Properties.Hot.Average.rho      =               PP.VapourDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
     220
     221"Hot Stream Inlet Mass Density"
    167222        Properties.Hot.Inlet.rho        =               PP.VapourDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
     223
     224"Hot Stream Outlet Mass Density"
    168225        Properties.Hot.Outlet.rho       =               PP.VapourDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    169226
    170 "Viscosity Hot Stream"
     227"Hot Stream Average Viscosity"
    171228        Properties.Hot.Average.Mu       =               PP.VapourViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);
     229
     230"Hot Stream Inlet Viscosity"
    172231        Properties.Hot.Inlet.Mu         =               PP.VapourViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);
     232
     233"Hot Stream Outlet Viscosity"
    173234        Properties.Hot.Outlet.Mu        =               PP.VapourViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);
    174235
    175 "Conductivity Hot Stream"
     236"Hot Stream Average Conductivity"
    176237        Properties.Hot.Average.K        =               PP.VapourThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);   
     238
     239"Hot Stream Inlet Conductivity"
    177240        Properties.Hot.Inlet.K          =               PP.VapourThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);     
     241       
     242"Hot Stream Outlet Conductivity"
    178243        Properties.Hot.Outlet.K         =               PP.VapourThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);   
    179244
    180 "Heat Capacity Hot Stream"
     245"Hot Stream Heat Capacity at Wall Temperature"
    181246        Properties.Hot.Wall.Cp          =               PP.VapourCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
    182247
    183 "Viscosity Hot Stream"
     248"Hot Stream Viscosity at Wall Temperature"
    184249        Properties.Hot.Wall.Mu          =               PP.VapourViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
    185250
    186 "Conductivity Hot Stream"
     251"Hot Stream Conductivity at Wall Temperature"
    187252        Properties.Hot.Wall.K           =               PP.VapourThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);   
    188253
    189254
    190255end
     256
    191257
    192258#=====================================================================
     
    417483
    418484Model HeatExchanger_NTU         as Heatex_Basic_NTU
    419 
    420 EQUATIONS
     485       
     486PARAMETERS
     487
     488Side as Integer (Brief="Flow Direction",Lower=0,Upper=1);
     489
     490SET
     491
     492Side = HE.FlowDir(); # Return Flow Direction
     493
     494EQUATIONS
     495
     496if Details.Cr equal 0
     497       
     498        then   
    421499"Effectiveness"
    422         Eft=HE.Effectiveness(Details.Cr,Details.NTU);
     500        Eft = 1-exp(-Details.NTU);
     501       
     502        else
     503
     504if 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
     512if 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       
     522end
     523
     524end
     525
     526
     527end
     528
    423529       
    424530end
  • mso/sample/heat_exchangers/Eshell_Detailed_LMTD.mso

    r45 r68  
    22#===============================================================
    33#   Heat Exchanger TEMA E Shell - LMTD Method
    4 #  converge após várias tentativas
    54#===============================================================
    65
  • mso/sample/heat_exchangers/Eshell_Detailed_NTU.mso

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

    r45 r68  
    100100#   Simulation Options
    101101#============================================
    102  mode                   = "steady";
    103 
     102 mode           = "steady";
     103guessFile   = "E_Shell_Discretized_NTU";
    104104end
  • mso/sample/heat_exchangers/Eshell_Discretized_NTU.mso

    r45 r68  
    9797#============================================
    9898 mode = "steady";
    99 
     99guessFile = "E_Shell_Discretized_NTU";
    100100end
  • mso/sample/heat_exchangers/Multipass_Detailed.mso

    r45 r68  
    9999#============================================
    100100 mode                   = "steady";
    101 
     101guessFile  = "Multipass_LMTD";
    102102end
  • mso/sample/heat_exchangers/NTU_Method.mso

    r45 r68  
    1111        streamcold      as streamTP;
    1212
    13 
    14        
    1513CONNECTIONS
    1614
     
    2018PARAMETERS
    2119
    22         PP                      as CalcObject   (File="vrpp.dll");
    23         HE                      as CalcObject   (File="heatex.dll");
     20        PP                      as CalcObject   (File="vrpp");
    2421        NComp           as Integer;
    2522       
     
    3330        PP.Components           = ["water"];
    3431        NComp                           = PP.NumberOfComponents;
    35         HE.HotSide          = "Shell";
     32        exchanger.HE.HotSide          = "Shell";
    3633       
    37 #       HE.TurbulentFlow     = "Petukhov";
    38         HE.TurbulentFlow     = "SiederTate";
     34        exchanger.HE.TurbulentFlow     = "Petukhov";
     35#       exchanger.HE.TurbulentFlow     = "SiederTate";
    3936       
    4037#=====================================================================
  • mso/sample/heat_exchangers/sampleEshell_LMTD.mso

    r45 r68  
    2323
    2424        exchanger.HE.LMTDcorrection     = "Bowmann";
    25 #       HE.LMTDcorrection       = "Fakeri";
     25#       exchanger.HE.LMTDcorrection     = "Fakeri";
    2626       
    2727
  • mso/sample/heat_exchangers/sampleNTU.mso

    r45 r68  
    99NComp as Integer;
    1010       
     11DEVICES
    1112
    12        
    13 DEVICES
    1413exchanger       as HeatExchanger_NTU;
    1514streamhot       as streamTP;
     
    2019PP.LiquidModel          = "PR";
    2120PP.VapourModel          = "PR";
    22 PP.Components   = ["water","n-butane", "benzene", "n-octane" ];
    23 NComp                                           = PP.NumberOfComponents;
    24 exchanger.HE.FlowDirection    = "Cocurrent";
     21PP.Components           = ["water","n-butane", "benzene", "n-octane" ];
     22NComp                           = PP.NumberOfComponents;
     23
     24exchanger.HE.FlowDirection = "Cocurrent";
    2525
    2626CONNECTIONS
     
    3232
    3333exchanger.Details.U                             = 210*"W/(m^2*K)";
    34 exchanger.PressureDrop.Hot.Pdrop        = 0.1*"kPa";
     34exchanger.PressureDrop.Hot.Pdrop        = 0.2*"kPa";
    3535exchanger.PressureDrop.Cold.Pdrop   = 0.2*"kPa";
     36
    3637streamhot.F             = 20    * "kmol/h";
    3738streamhot.T             = 450     * "K";
    3839streamhot.P             = 120    * "kPa";
    3940streamhot.z             = [1,0,0,0];
     41
    4042streamcold.F            = 10 * "kmol/h";
    4143streamcold.P            = 120 * "kPa";
    4244streamcold.T            = 300 * "K";   
    4345streamcold.z            = [0,0.5, 0.1, 0.4];
    44 exchanger.Outlet.Cold.T = 340*"K";
    45        
     46
     47exchanger.Details.A = 20*"m^2";
     48
    4649OPTIONS
    47  mode   = "steady";
    48  end
     50
     51 mode = "steady";
     52
     53end
  • mso/sample/heat_exchangers/samples1.mso

    r45 r68  
    2525SET
    2626
    27         PP.LiquidModel          = "RK";
    28         PP.VapourModel          = "RK";
     27        PP.LiquidModel          = "PR";
     28        PP.VapourModel          = "PR";
    2929        PP.Components           = ["water"];
    3030        NComp                           = PP.NumberOfComponents;
Note: See TracChangeset for help on using the changeset viewer.