Difference between revisions of "Control of Heat Equation with Actuator Placement"
FelixMueller (Talk | contribs) (→Mathematical formulation) |
FelixMueller (Talk | contribs) (→Discretization) |
||
Line 72: | Line 72: | ||
For the heat equation, we use a five-point-stencil in space and the implicit euler in time. | For the heat equation, we use a five-point-stencil in space and the implicit euler in time. | ||
For <math>i=1,\dots, N-1</math>, <math>j=1,\dots, M-1</math>, and <math>k=0,\dots, T_n-1</math>, this yields: | For <math>i=1,\dots, N-1</math>, <math>j=1,\dots, M-1</math>, and <math>k=0,\dots, T_n-1</math>, this yields: | ||
− | + | <p | |
<math> \frac{1}{h_t}(u_{i,j,k}-u_{i,j,k-1}) - \frac{\kappa_{i,j}}{h_x h_y}(u_{i-1,j,k}+u_{i+1,j,k}+u_{i,j-1,k}+u_{i,j+1,k}-4u_{i,j,k}) = \sum\limits_{l=1}^L(v_{k+1,l}+v_{k ,l} )f_l(ih_x,jh_y), </math> | <math> \frac{1}{h_t}(u_{i,j,k}-u_{i,j,k-1}) - \frac{\kappa_{i,j}}{h_x h_y}(u_{i-1,j,k}+u_{i+1,j,k}+u_{i,j-1,k}+u_{i,j+1,k}-4u_{i,j,k}) = \sum\limits_{l=1}^L(v_{k+1,l}+v_{k ,l} )f_l(ih_x,jh_y), </math> | ||
− | + | </p> | |
with <math>u_{i,j,k}=u(ih_x,jh_y,kh_t)</math>,<math>v_{k ,l}=v_l(kh_t)</math>, <math>\kappa_{i,j}=\kappa(ih_x,jh_y)</math>, the stepsizes in space <math>h_x,h_y</math>, and in time <math>h_t</math>, respectively. | with <math>u_{i,j,k}=u(ih_x,jh_y,kh_t)</math>,<math>v_{k ,l}=v_l(kh_t)</math>, <math>\kappa_{i,j}=\kappa(ih_x,jh_y)</math>, the stepsizes in space <math>h_x,h_y</math>, and in time <math>h_t</math>, respectively. | ||
Revision as of 23:03, 23 February 2016
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 , 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 by and adding the Big formulation.
Contents
Mathematical formulation
Parameters
We define the source term for all locations and a fix parameter : where is the coordinate of the mesh point of the th possible actuator location.
The parameters used are:
The parameter describes the thermal dissipativity of the material in the domain , it may vary in space: .
The parameter indicates the number of possible actuator locations. They are distributed as indicated in the picture.
The source budget is limited by and denotes the final time.
Discretization
To solve the problem we apply a "first discretize, then optimize" approach an discretize the components of the problem. For the heat equation, we use a five-point-stencil in space and the implicit euler in time. For , , and , this yields: <p
</p> with ,, , the stepsizes in space , and in time , respectively.
It holds for the source buget with the discretized binary controls for all :
This condition is called SOS- conditon.
We remark that the number of discretized binary variables does not depend on the space discretization but it depends on the time discretization. Thats why we taged this problem that it contains mesh-independend and mesh-dependend integer variables.
Reference Solution
Source Code
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 |