Ignore:
Timestamp:
Oct 17, 2007, 9:23:15 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

    r390 r391  
    44
    55        PARAMETERS
     6        esc as Real;
    67        fi as Real;
    78        beta as Real;
     
    2930
    3031        SET
    31         fi              = 0.04;
     32
     33        fi              = 0.133;
    3234        beta    = 8;
    33         delta   = 0.78;                                         
     35        delta   = 1;                                           
    3436        q               = 1;
    35         alfa    = 0.19;
    36         S               = 1.015;
    37         csi             = 0.32;
     37        alfa    = 1;
     38        S               = 0.01;
     39        csi             = 1;
    3840        delta1  = 10;
    39         delta2  = 0.952;
    40         gama    = 27.85;                       
    41         x1f             = 1;                   
     41        delta2  = 1;
     42        gama    = 1000;                 
     43        x1f             = 1;   
    4244        x2f             = 0;
    4345        x3f             = 0;
    4446        x4f             = -1;
    4547        qc              = 0;
    46 
     48        esc             = 1e6;
     49       
    4750        EQUATIONS
    4851 
    49         q*(x1f-u1)-u1*u5*fi - diff(u1)*'s' = 0;
    50         q*(x2f-u2)-u2*fi*S*u6+u1*fi*u5 - diff(u2)*'s' = 0;
    51         q*(x3f-u3)+delta*(u4-u3)+beta*fi*(u1*u5+alfa*u2*u6*S) - diff(u3)*'s' = 0;
     52
     53        q*(x1f-u1)-u1*u5*esc*fi - diff(u1)*'s' = 0;
     54
     55        q*(x2f-u2)-u2*fi*S*u6*esc+u1*fi*u5*esc - diff(u2)*'s' = 0;
     56
     57        q*(x3f-u3)+delta*(u4-u3)+beta*fi*(u1*u5*esc+alfa*u2*u6*esc*S) - diff(u3)*'s' = 0;
     58
    5259        delta1*(qc*(x4f-u4)+delta*delta2*(u3-u4)) - diff(u4)*'s' = 0;
    53         -u5 + exp(u3/(1+(u3/gama))) = 0;
    54         -u6 + exp(csi*u3/(1+(u3/gama))) = 0;
     60
     61        -u5*esc + exp(u3/(1+(u3/gama))) = 0;
     62
     63        -u6*esc + exp(csi*u3/(1+(u3/gama))) = 0;
    5564
    5665        GUESS
    57         u1 = 0.6419;
    58         u2 = 0.3270;
    59         u3 = 2.9101;
    60         u4 = 2.9101;
    61         u5 = 13.9403;
    62         u6 = 2.3236;                                     
     66        u1 = 0;
     67        u2 = 0;
     68        u3 = 16;
     69        u4 = 16;
     70        u5 = 3;
     71        u6 = 3;                                 
     72
    6373
    6474    INITIAL
     
    6878        u4 = 3;
    6979
     80
    7081        OPTIONS
    7182        TimeEnd = 10;
Note: See TracChangeset for help on using the changeset viewer.