Difference between revisions of "Direct Current Transmission Heating Problem"

From mintOC
Jump to: navigation, search
 
(16 intermediate revisions by 2 users not shown)
Line 1: Line 1:
 
{{Dimensions
 
{{Dimensions
|nd        = 1
+
|nd        = <math>|A|</math>
|nx        = 3
+
|nx        = <math>2 |A| + |P|</math>
|nw        = 1
+
|nw        = 0
 
|nre      = 3
 
|nre      = 3
 
}}
 
}}
  
The '''direct current transfer heating problem''' models a simplified flow of direct electrical current within an electrical network with defined producers and consumers of electrical power. The model includes heating and cooling that occurs in electrical conductors due to specific resistance and heat exchange with the surrounding air. An optimal strategy for regulating the power production of individual producers as well as voltages and currents on specific conductors with the goal of minimizing the net power loss due to heating on all power lines over a fixed time horizon while fully supplying all consumers.
+
The '''direct current transmission heating problem''' models a simplified flow of direct electrical current within an electrical network with defined producers and consumers of electrical power. The model includes heating and cooling that occurs in electrical conductors due to specific resistance and heat exchange with the surrounding air. The goal is finding an optimal strategy for regulating the power production of individual producers as well as voltages and currents on specific conductors while minimizing the net power loss due to heating on all power lines over a fixed time horizon while fully supplying all consumers.
  
 
The model is largely assembled from well-known basic descriptions of phenomena easily accessible to the audience of high-school level physics courses. Mainly, these are [https://en.wikipedia.org/wiki/Ohm's_law Ohm's law], [https://en.wikipedia.org/wiki/Electrical_resistivity_and_conductivity electrical resistivity], [https://en.wikipedia.org/wiki/Joule_heating Joule heating] and [https://en.wikipedia.org/wiki/Thermal_conductivity thermal conductivity]. Note that changes of resistivity with respect to conductor temperature are taken into account.
 
The model is largely assembled from well-known basic descriptions of phenomena easily accessible to the audience of high-school level physics courses. Mainly, these are [https://en.wikipedia.org/wiki/Ohm's_law Ohm's law], [https://en.wikipedia.org/wiki/Electrical_resistivity_and_conductivity electrical resistivity], [https://en.wikipedia.org/wiki/Joule_heating Joule heating] and [https://en.wikipedia.org/wiki/Thermal_conductivity thermal conductivity]. Note that changes of resistivity with respect to conductor temperature are taken into account.
Line 37: Line 37:
 
</math>
 
</math>
 
</p>
 
</p>
where <math>C</math> is the circumference of the cross-section of the power line. Assuming constant density, we can use the conducting materials volume-specific heat capacity which we will refer to as <math>c</math>, we can calculate the change in temperature from the change in heat:
+
where <math>C</math> is the circumference of the cross-section of the power line. Assuming constant density, we can use the conducting material's volume-specific heat capacity, which we will refer to as <math>c</math>, to calculate the change in temperature from the change in heat:
 
<p>
 
<p>
 
<math>
 
<math>
Line 44: Line 44:
 
</p>
 
</p>
  
For the examples in this article, we will assume the conductor to be solid copper surrounded by a layer of [https://en.wikipedia.org/wiki/Cross-linked_polyethylene cross-linked polyethylene] with a thermal conductivity of approximately <math>0.51 W / (m \cdot K)</math>. We will assume the insulator to be <math>1 cm</math> thick. Thus, <math>q = 51 W / (m^2 K)</math>.
+
For the examples in this article, we will assume the conductor to be solid copper surrounded by a layer of [https://en.wikipedia.org/wiki/Cross-linked_polyethylene cross-linked polyethylene (PEX)] with a thermal conductivity of approximately <math>0.51 \; W / (m \cdot K)</math>.
  
 
== Mathematical formulation ==
 
== Mathematical formulation ==
Line 54: Line 54:
 
<p>
 
<p>
 
<math>
 
<math>
\begin{array}{llcl}
+
\begin{array}{llcll}
  \displaystyle \min_{x, w} & \int_{t_0}^{t_f} \sum_{a \in A} P_1^a(t) \,\mathrm{d} t  \\[1.5ex]
+
  \displaystyle \min_{T, I, u} & \int_{t_0}^{t_f} \sum_{a \in A} P_{1,a}(t) \,\mathrm{d} t  \\[1.5ex]
  \mbox{s.t.} & \dot{T}^a(t) & = & x_0(t) - x_0(t) x_1(t) - \; c_0 x_0(t) \; w(t), \\
+
  \mbox{s.t.} & \dot{T}_a & = & \beta_{1,a} \cdot e^{\alpha_a T_a} \cdot I_a^2 - \beta_{2,a} \cdot T_a \qquad & \forall a \in A, \\
  & \dot{x}_1(t) & = & - x_1(t) + x_0(t) x_1(t) - \; c_1 x_1(t) \; w(t),  \\
+
  & d_v & = & \sum_{a \in \delta^{in}(v)} (U_a \cdot I_{a} - P_1^a) - \sum_{a \in \delta^{out}(v)} (U_a \cdot I_a) \qquad & \forall v \in V \setminus P,  \\
  & \dot{x}_2(t) & = & (x_0(t) - 1)^2 + (x_1(t) - 1)^2,  \\[1.5ex]
+
  & -u_v & = & \sum_{a \in \delta^{in}(v)} (U_a \cdot I_{a} - P_1^a) - \sum_{a \in \delta^{out}(v)} (U_a \cdot I_a) \qquad & \forall v \in P,  \\[1.5ex]
  & x(0) &=& (0.5, 0.7, 0)^T, \\
+
  & T_a(t_0) &=& 0.0 \qquad & \forall a \in A, \\
  & w(t) &\in&  \{0, 1\}.
+
& I_a(t) &\in& [0, I_{max, a}] \qquad & \forall a \in A, t \in [t_0, t_f], \\
 +
  & T_a(t) &\in& [0, T_{max, a}] \qquad & \forall a \in A, t \in [t_0, t_f], \\
 +
  & u_v(t) &\in& [0, u_{max, v}] \qquad & \forall v \in P, t \in [t_0, t_f].
 
\end{array}  
 
\end{array}  
 
</math>
 
</math>
 
</p>
 
</p>
  
Here the differential states <math>(x_0, x_1)</math> describe the biomasses of prey and predator, respectively. The third differential state is used here to transform the objective, an integrated deviation, into the Mayer formulation <math>\min \; x_2(t_f)</math>. The decision, whether the fishing fleet is actually fishing at time <math>t</math> is denoted by <math>w(t)</math>.
+
The differential state <math>T_a</math> describes the temperature of power line <math>a</math> while the controls <math>I_a</math> and <math>u_v</math> represent the current in a power line <math>a \in A</math> and the power produced in a production node <math>v \in P</math> respectively.
  
 
== Parameters ==
 
== Parameters ==
  
These fixed values are used within the model.
+
There are multiple test scenarios. For all scenarios, we allow three types of conductors, all of which are copper (<math>\alpha = 0.003862</math>, <math>\rho_0 = 1.68 \cdot 10^{-8} \; \Omega m</math>, <math>c = 3.446 \cdot 10^9 \; J / (K \cdot m^3)</math>) cylinders surrounded by an insulating layer of PEX:
 
+
{| class="wikitable"
<math>
+
|-
\begin{array}{rcl}
+
! Type
[t_0, t_f] &=& [0, 12],\\
+
! Voltage
(c_0, c_1) &=& (0.4, 0.2).
+
! Conductor diameter
\end{array}
+
! Insulator thickness
</math>
+
! <math>A</math>
 +
! <math>C</math>
 +
! <math>q</math>
 +
|-
 +
| A
 +
| <math>330 \;kV</math>
 +
| <math>7 \;cm</math>
 +
| <math>2 \;cm</math>
 +
| <math>0.003848 \;m^2</math>
 +
| <math>0.2199 \;m</math>
 +
| <math>25.5 \;W / (m^2 \cdot K)</math>
 +
|-
 +
| B
 +
| <math>220 \;kV</math>
 +
| <math>5 \;cm</math>
 +
| <math>1.5 \;cm</math>
 +
| <math>0.001963 \;m^2</math>
 +
| <math>0.1571 \;m</math>
 +
| <math>34 \;W / (m^2 \cdot K)</math>
 +
|-
 +
| C
 +
| <math>110 \;kV</math>
 +
| <math>3 \;cm</math>
 +
| <math>0.75 \;cm</math>
 +
| <math>0.0007069 \;m^2</math>
 +
| <math>0.09425 \;m</math>
 +
| <math>68 \;W / (m^2 \cdot K)</math>
 +
|}
 +
PEX looses structural integrity at operating temperatures roughly in excess of <math>120^{\circ} C</math> which means that <math>T_{max,a} = 100</math> for all three types of power lines.
  
 +
=== Delayed heating ===
 +
In this scenario, we consider a simple network consisting of two producers producing up to <math>1 \;MW</math> each and a single consumer with a demand of <math>1.25 \;MW</math>. Both producers have direct power links of type C and length <math>1 \;km</math> to the consumer. The time horizon is defined as <math>(t_0, t_f) = (0, 300)</math>. A JSON description of the network can be found under [[Direct Current Transmission Heating Problem: Delayed heating (Network)]].
 +
<!--
 
== Reference Solutions ==
 
== Reference Solutions ==
  
If the problem is relaxed, i.e., we demand that <math>w(t)</math> be in the continuous interval <math>[0, 1]</math> instead of the binary choice <math>\{0,1\}</math>, the optimal solution can be determined by means of [http://en.wikipedia.org/wiki/Pontryagin%27s_minimum_principle Pontryagins maximum principle]. The optimal solution contains a singular arc, as can be seen in the plot of the optimal control. The two differential states and corresponding adjoint variables in the indirect approach are also displayed. A different approach to solving the relaxed problem is by using a direct method such as collocation or Bock's direct multiple shooting method. Optimal solutions for different control discretizations are also plotted in the leftmost figure.
+
TODO
 
+
The optimal objective value of this relaxed problem is <math>x_2(t_f) = 1.34408</math>. As follows from MIOC theory<bibref>Sager2008</bibref> this is the best lower bound on the optimal value of the original problem with the integer restriction on the control function. In other words, this objective value can be approximated arbitrarily close, if the control only switches often enough between 0 and 1. As no optimal solution exists, two suboptimal ones are shown, one with only two switches and an objective function value of <math>x_2(t_f) = 1.38276</math>, and one with 56 switches and <math>x_2(t_f) = 1.34416</math>.
+
 
+
<gallery caption="Reference solution plots" widths="180px" heights="140px" perrow="2">
+
Image:lotkaRelaxedControls.png| Optimal relaxed control determined by an indirect approach and by a direct approach on different control discretization grids.
+
Image:lotkaindirektStates.png| Differential states and corresponding adjoint variables in the indirect approach.
+
Image:lotka2Switches.png| Control and differential states with only two switches.
+
Image:lotka56Switches.png| Control and differential states with 56 switches.
+
</gallery>
+
  
  
Line 96: Line 120:
 
Model descriptions are available in
 
Model descriptions are available in
  
* [[:Category:AMPL | AMPL code]] at [[Lotka Volterra fishing problem (AMPL)]]
+
* [[:Category: Casadi | Casadi code]] at [[Direct Current Transmission Heating Problem (Casadi)]]
* [[:Category:C | C code]] at [[Lotka Volterra fishing problem (C)]]
+
-->
* [[:Category:optimica | optimica code]] at [[Lotka Volterra fishing problem (optimica)]]
+
[[Category: MIOCP]]
* [[:Category: JuMP | JuMP code]] at [[Lotka Volterra fishing problem (JuMP)]]
+
* [[:Category: Muscod | Muscod code]] at [[Lotka Volterra fishing problem (Muscod)]]
+
* [[:Category: ACADO | ACADO code]] at [[Lotka Volterra fishing problem (ACADO)]]
+
* [[:Category: VPLAN | VPLAN code]] at [[Lotka Volterra fishing problem (VPLAN)]]
+
* [[:Category: Casadi | Casadi code]] at [[Lotka Volterra fishing problem (Casadi)]]
+
* [[:Category: GloOptCon | GloOptCon code]] at [[Lotka Volterra fishing problem (GloOptCon)]]
+
* [[:Category: switch | switch code]] at [[Lotka Volterra fishing problem (switch)]]
+
== Variants ==
+
 
+
There are several alternative formulations and variants of the above problem, in particular
+
 
+
* a prescribed time grid for the control function <bibref>Sager2006</bibref>, see also [[Lotka Volterra fishing problem (AMPL)]],
+
* a time-optimal formulation to get into a steady-state <bibref>Sager2005</bibref>,
+
* the usage of a different target steady-state, as the one corresponding to <math> w(t) = 1</math> which is <math>(1 + c_1, 1 - c_0)</math>,
+
* different fishing control functions for the two species,
+
* different parameters and start values.
+
 
+
== Miscellaneous and Further Reading ==
+
The Lotka Volterra fishing problem was introduced by Sebastian Sager in a proceedings paper <bibref>Sager2006</bibref> and revisited in his PhD thesis <bibref>Sager2005</bibref>. These are also the references to look for more details.
+
 
+
== References ==
+
<bibreferences/>
+
 
+
<!--List of all categories this page is part of. List characterization of solution behavior, model properties, ore presence of implementation details (e.g., AMPL for AMPL model) here -->
+
[[Category:MIOCP]]
+
 
[[Category:ODE model]]
 
[[Category:ODE model]]
[[Category:Tracking objective]]
 
[[Category:Chattering]]
 
[[Category:Sensitivity-seeking arcs]]
 
[[Category:Population dynamics]]
 
 
 
<!--Testing Graphviz
 
<graphviz border='frame' format='svg'>
 
digraph G {Hello->World!}
 
</graphviz>
 
-->
 

Latest revision as of 23:09, 27 June 2016

Direct Current Transmission Heating Problem
State dimension: |A|
Differential states: 2 |A| + |P|
Discrete control functions: 0
Interior point equalities: 3


The direct current transmission heating problem models a simplified flow of direct electrical current within an electrical network with defined producers and consumers of electrical power. The model includes heating and cooling that occurs in electrical conductors due to specific resistance and heat exchange with the surrounding air. The goal is finding an optimal strategy for regulating the power production of individual producers as well as voltages and currents on specific conductors while minimizing the net power loss due to heating on all power lines over a fixed time horizon while fully supplying all consumers.

The model is largely assembled from well-known basic descriptions of phenomena easily accessible to the audience of high-school level physics courses. Mainly, these are Ohm's law, electrical resistivity, Joule heating and thermal conductivity. Note that changes of resistivity with respect to conductor temperature are taken into account.

Physical background

The aim of the controller in this model is to provide all consumers with their required amount of electric power. Electric power is given by the product of voltage and current P = I \cdot U. The electric power is brought to the consumers from producers via a network of power lines. Producers are assumed to have the capability of accurately regulating their power output within a certain range. Power lines are assumed to be solid blocks of an electrically conductive material surrounded by a layer of insulating material. Their cross-section is assumed to be constant along the entire length of the power line. Since both heating and power loss increase with current but not with voltage, it is always preferable to transfer direct current at the highest possible voltage which is allowed by the design of the power line. Therefore, only the electrical current is considered as a control.

The power loss in a power line is given by W = I^2 \cdot R where R is the resistance of the power line which is derived from the resistivity of the conductor using the formula


R = \rho \cdot \frac{l}{A}

where \rho is the resistivity of the conductor, l is the length of the conductor and A is the the cross-sectional area of the conducting material within the power line. Resistivity is assumed to vary exponentially with the temperature of the conducting material, meaning that \rho = \rho_0 \cdot e^{\alpha T} where \rho_0 is the resistivity of the conducting material at room temperature (which is defined as T = 0) and \alpha is a material-specific dimensionless constant. This causes a steady power loss of


P_1 = \rho_0 \cdot e^{\alpha T} \cdot \frac{l}{A} \cdot I^2.

The lost energy is transformed into heat which raises the temperature of the conductor. As the conductor is now hotter than the air surrounding its layer of insulating material, heat is steadily lost at a rate determined by the temperature difference, the conductors surface area and the insulating material's heat conductivity. The heat loss is given by


P_2 = q \cdot S \cdot \Delta T

where q is the ratio between the isolating material's thermal conductivity and its thickness, S is the surface area of the power line and \Delta T is the difference between the conductor's temperature and the surrounding air's temperature (which we will assume to be the same as room temperature, i.e. T = 0). In summary, the rate at which the heat stored in a power line changes is given by


\dot{W} = P_1 - P_2 = \rho_0 \cdot e^{\alpha T} \cdot \frac{l}{A} \cdot I^2 - q \cdot l \cdot C \cdot T

where C is the circumference of the cross-section of the power line. Assuming constant density, we can use the conducting material's volume-specific heat capacity, which we will refer to as c, to calculate the change in temperature from the change in heat:


\dot{T} = \frac{\dot{W}}{c \cdot A \cdot l} = \frac{\rho_0}{c \cdot A^2} \cdot e^{\alpha T} \cdot I^2 - \frac{q \cdot C}{c \cdot A} \cdot T =: \beta_1 \cdot e^{\alpha T} \cdot I^2 - \beta_2 \cdot T.

For the examples in this article, we will assume the conductor to be solid copper surrounded by a layer of cross-linked polyethylene (PEX) with a thermal conductivity of approximately 0.51 \; W / (m \cdot K).

Mathematical formulation

Let G = (V, A) a directed graph with vertex set V and arc set A. Let P, C \subset V be disjoint subsets of the vertex set. P serves as the set of producer nodes while C serves as the set of consumers. For c \in C let d_c denote the power demand for c. For v \in V \setminus C, let d_v = 0.

The resulting optimal control problem is given by


\begin{array}{llcll}
 \displaystyle \min_{T, I, u} & \int_{t_0}^{t_f} \sum_{a \in A} P_{1,a}(t) \,\mathrm{d} t   \\[1.5ex]
 \mbox{s.t.} & \dot{T}_a & = & \beta_{1,a} \cdot e^{\alpha_a T_a} \cdot I_a^2 - \beta_{2,a} \cdot T_a \qquad & \forall a \in A, \\
 & d_v & = & \sum_{a \in \delta^{in}(v)} (U_a \cdot I_{a} - P_1^a) - \sum_{a \in \delta^{out}(v)} (U_a \cdot I_a) \qquad & \forall v \in V \setminus P,  \\
 & -u_v & = & \sum_{a \in \delta^{in}(v)} (U_a \cdot I_{a} - P_1^a) - \sum_{a \in \delta^{out}(v)} (U_a \cdot I_a) \qquad & \forall v \in P,  \\[1.5ex]
 & T_a(t_0) &=& 0.0 \qquad & \forall a \in A, \\
 & I_a(t) &\in& [0, I_{max, a}] \qquad & \forall a \in A, t \in [t_0, t_f], \\
 & T_a(t) &\in& [0, T_{max, a}] \qquad & \forall a \in A, t \in [t_0, t_f], \\
 & u_v(t) &\in& [0, u_{max, v}] \qquad & \forall v \in P, t \in [t_0, t_f].
\end{array}

The differential state T_a describes the temperature of power line a while the controls I_a and u_v represent the current in a power line a \in A and the power produced in a production node v \in P respectively.

Parameters

There are multiple test scenarios. For all scenarios, we allow three types of conductors, all of which are copper (\alpha = 0.003862, \rho_0 = 1.68 \cdot 10^{-8} \; \Omega m, c = 3.446 \cdot 10^9 \; J / (K \cdot m^3)) cylinders surrounded by an insulating layer of PEX:

Type Voltage Conductor diameter Insulator thickness A C q
A 330 \;kV 7 \;cm 2 \;cm 0.003848 \;m^2 0.2199 \;m 25.5 \;W / (m^2 \cdot K)
B 220 \;kV 5 \;cm 1.5 \;cm 0.001963 \;m^2 0.1571 \;m 34 \;W / (m^2 \cdot K)
C 110 \;kV 3 \;cm 0.75 \;cm 0.0007069 \;m^2 0.09425 \;m 68 \;W / (m^2 \cdot K)

PEX looses structural integrity at operating temperatures roughly in excess of 120^{\circ} C which means that T_{max,a} = 100 for all three types of power lines.

Delayed heating

In this scenario, we consider a simple network consisting of two producers producing up to 1 \;MW each and a single consumer with a demand of 1.25 \;MW. Both producers have direct power links of type C and length 1 \;km to the consumer. The time horizon is defined as (t_0, t_f) = (0, 300). A JSON description of the network can be found under Direct Current Transmission Heating Problem: Delayed heating (Network).