Ignore:
Timestamp:
Sep 21, 2007, 7:17:36 PM (15 years ago)
Author:
Argimiro Resende Secchi
Message:

Adding pressure drop variable (Pdrop) and pressure ratio (Pratio) in some models. Need to do the same for other models.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/eml/stage_separators/flash.mso

    r353 r372  
    6464        Across as area (Brief="Flash Cross section area");
    6565        vfrac as positive (Brief="Vapourization fraction", Symbol="\phi");
     66        Pratio as positive      (Brief = "Pressure Ratio", Symbol ="P_{ratio}");       
     67        Pdrop as press_delta (Brief = "Pressure Drop", DisplayUnit = 'kPa', Symbol ="\Delta P");
    6668
    6769        EQUATIONS
     
    102104        "Mechanical Equilibrium"
    103105        OutletV.P = OutletL.P;
    104        
     106
     107        "Pressure Drop"
     108        OutletL.P  = Inlet.P - Pdrop;
     109
     110        "Pressure Ratio"
     111        OutletL.P = Inlet.P * Pratio;
     112
    105113        "Geometry Constraint"
    106114        V = ML * vL + MV * vV;
     
    151159in      InletQ as energy_stream (Brief="Rate of heat supply", PosX=1, PosY=0.7559, Symbol="_{in}");
    152160        vfrac as fraction (Brief="Vapourization fraction", Symbol="\phi");
     161        Pratio as positive      (Brief = "Pressure Ratio", Symbol ="P_{ratio}");       
     162        Pdrop as press_delta (Brief = "Pressure Drop", DisplayUnit = 'kPa', Symbol ="\Delta P");
    153163
    154164        EQUATIONS
     
    170180        "Mechanical Equilibrium"
    171181        OutletV.P = OutletL.P;
     182
     183        "Pressure Drop"
     184        OutletL.P  = Inlet.P - Pdrop;
     185
     186        "Pressure Ratio"
     187        OutletL.P = Inlet.P * Pratio;
    172188end
    173189
     
    203219        vfrac as fraction(Brief="Vaporization fraction", Symbol="\phi");
    204220        h as enth_mol(Brief="Mixture enthalpy");
     221        Pratio as positive (Brief = "Pressure Ratio", Symbol ="P_{ratio}");     
     222        Pdrop as press_delta (Brief = "Pressure Drop", DisplayUnit = 'kPa', Symbol ="\Delta P");
    205223
    206224        EQUATIONS
     
    222240        "Mechanical Equilibrium"
    223241        OutletV.P = OutletL.P;
     242
     243        "Pressure Drop"
     244        OutletL.P  = Inlet.P - Pdrop;
     245
     246        "Pressure Ratio"
     247        OutletL.P = Inlet.P * Pratio;
    224248end
    225249
     
    264288        xsat(NComp) as Real(Lower=0, Upper=1, Brief="Liquid composition if saturated");
    265289        ysat(NComp) as Real(Lower=0, Upper=1, Brief="Vapour composition if saturated");
     290        Pratio as positive (Brief = "Pressure Ratio", Symbol ="P_{ratio}");     
     291        Pdrop as press_delta (Brief = "Pressure Drop", DisplayUnit = 'kPa', Symbol ="\Delta P");
    266292       
    267293        zero_one as fraction(Brief="Regularization Variable");
     
    295321        "Mechanical Equilibrium"
    296322        OutletV.P = OutletL.P;
    297        
     323
     324        "Pressure Drop"
     325        OutletL.P  = Inlet.P - Pdrop;
     326
     327        "Pressure Ratio"
     328        OutletL.P = Inlet.P * Pratio;
     329
    298330        # regularization functions
    299331        zero_one = (1 + tanh(B * vsat))/2;
Note: See TracChangeset for help on using the changeset viewer.