Batch distillation problem (TACO) - Revision history

JonasSchulze: Text replacement - "\(.*)\<\/bibref\>" to ""

2015-12-30T20:32:10Z

Text replacement - "\<bibref\>(.*)\<\/bibref\>" to "<bib id="$1" />"

← Older revision		Revision as of 20:32, 30 December 2015
Line 1:		Line 1:
	This page contains a model of the [[Batch distillation problem]] in [http://www.ampl.org AMPL] format, making use of the TACO toolkit for AMPL control optimization extensions. The original model can be found in <~~bibref>~~Diehl2006c</~~bibref~~>.		This page contains a model of the [[Batch distillation problem]] in [http://www.ampl.org AMPL] format, making use of the TACO toolkit for AMPL control optimization extensions. The original model can be found in <bib id="Diehl2006c" />.
	Note that you will need to include a generic [[support AMPL files\|AMPL/TACO support file]], OptimalControl.mod.		Note that you will need to include a generic [[support AMPL files\|AMPL/TACO support file]], OptimalControl.mod.
	To solve this model, you require an optimal control or NLP code that uses the TACO toolkit to support the AMPL optimal control extensions.		To solve this model, you require an optimal control or NLP code that uses the TACO toolkit to support the AMPL optimal control extensions.

JonasSchulze: Text replacement - "" to ""

2015-12-30T20:28:26Z

Text replacement - "<bibreferences/>" to "<biblist />"

← Older revision		Revision as of 20:28, 30 December 2015
Line 124:		Line 124:
	== References ==		== References ==

	<~~bibreferences~~/>		<biblist />

	[[Category:AMPL/TACO]]		[[Category:AMPL/TACO]]

Ckirches at 20:29, 3 October 2011

2011-10-03T20:29:03Z

← Older revision		Revision as of 20:29, 3 October 2011
Line 1:		Line 1:
	This page contains a model of the [[Batch distillation problem]] in [http://www.ampl.org AMPL] format, making use of the TACO toolkit for AMPL control optimization extensions. The original model can be found in <bibref>~~Diehl2006cf~~</bibref>.		This page contains a model of the [[Batch distillation problem]] in [http://www.ampl.org AMPL] format, making use of the TACO toolkit for AMPL control optimization extensions. The original model can be found in <bibref>Diehl2006c</bibref>.
	Note that you will need to include a generic [[support AMPL files\|AMPL/TACO support file]], OptimalControl.mod.		Note that you will need to include a generic [[support AMPL files\|AMPL/TACO support file]], OptimalControl.mod.
	To solve this model, you require an optimal control or NLP code that uses the TACO toolkit to support the AMPL optimal control extensions.		To solve this model, you require an optimal control or NLP code that uses the TACO toolkit to support the AMPL optimal control extensions.

Ckirches: fixed a reference

2011-10-03T20:27:51Z

fixed a reference

← Older revision		Revision as of 20:27, 3 October 2011
Line 1:		Line 1:
	This page contains a model of the [[Batch distillation problem]] in [http://www.ampl.org AMPL] format, making use of the TACO toolkit for AMPL control optimization extensions. The original model can be found in <bibref>~~Diehl2006b~~</bibref>.		This page contains a model of the [[Batch distillation problem]] in [http://www.ampl.org AMPL] format, making use of the TACO toolkit for AMPL control optimization extensions. The original model can be found in <bibref>Diehl2006cf</bibref>.
	Note that you will need to include a generic [[support AMPL files\|AMPL/TACO support file]], OptimalControl.mod.		Note that you will need to include a generic [[support AMPL files\|AMPL/TACO support file]], OptimalControl.mod.
	To solve this model, you require an optimal control or NLP code that uses the TACO toolkit to support the AMPL optimal control extensions.		To solve this model, you require an optimal control or NLP code that uses the TACO toolkit to support the AMPL optimal control extensions.

Ckirches: Batch distillation problem (TACO)

2011-10-03T20:21:50Z

Batch distillation problem (TACO)

New page

This page contains a model of the [[Batch distillation problem]] in [http://www.ampl.org AMPL] format, making use of the TACO toolkit for AMPL control optimization extensions. The original model can be found in <bibref>Diehl2006b</bibref>.
Note that you will need to include a generic [[support AMPL files|AMPL/TACO support file]], OptimalControl.mod.
To solve this model, you require an optimal control or NLP code that uses the TACO toolkit to support the AMPL optimal control extensions.

=== AMPL ===

This is the source file batchdist_taco.mod
<source lang="AMPL">
# ----------------------------------------------------------------
# Batch distillation problem using AMPL and TACO
# (c) Christian Kirches, Sven Leyffer
#
# Source: M.Diehl/H.G.Bock/E.Kostina'06
# ----------------------------------------------------------------
include OptimalControl.mod;

# time and free end-time

var t;
var tf := 2.5, >= 0.5, <= 10.0;

# constant parameters

param Pur := 0.99; # percent
param V := 100.0; # mol/h
param m := 0.1; # mol
param mC := 0.1; # mol

# control

var R := 8.0, >= 0.0, <= 15.0;
let R.type := "u1";
let R.scale := 0.1;

# differential states

param NDIS := 5; # PDE discretization points

var M0;
var x{0..NDIS+1};
var MD;
var xD;
var alpha;

# algebraic expressions eliminated by AMPL's presolve

var L = R/(1+R)*V;
var y{i in 0..NDIS} = (1+alpha)*x[i]/(alpha+x[i]);

var dot0 = ( L*x[1] - V*y[0] + (V-L)*x[0] ) / M0;
var dot{i in 1..NDIS} = ( L*x[i+1] - V*y[i] + V*y[i-1] - L*x[i] )/m;
var dotNDISp1 = V/mC * (-x[NDIS+1] + y[NDIS]);

# objective function

minimize Compromise:
eval (t - MD, tf);

# terminal constraint

subject to Purity_Constraint:
eval(xD, tf) >= Pur;

# ODE system

subject to ODE_M0:
diff(M0, t) = -V+L;

subject to ODE_x_0:
diff(x[0], t) = dot0;

subject to ODE_x{i in 1..NDIS}:
diff(x[i], t) = dot[i];

subject to ODE_x_NDISp1:
diff(x[NDIS+1], t) = dotNDISp1;

subject to ODE_MD:
diff(MD, t) = V-L;

subject to ODE_xD:
diff(xD, t) = (V-L) * (x[NDIS] - xD)/MD;

subject to ODE_alpha:
diff(alpha, t) = 0.0;

# Initial value constraints

subject to IVC_M0:
eval(M0, 0) = 100.0;

subject to IVC_x_0:
eval(x[0], 0) = 0.5;

subject to IVC_x{i in 1..NDIS+1}:
eval(x[i], 0) = 1;

subject to IVC_MD:
eval(MD, 0) = 0.1;

subject to IVC_xD:
eval(xD, 0) = 1;

subject to IVC_alpha:
eval(alpha, 0) = 0.2;

option solver ...;

solve;

</source>

== Other Descriptions ==

Other descriptions of this problem are available in

* Mathematical notation at [[Batch distillation problem]]

== References ==

<bibreferences/>

[[Category:AMPL/TACO]]