Ignore:
Timestamp:
Jun 13, 2016, 2:21:42 AM (6 years ago)
Author:
Argimiro Resende Secchi
Message:

fixing some examples

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/sample/heat_exchangers/PHE_diagram.mso

    r749 r992  
    33*----------------------------------------------*#
    44
    5 using "streams";
    6 using "heat_exchangers/PHE";
     5using "streams.mso";
     6using "heat_exchangers\PHE.mso";
    77
    88FlowSheet PHE_diagram
    99        PARAMETERS
    1010        PP as Plugin(Brief="Physical Properties",
    11                 Type="PP",
    12                 Components = ["toluene","benzene"],
    13                 LiquidModel = "PR",
    14                 VapourModel = "PR"
     11                Type="PP"
     12                ,Components = ["toluene","benzene"]
     13                ,LiquidModel = "PR"
     14                ,VapourModel = "PR"
    1515        );
    1616        NComp as Integer;
    … …  
    2020        DEVICES
    2121        simple_sink_1 as simple_sink;
    22 
    23         SET
    24 
    25         SPECIFY
    26 
    27         INITIAL
    28 
    29         GUESS
    30         simple_sink_1.Inlet.T = 350 * 'K';
    31         simple_sink_1.Inlet.P = 2 * 'atm';
    32 
    33         DEVICES
    3422        simple_sink_2 as simple_sink;
    35 
    36         SET
    37 
    38         SPECIFY
    39 
    40         INITIAL
    41 
    42         GUESS
    43         simple_sink_2.Inlet.T = 300 * 'K';
    44         simple_sink_2.Inlet.P = 2 * 'atm';
    45 
    46         DEVICES
    4723        PHE_1 as PHE;
    48 
    49         SET
    50         PHE_1.ChevronAngle = "A30_Deg";
    51         PHE_1.SideOne = "cold";
    52         PHE_1.Geometry.Lv = 155 * 'cm';
    53         PHE_1.Geometry.Nplates = 129;
    54         PHE_1.Geometry.NpassHot = 2;
    55         PHE_1.Geometry.NpassCold = 2;
    56         PHE_1.Geometry.Dports = 15 * 'cm';
    57         PHE_1.Geometry.Lw = 70 * 'cm';
    58         PHE_1.Geometry.pitch = 3.7 * 'mm';
    59         PHE_1.Geometry.pt = 0.70 * 'mm';
    60         PHE_1.Geometry.Kwall = 17 * 'W/m/K';
    61         PHE_1.Geometry.Rfh = 0 * 'm^2*K/kW';
    62         PHE_1.Geometry.Rfc = 0 * 'm^2*K/kW';
    63         PHE_1.Geometry.PhiFactor = 1.17;
    64 
    65         SPECIFY
    66 
    67         INITIAL
    68 
    69         GUESS
    70         PHE_1.InletHot.z(1) = 1;
    71         PHE_1.InletHot.z(2) = 0;
    72         PHE_1.InletCold.z(1) = 0;
    73         PHE_1.InletCold.z(2) = 1;
    74 
    75         DEVICES
    76         source_1 as source;
    77 
    78         SET
    79         source_1.ValidPhases = "Vapour-Liquid";
    80 
    81         SPECIFY
    82         source_1.Composition(1) = 0;
    83         source_1.Composition(2) = 1;
    84         source_1.F = 290 * 'kmol/h';
    85         source_1.T = 293 * 'K';
    86         source_1.P = 2 * 'atm';
    87 
    88         INITIAL
    89 
    90         GUESS
    91 
    92         DEVICES
    93         source_2 as source;
    94 
    95         SET
    96         source_2.ValidPhases = "Vapour-Liquid";
    97 
    98         SPECIFY
    99         source_2.Composition(1) = 1;
    100         source_2.Composition(2) = 0;
    101         source_2.F = 82.5 * 'kmol/h';
    102         source_2.T = 360 * 'K';
    103         source_2.P = 2 * 'atm';
    104 
    105         INITIAL
    106 
    107         GUESS
     24        source_1 as simple_source;
     25        source_2 as simple_source;
    10826
    10927        CONNECTIONS
    … …  
    11230        source_1.Outlet to PHE_1.InletCold;
    11331        source_2.Outlet to PHE_1.InletHot;
     32#       source_1.Outlet to simple_sink_2.Inlet;
     33#       source_2.Outlet to simple_sink_1.Inlet;
     34
     35        SET
     36        PHE_1.Geometry.Lv = 155 * 'cm' ;
     37        PHE_1.Geometry.Nplates = 129 ;
     38        PHE_1.Geometry.NpassHot = 2 ;
     39        PHE_1.Geometry.NpassCold = 2 ;
     40        PHE_1.Geometry.Dports = 15 * 'cm' ;
     41        PHE_1.Geometry.Lw = 70 * 'cm' ;
     42        PHE_1.Geometry.pitch = 3.7 * 'mm' ;
     43        PHE_1.Geometry.pt = 0.70 * 'mm' ;
     44        PHE_1.Geometry.Kwall = 17 * 'W/m/K' ;
     45        PHE_1.Geometry.Rfh = 0 * 'm^2*K/kW' ;
     46        PHE_1.Geometry.Rfc = 0 * 'm^2*K/kW' ;
     47        PHE_1.Geometry.PhiFactor = 1.17 ;
     48        PHE_1.ChevronAngle = "A30_Deg" ;
     49        PHE_1.SideOne = "cold" ;
     50        source_1.ValidPhases = "Liquid-Only" ;
     51        source_2.ValidPhases = "Liquid-Only" ;
     52
     53        SPECIFY
     54#       source_1.Composition(1) = 0 ;
     55#       source_1.Composition(2) = 1 ;
     56        source_1.MolarComposition(1) = 0 ;
     57        source_1.MolarComposition(2) = 1 ;     
     58        source_1.F = 290 * 'kmol/h' ;
     59        source_1.T = 293 * 'K' ;
     60        source_1.P = 2 * 'atm' ;
     61#       source_2.Composition(1) = 1 ;
     62#       source_2.Composition(2) = 0 ;
     63        source_2.MolarComposition(1) = 1 ;
     64        source_2.MolarComposition(2) = 0 ;             
     65        source_2.F = 82.5 * 'kmol/h' ;
     66        source_2.T = 360 * 'K' ;
     67        source_2.P = 2 * 'atm' ;
     68
     69        INITIAL
     70
     71        GUESS
     72        simple_sink_1.Inlet.T = 350 * 'K' ;
     73        simple_sink_1.Inlet.P = 2 * 'atm' ;
     74        simple_sink_2.Inlet.T = 300 * 'K' ;
     75        simple_sink_2.Inlet.P = 2 * 'atm' ;
     76        PHE_1.InletHot.z(1) = 1 ;
     77        PHE_1.InletHot.z(2) = 0 ;
     78        PHE_1.InletCold.z(1) = 0 ;
     79        PHE_1.InletCold.z(2) = 1 ;
     80        PHE_1.OutletCold.T = 300 *'K';
     81        PHE_1.OutletHot.T = 350 *'K';
     82        PHE_1.OutletCold.P = 2 *'atm';
     83        PHE_1.OutletHot.P = 2 *'atm';
     84        PHE_1.OutletHot.z(1) = 1 ;
     85        PHE_1.OutletHot.z(2) = 0 ;
     86        PHE_1.OutletCold.z(1) = 0 ;
     87        PHE_1.OutletCold.z(2) = 1 ;
     88        PHE_1.ColdSide.Properties.Outlet.rho=900*'kg/m^3';
     89        PHE_1.HotSide.Properties.Outlet.rho=800*'kg/m^3';
     90        PHE_1.ColdSide.Properties.Average.rho=900*'kg/m^3';
     91        PHE_1.HotSide.Properties.Average.rho=800*'kg/m^3';
    11492
    11593        OPTIONS
    11694        Dynamic = false;
    117         Integration = "index1";
     95        Integration = "original";
     96        GuessFile = "PHE_diagram";
    11897        NLASolver(
    11998                File = "sundials",
    … …  
    123102        );
    124103        DAESolver(
    125                 File = "",
     104                File = "dasslc",
    126105                RelativeAccuracy = 1e-3,
    127106                AbsoluteAccuracy = 1e-6,
    128107                EventAccuracy = 1e-2
    129108        );
     109
    130110end
Note: See TracChangeset for help on using the changeset viewer.