Changeset 708 for trunk/sample


Ignore:
Timestamp:
Feb 2, 2009, 10:50:20 PM (14 years ago)
Author:
Argimiro Resende Secchi
Message:

Changing name of reverse integration.

Location:
trunk/sample/miscellaneous
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/sample/miscellaneous/bio.rlt

    r707 r708  
    440 0
    55
    6 9
     610
    77miMax
    881 1 21 0.53 0.53 0.53 0.53 0.53 0.53 0.53 0.53 0.53 0.53 0.53 0.53 0.53 0.53 0.53 0.53 0.53 0.53 0.53 0.53 0.53 0
    9 
    10 0
    11 D
    12 1 1 21 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0
    139
    14100
     
    30260
    3127biomass
    32 1 1 21 1 1.03103 1.06258 1.09461 1.12729 1.16095 1.19594 1.23249 1.27064 1.31005 1.34963 1.38729 1.42021 1.44596 1.46435 1.47707 1.48614 1.49303 1.49854 1.5031 1.50695 0
     281 1 21 1 1.02973 1.05276 1.06421 1.06479 1.06177 1.06715 1.09438 1.15516 1.256 1.38697 1.47437 1.48025 1.46135 1.41565 1.33337 1.22228 1.11148 1.03395 1.01347 1.06688 0
    3329
    34300
    3531substrate
    36 1 1 21 0.5 0.561773 0.602834 0.625941 0.633031 0.625316 0.603636 0.568876 0.522242 0.465661 0.402831 0.339745 0.284194 0.242835 0.216669 0.201929 0.193944 0.189352 0.186445 0.184402 0.18284 0
     321 1 21 0.5 0.566 0.6339 0.71861 0.817116 0.908706 0.961833 0.943574 0.8273 0.600602 0.293024 0.0936341 0.101276 0.175659 0.320251 0.55582 0.859213 1.15584 1.36332 1.42252 1.29324 0
    3733
    38340
    3935mi
    40 1 1 21 0.36122 0.360806 0.359798 0.359024 0.358758 0.35906 0.359768 0.360659 0.361251 0.360705 0.357875 0.351543 0.341629 0.330315 0.320759 0.31436 0.310515 0.308177 0.30664 0.305576 0.304663 0
     361 1 21 0.36122 0.360727 0.358735 0.354846 0.34909 0.34303 0.339333 0.340614 0.348445 0.359859 0.343553 0.228047 0.23757 0.300631 0.34863 0.360933 0.346411 0.325332 0.310341 0.306177 0.315312 0
    4137
    42380
     39D
     401 1 21 0.3 0.307191 0.325243 0.344888 0.354562 0.344878 0.312698 0.263168 0.209187 0.168036 0.156161 0.183586 0.249705 0.341948 0.437967 0.51105 0.537446 0.503614 0.411272 0.278581 0.13679 0
     41
     420
     43tempo
     441 1 21 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 8.5 9 9.5 10 0
     45
     460
  • trunk/sample/miscellaneous/sample_reverse.mso

    r707 r708  
    2525        PARAMETERS
    2626        miMax   as frequency(Brief="Monod parameter", DisplayUnit='1/h');
    27         D       as frequency(Brief="Dilution rate F/V", DisplayUnit='1/h');
    2827        Km      as conc_mass(Brief="Monod parameter", Default=0.12);
    2928        K1      as inv_conc_mass(Brief="Monod parameter", Default=0.4545);
     
    3534        substrate as conc_mass (Brief="Substrate concentration");
    3635        mi        as frequency (Brief="Specific growth rate", DisplayUnit='1/h');
     36        D         as frequency(Brief="Dilution rate F/V", DisplayUnit='1/h');
     37        tempo     as time_h (Brief="Independent variable");
    3738       
    3839        EQUATIONS
     
    4647        diff(substrate) = D*(x2f-substrate)-(biomass*mi/Y);
    4748       
     49        "Independent variable to use in explicit functions"
     50        diff(tempo) = 1;
     51       
     52        "Manipulated variable"
     53        D = 0.3 * (1 + tempo/'10*h'*sin(tempo*'rad/h'))/'h';
     54       
    4855        SET
    4956        miMax = 0.53  *  '1/h';
    50         D      = 0.3  *  '1/h';
    5157       
    5258        INITIAL
    5359        biomass = 1  *  'kg/m^3';
    5460        substrate = 0.5  *  'kg/m^3';
     61        tempo = 0 * 'h';
    5562 
    5663        OPTIONS
     
    5966        TimeEnd = 10;
    6067        TimeUnit = 'h';
    61         InitialFile = "bio"; # use with adjoint=1 and compare with adjoint=0
    62 #       Linearize = true;
    63         SparseAlgebra = true;
    64     DAESolver(File = "adjoint", adjoint=1);
     68#       InitialFile = "bio"; # use with reverse=1 and compare with reverse=0
     69    DAESolver(File = "reverse", reverse=0);
    6570end
Note: See TracChangeset for help on using the changeset viewer.