Jump to content

Methanol to Hydrocarbons problem: Difference between revisions

From mintOC
No edit summary
No edit summary
 
(4 intermediate revisions by the same user not shown)
Line 2: Line 2:
|nz        = 3
|nz        = 3
|np        = 5
|np        = 5
|nc        = 5
}}<!-- Do not insert line break here or Dimensions Box moves up in the layout...
}}<!-- Do not insert line break here or Dimensions Box moves up in the layout...


Line 17: Line 18:
  \displaystyle \min_{\theta} &\sum\limits_{j=1}^{16} &&||y(\tau_j; \theta) - z_j||^2  \\[1.5ex]
  \displaystyle \min_{\theta} &\sum\limits_{j=1}^{16} &&||y(\tau_j; \theta) - z_j||^2  \\[1.5ex]
  \mbox{s.t.}  
  \mbox{s.t.}  
  & \dot{y}_1 & = &  -( 2 \theta_2 - \frac{\theta_1 y_2}{(\theta_2 + \theta_5) y_1 + y_2} + \theta_3 + \theta_4) y_1, \\
  & \dot{y}_1 & = &  -( 2 \theta_2 - \frac{\theta_1 y_2}{(\theta_2 + \theta_5) y_1 + y_2} + \theta_3 + \theta_4) y_1, \\[0.3cm]
  & \dot{y}_2 & = & \frac{\theta_1 y_1 (\theta_2 y_1 - y_2)}{(\theta_2 + \theta_5) y_1 + y_2} + \theta_3 y_1,  \\
  & \dot{y}_2 & = & \frac{\theta_1 y_1 (\theta_2 y_1 - y_2)}{(\theta_2 + \theta_5) y_1 + y_2} + \theta_3 y_1,  \\[0.3cm]
  & \dot{y}_2 & = & \frac{\theta_1 y_1 (y_2 + \theta_5 y_1}{(\theta_2 + \theta_5) y_1 + y_2} + \theta_4 y_1,  \\
  & \dot{y}_3 & = & \frac{\theta_1 y_1 (y_2 + \theta_5 y_1}{(\theta_2 + \theta_5) y_1 + y_2} + \theta_4 y_1,  \\[0.3cm]
  & \theta_i & \geq & 0
  & \theta_i & \geq & 0 \quad i = 1,...,5.
\end{array}  
\end{array}  
</math>
</math>
Line 38: Line 39:
[[Category:MIOCP]]
[[Category:MIOCP]]
[[Category:ODE model]]
[[Category:ODE model]]
[[Category:Chemical engineering]]

Latest revision as of 08:32, 27 July 2016

Methanol to Hydrocarbons problem
Algebraic states: 3
Continuous control values: 5
Path constraints: 5

The Methanol to Hydrocarbons problem tries to determine "reaction coefficients for the conversion of methanol into various hydrocarbons." (Cite and problem taken from the COPS library)


Mathematical formulation

The problem is given by

minθj=116||y(τj;θ)zj||2s.t.y˙1=(2θ2θ1y2(θ2+θ5)y1+y2+θ3+θ4)y1,y˙2=θ1y1(θ2y1y2)(θ2+θ5)y1+y2+θ3y1,y˙3=θ1y1(y2+θ5y1(θ2+θ5)y1+y2+θ4y1,θi0i=1,...,5.

Parameters

The values zj are measurements for the concentration for y at time points τ1,...,τ16 and initial conditions are known.

Source Code

Model descriptions are available in