Changeset 523 for trunk/eml


Ignore:
Timestamp:
May 23, 2008, 4:44:08 PM (14 years ago)
Author:
gerson bicca
Message:

improved the plate heat exchanger model

Location:
trunk/eml
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/eml/heat_exchangers/HEX_Engine.mso

    r492 r523  
    416416
    417417end
    418 
    419 Model PHE_HeatTransfer
    420        
    421 ATTRIBUTES
    422         Pallete = false;
    423         Brief = "to be documented";
    424         Info =
    425         "to be documented";
    426        
    427 VARIABLES
    428 
    429 Re                              as positive                                     (Brief="Reynolds Number",Default=100,Lower=1);
    430 PR                              as positive                                     (Brief="Prandtl Number",Default=0.5,Lower=1e-8);
    431 NTU                     as positive                                     (Brief="Number of Units Transference",Default=0.05,Lower=1e-10);
    432 WCp                     as positive                                     (Brief="Stream Heat Capacity",Lower=1e-3,Default=1e3,Unit='W/K');
    433 hcoeff          as heat_trans_coeff     (Brief="Film Coefficient",Default=1,Lower=1e-12, Upper=1e6);
    434 Gchannel        as flux_mass                            (Brief ="Channel Mass Flux", Default=1, Lower=1e-6, Symbol ="G^{channel}");
    435 Gports          as flux_mass                            (Brief ="Ports Mass Flux", Default=1, Lower=1e-6, Symbol ="G^{ports}");
    436 Phi                     as positive                                     (Brief="Viscosity Correction",Default=1,Lower=1e-6, Symbol="\phi");
    437 
    438 end
    439 
    440 Model PHE_PressureDrop
    441        
    442 ATTRIBUTES
    443         Pallete = false;
    444         Brief = "to be documented";
    445         Info =
    446         "to be documented";
    447        
    448 VARIABLES
    449 
    450 DPchannel                       as press_delta  (Brief="Channel Pressure Drop",Default=0.01, Lower=1e-10,DisplayUnit='kPa', Symbol ="\Delta P^{channel}");
    451 DPports                         as press_delta  (Brief="Ports Pressure Drop",Default=0.01, Lower=1e-10,DisplayUnit='kPa', Symbol ="\Delta P^{ports}");
    452 Pdrop                           as press_delta  (Brief="Total Pressure Drop",Default=0.01, Lower=1e-10,DisplayUnit='kPa', Symbol ="\Delta P");
    453 fi                                      as fricfactor           (Brief="Friction Factor", Default=0.05, Lower=1e-10, Upper=2000);
    454 Vchannel                        as velocity             (Brief="Stream Velocity in Channel",Lower=1e-8, Symbol ="V^{channel}");
    455 Vports                          as velocity             (Brief="Stream Velocity in Ports",Lower=1e-8, Symbol ="V^{ports}");
    456 Npassage                        as Real                                 (Brief="Number of  Channels per Pass", Symbol ="N^{passage}");
    457 end
    458 
    459 Model Main_PHE
    460        
    461 ATTRIBUTES
    462         Pallete = false;
    463         Brief = "to be documented";
    464         Info =
    465         "to be documented";
    466        
    467 VARIABLES
    468 
    469 HeatTransfer            as PHE_HeatTransfer             (Brief="PHE Heat Transfer", Symbol = " ");
    470 PressureDrop    as PHE_PressureDrop             (Brief="PHE Pressure Drop", Symbol = " ");
    471 Properties              as Physical_Properties          (Brief="PHE Properties", Symbol = " ");
    472 
    473 end
    474 
    475 Model Thermal_PHE
    476 
    477 ATTRIBUTES
    478         Pallete = false;
    479         Brief = "to be documented";
    480         Info =
    481         "to be documented";
    482        
    483 VARIABLES
    484 Cr      as positive                                     (Brief="Heat Capacity Ratio",Default=0.5,Lower=1e-6);
    485 Cmin    as positive                                     (Brief="Minimum Heat Capacity",Lower=1e-10,Default=1e3,Unit='W/K');
    486 Cmax    as positive                                     (Brief="Maximum Heat Capacity",Lower=1e-10,Default=1e3,Unit='W/K');
    487 NTU     as positive                                     (Brief="Number of Units Transference",Default=0.05,Lower=1e-10);
    488 Eft             as positive                             (Brief="Effectiveness",Default=0.5,Lower=0.4,Upper=1, Symbol = "\varepsilon");
    489 Q                       as power                                        (Brief="Heat Transfer", Default=7000, Lower=1e-6, Upper=1e10);
    490 Uc              as heat_trans_coeff     (Brief="Overall Heat Transfer Coefficient Clean",Default=1,Lower=1e-6,Upper=1e10);
    491 Ud              as heat_trans_coeff     (Brief="Overall Heat Transfer Coefficient Dirty",Default=1,Lower=1e-6,Upper=1e10);
    492 end
  • trunk/eml/heat_exchangers/PHE.mso

    r420 r523  
    1818using "HEX_Engine";
    1919
     20Model PHE_PressureDrop
     21       
     22ATTRIBUTES
     23        Pallete = false;
     24        Brief = "to be documented";
     25        Info =
     26        "to be documented";
     27       
     28VARIABLES
     29
     30DPchannel               as press_delta  (Brief="Channel Pressure Drop",Default=0.01, Lower=1E10,DisplayUnit='kPa', Symbol ="\Delta P^{channel}");
     31DPports                 as press_delta  (Brief="Ports Pressure Drop",Default=0.01, Lower=1E-10,DisplayUnit='kPa', Symbol ="\Delta P^{ports}");
     32Pdrop                           as press_delta  (Brief="Total Pressure Drop",Default=0.01, Lower=1E-10,DisplayUnit='kPa', Symbol ="\Delta P");
     33fi                                      as fricfactor           (Brief="Friction Factor", Default=0.05, Lower=1E-10, Upper=2000);
     34Vchannel                        as velocity             (Brief="Stream Velocity in Channel",Lower=1E-8, Symbol ="V^{channel}");
     35Vports                          as velocity             (Brief="Stream Velocity in Ports",Lower=1E-8, Symbol ="V^{ports}");
     36Npassage                        as positive             (Brief="Number of  Channels per Pass", Symbol ="N^{passage}");
     37
     38end
     39
     40Model PHE_HeatTransfer
     41       
     42ATTRIBUTES
     43        Pallete = false;
     44        Brief = "to be documented";
     45        Info =
     46        "to be documented";
     47       
     48VARIABLES
     49
     50Re                              as positive                                     (Brief="Reynolds Number",Default=100,Lower=1);
     51PR                              as positive                                     (Brief="Prandtl Number",Default=0.5,Lower=1e-8);
     52NTU                             as positive                                     (Brief="Number of Units Transference",Default=0.05,Lower=1E-10);
     53WCp                     as positive                                     (Brief="Stream Heat Capacity",Lower=1E-3,Default=1E3,Unit='W/K');
     54hcoeff                  as heat_trans_coeff             (Brief="Film Coefficient",Default=1,Lower=1E-12, Upper=1E6);
     55Gchannel                as flux_mass                            (Brief ="Channel Mass Flux", Default=1, Lower=1E-6, Symbol ="G^{channel}");
     56Gports                  as flux_mass                            (Brief ="Ports Mass Flux", Default=1, Lower=1E-6, Symbol ="G^{ports}");
     57Phi                     as positive                                     (Brief="Viscosity Correction",Default=1,Lower=1E-6, Symbol="\phi");
     58
     59end
     60
     61Model Main_PHE
     62       
     63ATTRIBUTES
     64        Pallete = false;
     65        Brief = "to be documented";
     66        Info =
     67        "to be documented";
     68       
     69VARIABLES
     70
     71HeatTransfer    as PHE_HeatTransfer             (Brief="PHE Heat Transfer", Symbol = " ");
     72PressureDrop    as PHE_PressureDrop     (Brief="PHE Pressure Drop", Symbol = " ");
     73Properties              as Physical_Properties          (Brief="PHE Properties", Symbol = " ");
     74
     75end
     76
     77Model Thermal_PHE
     78
     79ATTRIBUTES
     80        Pallete = false;
     81        Brief = "to be documented";
     82        Info =
     83        "to be documented";
     84       
     85VARIABLES
     86Cr      as positive                             (Brief="Heat Capacity Ratio",Default=0.5,Lower=1E-6);
     87Cmin    as positive                             (Brief="Minimum Heat Capacity",Lower=1E-10,Default=1E3,Unit='W/K');
     88Cmax    as positive                             (Brief="Maximum Heat Capacity",Lower=1E-10,Default=1E3,Unit='W/K');
     89NTU             as positive                             (Brief="Number of Units Transference",Default=0.05,Lower=1E-10);
     90Eft             as positive                             (Brief="Effectiveness",Default=0.5,Lower=0.1,Upper=1.1, Symbol = "\varepsilon");
     91Q                       as power                                        (Brief="Heat Transfer", Default=7000, Lower=1E-6, Upper=1E10);
     92Uc              as heat_trans_coeff     (Brief="Overall Heat Transfer Coefficient Clean",Default=1,Lower=1E-6,Upper=1E10);
     93Ud              as heat_trans_coeff     (Brief="Overall Heat Transfer Coefficient Dirty",Default=1,Lower=1E-6,Upper=1E10);
     94
     95end
     96
     97Model PHE_Geometry
     98
     99 ATTRIBUTES
     100        Pallete         = false;
     101        Brief           = "Parameters for a gasketed plate heat exchanger.";
     102
     103 PARAMETERS
     104
     105outer PP                as Plugin               (Brief="External Physical Properties", Type="PP");
     106outer NComp   as Integer        (Brief="Number of Chemical Components",Hidden=true);
     107       
     108        Pi                                              as constant             (Brief="Pi Number",Default=3.14159265, Hidden=true,Symbol = "\pi");
     109        N1                                      as Integer              (Brief="Auxiliar Constant", Hidden=true,Default = 15);
     110        N2                                      as Integer              (Brief="Auxiliar Constant",Hidden=true,Default = 14);
     111        Kp1(N1)                 as constant             (Brief="First constant in Kumar calculation for Pressure Drop", Hidden=true);
     112        Kp2(N1)                 as constant             (Brief="Second constant in Kumar calculation for Pressure Drop", Hidden=true);
     113        Kc1(N2)                 as constant             (Brief="First constant in Kumar calculation for Heat Transfer", Hidden=true);
     114        Kc2(N2)                 as constant             (Brief="Second constant Kumar calculation for Heat Transfer", Hidden=true);
     115        M(NComp)        as molweight    (Brief="Component Mol Weight", Hidden=true);
     116       
     117       
     118        Lv                                      as length                               (Brief="Vertical Ports Distance",Lower=0.1);
     119        Nplates                 as Integer                      (Brief="Total Number of Plates in The Whole Heat Exchanger",Default=25, Symbol ="N_{plates}");
     120        NpassHot                as Integer                      (Brief="Number of Passes for Hot Side", Symbol ="Npasshot");
     121        NpassCold               as Integer                      (Brief="Number of Passes for Cold Side", Symbol ="Npasscold");
     122        Dports                          as length                               (Brief="Ports Diameter",Lower=1e-6, Symbol ="D_{ports}");
     123        Lw                                      as length                               (Brief="Plate Width",Lower=0.1);
     124        pitch                           as length                               (Brief="Plate Pitch",Lower=0.1);
     125        pt                                              as length                               (Brief="Plate Thickness",Lower=0.1);
     126        Kwall                           as conductivity         (Brief="Plate Thermal Conductivity",Default=1.0, Symbol ="K_{wall}");
     127        Rfh                                     as positive                     (Brief="Hot Side Fouling Resistance",Unit='m^2*K/kW',Default=1e-6,Lower=0);
     128        Rfc                                     as positive                     (Brief="Cold Side Fouling Resistance",Unit='m^2*K/kW',Default=1e-6,Lower=0);
     129        PhiFactor                       as Real                                 (Brief="Enlargement Factor",Lower=1e-6, Symbol ="\phi");
     130       
     131        Atotal                                  as area                                 (Brief="Total Effective  Area",Lower=1e-6, Symbol ="A_{total}", Protected=true);
     132        Aports                          as area                                 (Brief="Port Opening  Area of Plate",Lower=1e-6, Symbol ="A_{ports}", Protected=true);
     133        Achannel                as area                                 (Brief="Cross-Sectional Area for Channel Flow",Lower=1e-6, Symbol ="A_{channel}", Protected=true);
     134        Dh                                      as length                               (Brief="Equivalent Diameter of Channel",Lower=1e-6, Protected=true);
     135        Depth                   as length                               (Brief="Corrugation Depth",Lower=1e-6, Protected=true);
     136        Nchannels               as Integer                      (Brief="Total Number of Channels in The Whole Heat Exchanger", Protected=true);
     137        Lp                                      as length                               (Brief="Plate Vertical Distance between Port Centers",Lower=0.1, Protected=true);
     138        Lpack                           as length                               (Brief="Compact Plate Pack Length",Lower=0.1, Protected=true);
     139        Lh                                      as length                               (Brief="Plate Horizontal Distance between Port Centers",Lower=0.1, Protected=true);
     140
     141SET
     142
     143#"Vector Length of constants for Kumar's calculating Pressure Drop"
     144        N1 = 15;
     145
     146#"Vector Length of constants for Kumar's calculating Heat Transfer"
     147        N2 = 14;
     148
     149#"First constant for Kumar's calculating Pressure Drop"
     150        Kp1 = [50,19.40,2.990,47,18.290,1.441,34,11.250,0.772,24,3.240,0.760,24,2.80,0.639];
     151
     152#"Second constant for Kumar's calculating Pressure Drop"
     153        Kp2 = [1,0.589,0.183,1,0.652,0.206,1,0.631,0.161,1,0.457,0.215,1,0.451,0.213];
     154
     155#"First constant for Kumar's calculating Heat Transfer"
     156        Kc1 = [0.718,0.348,0.718,0.400,0.300,0.630,0.291,0.130,0.562,0.306,0.108,0.562,0.331,0.087];
     157
     158#"Second constant for Kumar's calculating Heat Transfer"
     159        Kc2 = [0.349,0.663,0.349,0.598,0.663,0.333,0.591,0.732,0.326,0.529,0.703,0.326,0.503,0.718];
     160       
     161#"Component Molecular Weight"
     162        M  = PP.MolecularWeight();
     163       
     164#"Pi Number"
     165        Pi      = 3.14159265;
     166       
     167#"Plate Vertical Distance between Port Centers"
     168        Lp = Lv - Dports;
     169       
     170#"Corrugation Depth"
     171        Depth=pitch-pt;
     172       
     173#"Plate Horizontal Distance between Port Centers"
     174        Lh=Lw-Dports;
     175       
     176#"Hydraulic Diameter"
     177        Dh=2*Depth/PhiFactor;
     178
     179#"Ports Area"
     180        Aports=0.25*Pi*Dports*Dports;
     181       
     182#"Channel Area"
     183        Achannel=Depth*Lw;
     184       
     185#"Pack Length"
     186        Lpack=Depth*(Nplates-1)+Nplates*pt;
     187       
     188#"Total Number of  Channels"
     189        Nchannels = Nplates -1;
     190       
     191#"Exchange Surface Area"
     192        Atotal =(Nplates-2)*Lw*Lp*PhiFactor;
     193
     194end
     195
    20196Model PHE
    21197
     
    23199        Icon            = "icon/phe";
    24200        Pallete         = true;
    25         Brief           = "Shortcut model  for plate and Frame heat exchanger.";
     201        Brief           = "Shortcut model  for Plate and Frame heat exchanger.";
    26202        Info            =
    27203"Model of a gasketed plate heat exchanger.
     
    69245
    70246outer PP                as Plugin               (Brief="External Physical Properties", Type="PP");
    71 outer NComp             as Integer      (Brief="Number of Chemical Components");
    72         Pi                              as constant     (Brief="Pi Number",Default=3.14159265, Symbol = "\pi");
    73         N1                      as Integer      (Brief="Auxiliar Constant",Default = 15);
    74         N2                      as Integer      (Brief="Auxiliar Constant",Default = 14);
    75         Kp1(N1)                 as constant     (Brief="First constant in Kumar calculation for Pressure Drop");
    76         Kp2(N1)                 as constant     (Brief="Second constant in Kumar calculation for Pressure Drop");
    77         Kc1(N2)                 as constant     (Brief="First constant in Kumar calculation for Heat Transfer");
    78         Kc2(N2)                 as constant     (Brief="Second constant Kumar calculation for Heat Transfer");
    79         M(NComp)                as molweight    (Brief="Component Mol Weight");
    80        
    81         ChevronAngle    as Switcher             (Brief="Chevron Corrugation Inclination Angle in Degrees ",Valid=["A30_Deg","A45_Deg","A50_Deg","A60_Deg","A65_Deg"],Default="A30_Deg");
    82         SideOne                 as Switcher             (Brief="Fluid Alocation in the Side I - (The odd channels)",Valid=["hot","cold"],Default="hot");
    83         Nchannels               as Integer              (Brief="Total Number of Channels in The Whole Heat Exchanger");
    84         Nplates                 as Integer              (Brief="Total Number of Plates in The Whole Heat Exchanger",Default=25, Symbol ="N_{plates}");
    85         NpassHot                as Integer              (Brief="Number of Passes for Hot Side", Symbol ="Npasshot");
    86         NpassCold               as Integer              (Brief="Number of Passes for Cold Side", Symbol ="Npasscold");
    87         Dports                  as length               (Brief="Ports Diameter",Lower=1e-6, Symbol ="D_{ports}");
    88         Atotal                  as area                 (Brief="Total Effective  Area",Lower=1e-6, Symbol ="A_{total}");
    89         Aports                  as area                 (Brief="Port Opening  Area of Plate",Lower=1e-6, Symbol ="A_{ports}");
    90         Achannel                as area                 (Brief="Cross-Sectional Area for Channel Flow",Lower=1e-6, Symbol ="A_{channel}");
    91         Dh                      as length               (Brief="Equivalent Diameter of Channel",Lower=1e-6);
    92         Depth                   as length               (Brief="Corrugation Depth",Lower=1e-6);
    93         PhiFactor               as Real                 (Brief="Enlargement Factor",Lower=1e-6, Symbol ="\phi");
    94         Lp                              as length               (Brief="Plate Vertical Distance between Port Centers",Lower=0.1);
    95         Lpack                   as length               (Brief="Compact Plate Pack Length",Lower=0.1);
    96         Lv                              as length               (Brief="Vertical Ports Distance",Lower=0.1);
    97         Lh                              as length               (Brief="Plate Horizontal Distance between Port Centers",Lower=0.1);
    98         Lw                              as length               (Brief="Plate Width",Lower=0.1);
    99         pitch                   as length               (Brief="Plate Pitch",Lower=0.1);
    100         pt                              as length               (Brief="Plate Thickness",Lower=0.1);
    101         Kwall                   as conductivity (Brief="Plate Thermal Conductivity",Default=1.0, Symbol ="K_{wall}");
    102         Rfh                             as positive             (Brief="Hot Side Fouling Resistance",Unit='m^2*K/kW',Default=1e-6,Lower=0);
    103         Rfc                             as positive             (Brief="Cold Side Fouling Resistance",Unit='m^2*K/kW',Default=1e-6,Lower=0);
     247outer NComp   as Integer        (Brief="Number of Chemical Components");
     248       
     249        ChevronAngle    as Switcher                     (Brief="Chevron Corrugation Inclination Angle in Degrees ",Valid=["A30_Deg","A45_Deg","A50_Deg","A60_Deg","A65_Deg"],Default="A30_Deg");
     250        SideOne                         as Switcher                     (Brief="Fluid Alocation in the Side I - (The odd channels)",Valid=["hot","cold"],Default="hot");
    104251
    105252 VARIABLES
    106253
    107 in  InletHot        as stream           (Brief="Inlet Hot Stream", PosX=0, PosY=0.75, Symbol="^{inHot}");       
    108 in  InletCold       as stream           (Brief="Inlet Cold Stream", PosX=0, PosY=0.25, Symbol="^{inCold}");
    109 out OutletHot           as streamPH     (Brief="Outlet Hot Stream", PosX=1, PosY=0.25, Symbol="^{outHot}");
    110 out OutletCold          as streamPH     (Brief="Outlet Cold Stream", PosX=1, PosY=0.75, Symbol="^{outCold}");
    111 
    112         HotSide                 as Main_PHE             (Brief="Plate Heat Exchanger Hot Side", Symbol="_{hot}");
    113         ColdSide                as Main_PHE             (Brief="Plate Heat Exchanger Cold Side", Symbol="_{cold}");
    114         Thermal                 as Thermal_PHE  (Brief="Thermal Results", Symbol = " ");
    115 
    116  SET
    117 #"Vector Length of constants for Kumar's calculating Pressure Drop"
    118         N1 = 15;
    119 
    120 #"Vector Length of constants for Kumar's calculating Heat Transfer"
    121         N2 = 14;
    122 
    123 #"First constant for Kumar's calculating Pressure Drop"
    124         Kp1 = [50,19.40,2.990,47,18.290,1.441,34,11.250,0.772,24,3.240,0.760,24,2.80,0.639];
    125 
    126 #"Second constant for Kumar's calculating Pressure Drop"
    127         Kp2 = [1,0.589,0.183,1,0.652,0.206,1,0.631,0.161,1,0.457,0.215,1,0.451,0.213];
    128 
    129 #"First constant for Kumar's calculating Heat Transfer"
    130         Kc1 = [0.718,0.348,0.718,0.400,0.300,0.630,0.291,0.130,0.562,0.306,0.108,0.562,0.331,0.087];
    131 
    132 #"Second constant for Kumar's calculating Heat Transfer"
    133         Kc2 = [0.349,0.663,0.349,0.598,0.663,0.333,0.591,0.732,0.326,0.529,0.703,0.326,0.503,0.718];
    134        
    135 #"Component Molecular Weight"
    136         M  = PP.MolecularWeight();
    137        
    138 #"Pi Number"
    139         Pi      = 3.14159265;
    140        
    141 #"Plate Vertical Distance between Port Centers"
    142         Lp = Lv - Dports;
    143        
    144 #"Corrugation Depth"
    145         Depth=pitch-pt;
    146        
    147 #"Plate Horizontal Distance between Port Centers"
    148         Lh=Lw-Dports;
    149        
    150 #"Hydraulic Diameter"
    151         Dh=2*Depth/PhiFactor;
    152 
    153 #"Ports Area"
    154         Aports=0.25*Pi*Dports*Dports;
    155        
    156 #"Channel Area"
    157         Achannel=Depth*Lw;
    158        
    159 #"Pack Length"
    160         Lpack=Depth*(Nplates-1)+Nplates*pt;
    161        
    162 #"Total Number of  Channels"
    163         Nchannels = Nplates -1;
    164        
    165 #"Exchange Surface Area"
    166         Atotal =(Nplates-2)*Lw*Lp*PhiFactor;
    167        
     254Geometry                        as PHE_Geometry         (Brief="Plate Heat Exchanger Geometrical Parameters", Symbol=" ");
     255in  InletHot            as stream                                       (Brief="Inlet Hot Stream", PosX=0, PosY=0.75, Symbol="^{inHot}");       
     256in  InletCold           as stream                                       (Brief="Inlet Cold Stream", PosX=0, PosY=0.25, Symbol="^{inCold}");
     257out OutletHot           as streamPH                             (Brief="Outlet Hot Stream", PosX=1, PosY=0.25, Symbol="^{outHot}");
     258out OutletCold          as streamPH                             (Brief="Outlet Cold Stream", PosX=1, PosY=0.75, Symbol="^{outCold}");
     259
     260       
     261        HotSide                 as Main_PHE                             (Brief="Plate Heat Exchanger Hot Side", Symbol="_{hot}");
     262        ColdSide        as Main_PHE                             (Brief="Plate Heat Exchanger Cold Side", Symbol="_{cold}");
     263        Thermal                 as Thermal_PHE                  (Brief="Thermal Results", Symbol = " ");
     264
    168265 EQUATIONS
    169266
     
    187284
    188285"Hot Stream Average Molecular Weight"
    189         HotSide.Properties.Average.Mw = sum(M*InletHot.z);
     286        HotSide.Properties.Average.Mw = sum(Geometry.M*InletHot.z);
    190287
    191288"Cold Stream Average Molecular Weight"
    192         ColdSide.Properties.Average.Mw = sum(M*InletCold.z);
     289        ColdSide.Properties.Average.Mw = sum(Geometry.M*InletCold.z);
    193290
    194291 if InletCold.v equal 0
     
    300397
    301398"Flow Mass Inlet Cold Stream"
    302         ColdSide.Properties.Inlet.Fw    =  sum(M*InletCold.z)*InletCold.F;
     399        ColdSide.Properties.Inlet.Fw    =  sum(Geometry.M*InletCold.z)*InletCold.F;
    303400
    304401"Flow Mass Outlet Cold Stream"
    305         ColdSide.Properties.Outlet.Fw   =  sum(M*OutletCold.z)*OutletCold.F;
     402        ColdSide.Properties.Outlet.Fw   =  sum(Geometry.M*OutletCold.z)*OutletCold.F;
    306403
    307404"Flow Mass Inlet Hot Stream"
    308         HotSide.Properties.Inlet.Fw             =  sum(M*InletHot.z)*InletHot.F;
     405        HotSide.Properties.Inlet.Fw             =  sum(Geometry.M*InletHot.z)*InletHot.F;
    309406
    310407"Flow Mass Outlet Hot Stream"   
    311         HotSide.Properties.Outlet.Fw    =  sum(M*OutletHot.z)*OutletHot.F;
     408        HotSide.Properties.Outlet.Fw    =  sum(Geometry.M*OutletHot.z)*OutletHot.F;
    312409
    313410"Molar Balance Hot Stream"
     
    328425       
    329426"Total Number of  Passages Cold Side"
    330         ColdSide.PressureDrop.Npassage = (2*Nchannels+1+(-1)^(Nchannels+1))/(4*NpassCold);
     427        ColdSide.PressureDrop.Npassage = (2*Geometry.Nchannels+1+(-1)^(Geometry.Nchannels+1))/(4*Geometry.NpassCold);
    331428
    332429"Total Number of  Passages Hot Side"
    333         HotSide.PressureDrop.Npassage = (2*Nchannels-1+(-1)^(Nchannels))/(4*NpassHot);
     430        HotSide.PressureDrop.Npassage = (2*Geometry.Nchannels-1+(-1)^(Geometry.Nchannels))/(4*Geometry.NpassHot);
    334431       
    335432        case "hot":
    336433       
    337434"Total Number of  Passages Cold Side"
    338         HotSide.PressureDrop.Npassage = (2*Nchannels+1+(-1)^(Nchannels+1))/(4*NpassHot);
     435        HotSide.PressureDrop.Npassage = (2*Geometry.Nchannels+1+(-1)^(Geometry.Nchannels+1))/(4*Geometry.NpassHot);
    339436
    340437"Total Number of  Passages Hot Side"
    341         ColdSide.PressureDrop.Npassage = (2*Nchannels-1+(-1)^(Nchannels))/(4*NpassCold);
     438        ColdSide.PressureDrop.Npassage = (2*Geometry.Nchannels-1+(-1)^(Geometry.Nchannels))/(4*Geometry.NpassCold);
    342439
    343440 end
    344441
    345442"Hot Stream Mass Flux in the Channel"
    346         HotSide.HeatTransfer.Gchannel=HotSide.Properties.Inlet.Fw/(HotSide.PressureDrop.Npassage*Achannel);
     443        HotSide.HeatTransfer.Gchannel=HotSide.Properties.Inlet.Fw/(HotSide.PressureDrop.Npassage*Geometry.Achannel);
    347444
    348445"Hot Stream Mass Flux in the Ports"
    349         HotSide.HeatTransfer.Gports=HotSide.Properties.Inlet.Fw/Aports;
     446        HotSide.HeatTransfer.Gports=HotSide.Properties.Inlet.Fw/Geometry.Aports;
    350447
    351448"Cold Stream Mass Flux in the Ports"
    352         ColdSide.HeatTransfer.Gports=ColdSide.Properties.Inlet.Fw/Aports;
     449        ColdSide.HeatTransfer.Gports=ColdSide.Properties.Inlet.Fw/Geometry.Aports;
    353450
    354451"Cold Stream Mass Flux in the Channel"
    355         ColdSide.HeatTransfer.Gchannel=ColdSide.Properties.Inlet.Fw/(ColdSide.PressureDrop.Npassage*Achannel);
     452        ColdSide.HeatTransfer.Gchannel=ColdSide.Properties.Inlet.Fw/(ColdSide.PressureDrop.Npassage*Geometry.Achannel);
    356453
    357454"Hot Stream Pressure Drop in Ports"
    358         HotSide.PressureDrop.DPports =1.5*NpassHot*HotSide.HeatTransfer.Gports^2/(2*HotSide.Properties.Average.rho);
     455        HotSide.PressureDrop.DPports =1.5*Geometry.NpassHot*HotSide.HeatTransfer.Gports^2/(2*HotSide.Properties.Average.rho);
    359456
    360457"Cold Stream Pressure Drop in Ports"
    361         ColdSide.PressureDrop.DPports =1.5*NpassCold*ColdSide.HeatTransfer.Gports^2/(2*ColdSide.Properties.Average.rho);
     458        ColdSide.PressureDrop.DPports =1.5*Geometry.NpassCold*ColdSide.HeatTransfer.Gports^2/(2*ColdSide.Properties.Average.rho);
    362459
    363460"Hot Stream Pressure Drop in Channels"
    364         HotSide.PressureDrop.DPchannel =2*HotSide.PressureDrop.fi*NpassHot*Lv*HotSide.HeatTransfer.Gchannel^2/(HotSide.Properties.Average.rho*Dh*HotSide.HeatTransfer.Phi^0.17);
     461        HotSide.PressureDrop.DPchannel =2*HotSide.PressureDrop.fi*Geometry.NpassHot*Geometry.Lv*HotSide.HeatTransfer.Gchannel^2/(HotSide.Properties.Average.rho*Geometry.Dh*HotSide.HeatTransfer.Phi^0.17);
    365462
    366463"Cold Stream Pressure Drop in Channels"
    367         ColdSide.PressureDrop.DPchannel =2*ColdSide.PressureDrop.fi*NpassCold*Lv*ColdSide.HeatTransfer.Gchannel^2/(ColdSide.Properties.Average.rho*Dh*ColdSide.HeatTransfer.Phi^0.17);
     464        ColdSide.PressureDrop.DPchannel =2*ColdSide.PressureDrop.fi*Geometry.NpassCold*Geometry.Lv*ColdSide.HeatTransfer.Gchannel^2/(ColdSide.Properties.Average.rho*Geometry.Dh*ColdSide.HeatTransfer.Phi^0.17);
    368465
    369466"Hot Stream Total Pressure Drop"
     
    379476        if      HotSide.HeatTransfer.Re < 10
    380477                then
    381                                 HotSide.PressureDrop.fi                 = Kp1(1)/HotSide.HeatTransfer.Re^Kp2(1);
    382                                 ColdSide.PressureDrop.fi        =       Kp1(1)/ColdSide.HeatTransfer.Re^Kp2(1);
     478                                HotSide.PressureDrop.fi                 = Geometry.Kp1(1)/HotSide.HeatTransfer.Re^Geometry.Kp2(1);
     479                                ColdSide.PressureDrop.fi        =       Geometry.Kp1(1)/ColdSide.HeatTransfer.Re^Geometry.Kp2(1);
    383480                else
    384481                                if      HotSide.HeatTransfer.Re < 100
    385482                                        then
    386                                                         HotSide.PressureDrop.fi                 = Kp1(2)/HotSide.HeatTransfer.Re^Kp2(2);
    387                                                         ColdSide.PressureDrop.fi        =       Kp1(2)/ColdSide.HeatTransfer.Re^Kp2(2);
    388                                         else
    389                                                         HotSide.PressureDrop.fi                 = Kp1(3)/HotSide.HeatTransfer.Re^Kp2(3);
    390                                                         ColdSide.PressureDrop.fi        =       Kp1(3)/ColdSide.HeatTransfer.Re^Kp2(3);
     483                                                        HotSide.PressureDrop.fi                 = Geometry.Kp1(2)/HotSide.HeatTransfer.Re^Geometry.Kp2(2);
     484                                                        ColdSide.PressureDrop.fi        =       Geometry.Kp1(2)/ColdSide.HeatTransfer.Re^Geometry.Kp2(2);
     485                                        else
     486                                                        HotSide.PressureDrop.fi                 = Geometry.Kp1(3)/HotSide.HeatTransfer.Re^Geometry.Kp2(3);
     487                                                        ColdSide.PressureDrop.fi        =       Geometry.Kp1(3)/ColdSide.HeatTransfer.Re^Geometry.Kp2(3);
    391488                                end
    392489                       
     
    397494        if      HotSide.HeatTransfer.Re < 15
    398495                then
    399                                 HotSide.PressureDrop.fi                 = Kp1(4)/HotSide.HeatTransfer.Re^Kp2(4);
    400                                 ColdSide.PressureDrop.fi        =       Kp1(4)/ColdSide.HeatTransfer.Re^Kp2(4);
     496                                HotSide.PressureDrop.fi                 = Geometry.Kp1(4)/HotSide.HeatTransfer.Re^Geometry.Kp2(4);
     497                                ColdSide.PressureDrop.fi        =       Geometry.Kp1(4)/ColdSide.HeatTransfer.Re^Geometry.Kp2(4);
    401498                else
    402499                                if      HotSide.HeatTransfer.Re < 300
    403500                                        then
    404                                                         HotSide.PressureDrop.fi                 = Kp1(5)/HotSide.HeatTransfer.Re^Kp2(5);
    405                                                         ColdSide.PressureDrop.fi        =       Kp1(5)/ColdSide.HeatTransfer.Re^Kp2(5);
    406                                         else
    407                                                         HotSide.PressureDrop.fi                 = Kp1(6)/HotSide.HeatTransfer.Re^Kp2(6);
    408                                                         ColdSide.PressureDrop.fi        =       Kp1(6)/ColdSide.HeatTransfer.Re^Kp2(6);
     501                                                        HotSide.PressureDrop.fi                 = Geometry.Kp1(5)/HotSide.HeatTransfer.Re^Geometry.Kp2(5);
     502                                                        ColdSide.PressureDrop.fi        =       Geometry.Kp1(5)/ColdSide.HeatTransfer.Re^Geometry.Kp2(5);
     503                                        else
     504                                                        HotSide.PressureDrop.fi                 = Geometry.Kp1(6)/HotSide.HeatTransfer.Re^Geometry.Kp2(6);
     505                                                        ColdSide.PressureDrop.fi        =       Geometry.Kp1(6)/ColdSide.HeatTransfer.Re^Geometry.Kp2(6);
    409506                                end
    410507                       
     
    415512        if      HotSide.HeatTransfer.Re < 20
    416513                then
    417                                 HotSide.PressureDrop.fi                 = Kp1(7)/HotSide.HeatTransfer.Re^Kp2(7);
    418                                 ColdSide.PressureDrop.fi        =       Kp1(7)/ColdSide.HeatTransfer.Re^Kp2(7);
     514                                HotSide.PressureDrop.fi                 = Geometry.Kp1(7)/HotSide.HeatTransfer.Re^Geometry.Kp2(7);
     515                                ColdSide.PressureDrop.fi        =       Geometry.Kp1(7)/ColdSide.HeatTransfer.Re^Geometry.Kp2(7);
    419516                else
    420517                                if      HotSide.HeatTransfer.Re < 300
    421518                                        then
    422                                                         HotSide.PressureDrop.fi                 = Kp1(8)/HotSide.HeatTransfer.Re^Kp2(8);
    423                                                         ColdSide.PressureDrop.fi        =       Kp1(8)/ColdSide.HeatTransfer.Re^Kp2(8);
    424                                         else
    425                                                         HotSide.PressureDrop.fi                 = Kp1(9)/HotSide.HeatTransfer.Re^Kp2(9);
    426                                                         ColdSide.PressureDrop.fi        =       Kp1(9)/ColdSide.HeatTransfer.Re^Kp2(9);
     519                                                        HotSide.PressureDrop.fi                 = Geometry.Kp1(8)/HotSide.HeatTransfer.Re^Geometry.Kp2(8);
     520                                                        ColdSide.PressureDrop.fi        =       Geometry.Kp1(8)/ColdSide.HeatTransfer.Re^Geometry.Kp2(8);
     521                                        else
     522                                                        HotSide.PressureDrop.fi                 = Geometry.Kp1(9)/HotSide.HeatTransfer.Re^Geometry.Kp2(9);
     523                                                        ColdSide.PressureDrop.fi        =       Geometry.Kp1(9)/ColdSide.HeatTransfer.Re^Geometry.Kp2(9);
    427524                                end
    428525                       
     
    433530        if      HotSide.HeatTransfer.Re < 40
    434531                then
    435                                 HotSide.PressureDrop.fi                 = Kp1(10)/HotSide.HeatTransfer.Re^Kp2(10);
    436                                 ColdSide.PressureDrop.fi        =       Kp1(10)/ColdSide.HeatTransfer.Re^Kp2(10);
     532                                HotSide.PressureDrop.fi                 = Geometry.Kp1(10)/HotSide.HeatTransfer.Re^Geometry.Kp2(10);
     533                                ColdSide.PressureDrop.fi        =       Geometry.Kp1(10)/ColdSide.HeatTransfer.Re^Geometry.Kp2(10);
    437534                else
    438535                                if      HotSide.HeatTransfer.Re < 400
    439536                                        then
    440                                                         HotSide.PressureDrop.fi                 = Kp1(11)/HotSide.HeatTransfer.Re^Kp2(11);
    441                                                         ColdSide.PressureDrop.fi        =       Kp1(11)/ColdSide.HeatTransfer.Re^Kp2(11);
    442                                         else
    443                                                         HotSide.PressureDrop.fi                 = Kp1(12)/HotSide.HeatTransfer.Re^Kp2(12);
    444                                                         ColdSide.PressureDrop.fi        =       Kp1(12)/ColdSide.HeatTransfer.Re^Kp2(12);
     537                                                        HotSide.PressureDrop.fi                 = Geometry.Kp1(11)/HotSide.HeatTransfer.Re^Geometry.Kp2(11);
     538                                                        ColdSide.PressureDrop.fi        =       Geometry.Kp1(11)/ColdSide.HeatTransfer.Re^Geometry.Kp2(11);
     539                                        else
     540                                                        HotSide.PressureDrop.fi                 = Geometry.Kp1(12)/HotSide.HeatTransfer.Re^Geometry.Kp2(12);
     541                                                        ColdSide.PressureDrop.fi        =       Geometry.Kp1(12)/ColdSide.HeatTransfer.Re^Geometry.Kp2(12);
    445542                                end
    446543                       
     
    451548        if      HotSide.HeatTransfer.Re < 50
    452549                then
    453                                 HotSide.PressureDrop.fi                 = Kp1(13)/HotSide.HeatTransfer.Re^Kp2(13);
    454                                 ColdSide.PressureDrop.fi        =       Kp1(13)/ColdSide.HeatTransfer.Re^Kp2(13);
     550                                HotSide.PressureDrop.fi                 = Geometry.Kp1(13)/HotSide.HeatTransfer.Re^Geometry.Kp2(13);
     551                                ColdSide.PressureDrop.fi        =       Geometry.Kp1(13)/ColdSide.HeatTransfer.Re^Geometry.Kp2(13);
    455552                else
    456553                                if      HotSide.HeatTransfer.Re < 500
    457554                                        then
    458                                                         HotSide.PressureDrop.fi                 = Kp1(14)/HotSide.HeatTransfer.Re^Kp2(14);
    459                                                         ColdSide.PressureDrop.fi        =       Kp1(14)/ColdSide.HeatTransfer.Re^Kp2(14);
    460                                         else
    461                                                         HotSide.PressureDrop.fi                 = Kp1(15)/HotSide.HeatTransfer.Re^Kp2(15);
    462                                                         ColdSide.PressureDrop.fi        =       Kp1(15)/ColdSide.HeatTransfer.Re^Kp2(15);
     555                                                        HotSide.PressureDrop.fi                 = Geometry.Kp1(14)/HotSide.HeatTransfer.Re^Geometry.Kp2(14);
     556                                                        ColdSide.PressureDrop.fi        =       Geometry.Kp1(14)/ColdSide.HeatTransfer.Re^Geometry.Kp2(14);
     557                                        else
     558                                                        HotSide.PressureDrop.fi                 = Geometry.Kp1(15)/HotSide.HeatTransfer.Re^Geometry.Kp2(15);
     559                                                        ColdSide.PressureDrop.fi        =       Geometry.Kp1(15)/ColdSide.HeatTransfer.Re^Geometry.Kp2(15);
    463560                                end
    464561                       
     
    473570        if      HotSide.HeatTransfer.Re < 10
    474571                then
    475                                 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Kc1(1)*HotSide.HeatTransfer.Re^Kc2(1))/Dh;
    476                                 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Kc1(1)*ColdSide.HeatTransfer.Re^Kc2(1))/Dh;
    477                 else
    478                                 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Kc1(2)*HotSide.HeatTransfer.Re^Kc2(2))/Dh;
    479                                 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Kc1(2)*ColdSide.HeatTransfer.Re^Kc2(2))/Dh;
     572                                HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(1)*HotSide.HeatTransfer.Re^Geometry.Kc2(1))/Geometry.Dh;
     573                                ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(1)*ColdSide.HeatTransfer.Re^Geometry.Kc2(1))/Geometry.Dh;
     574                else
     575                                HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(2)*HotSide.HeatTransfer.Re^Geometry.Kc2(2))/Geometry.Dh;
     576                                ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(2)*ColdSide.HeatTransfer.Re^Geometry.Kc2(2))/Geometry.Dh;
    480577        end     
    481578       
     
    484581        if      HotSide.HeatTransfer.Re < 10
    485582                then
    486                                 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Kc1(3)*HotSide.HeatTransfer.Re^Kc2(3))/Dh;
    487                                 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Kc1(3)*ColdSide.HeatTransfer.Re^Kc2(3))/Dh;
     583                                HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(3)*HotSide.HeatTransfer.Re^Geometry.Kc2(3))/Geometry.Dh;
     584                                ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(3)*ColdSide.HeatTransfer.Re^Geometry.Kc2(3))/Geometry.Dh;
    488585                else
    489586                                if      HotSide.HeatTransfer.Re < 100
    490587                                        then
    491                                                         HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Kc1(4)*HotSide.HeatTransfer.Re^Kc2(4))/Dh;
    492                                                         ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Kc1(4)*ColdSide.HeatTransfer.Re^Kc2(4))/Dh;
    493                                         else
    494                                                         HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Kc1(5)*HotSide.HeatTransfer.Re^Kc2(5))/Dh;
    495                                                         ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Kc1(5)*ColdSide.HeatTransfer.Re^Kc2(5))/Dh;
     588                                                        HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(4)*HotSide.HeatTransfer.Re^Geometry.Kc2(4))/Geometry.Dh;
     589                                                        ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(4)*ColdSide.HeatTransfer.Re^Geometry.Kc2(4))/Geometry.Dh;
     590                                        else
     591                                                        HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(5)*HotSide.HeatTransfer.Re^Geometry.Kc2(5))/Geometry.Dh;
     592                                                        ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(5)*ColdSide.HeatTransfer.Re^Geometry.Kc2(5))/Geometry.Dh;
    496593                                end
    497594        end
     
    501598        if      HotSide.HeatTransfer.Re < 20
    502599                then
    503                                 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Kc1(6)*HotSide.HeatTransfer.Re^Kc2(6))/Dh;
    504                                 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Kc1(6)*ColdSide.HeatTransfer.Re^Kc2(6))/Dh;
     600                                HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(6)*HotSide.HeatTransfer.Re^Geometry.Kc2(6))/Geometry.Dh;
     601                                ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(6)*ColdSide.HeatTransfer.Re^Geometry.Kc2(6))/Geometry.Dh;
    505602                else
    506603                                if      HotSide.HeatTransfer.Re < 300
    507604                                        then
    508                                                         HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Kc1(7)*HotSide.HeatTransfer.Re^Kc2(7))/Dh;
    509                                                         ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Kc1(7)*ColdSide.HeatTransfer.Re^Kc2(7))/Dh;
    510                                         else
    511                                                         HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Kc1(8)*HotSide.HeatTransfer.Re^Kc2(8))/Dh;
    512                                                         ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Kc1(8)*ColdSide.HeatTransfer.Re^Kc2(8))/Dh;
     605                                                        HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(7)*HotSide.HeatTransfer.Re^Geometry.Kc2(7))/Geometry.Dh;
     606                                                        ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(7)*ColdSide.HeatTransfer.Re^Geometry.Kc2(7))/Geometry.Dh;
     607                                        else
     608                                                        HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(8)*HotSide.HeatTransfer.Re^Geometry.Kc2(8))/Geometry.Dh;
     609                                                        ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(8)*ColdSide.HeatTransfer.Re^Geometry.Kc2(8))/Geometry.Dh;
    513610                                end
    514611        end
     
    518615        if      HotSide.HeatTransfer.Re < 20
    519616                then
    520                                 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Kc1(9)*HotSide.HeatTransfer.Re^Kc2(9))/Dh;
    521                                 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Kc1(9)*ColdSide.HeatTransfer.Re^Kc2(9))/Dh;
     617                                HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(9)*HotSide.HeatTransfer.Re^Geometry.Kc2(9))/Geometry.Dh;
     618                                ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(9)*ColdSide.HeatTransfer.Re^Geometry.Kc2(9))/Geometry.Dh;
    522619                else
    523620                                if      HotSide.HeatTransfer.Re < 400
    524621                                        then
    525                                                         HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Kc1(10)*HotSide.HeatTransfer.Re^Kc2(10))/Dh;
    526                                                         ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Kc1(10)*ColdSide.HeatTransfer.Re^Kc2(10))/Dh;
    527                                         else
    528                                                         HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Kc1(11)*HotSide.HeatTransfer.Re^Kc2(11))/Dh;
    529                                                         ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Kc1(11)*ColdSide.HeatTransfer.Re^Kc2(11))/Dh;
     622                                                        HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(10)*HotSide.HeatTransfer.Re^Geometry.Kc2(10))/Geometry.Dh;
     623                                                        ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(10)*ColdSide.HeatTransfer.Re^Geometry.Kc2(10))/Geometry.Dh;
     624                                        else
     625                                                        HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(11)*HotSide.HeatTransfer.Re^Geometry.Kc2(11))/Geometry.Dh;
     626                                                        ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(11)*ColdSide.HeatTransfer.Re^Geometry.Kc2(11))/Geometry.Dh;
    530627                                end
    531628        end
     
    535632        if      HotSide.HeatTransfer.Re < 20
    536633                then
    537                                 HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Kc1(12)*HotSide.HeatTransfer.Re^Kc2(12))/Dh;
    538                                 ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Kc1(12)*ColdSide.HeatTransfer.Re^Kc2(12))/Dh;
     634                                HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(12)*HotSide.HeatTransfer.Re^Geometry.Kc2(12))/Geometry.Dh;
     635                                ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(12)*ColdSide.HeatTransfer.Re^Geometry.Kc2(12))/Geometry.Dh;
    539636                else
    540637                                if      HotSide.HeatTransfer.Re < 500
    541638                                        then
    542                                                         HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Kc1(13)*HotSide.HeatTransfer.Re^Kc2(13))/Dh;
    543                                                         ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Kc1(13)*ColdSide.HeatTransfer.Re^Kc2(13))/Dh;
    544                                         else
    545                                                         HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Kc1(14)*HotSide.HeatTransfer.Re^Kc2(14))/Dh;
    546                                                         ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Kc1(14)*ColdSide.HeatTransfer.Re^Kc2(14))/Dh;
     639                                                        HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(13)*HotSide.HeatTransfer.Re^Geometry.Kc2(13))/Geometry.Dh;
     640                                                        ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(13)*ColdSide.HeatTransfer.Re^Geometry.Kc2(13))/Geometry.Dh;
     641                                        else
     642                                                        HotSide.HeatTransfer.hcoeff=(HotSide.Properties.Average.K*HotSide.HeatTransfer.PR^(1/3)*HotSide.HeatTransfer.Phi^0.17*Geometry.Kc1(14)*HotSide.HeatTransfer.Re^Geometry.Kc2(14))/Geometry.Dh;
     643                                                        ColdSide.HeatTransfer.hcoeff =(ColdSide.Properties.Average.K*ColdSide.HeatTransfer.PR^(1/3)*ColdSide.HeatTransfer.Phi^0.17*Geometry.Kc1(14)*ColdSide.HeatTransfer.Re^Geometry.Kc2(14))/Geometry.Dh;
    547644                                end
    548645        end
     
    557654
    558655"Hot Stream Velocity in Ports"
    559         HotSide.PressureDrop.Vports =HotSide.Properties.Inlet.Fw/(Aports*HotSide.Properties.Inlet.rho);
     656        HotSide.PressureDrop.Vports =HotSide.Properties.Inlet.Fw/(Geometry.Aports*HotSide.Properties.Inlet.rho);
    560657
    561658"Cold Stream Velocity in Ports"
    562         ColdSide.PressureDrop.Vports =ColdSide.Properties.Inlet.Fw/(Aports*ColdSide.Properties.Inlet.rho);
     659        ColdSide.PressureDrop.Vports =ColdSide.Properties.Inlet.Fw/(Geometry.Aports*ColdSide.Properties.Inlet.rho);
    563660
    564661"Hot Stream Reynolds Number"
    565         HotSide.HeatTransfer.Re =Dh*HotSide.HeatTransfer.Gchannel/HotSide.Properties.Average.Mu;
     662        HotSide.HeatTransfer.Re =Geometry.Dh*HotSide.HeatTransfer.Gchannel/HotSide.Properties.Average.Mu;
    566663
    567664"Cold Stream Reynolds Number"
    568         ColdSide.HeatTransfer.Re =Dh*ColdSide.HeatTransfer.Gchannel/ColdSide.Properties.Average.Mu;
     665        ColdSide.HeatTransfer.Re =Geometry.Dh*ColdSide.HeatTransfer.Gchannel/ColdSide.Properties.Average.Mu;
    569666
    570667"Hot Stream Prandtl Number"
     
    587684       
    588685"Overall Heat Transfer Coefficient Clean"
    589         Thermal.Uc/HotSide.HeatTransfer.hcoeff +Thermal.Uc*pt/Kwall+Thermal.Uc/ColdSide.HeatTransfer.hcoeff=1;
     686        Thermal.Uc/HotSide.HeatTransfer.hcoeff +Thermal.Uc*Geometry.pt/Geometry.Kwall+Thermal.Uc/ColdSide.HeatTransfer.hcoeff=1;
    590687
    591688"Overall Heat Transfer Coefficient Dirty"
    592         Thermal.Ud*(1/HotSide.HeatTransfer.hcoeff +pt/Kwall+1/ColdSide.HeatTransfer.hcoeff + Rfc + Rfh)=1;
     689        Thermal.Ud*(1/HotSide.HeatTransfer.hcoeff +Geometry.pt/Geometry.Kwall+1/ColdSide.HeatTransfer.hcoeff + Geometry.Rfc + Geometry.Rfh)=1;
    593690
    594691 "Duty"
     
    614711
    615712"Number of Units Transference for Hot Side"
    616         HotSide.HeatTransfer.NTU*HotSide.HeatTransfer.WCp = Thermal.Ud*Atotal;
     713        HotSide.HeatTransfer.NTU*HotSide.HeatTransfer.WCp = Thermal.Ud*Geometry.Atotal;
    617714
    618715"Number of Units Transference for Cold Side"
    619         ColdSide.HeatTransfer.NTU*ColdSide.HeatTransfer.WCp = Thermal.Ud*Atotal;
     716        ColdSide.HeatTransfer.NTU*ColdSide.HeatTransfer.WCp = Thermal.Ud*Geometry.Atotal;
    620717
    621718if  Thermal.Eft >= 1  #To be Fixed: Effectiveness in true counter flow !
  • trunk/eml/streams.mso

    r501 r523  
    3838        T as temperature                (Brief = "Stream Temperature");
    3939        P as pressure                   (Brief = "Stream Pressure");
    40         z(NComp) as fraction    (Brief = "Stream Molar Fraction");
    4140        h as enth_mol                   (Brief = "Stream Enthalpy");
    4241        v as fraction                   (Brief = "Vapourization fraction");
     42        z(NComp) as fraction    (Brief = "Stream Molar Fraction");
    4343end
    4444
Note: See TracChangeset for help on using the changeset viewer.