Difference between revisions of "Catalyst mixing problem (TACO)"

# ----------------------------------------------------------------
# Catalyst mixing problem using AMPL and TACO
# (c) Christian Kirches, Sven Leyffer
#
# Source: COPS 3.1 collocation formulation - March 2004
# ----------------------------------------------------------------
include OptimalControl.mod;
 
param ne := 2;    	  	# number of differential equations
 
var tf := 1;      		# Final time
var t;
 
param bc {1..ne};    		# Boundary conditions for x
 
var u;
let u.type := "u1";
 
var v {1..ne};
 
minimize objective: eval (-1 + v[1] + v[2], tf);
let objective.scale := 0.01;
 
subject to u_bounds: 0.0 <= u <= 1.0;
 
subject to de1:
  diff(v[1],t) = u*(10*v[2] - v[1]);
 
subject to de2:
  diff(v[2],t) = u*(v[1] - 10*v[2]) - (1 - u)*v[2];
 
subject to b_eqn {s in 1..ne}: eval(v[s],0) = bc[s];
 
data catmix_taco.dat;
 
option solver ...;
 
solve;

# Set the design parameters
 
param bc :=
  1   1
  2   0;