Ignore:
Timestamp:
Oct 17, 2007, 9:08:01 AM (15 years ago)
Author:
Argimiro Resende Secchi
Message:

Adding AUTO DAE examples.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/sample/auto/cstr_dae.mso

    r388 r390  
    11using "types";
    22
    3 FlowSheet cstr2_dae
     3FlowSheet cstr_dae
    44
    55        PARAMETERS
    … …  
    3030        SET
    3131        fi              = 0.04;
    32 
    3332        beta    = 8;
    34 
    3533        delta   = 0.78;                                         
    36 
    3734        q               = 1;
    38 
    3935        alfa    = 0.19;
    40 
    4136        S               = 1.015;
    42 
    4337        csi             = 0.32;
    44 
    4538        delta1  = 10;
    46 
    4739        delta2  = 0.952;
    48 
    4940        gama    = 27.85;                       
    50 
    5141        x1f             = 1;                   
    52 
    5342        x2f             = 0;
    54 
    5543        x3f             = 0;
    56 
    5744        x4f             = -1;
    58 
    5945        qc              = 0;
    6046
    6147        EQUATIONS
    6248 
    63 
    6449        q*(x1f-u1)-u1*u5*fi - diff(u1)*'s' = 0;
    65 
    6650        q*(x2f-u2)-u2*fi*S*u6+u1*fi*u5 - diff(u2)*'s' = 0;
    67 
    6851        q*(x3f-u3)+delta*(u4-u3)+beta*fi*(u1*u5+alfa*u2*u6*S) - diff(u3)*'s' = 0;
    69 
    7052        delta1*(qc*(x4f-u4)+delta*delta2*(u3-u4)) - diff(u4)*'s' = 0;
    71 
    7253        -u5 + exp(u3/(1+(u3/gama))) = 0;
    73 
    7454        -u6 + exp(csi*u3/(1+(u3/gama))) = 0;
    7555
    7656        GUESS
    7757        u1 = 0.6419;
    78 
    7958        u2 = 0.3270;
    80 
    8159        u3 = 2.9101;
    82 
    8360        u4 = 2.9101;
    84 
    8561        u5 = 13.9403;
    86 
    8762        u6 = 2.3236;                                     
    88 
    8963
    9064    INITIAL
    9165        u1 = 0.6;
    92 
    9366        u2 = 0.3;
    94 
    9567        u3 = 3;
    96 
    9768        u4 = 3;
    98 
    9969
    10070        OPTIONS
Note: See TracChangeset for help on using the changeset viewer.