Changeset 491 for trunk/eml


Ignore:
Timestamp:
Apr 3, 2008, 9:22:56 PM (14 years ago)
Author:
gerson bicca
Message:

fixed heat exchanger cost model

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/eml/costs/HeatExchangerSimplified_cost.mso

    r331 r491  
    1818*------------------------------------------------------------------*#
    1919
    20 using "heat_exchangers/HeatExchangerSimplified";
     20using "heat_exchangers/Heatex";
    2121
    2222# This model is valid if 14m^2 < A < 1100 m^2
    2323
    24 Model HeatExchanger_LMTD_cost as HeatExchanger_LMTD
     24Model HeatExchanger_LMTD_cost as Heatex_LMTD
    2525        ATTRIBUTES
    2626        Pallete         = true;
     
    7777end
    7878
    79 Model Shell_and_tubes_LMTD_cost as Shell_and_Tubes_LMTD
     79Model Shell_and_tubes_LMTD_cost as Heatex_LMTD
    8080        ATTRIBUTES
    8181        Pallete         = true;
     
    9898        Pmax  as pressure       (Brief="Average  Pressure");
    9999
    100         EQUATIONS
    101         "Average pressure"
    102         Pmax = max( [HotSide.Properties.Average.P , ColdSide.Properties.Average.P] );
    103        
    104         "Capital Cost"
    105         Ce = Cb*Fd*Fp*Fm;
    106        
    107         "Basic Cost"
    108         Cb = 'US$'*exp(Cost(1,1) + Cost(1,2)*ln(Details.A/'m^2') + Cost(1,3)*(ln(Details.A/'m^2'))^2);
    109 
    110         "Cost Factor based on the type of the heat exchanger"
    111         Fd = exp(Cost(2,1) + Cost(2,2)*ln(Details.A/'m^2') + Cost(2,3)*ln(Details.A/'m^2'));
    112 
    113         if Pmax <= 700 * 'kPa' then
    114                 "Cost Factor based on the project pressure" # verificar
    115                 Fp = Cost(5,1) + Cost(5,2)*ln(Details.A/'m^2');
    116         else if 700 * 'kPa' < Pmax and Pmax < 2100 * 'kPa' then
    117                 Fp = Cost(5,1) + Cost(5,2)*ln(Details.A/'m^2');
    118         else if 2100 * 'kPa' < Pmax and Pmax < 4200 * 'kPa' then
    119                 Fp = Cost(6,1) + Cost(6,2)*ln(Details.A/'m^2');
    120         else if 4200 * 'kPa' < Pmax and Pmax < 6200 * 'kPa' then
    121                 Fp = Cost(7,1) + Cost(7,2)*ln(Details.A/'m^2');
    122         else
    123                 Fp = Cost(7,1) + Cost(7,2)*ln(Details.A/'m^2'); # verificar
    124         end
    125         end
    126         end
    127         end
    128        
    129         "Cost Factor based on the construction material"
    130         Fm = Cost(8,1) + Cost(8,2)*ln(Details.A/'m^2');
    131        
    132 end
    133 
    134 Model HeatExchanger_NTU_cost as HeatExchanger_NTU
     100SET
     101
     102ExchangerType = "Shell and Tube";
     103
     104        EQUATIONS
     105        "Average pressure"
     106        Pmax = max( [HotSide.Properties.Average.P , ColdSide.Properties.Average.P] );
     107       
     108        "Capital Cost"
     109        Ce = Cb*Fd*Fp*Fm;
     110       
     111        "Basic Cost"
     112        Cb = 'US$'*exp(Cost(1,1) + Cost(1,2)*ln(Details.A/'m^2') + Cost(1,3)*(ln(Details.A/'m^2'))^2);
     113
     114        "Cost Factor based on the type of the heat exchanger"
     115        Fd = exp(Cost(2,1) + Cost(2,2)*ln(Details.A/'m^2') + Cost(2,3)*ln(Details.A/'m^2'));
     116
     117        if Pmax <= 700 * 'kPa' then
     118                "Cost Factor based on the project pressure" # verificar
     119                Fp = Cost(5,1) + Cost(5,2)*ln(Details.A/'m^2');
     120        else if 700 * 'kPa' < Pmax and Pmax < 2100 * 'kPa' then
     121                Fp = Cost(5,1) + Cost(5,2)*ln(Details.A/'m^2');
     122        else if 2100 * 'kPa' < Pmax and Pmax < 4200 * 'kPa' then
     123                Fp = Cost(6,1) + Cost(6,2)*ln(Details.A/'m^2');
     124        else if 4200 * 'kPa' < Pmax and Pmax < 6200 * 'kPa' then
     125                Fp = Cost(7,1) + Cost(7,2)*ln(Details.A/'m^2');
     126        else
     127                Fp = Cost(7,1) + Cost(7,2)*ln(Details.A/'m^2'); # verificar
     128        end
     129        end
     130        end
     131        end
     132       
     133        "Cost Factor based on the construction material"
     134        Fm = Cost(8,1) + Cost(8,2)*ln(Details.A/'m^2');
     135       
     136end
     137
     138Model HeatExchanger_NTU_cost as Heatex_NTU
    135139        ATTRIBUTES
    136140        Pallete         = true;
     
    187191end
    188192
    189 Model Shell_and_tubes_NTU_cost as Shell_and_Tubes_NTU
     193Model Shell_and_tubes_NTU_cost as Heatex_NTU
    190194        ATTRIBUTES
    191195        Pallete         = true;
     
    208212        Pmax  as pressure       (Brief="Average  Pressure");
    209213
    210         EQUATIONS
    211         "Average pressure"
    212         Pmax = max( [HotSide.Properties.Average.P , ColdSide.Properties.Average.P] );
    213        
    214         "Capital Cost"
    215         Ce = Cb*Fd*Fp*Fm;
    216        
    217         "Basic Cost"
    218         Cb = 'US$'*exp(Cost(1,1) + Cost(1,2)*ln(Details.A/'m^2') + Cost(1,3)*(ln(Details.A/'m^2'))^2);
    219 
    220         "Cost Factor based on the type of the heat exchanger"
    221         Fd = exp(Cost(2,1) + Cost(2,2)*ln(Details.A/'m^2') + Cost(2,3)*ln(Details.A/'m^2'));
    222 
    223         if Pmax <= 700 * 'kPa' then
    224                 "Cost Factor based on the project pressure" # verificar
    225                 Fp = Cost(5,1) + Cost(5,2)*ln(Details.A/'m^2');
    226         else if 700 * 'kPa' < Pmax and Pmax < 2100 * 'kPa' then
    227                 Fp = Cost(5,1) + Cost(5,2)*ln(Details.A/'m^2');
    228         else if 2100 * 'kPa' < Pmax and Pmax < 4200 * 'kPa' then
    229                 Fp = Cost(6,1) + Cost(6,2)*ln(Details.A/'m^2');
    230         else if 4200 * 'kPa' < Pmax and Pmax < 6200 * 'kPa' then
    231                 Fp = Cost(7,1) + Cost(7,2)*ln(Details.A/'m^2');
    232         else
    233                 Fp = Cost(7,1) + Cost(7,2)*ln(Details.A/'m^2'); # verificar
    234         end
    235         end
    236         end
    237         end
    238        
    239         "Cost Factor based on the construction material"
    240         Fm = Cost(8,1) + Cost(8,2)*ln(Details.A/'m^2');
    241        
    242 end
    243 
     214SET
     215
     216ExchangerType = "Shell and Tube";
     217
     218        EQUATIONS
     219        "Average pressure"
     220        Pmax = max( [HotSide.Properties.Average.P , ColdSide.Properties.Average.P] );
     221       
     222        "Capital Cost"
     223        Ce = Cb*Fd*Fp*Fm;
     224       
     225        "Basic Cost"
     226        Cb = 'US$'*exp(Cost(1,1) + Cost(1,2)*ln(Details.A/'m^2') + Cost(1,3)*(ln(Details.A/'m^2'))^2);
     227
     228        "Cost Factor based on the type of the heat exchanger"
     229        Fd = exp(Cost(2,1) + Cost(2,2)*ln(Details.A/'m^2') + Cost(2,3)*ln(Details.A/'m^2'));
     230
     231        if Pmax <= 700 * 'kPa' then
     232                "Cost Factor based on the project pressure" # verificar
     233                Fp = Cost(5,1) + Cost(5,2)*ln(Details.A/'m^2');
     234        else if 700 * 'kPa' < Pmax and Pmax < 2100 * 'kPa' then
     235                Fp = Cost(5,1) + Cost(5,2)*ln(Details.A/'m^2');
     236        else if 2100 * 'kPa' < Pmax and Pmax < 4200 * 'kPa' then
     237                Fp = Cost(6,1) + Cost(6,2)*ln(Details.A/'m^2');
     238        else if 4200 * 'kPa' < Pmax and Pmax < 6200 * 'kPa' then
     239                Fp = Cost(7,1) + Cost(7,2)*ln(Details.A/'m^2');
     240        else
     241                Fp = Cost(7,1) + Cost(7,2)*ln(Details.A/'m^2'); # verificar
     242        end
     243        end
     244        end
     245        end
     246       
     247        "Cost Factor based on the construction material"
     248        Fm = Cost(8,1) + Cost(8,2)*ln(Details.A/'m^2');
     249       
     250end
     251
Note: See TracChangeset for help on using the changeset viewer.