Ignore:
Timestamp:
Mar 4, 2007, 1:56:54 AM (15 years ago)
Author:
Argimiro Resende Secchi
Message:

Fix some new language syntax.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/newlanguage/eml/stage_separators/column.mso

    r125 r176  
    6464        trays([top:topdown:bot]).OutletV.F = (1 + tanh(1 *
    6565            (trays([top:topdown:bot]).OutletV.P -
    66             trays([top:topdown:bot]).InletL.P)/"Pa"))/2 *
     66            trays([top:topdown:bot]).InletL.P)/'Pa'))/2 *
    6767            trays([top:topdown:bot]).Ah/trays([top:topdown:bot]).vV *
    6868            sqrt(2*(trays([top:topdown:bot]).OutletV.P -
    69                 trays([top:topdown:bot]).InletL.P + 1e-8 * "atm") /
     69                trays([top:topdown:bot]).InletL.P + 1e-8 * 'atm') /
    7070            (trays([top:topdown:bot]).alfa*trays([top:topdown:bot]).rhoV));
    7171
     
    111111        else
    112112                "No flow in reboiler"
    113                 reb.OutletV.F = 0.0 * "mol/s";
     113                reb.OutletV.F = 0.0 * 'mol/s';
    114114        end
    115115
    116116        "Pressure Drop through the tray"
    117         trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - cond.OutletL.P)/"Pa"))/2 *
     117        trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - cond.OutletL.P)/'Pa'))/2 *
    118118                trays(top).Ah/trays(top).vV * sqrt(2*(trays(top).OutletV.P -
    119                 cond.OutletL.P + 1e-8 * "atm") / (trays(top).alfa*trays(top).rhoV));
     119                cond.OutletL.P + 1e-8 * 'atm') / (trays(top).alfa*trays(top).rhoV));
    120120               
    121121        trays([top+topdown:topdown:bot]).OutletV.F = (1 + tanh(1 *
    122122            (trays([top+topdown:topdown:bot]).OutletV.P -
    123             trays([top+topdown:topdown:bot]).InletL.P)/"Pa"))/2 *
     123            trays([top+topdown:topdown:bot]).InletL.P)/'Pa'))/2 *
    124124            trays([top+topdown:topdown:bot]).Ah/trays([top+topdown:topdown:bot]).vV *
    125125            sqrt(2*(trays([top+topdown:topdown:bot]).OutletV.P -
    126                 trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * "atm") /
     126                trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * 'atm') /
    127127            (trays([top+topdown:topdown:bot]).alfa*trays([top+topdown:topdown:bot]).rhoV));
    128128
     
    185185        else
    186186                "No flow in reboiler"
    187                 reb.OutletV.F = 0.0 * "mol/s";
     187                reb.OutletV.F = 0.0 * 'mol/s';
    188188        end     
    189189       
    190190        "Pressure Drop through the tray"
    191         trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - ttop.Outlet.P)/"Pa"))/2 *
     191        trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - ttop.Outlet.P)/'Pa'))/2 *
    192192                trays(top).Ah/trays(top).vV /2* sqrt(2*(trays(top).OutletV.P -
    193                 ttop.Outlet.P + 1e-8 * "atm") / (trays(top).alfa*trays(top).rhoV));
     193                ttop.Outlet.P + 1e-8 * 'atm') / (trays(top).alfa*trays(top).rhoV));
    194194               
    195195        trays([top+topdown:topdown:bot]).OutletV.F = (1 + tanh(1 *
    196196            (trays([top+topdown:topdown:bot]).OutletV.P -
    197             trays([top+topdown:topdown:bot]).InletL.P)/"Pa"))/2 *
     197            trays([top+topdown:topdown:bot]).InletL.P)/'Pa'))/2 *
    198198            trays([top+topdown:topdown:bot]).Ah/trays([top+topdown:topdown:bot]).vV /2*
    199199            sqrt(2*(trays([top+topdown:topdown:bot]).OutletV.P -
    200                 trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * "atm") /
     200                trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * 'atm') /
    201201            (trays([top+topdown:topdown:bot]).alfa*trays([top+topdown:topdown:bot]).rhoV));
    202202       
     
    260260        else
    261261                "No flow in reboiler"
    262                 reb.OutletV.F = 0.0 * "mol/s";
     262                reb.OutletV.F = 0.0 * 'mol/s';
    263263        end
    264264
    265265        "Pressure Drop through the tray"
    266         trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - cond.OutletL.P)/"Pa"))/2 *
     266        trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - cond.OutletL.P)/'Pa'))/2 *
    267267                trays(top).Ah/trays(top).vV * sqrt(2*(trays(top).OutletV.P -
    268                 cond.OutletL.P + 1e-8 * "atm") / (trays(top).alfa*trays(top).rhoV));
     268                cond.OutletL.P + 1e-8 * 'atm') / (trays(top).alfa*trays(top).rhoV));
    269269               
    270270        trays([top+topdown:topdown:bot]).OutletV.F = (1 + tanh(1 *
    271271            (trays([top+topdown:topdown:bot]).OutletV.P -
    272             trays([top+topdown:topdown:bot]).InletL.P)/"Pa"))/2 *
     272            trays([top+topdown:topdown:bot]).InletL.P)/'Pa'))/2 *
    273273            trays([top+topdown:topdown:bot]).Ah/trays([top+topdown:topdown:bot]).vV *
    274274            sqrt(2*(trays([top+topdown:topdown:bot]).OutletV.P -
    275                 trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * "atm") /
     275                trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * 'atm') /
    276276            (trays([top+topdown:topdown:bot]).alfa*trays([top+topdown:topdown:bot]).rhoV));
    277277                       
     
    331331        else
    332332                "No flow in reboiler"
    333                 reb.OutletV.F = 0.0 * "mol/s";
     333                reb.OutletV.F = 0.0 * 'mol/s';
    334334        end
    335335
    336336        "Pressure Drop through the tray"
    337         trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - ttop.Outlet.P)/"Pa"))/2 *
     337        trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - ttop.Outlet.P)/'Pa'))/2 *
    338338                trays(top).Ah/trays(top).vV * sqrt(2*(trays(top).OutletV.P -
    339                 ttop.Outlet.P + 1e-8 * "atm") / (trays(top).alfa*trays(top).rhoV));
     339                ttop.Outlet.P + 1e-8 * 'atm') / (trays(top).alfa*trays(top).rhoV));
    340340               
    341341        trays([top+topdown:topdown:bot]).OutletV.F = (1 + tanh(1 *
    342342            (trays([top+topdown:topdown:bot]).OutletV.P -
    343             trays([top+topdown:topdown:bot]).InletL.P)/"Pa"))/2 *
     343            trays([top+topdown:topdown:bot]).InletL.P)/'Pa'))/2 *
    344344            trays([top+topdown:topdown:bot]).Ah/trays([top+topdown:topdown:bot]).vV *
    345345            sqrt(2*(trays([top+topdown:topdown:bot]).OutletV.P -
    346                 trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * "atm") /
     346                trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * 'atm') /
    347347            (trays([top+topdown:topdown:bot]).alfa*trays([top+topdown:topdown:bot]).rhoV));
    348348       
     
    391391        EQUATIONS
    392392        "Pressure Drop through the tray"
    393         trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - cond.OutletL.P)/"Pa"))/2 *
     393        trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - cond.OutletL.P)/'Pa'))/2 *
    394394                trays(top).Ah/trays(top).vV * sqrt(2*(trays(top).OutletV.P -
    395                 cond.OutletL.P + 1e-8 * "atm") / (trays(top).alfa*trays(top).rhoV));
     395                cond.OutletL.P + 1e-8 * 'atm') / (trays(top).alfa*trays(top).rhoV));
    396396               
    397397        trays([top+topdown:topdown:bot]).OutletV.F = (1 + tanh(1 *
    398398            (trays([top+topdown:topdown:bot]).OutletV.P -
    399             trays([top+topdown:topdown:bot]).InletL.P)/"Pa"))/2 *
     399            trays([top+topdown:topdown:bot]).InletL.P)/'Pa'))/2 *
    400400            trays([top+topdown:topdown:bot]).Ah/trays([top+topdown:topdown:bot]).vV *
    401401            sqrt(2*(trays([top+topdown:topdown:bot]).OutletV.P -
    402                 trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * "atm") /
     402                trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * 'atm') /
    403403            (trays([top+topdown:topdown:bot]).alfa*trays([top+topdown:topdown:bot]).rhoV));
    404404       
     
    448448        EQUATIONS
    449449        "Pressure Drop through the tray"
    450         trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - ttop.Outlet.P)/"Pa"))/2 *
     450        trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - ttop.Outlet.P)/'Pa'))/2 *
    451451                trays(top).Ah/trays(top).vV * sqrt(2*(trays(top).OutletV.P -
    452                 ttop.Outlet.P + 1e-8 * "atm") / (trays(top).alfa*trays(top).rhoV));
     452                ttop.Outlet.P + 1e-8 * 'atm') / (trays(top).alfa*trays(top).rhoV));
    453453               
    454454        trays([top+topdown:topdown:bot]).OutletV.F = (1 + tanh(1 *
    455455            (trays([top+topdown:topdown:bot]).OutletV.P -
    456             trays([top+topdown:topdown:bot]).InletL.P)/"Pa"))/2 *
     456            trays([top+topdown:topdown:bot]).InletL.P)/'Pa'))/2 *
    457457            trays([top+topdown:topdown:bot]).Ah/trays([top+topdown:topdown:bot]).vV *
    458458            sqrt(2*(trays([top+topdown:topdown:bot]).OutletV.P -
    459                 trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * "atm") /
     459                trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * 'atm') /
    460460            (trays([top+topdown:topdown:bot]).alfa*trays([top+topdown:topdown:bot]).rhoV));
    461461       
     
    502502        EQUATIONS
    503503        "Pressure Drop through the tray"
    504         trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - cond.OutletL.P)/"Pa"))/2 *
     504        trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - cond.OutletL.P)/'Pa'))/2 *
    505505                trays(top).Ah/trays(top).vV * sqrt(2*(trays(top).OutletV.P -
    506                 cond.OutletL.P + 1e-8 * "atm") / (trays(top).alfa*trays(top).rhoV));
     506                cond.OutletL.P + 1e-8 * 'atm') / (trays(top).alfa*trays(top).rhoV));
    507507               
    508508        trays([top+topdown:topdown:bot]).OutletV.F = (1 + tanh(1 *
    509509            (trays([top+topdown:topdown:bot]).OutletV.P -
    510             trays([top+topdown:topdown:bot]).InletL.P)/"Pa"))/2 *
     510            trays([top+topdown:topdown:bot]).InletL.P)/'Pa'))/2 *
    511511            trays([top+topdown:topdown:bot]).Ah/trays([top+topdown:topdown:bot]).vV *
    512512            sqrt(2*(trays([top+topdown:topdown:bot]).OutletV.P -
    513                 trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * "atm") /
     513                trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * 'atm') /
    514514            (trays([top+topdown:topdown:bot]).alfa*trays([top+topdown:topdown:bot]).rhoV));
    515515       
     
    559559        EQUATIONS
    560560        "Pressure Drop through the tray"
    561         trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - ttop.Outlet.P)/"Pa"))/2 *
     561        trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - ttop.Outlet.P)/'Pa'))/2 *
    562562                trays(top).Ah/trays(top).vV * sqrt(2*(trays(top).OutletV.P -
    563                 ttop.Outlet.P + 1e-8 * "atm") / (trays(top).alfa*trays(top).rhoV));
     563                ttop.Outlet.P + 1e-8 * 'atm') / (trays(top).alfa*trays(top).rhoV));
    564564               
    565565        trays([top+topdown:topdown:bot]).OutletV.F = (1 + tanh(1 *
    566566            (trays([top+topdown:topdown:bot]).OutletV.P -
    567             trays([top+topdown:topdown:bot]).InletL.P)/"Pa"))/2 *
     567            trays([top+topdown:topdown:bot]).InletL.P)/'Pa'))/2 *
    568568            trays([top+topdown:topdown:bot]).Ah/trays([top+topdown:topdown:bot]).vV *
    569569            sqrt(2*(trays([top+topdown:topdown:bot]).OutletV.P -
    570                 trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * "atm") /
     570                trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * 'atm') /
    571571            (trays([top+topdown:topdown:bot]).alfa*trays([top+topdown:topdown:bot]).rhoV));
    572572       
     
    613613        EQUATIONS
    614614        "Pressure Drop through the tray"
    615         trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - cond.OutletL.P)/"Pa"))/2 *
     615        trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - cond.OutletL.P)/'Pa'))/2 *
    616616                trays(top).Ah/trays(top).vV * sqrt(2*(trays(top).OutletV.P -
    617                 cond.OutletL.P + 1e-8 * "atm") / (trays(top).alfa*trays(top).rhoV));
     617                cond.OutletL.P + 1e-8 * 'atm') / (trays(top).alfa*trays(top).rhoV));
    618618               
    619619        trays([top+topdown:topdown:bot]).OutletV.F = (1 + tanh(1 *
    620620            (trays([top+topdown:topdown:bot]).OutletV.P -
    621             trays([top+topdown:topdown:bot]).InletL.P)/"Pa"))/2 *
     621            trays([top+topdown:topdown:bot]).InletL.P)/'Pa'))/2 *
    622622            trays([top+topdown:topdown:bot]).Ah/trays([top+topdown:topdown:bot]).vV *
    623623            sqrt(2*(trays([top+topdown:topdown:bot]).OutletV.P -
    624                 trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * "atm") /
     624                trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * 'atm') /
    625625            (trays([top+topdown:topdown:bot]).alfa*trays([top+topdown:topdown:bot]).rhoV));
    626626       
     
    671671        EQUATIONS
    672672        "Pressure Drop through the tray"
    673         trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - ttop.Outlet.P)/"Pa"))/2 *
     673        trays(top).OutletV.F = (1 + tanh(1 * (trays(top).OutletV.P - ttop.Outlet.P)/'Pa'))/2 *
    674674                trays(top).Ah/trays(top).vV * sqrt(2*(trays(top).OutletV.P -
    675                 ttop.Outlet.P + 1e-8 * "atm") / (trays(top).alfa*trays(top).rhoV));
     675                ttop.Outlet.P + 1e-8 * 'atm') / (trays(top).alfa*trays(top).rhoV));
    676676               
    677677        trays([top+topdown:topdown:bot]).OutletV.F = (1 + tanh(1 *
    678678            (trays([top+topdown:topdown:bot]).OutletV.P -
    679             trays([top+topdown:topdown:bot]).InletL.P)/"Pa"))/2 *
     679            trays([top+topdown:topdown:bot]).InletL.P)/'Pa'))/2 *
    680680            trays([top+topdown:topdown:bot]).Ah/trays([top+topdown:topdown:bot]).vV *
    681681            sqrt(2*(trays([top+topdown:topdown:bot]).OutletV.P -
    682                 trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * "atm") /
     682                trays([top+topdown:topdown:bot]).InletL.P + 1e-8 * 'atm') /
    683683            (trays([top+topdown:topdown:bot]).alfa*trays([top+topdown:topdown:bot]).rhoV));
    684684       
     
    728728        else
    729729                "Prato selado"
    730                 reb.OutletV.F = 0.0 * "mol/s";
     730                reb.OutletV.F = 0.0 * 'mol/s';
    731731        end
    732732
     
    734734        trays([top:topdown:bot]).OutletV.F = (1 + tanh(1 *
    735735            (trays([top:topdown:bot]).OutletV.P -
    736             trays([top:topdown:bot]).InletL.P)/"Pa"))/2 *
     736            trays([top:topdown:bot]).InletL.P)/'Pa'))/2 *
    737737            trays([top:topdown:bot]).Ah/trays([top:topdown:bot]).vV *
    738738            sqrt(2*(trays([top:topdown:bot]).OutletV.P -
    739                 trays([top:topdown:bot]).InletL.P + 1e-8 * "atm") /
     739                trays([top:topdown:bot]).InletL.P + 1e-8 * 'atm') /
    740740            (trays([top:topdown:bot]).alfa*trays([top:topdown:bot]).rhoV));
    741741       
     
    786786        else
    787787                "Prato selado"
    788                 reb.OutletV.F = 0.0 * "mol/s";
     788                reb.OutletV.F = 0.0 * 'mol/s';
    789789        end
    790790
     
    792792        trays([top:topdown:bot]).OutletV.F = (1 + tanh(1 *
    793793            (trays([top:topdown:bot]).OutletV.P -
    794             trays([top:topdown:bot]).InletL.P)/"Pa"))/2 *
     794            trays([top:topdown:bot]).InletL.P)/'Pa'))/2 *
    795795            trays([top:topdown:bot]).Ah/trays([top:topdown:bot]).vV *
    796796            sqrt(2*(trays([top:topdown:bot]).OutletV.P -
    797                 trays([top:topdown:bot]).InletL.P + 1e-8 * "atm") /
     797                trays([top:topdown:bot]).InletL.P + 1e-8 * 'atm') /
    798798            (trays([top:topdown:bot]).alfa*trays([top:topdown:bot]).rhoV));
    799799       
     
    842842        else
    843843                "Prato selado"
    844                 reb.OutletV.F = 0.0 * "mol/s";
     844                reb.OutletV.F = 0.0 * 'mol/s';
    845845        end
    846846
     
    848848        trays([top:topdown:bot]).OutletV.F = (1 + tanh(1 *
    849849            (trays([top:topdown:bot]).OutletV.P -
    850             trays([top:topdown:bot]).InletL.P)/"Pa"))/2 *
     850            trays([top:topdown:bot]).InletL.P)/'Pa'))/2 *
    851851            trays([top:topdown:bot]).Ah/trays([top:topdown:bot]).vV *
    852852            sqrt(2*(trays([top:topdown:bot]).OutletV.P -
    853                 trays([top:topdown:bot]).InletL.P + 1e-8 * "atm") /
     853                trays([top:topdown:bot]).InletL.P + 1e-8 * 'atm') /
    854854            (trays([top:topdown:bot]).alfa*trays([top:topdown:bot]).rhoV));
    855855       
     
    900900        else
    901901                "Prato selado"
    902                 reb.OutletV.F = 0.0 * "mol/s";
     902                reb.OutletV.F = 0.0 * 'mol/s';
    903903        end
    904904
     
    906906        trays([top:topdown:bot]).OutletV.F = (1 + tanh(1 *
    907907            (trays([top:topdown:bot]).OutletV.P -
    908             trays([top:topdown:bot]).InletL.P)/"Pa"))/2 *
     908            trays([top:topdown:bot]).InletL.P)/'Pa'))/2 *
    909909            trays([top:topdown:bot]).Ah/trays([top:topdown:bot]).vV *
    910910            sqrt(2*(trays([top:topdown:bot]).OutletV.P -
    911                 trays([top:topdown:bot]).InletL.P + 1e-8 * "atm") /
     911                trays([top:topdown:bot]).InletL.P + 1e-8 * 'atm') /
    912912            (trays([top:topdown:bot]).alfa*trays([top:topdown:bot]).rhoV));
    913913       
     
    942942       
    943943        EQUATIONS
    944         if ( reb.OutletV.P > 1 * "atm" ) then
     944        if ( reb.OutletV.P > 1 * 'atm' ) then
    945945                "Pressure Drop through the tray"
    946                 reb.OutletV.F = trays(1).Ah/reb.vV * sqrt((reb.OutletV.P - 1*"atm") / (0.15*reb.rhoV) );
     946                reb.OutletV.F = trays(1).Ah/reb.vV * sqrt((reb.OutletV.P - 1*'atm') / (0.15*reb.rhoV) );
    947947        else
    948948                "Prato selado"
    949                 reb.OutletV.F = 0.0 * "mol/s";
     949                reb.OutletV.F = 0.0 * 'mol/s';
    950950        end
    951951       
Note: See TracChangeset for help on using the changeset viewer.