Changeset 391 for trunk/sample/auto


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

Adding AUTO DAE examples.

Location:
trunk/sample/auto
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • trunk/sample/auto/ab_dae.dat

    r388 r391  
    55p2
    66p3
    7 1
    8 /home/aps/winapps/EMSO/mso/eml
    9 1
    10 /home/aps/winapps/EMSO/interface
     70
     80

  • trunk/sample/auto/ab_ode.dat

    r389 r391  
    55p2
    66p3
    7 1
    8 /home/aps/winapps/EMSO/mso/eml
    9 1
    10 /home/aps/winapps/EMSO/interface
     70
     80
     9

  • trunk/sample/auto/c.cstr_dae

    r388 r391  
    33100 4 3 1 1 3 0 0           NTST,NCOL,IAD,ISP,ISW,IPLT,NBC,NINT
    44100000 -1 20 0.0 100.0           NMX,RL0,RL1,A0,A1
    5 10000 10 2 8 5 3 0           NPR,MXBF,IID,ITMX,ITNW,NWTN,JAC
     5100 10 2 8 5 3 0           NPR,MXBF,IID,ITMX,ITNW,NWTN,JAC
    661e-06 1e-06 0.0001           EPSL,EPSU,EPSS
    770.01 0.005 0.05 1           DS,DSMIN,DSMAX,IADS

  • trunk/sample/auto/cstr_dae.dat

    r388 r391  
    331
    44qc
    5 1
    6 /home/aps/winapps/EMSO/mso/eml
    7 1
    8 /home/aps/winapps/EMSO/interface
     50
     60
     7

  • 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.