Ignore:
Timestamp:
Sep 27, 2008, 1:34:26 AM (15 years ago)
Author:
Argimiro Resende Secchi
Message:

Fix cost model due to Heatex changes.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/eml/costs/HeatExchangerSimplified_cost.mso

    r491 r639  
    4545        EQUATIONS
    4646        "Average pressure"
    47         Pmax = max( [HotSide.Properties.Average.P , ColdSide.Properties.Average.P] );
    48        
    49         "Capital Cost"
    50         Ce = Cb*Fd*Fp*Fm;
    51        
    52         "Basic Cost"
    53         Cb = 'US$'*exp(Cost(1,1) + Cost(1,2)*ln(Details.A/'m^2') + Cost(1,3)*(ln(Details.A/'m^2'))^2);
    54 
    55         "Cost Factor based on the type of the heat exchanger"
    56         Fd = exp(Cost(2,1) + Cost(2,2)*ln(Details.A/'m^2') + Cost(2,3)*ln(Details.A/'m^2'));
    57 
    58         if Pmax <= 700 * 'kPa' then
    59                 "Cost Factor based on the project pressure" # verificar
    60                 Fp = Cost(5,1) + Cost(5,2)*ln(Details.A/'m^2');
    61         else if 700 * 'kPa' < Pmax and Pmax < 2100 * 'kPa' then
    62                 Fp = Cost(5,1) + Cost(5,2)*ln(Details.A/'m^2');
    63         else if 2100 * 'kPa' < Pmax and Pmax < 4200 * 'kPa' then
    64                 Fp = Cost(6,1) + Cost(6,2)*ln(Details.A/'m^2');
    65         else if 4200 * 'kPa' < Pmax and Pmax < 6200 * 'kPa' then
    66                 Fp = Cost(7,1) + Cost(7,2)*ln(Details.A/'m^2');
    67         else
    68                 Fp = Cost(7,1) + Cost(7,2)*ln(Details.A/'m^2'); # verificar
    69         end
    70         end
    71         end
    72         end
    73        
    74         "Cost Factor based on the construction material"
    75         Fm = Cost(8,1) + Cost(8,2)*ln(Details.A/'m^2');
     47        Pmax = max( [OutletHot.P , OutletCold.P] );
     48       
     49        "Capital Cost"
     50        Ce = Cb*Fd*Fp*Fm;
     51       
     52        "Basic Cost"
     53        Cb = 'US$'*exp(Cost(1,1) + Cost(1,2)*ln(A/'m^2') + Cost(1,3)*(ln(A/'m^2'))^2);
     54
     55        "Cost Factor based on the type of the heat exchanger"
     56        Fd = exp(Cost(2,1) + Cost(2,2)*ln(A/'m^2') + Cost(2,3)*ln(A/'m^2'));
     57
     58        if Pmax <= 700 * 'kPa' then
     59                "Cost Factor based on the project pressure" # verificar
     60                Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
     61        else if 700 * 'kPa' < Pmax and Pmax < 2100 * 'kPa' then
     62                Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
     63        else if 2100 * 'kPa' < Pmax and Pmax < 4200 * 'kPa' then
     64                Fp = Cost(6,1) + Cost(6,2)*ln(A/'m^2');
     65        else if 4200 * 'kPa' < Pmax and Pmax < 6200 * 'kPa' then
     66                Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2');
     67        else
     68                Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2'); # verificar
     69        end
     70        end
     71        end
     72        end
     73       
     74        "Cost Factor based on the construction material"
     75        Fm = Cost(8,1) + Cost(8,2)*ln(A/'m^2');
    7676       
    7777end
    … …  
    104104        EQUATIONS
    105105        "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');
     106        Pmax = max( [OutletHot.P , OutletCold.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(A/'m^2') + Cost(1,3)*(ln(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(A/'m^2') + Cost(2,3)*ln(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(A/'m^2');
     120        else if 700 * 'kPa' < Pmax and Pmax < 2100 * 'kPa' then
     121                Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
     122        else if 2100 * 'kPa' < Pmax and Pmax < 4200 * 'kPa' then
     123                Fp = Cost(6,1) + Cost(6,2)*ln(A/'m^2');
     124        else if 4200 * 'kPa' < Pmax and Pmax < 6200 * 'kPa' then
     125                Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2');
     126        else
     127                Fp = Cost(7,1) + Cost(7,2)*ln(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(A/'m^2');
    135135       
    136136end
    … …  
    159159        EQUATIONS
    160160        "Average pressure"
    161         Pmax = max( [HotSide.Properties.Average.P , ColdSide.Properties.Average.P] );
    162        
    163         "Capital Cost"
    164         Ce = Cb*Fd*Fp*Fm;
    165        
    166         "Basic Cost"
    167         Cb = 'US$'*exp(Cost(1,1) + Cost(1,2)*ln(Details.A/'m^2') + Cost(1,3)*(ln(Details.A/'m^2'))^2);
    168 
    169         "Cost Factor based on the type of the heat exchanger"
    170         Fd = exp(Cost(2,1) + Cost(2,2)*ln(Details.A/'m^2') + Cost(2,3)*ln(Details.A/'m^2'));
    171 
    172         if Pmax <= 700 * 'kPa' then
    173                 "Cost Factor based on the project pressure" # verificar
    174                 Fp = Cost(5,1) + Cost(5,2)*ln(Details.A/'m^2');
    175         else if 700 * 'kPa' < Pmax and Pmax < 2100 * 'kPa' then
    176                 Fp = Cost(5,1) + Cost(5,2)*ln(Details.A/'m^2');
    177         else if 2100 * 'kPa' < Pmax and Pmax < 4200 * 'kPa' then
    178                 Fp = Cost(6,1) + Cost(6,2)*ln(Details.A/'m^2');
    179         else if 4200 * 'kPa' < Pmax and Pmax < 6200 * 'kPa' then
    180                 Fp = Cost(7,1) + Cost(7,2)*ln(Details.A/'m^2');
    181         else
    182                 Fp = Cost(7,1) + Cost(7,2)*ln(Details.A/'m^2'); # verificar
    183         end
    184         end
    185         end
    186         end
    187        
    188         "Cost Factor based on the construction material"
    189         Fm = Cost(8,1) + Cost(8,2)*ln(Details.A/'m^2');
     161        Pmax = max( [OutletHot.P , OutletCold.P] );
     162       
     163        "Capital Cost"
     164        Ce = Cb*Fd*Fp*Fm;
     165       
     166        "Basic Cost"
     167        Cb = 'US$'*exp(Cost(1,1) + Cost(1,2)*ln(A/'m^2') + Cost(1,3)*(ln(A/'m^2'))^2);
     168
     169        "Cost Factor based on the type of the heat exchanger"
     170        Fd = exp(Cost(2,1) + Cost(2,2)*ln(A/'m^2') + Cost(2,3)*ln(A/'m^2'));
     171
     172        if Pmax <= 700 * 'kPa' then
     173                "Cost Factor based on the project pressure" # verificar
     174                Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
     175        else if 700 * 'kPa' < Pmax and Pmax < 2100 * 'kPa' then
     176                Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
     177        else if 2100 * 'kPa' < Pmax and Pmax < 4200 * 'kPa' then
     178                Fp = Cost(6,1) + Cost(6,2)*ln(A/'m^2');
     179        else if 4200 * 'kPa' < Pmax and Pmax < 6200 * 'kPa' then
     180                Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2');
     181        else
     182                Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2'); # verificar
     183        end
     184        end
     185        end
     186        end
     187       
     188        "Cost Factor based on the construction material"
     189        Fm = Cost(8,1) + Cost(8,2)*ln(A/'m^2');
    190190       
    191191end
    … …  
    218218        EQUATIONS
    219219        "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        
    250 end
    251 
     220        Pmax = max( [OutletHot.P , OutletCold.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(A/'m^2') + Cost(1,3)*(ln(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(A/'m^2') + Cost(2,3)*ln(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(A/'m^2');
     234        else if 700 * 'kPa' < Pmax and Pmax < 2100 * 'kPa' then
     235                Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
     236        else if 2100 * 'kPa' < Pmax and Pmax < 4200 * 'kPa' then
     237                Fp = Cost(6,1) + Cost(6,2)*ln(A/'m^2');
     238        else if 4200 * 'kPa' < Pmax and Pmax < 6200 * 'kPa' then
     239                Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2');
     240        else
     241                Fp = Cost(7,1) + Cost(7,2)*ln(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(A/'m^2');
     249       
     250end
     251
Note: See TracChangeset for help on using the changeset viewer.