Direct Current Transmission Heating Problem
|Direct Current Transmission Heating Problem
|Discrete control functions:
|Interior point equalities:
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.
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 . 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 where is the resistance of the power line which is derived from the resistivity of the conductor using the formula
where is the resistivity of the conductor, is the length of the conductor and 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 where is the resistivity of the conducting material at room temperature (which is defined as ) and is a material-specific dimensionless constant. This causes a steady power loss of
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
where is the ratio between the isolating material's thermal conductivity and its thickness, is the surface area of the power line and 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. ). In summary, the rate at which the heat stored in a power line changes is given by
where 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 , to calculate the change in temperature from the change in heat:
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 .
Let a directed graph with vertex set and arc set . Let be disjoint subsets of the vertex set. serves as the set of producer nodes while serves as the set of consumers. For let denote the power demand for . For , let .
The resulting optimal control problem is given by
The differential state describes the temperature of power line while the controls and represent the current in a power line and the power produced in a production node respectively.
There are multiple test scenarios. For all scenarios, we allow three types of conductors, all of which are copper (, , ) cylinders surrounded by an insulating layer of PEX:
PEX looses structural integrity at operating temperatures roughly in excess of which means that for all three types of power lines.
In this scenario, we consider a simple network consisting of two producers producing up to each and a single consumer with a demand of . Both producers have direct power links of type C and length to the consumer. The time horizon is defined as . A JSON description of the network can be found under Direct Current Transmission Heating Problem: Delayed heating (Network).