Changeset 565


Ignore:
Timestamp:
Jul 23, 2008, 8:30:02 PM (15 years ago)
Author:
gerson bicca
Message:

updated

Location:
branches/gui
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • branches/gui/eml/costs/HeatExchangerSimplified_cost.mso

    r491 r565  
    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');
    76        
     47        Pmax = max( [(0.5*InletHot.P+0.5*OutletHot.P) , (0.5*InletCold.P+0.5*OutletCold.P)] );
     48       
     49
     50        "Capital Cost"
     51        Ce = Cb*Fd*Fp*Fm;
     52       
     53        "Basic Cost"
     54        Cb = 'US$'*exp(Cost(1,1) + Cost(1,2)*ln(A/'m^2') + Cost(1,3)*(ln(A/'m^2'))^2);
     55
     56        "Cost Factor based on the type of the heat exchanger"
     57        Fd = exp(Cost(2,1) + Cost(2,2)*ln(A/'m^2') + Cost(2,3)*ln(A/'m^2'));
     58
     59        if Pmax <= 700 * 'kPa' then
     60                "Cost Factor based on the project pressure" # verificar
     61                Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
     62        else if 700 * 'kPa' < Pmax and Pmax < 2100 * 'kPa' then
     63                Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
     64        else if 2100 * 'kPa' < Pmax and Pmax < 4200 * 'kPa' then
     65                Fp = Cost(6,1) + Cost(6,2)*ln(A/'m^2');
     66        else if 4200 * 'kPa' < Pmax and Pmax < 6200 * 'kPa' then
     67                Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2');
     68        else
     69                Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2'); # verificar
     70        end
     71        end
     72        end
     73        end
     74       
     75        "Cost Factor based on the construction material"
     76        Fm = Cost(8,1) + Cost(8,2)*ln(A/'m^2');
     77
    7778end
    7879
    … …  
    104105        EQUATIONS
    105106        "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');
     107        Pmax = max( [(0.5*InletHot.P+0.5*OutletHot.P) , (0.5*InletCold.P+0.5*OutletCold.P)] );
     108       
     109        "Capital Cost"
     110        Ce = Cb*Fd*Fp*Fm;
     111       
     112        "Basic Cost"
     113        Cb = 'US$'*exp(Cost(1,1) + Cost(1,2)*ln(A/'m^2') + Cost(1,3)*(ln(A/'m^2'))^2);
     114
     115        "Cost Factor based on the type of the heat exchanger"
     116        Fd = exp(Cost(2,1) + Cost(2,2)*ln(A/'m^2') + Cost(2,3)*ln(A/'m^2'));
     117
     118        if Pmax <= 700 * 'kPa' then
     119                "Cost Factor based on the project pressure" # verificar
     120                Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
     121        else if 700 * 'kPa' < Pmax and Pmax < 2100 * 'kPa' then
     122                Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
     123        else if 2100 * 'kPa' < Pmax and Pmax < 4200 * 'kPa' then
     124                Fp = Cost(6,1) + Cost(6,2)*ln(A/'m^2');
     125        else if 4200 * 'kPa' < Pmax and Pmax < 6200 * 'kPa' then
     126                Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2');
     127        else
     128                Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2'); # verificar
     129        end
     130        end
     131        end
     132        end
     133       
     134        "Cost Factor based on the construction material"
     135        Fm = Cost(8,1) + Cost(8,2)*ln(A/'m^2');
    135136       
    136137end
    … …  
    159160        EQUATIONS
    160161        "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');
     162        Pmax = max( [(0.5*InletHot.P+0.5*OutletHot.P) , (0.5*InletCold.P+0.5*OutletCold.P)] );
     163       
     164        "Capital Cost"
     165        Ce = Cb*Fd*Fp*Fm;
     166       
     167        "Basic Cost"
     168        Cb = 'US$'*exp(Cost(1,1) + Cost(1,2)*ln(A/'m^2') + Cost(1,3)*(ln(A/'m^2'))^2);
     169
     170        "Cost Factor based on the type of the heat exchanger"
     171        Fd = exp(Cost(2,1) + Cost(2,2)*ln(A/'m^2') + Cost(2,3)*ln(A/'m^2'));
     172
     173        if Pmax <= 700 * 'kPa' then
     174                "Cost Factor based on the project pressure" # verificar
     175                Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
     176        else if 700 * 'kPa' < Pmax and Pmax < 2100 * 'kPa' then
     177                Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
     178        else if 2100 * 'kPa' < Pmax and Pmax < 4200 * 'kPa' then
     179                Fp = Cost(6,1) + Cost(6,2)*ln(A/'m^2');
     180        else if 4200 * 'kPa' < Pmax and Pmax < 6200 * 'kPa' then
     181                Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2');
     182        else
     183                Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2'); # verificar
     184        end
     185        end
     186        end
     187        end
     188       
     189        "Cost Factor based on the construction material"
     190        Fm = Cost(8,1) + Cost(8,2)*ln(A/'m^2');
    190191       
    191192end
    … …  
    218219        EQUATIONS
    219220        "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 
     221                Pmax = max( [(0.5*InletHot.P+0.5*OutletHot.P) , (0.5*InletCold.P+0.5*OutletCold.P)] );
     222       
     223        "Capital Cost"
     224        Ce = Cb*Fd*Fp*Fm;
     225       
     226        "Basic Cost"
     227        Cb = 'US$'*exp(Cost(1,1) + Cost(1,2)*ln(A/'m^2') + Cost(1,3)*(ln(A/'m^2'))^2);
     228
     229        "Cost Factor based on the type of the heat exchanger"
     230        Fd = exp(Cost(2,1) + Cost(2,2)*ln(A/'m^2') + Cost(2,3)*ln(A/'m^2'));
     231
     232        if Pmax <= 700 * 'kPa' then
     233                "Cost Factor based on the project pressure" # verificar
     234                Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
     235        else if 700 * 'kPa' < Pmax and Pmax < 2100 * 'kPa' then
     236                Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
     237        else if 2100 * 'kPa' < Pmax and Pmax < 4200 * 'kPa' then
     238                Fp = Cost(6,1) + Cost(6,2)*ln(A/'m^2');
     239        else if 4200 * 'kPa' < Pmax and Pmax < 6200 * 'kPa' then
     240                Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2');
     241        else
     242                Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2'); # verificar
     243        end
     244        end
     245        end
     246        end
     247       
     248        "Cost Factor based on the construction material"
     249        Fm = Cost(8,1) + Cost(8,2)*ln(A/'m^2');
     250
     251end
     252

  • branches/gui/sample/costs/sample_HeatExchangerSimplified_cost.mso

    r491 r565  
    6767SPECIFY
    6868
    69         exchanger.Details.A                                      = 139.35*'m^2';
    70         exchanger.Details.Ud                                     = 0.75*'W/(m^2*K)';
    71         exchanger.Details.Uc                                     = 0.75*'W/(m^2*K)';
    72        
    73         exchanger.HotSide.PressureDrop.Pdrop  = 0*'atm';
    74         exchanger.ColdSide.PressureDrop.Pdrop = 0*'atm';
     69        exchanger.A                                      = 139.35*'m^2';
     70        exchanger.U                                      = 0.75*'W/(m^2*K)';
     71       
     72        exchanger.PdropHotSide  = 0*'atm';
     73        exchanger.PdropColdSide = 0*'atm';
    7574       
    7675        streamhot.Outlet.F                                      = 1000    * 'kmol/h';
    … …  
    131130SPECIFY
    132131
    133         exchanger.Details.A                                      = 65.031*'m^2';
    134         exchanger.Details.Ud                                     = 0.75*'W/(m^2*K)';
    135         exchanger.Details.Uc                                     = 0.75*'W/(m^2*K)';
    136        
    137         exchanger.HotSide.PressureDrop.Pdrop  = 0*'atm';
    138         exchanger.ColdSide.PressureDrop.Pdrop = 0*'atm';
     132        exchanger.A                                      = 65.031*'m^2';
     133        exchanger.U                                      = 0.75*'W/(m^2*K)';
     134       
     135        exchanger.PdropHotSide  = 0*'atm';
     136        exchanger.PdropColdSide = 0*'atm';
    139137       
    140138        streamhot.Outlet.F                                      = 1000    * 'kmol/h';
    … …  
    201199SPECIFY
    202200
    203 exchanger.Details.Ud                                            = 210*'W/(m^2*K)';
    204 exchanger.Details.Uc                                            = 210*'W/(m^2*K)';
    205 
    206 exchanger.HotSide.PressureDrop.Pdrop    = 0.1*'kPa';
    207 exchanger.ColdSide.PressureDrop.Pdrop           = 0.1*'kPa';
     201exchanger.U                                             = 210*'W/(m^2*K)';
     202
     203exchanger.PdropHotSide  = 0.1*'kPa';
     204exchanger.PdropColdSide         = 0.1*'kPa';
    208205
    209206streamhot.Outlet.F              = 36    * 'kmol/h';
    … …  
    217214streamcold.Outlet.z             = [0.1,0.7, 0.2];
    218215
    219 exchanger.Details.A     = 15*'m^2';
     216exchanger.A     = 15*'m^2';
    220217       
    221218OPTIONS
    222219
    223220 Dynamic        = false;
    224  GuessFile   = "GuessShellandTubes";
    225221
    226222end
    … …  
    269265SPECIFY
    270266
    271 exchanger.Details.Ud                                            = 210*'W/(m^2*K)';
    272 exchanger.Details.Uc                                            = 210*'W/(m^2*K)';
    273 
    274 exchanger.HotSide.PressureDrop.Pdrop    = 0.1*'kPa';
    275 exchanger.ColdSide.PressureDrop.Pdrop           = 0.1*'kPa';
     267exchanger.U                                             = 210*'W/(m^2*K)';
     268
     269
     270exchanger.PdropHotSide  = 0.1*'kPa';
     271exchanger.PdropColdSide         = 0.1*'kPa';
    276272
    277273streamhot.Outlet.F              = 36    * 'kmol/h';
    … …  
    285281streamcold.Outlet.z             = [0.1,0.7, 0.2];
    286282
    287 exchanger.Details.A     = 15*'m^2';
     283exchanger.A     = 15*'m^2';
    288284       
    289285OPTIONS
Note: See TracChangeset for help on using the changeset viewer.