Control of Heat Equation with Actuator Placement

Control of Heat Equation with Actuator Placement
State dimension:	1
Differential states:	1
Continuous control functions:	9
Discrete control functions:	9
Path constraints:	3
Interior point equalities:	2

This problem is governed by the heat equation and is adapted from Iftime and Demetriou ([Iftime2009]Author: Orest V. Iftime; Michael A. Demetriou
Journal: {A}utomatica
Number: 2
Pages: 312--323
Title: {O}ptimal control of switched distributed parameter systems with spatially scheduled actuators
Volume: 45
Year: 2009
). Its goal is to choose a place to apply an actuator in a given area depending on time. The objective function is quadratic, its first term captures the desired final state $\overline{u}\equiv 0$, the second term regularize the state over time and the third term regularize the continuous controls. The constraints are a source budget, which limits the quantity of placed actuators, and the two-dimensional heat equation with some source function. Additionally, we assume Dirichlet boundary conditions and initial conditions. Originally, the problem formulation was non-convex. We overcome this issue by substitution of $v(t)w_l(t)$ by $v_l(t)$ and adding the Big $M$ formulation \eqref{Eq:BigM}.

Mathematical formulation

$\begin{array}{llc}\min {}_{u,v,w}& J(u,v)=||u(\cdot ,\cdot ,t_f)|{|}_{2,\mathrm{\Omega }}^2+2||u(\cdot ,\cdot ,\cdot )|{|}_{2,\mathrm{\Omega }\times T}^2+\frac{1}{500}\sum _{l=1}^L||v_l(\cdot )|{|}_{2,T}^2& \\ \text{ s.t.}& \frac{\partial u}{\partial t}(x,y,t)-\kappa \mathrm{\Delta }u(x,y,t)=\sum _{l=1}^9{v}_l(t)f_l(x,y)& \text{ in }& \mathrm{\Omega }\times T\\ & u(x,y,t)=0& \text{ on }& \partial \mathrm{\Omega }\times T\\ & u(x,y,0)=100\sin (\pi x)\sin (\pi y)& \text{ in }& \mathrm{\Omega }\\ & -M{w}_l(t)\le v_l(t)\le M{w}_l(t)\text{ for all }l\in \{1,\dots ,L\}& \text{ in }& T\\ & \sum _{l=1}^L{w}_l(t)=W& \text{ in }& T\\ & w_l(t)\in \{0,1\}\text{ for all }l\in \{1,\dots ,L\}& \text{ in }& T.\end{array}$

Parameters

We define the source term for all locations $l\in \{1,\dots ,L\}$ and a fix parameter $\epsilon \in {\mathbb{ℝ}}_{+}$: $f_l(x,y)=\frac{1}{\sqrt{2\pi \epsilon }}{e}^{\frac{-((x_l-x{)}^2+(y_l-y{)}^2)}{2\epsilon }}$ where $(x_l,y_l)$ is the coordinate of the mesh point of the $l$th possible actuator location.

The parameters used are:

$\begin{array}{rcl}\mathrm{\Omega }& =& [0,1]\times [0,2],\\ L& =& 9,\\ \kappa & =& 0.01,\\ t_f& =& 10,\\ W& =& 1,\\ \epsilon & =& 0.01.\end{array}$

The parameter $\kappa$ describes the thermal dissipativity of the material in the domain $\mathrm{\Omega }$, it may vary in space: $\kappa (x,y)$. The parameter $L$ indicates the number of possible actuator locations. They are distributed as indicated in the picture. The source budget is limited by $W$.

References

[Iftime2009]

Orest V. Iftime; Michael A. Demetriou (2009): {O}ptimal control of switched distributed parameter systems with spatially scheduled actuators . {A}utomatica, 45, 312--323