<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">EPE</journal-id><journal-title-group><journal-title>Energy and Power Engineering</journal-title></journal-title-group><issn pub-type="epub">1949-243X</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/epe.2017.94B014</article-id><article-id pub-id-type="publisher-id">EPE-75238</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject></subj-group></article-categories><title-group><article-title>
 
 
  Economic Model Predictive Control for Hot Water Based Heating Systems in Smart Buildings
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>M.</surname><given-names>A. Ahmed Awadelrahman</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Yi</surname><given-names>Zong</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Hongwei</surname><given-names>Li</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Carsten</surname><given-names>Agert</given-names></name><xref ref-type="aff" rid="aff4"><sup>4</sup></xref></contrib></contrib-group><aff id="aff4"><addr-line>NEXT ENERGY, EWE Research Centre for Energy Technology at the University of Oldenburg, Oldenburg, Germany</addr-line></aff><aff id="aff2"><addr-line>Centre for Electric Power and Energy, Department of Electrical Engineering, Technical University of Denmark, Roskilde, 
Denmark</addr-line></aff><aff id="aff1"><addr-line>Institute of Physics, University of Oldenburg, Oldenburg, Germany</addr-line></aff><aff id="aff3"><addr-line>Section of Building Energy, Department of Civil Engineering, Technical University of Denmark, Lyngby, Denmark</addr-line></aff><pub-date pub-type="epub"><day>06</day><month>04</month><year>2017</year></pub-date><volume>09</volume><issue>04</issue><fpage>112</fpage><lpage>119</lpage><history><date date-type="received"><day>December</day>	<month>9,</month>	<year>2016</year></date><date date-type="rev-recd"><day>Accepted:</day>	<month>March</month>	<year>30,</year>	</date><date date-type="accepted"><day>April</day>	<month>6,</month>	<year>2017</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
   
   This paper presents a study to optimize the heating energy costs in a residential building with varying electricity price signals based on an Economic Model Predictive Controller (EMPC). The investigated heating system consists of an air source heat pump (ASHP) incorporated with a hot water tank as active Thermal Energy Storage (TES), where two optimization problems are integrated together to optimize both the ASHP electricity consumption and the building heating consumption utilizing a heat dynamic model of the building. The results show that the proposed EMPC can save the energy cost by load shifting compared with some reference cases. 
  
 
</p></abstract><kwd-group><kwd>Building Energy Management System</kwd><kwd> Demand Response</kwd><kwd> Economic Model Predictive Control</kwd><kwd> Heat Pumps</kwd><kwd> Smart Buildings</kwd><kwd> Thermal Energy Storage</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Setting a target of 100% renewable energy including transport sector by 2050 and 50% of electricity consumption supplied by wind in 2020 in Denmark, under the European 20-20-20 target regarding renewable energy policies [<xref ref-type="bibr" rid="scirp.75238-ref1">1</xref>], provides power systems with new challenges due to renewable energy sources fluctuating nature which calls for additional flexibility on the demand side. In the European Union, buildings are responsible for 40% of total energy consumption including approximately 20% absorbed in heating [<xref ref-type="bibr" rid="scirp.75238-ref2">2</xref>], which can be effectively used in the demand side management (DSM) strategy as a shiftable load.The electricity prices can reflect surplus or deficiency in electricity supply and also the low prices can indicate the high renewable energy penetration in Denmark.</p><p>The EMPC is a validated strategy in designing smart buildings’ control and has been studied thoroughly in the recent years. An MPC-based scheduling algorithms developed in [<xref ref-type="bibr" rid="scirp.75238-ref3">3</xref>] for the building energy management system to exploit the thermal mass of the building as well as the non-thermal appliances flexibility to incur energy savings. A pilot study in [<xref ref-type="bibr" rid="scirp.75238-ref4">4</xref>] showed that energy savings and load shifting can be achieved by applying EMPC for electric heaters consumption utilizing a residential building mass as heat storage.</p><p>Although the Thermal Energy Storage (TES) is a very potential powerful instrument in DSM programs [<xref ref-type="bibr" rid="scirp.75238-ref5">5</xref>], and the efficiency of the whole system can be extremely enhanced by using heat pumps as energy efficient heating equipment; but this combination has not been studied intensively in the DSM strategies. A study in [<xref ref-type="bibr" rid="scirp.75238-ref6">6</xref>] analyzed heat pump coupled with TES under a DSM strategy designed to flatten the electricity load curve by switching off the heat pump during specific peak hours.</p><p>Particularly in this paper we studied a system which consists of an ASHP incorporated with a hot water tank as active TES connected to hot-water radiators utilizing the building mass as passive energy storage. The building model used in this paper is discussed in [<xref ref-type="bibr" rid="scirp.75238-ref4">4</xref>], and the stratification phenomenon is considered in the TES model. Based on the dynamic power price signals and the weather forecast, EMPC optimizes both the ASHP electricity consumption and the building heating consumption; the system’s energy consumption is compared with some reference cases.</p><p>The rest of this paper is organized as follows: The system models are introduced in Section 2; the EMPC algorithm for the studied system is developed in Section 3; the corresponding results are shown in Section 4. Finally, the conclusion is drawn in Section 5.</p></sec><sec id="s2"><title>2. Modeling</title><p>The heating units in the building are hot-water radiators distributed among the rooms and fed by hot water via inlet electro-valves which can control the water flow rate. The building model was developed in [<xref ref-type="bibr" rid="scirp.75238-ref4">4</xref>], what-so-called T<sub>i</sub>T<sub>e</sub> model as a state space model that considers the building envelope temperature as well as the indoor air temperature.</p><sec id="s2_1"><title>2.1. ASHP Model</title><p>The ASHP heat output rate Q<sub>h</sub> can be calculated by Equation (1) where P<sub>el</sub> is the ASHP electric power consumption.</p><disp-formula id="scirp.75238-formula87"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75238x2.png"  xlink:type="simple"/></disp-formula><p>The model considers the air source temperature dependency in a linearized relation. Then the COP and Q<sub>h</sub> can be given by Equations (2) and (3), where COP<sub>0</sub>, Q<sub>0</sub>, a<sub>1</sub> and b<sub>1</sub> are identified according to manufacturer’s specifications [<xref ref-type="bibr" rid="scirp.75238-ref7">7</xref>].</p><disp-formula id="scirp.75238-formula88"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75238x3.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75238-formula89"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75238x4.png"  xlink:type="simple"/></disp-formula></sec><sec id="s2_2"><title>2.2. TES Model</title><p>Thermal stratification is an important concept when dealing with hot water storage, which results from the density difference between the warm water and the cold water causing the hot water to ascend to the top. Studies as in [<xref ref-type="bibr" rid="scirp.75238-ref8">8</xref>] showed that improving the thermal stratification can significantly improve the system efficiency. The TES model is developed based on the concept illustrated in [<xref ref-type="bibr" rid="scirp.75238-ref9">9</xref>]. <xref ref-type="fig" rid="fig1">Figure 1</xref> shows a stratified tank with L meters height and N layers. The hot water is fed from the production side to the top layer with <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x5.png" xlink:type="simple"/></inline-formula> mass flow rate at temperature T<sub>hp</sub><sub>,out</sub> ˚C. The heat load is supplied from the top layer with flow rate of and the return is to the bottom layer to avoid thermal mixing.</p><p>The heat balance can be given by Equation (4). Where m is the mass of the hot water; C<sub>p</sub> is the water thermal capacity, and T<sub>w,i</sub> is the water temperature in the i<sup>th</sup> layer.</p><disp-formula id="scirp.75238-formula90"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75238x6.png"  xlink:type="simple"/></disp-formula><p>Beside the ASHP heat input<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x7.png" xlink:type="simple"/></inline-formula>, the load consumption <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x8.png" xlink:type="simple"/></inline-formula> and the losses to the environment <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x9.png" xlink:type="simple"/></inline-formula> the term <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x10.png" xlink:type="simple"/></inline-formula> represents the amount if heat exchanged in layer i with the surrounding layers i − 1 and i + 1 by natural convection and thermal conduction. In the simplest approximation, this can be summarized in effective vertical heat conductivity<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x11.png" xlink:type="simple"/></inline-formula>, and the heat flow between layer i − 1 and i, or from i to i + 1. The net heat flow for layer i of height z and cross section A<sub>q</sub> results in the difference of the two heat flows. The effective heat conductivity <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x12.png" xlink:type="simple"/></inline-formula> can be between 1 - 1.5 W/mK. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x13.png" xlink:type="simple"/></inline-formula>is calculated by Equation (5). The forced convection is represented by <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x14.png" xlink:type="simple"/></inline-formula> and can be calculated by Equation (6), where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x15.png" xlink:type="simple"/></inline-formula> is the mass flow rate.</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Stratified hot water tank with N layers</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/75238x16.png"/></fig><disp-formula id="scirp.75238-formula91"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75238x17.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75238-formula92"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75238x18.png"  xlink:type="simple"/></disp-formula><p>So T<sub>w,i</sub> can be calculated downward successively according the Equation (7).</p><disp-formula id="scirp.75238-formula93"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75238x19.png"  xlink:type="simple"/></disp-formula><p>where:</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x20.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x21.png" xlink:type="simple"/></inline-formula></p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x22.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x23.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x24.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x25.png" xlink:type="simple"/></inline-formula> are control parameters as follows:</p><p>ASHP is connected to the top layer, and the heat load is returned to the bottom layer:</p><disp-formula id="scirp.75238-formula94"><graphic  xlink:href="http://html.scirp.org/file/75238x26.png"  xlink:type="simple"/></disp-formula><p>Energy input from layer i − 1 to layer i; and the energy input from layer i + 1 to i:</p><disp-formula id="scirp.75238-formula95"><graphic  xlink:href="http://html.scirp.org/file/75238x27.png"  xlink:type="simple"/></disp-formula></sec><sec id="s2_3"><title>2.3. Hot-Water Radiator Model</title><p>The hot-water radiator heat transfer rate and the hot water mass flow rate are calculated based on Equations (8) and (9) respectively. The design values are based on the EN 442 [<xref ref-type="bibr" rid="scirp.75238-ref10">10</xref>] and the manufacturer’s technical manual. Where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x28.png" xlink:type="simple"/></inline-formula> is the radiator design heat transfer rate; <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x29.png" xlink:type="simple"/></inline-formula>is the design surface temperature; <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x30.png" xlink:type="simple"/></inline-formula>is the design indoor temperature; <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x31.png" xlink:type="simple"/></inline-formula>is the radiator return temperature; <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x32.png" xlink:type="simple"/></inline-formula>is the flow temperature, n is the correction factor and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x33.png" xlink:type="simple"/></inline-formula> is the mass flow rate.</p><disp-formula id="scirp.75238-formula96"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75238x34.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75238-formula97"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75238x35.png"  xlink:type="simple"/></disp-formula></sec></sec><sec id="s3"><title>3. Case Study</title><p>The studied system is shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>; the ASHP extracts the heat from the air and supplies the TES with hot water. Upon request, the hot water flows through the radiators transferring the required heat by convection to the ambient in each floor in the building, then the cold water flows back to the TES.</p><p>In order to utilize both storage capacities, i.e. the building mass and TES, the</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> A generic description for the studied system</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/75238x36.png"/></fig><p>EMPC algorithm which is shown in <xref ref-type="fig" rid="fig3">Figure 3</xref> solves two sequential optimization problems for every sampling time within the prediction horizon. The PFH3 algorithm determines the radiators’ output to supply the building with the required heat; it is the same as developed in [<xref ref-type="bibr" rid="scirp.75238-ref4">4</xref>]. The second algorithm is to optimize the ASHP electric power consumption to set the TES temperature to the optimal temperature within the operating boundaries.</p><p>The objective function of the ASHP optimization is given by Equation (10).</p><disp-formula id="scirp.75238-formula98"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75238x37.png"  xlink:type="simple"/></disp-formula><p>where: <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x38.png" xlink:type="simple"/></inline-formula>is the price vector for the prediction horizon N sampled every k period of time. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x39.png" xlink:type="simple"/></inline-formula>is the ASHP electric consumption, which is the manipulated input within the operational limits <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x40.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x41.png" xlink:type="simple"/></inline-formula>. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x41.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x42.png" xlink:type="simple"/></inline-formula>is a nonlinear function which represents the temperature of the stored water for the TES top layer, and it is calculated using the TES in Equation (7).</p><p>Particularly, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x43.png" xlink:type="simple"/></inline-formula>is used as the manipulated variable. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x44.png" xlink:type="simple"/></inline-formula>is the minimum stored temperature, in this study the boundary is optimized to meet the demand at a every certain time and not fixed to a constant value; this minimum value is calculated by utilizing the radiator model described in Equation (8) to calculate the optimal <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x44.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x45.png" xlink:type="simple"/></inline-formula> and it should be sufficient to the maximum required load among the floors. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x44.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75238x46.png" xlink:type="simple"/></inline-formula>is the maximum boundary and set to 100˚C temperature to keep the water from being evaporated.</p></sec><sec id="s4"><title>4. Simulation Results</title><p>The system is simulated using day-ahead electricity market price signals and weather data from the local weather station records for 5 days in winter; a reference case is a traditional PID system realized by setting the low and high temperature constraints at the same values to follow the lower comfort band temperature in the building using electric heaters without TES. Also the system is compared to another case in which the ASHP is directly connected to the radiators without the TES, this case is simulated using the PFH3 developed in [<xref ref-type="bibr" rid="scirp.75238-ref4">4</xref>] by utilizing the ASHP model.</p><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Functional block diagram for EMPC system</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/75238x47.png"/></fig><p><xref ref-type="fig" rid="fig4">Figure 4</xref> shows the radiators optimal power consumption, the price signal (bottom plot) and indoor temperature (top plot); it is obvious that the radiators consumption is shifted to the low price periods while the indoor temperature is controlled within the comfort constraints (dotted lines). In <xref ref-type="fig" rid="fig5">Figure 5</xref> the bottom plot shows the ASHP optimal electric power consumption shifted to the low price periods; and the temperature of the top layer in the TES (top plot), the lower red dotted curve is the optimal flow temperature calculated from the optimal heat output resulted from PFH3 for every sample time; we can observe the tank top layer temperature is kept within these varying constraints, even though the maximum temperature in the tank is limited to the ASHP maximum flow temperature (55˚C).</p><p>The savings are evaluated based on the reference case (the traditional PID) for 10 days simulation and shown in <xref ref-type="table" rid="table1">Table 1</xref>. Using the EMPC with electric heaters achieved 4.3% savings over the traditional PID control which reflects the effectiveness of shifting the loads intelligently to the low price periods. Using ASHP with hot-water radiators; the savings jumped significantly to 73.6% over the reference case when using ASHP with hot-water radiators because of the ASHP principle of operation. The studied system in this paper; incorporating the TES with the ASHP saved 7.9% more where the total saving incurred is 81.5% over the reference case.</p></sec><sec id="s5"><title>5. Conclusion</title><p>The simulation results of the EMPC algorithms effectively showed a successful shift for the electricity demand based on the price where the power consumption is shifted to the low price periods in the study case while the systems states are kept within the limits. In addition, for households, the monetary value of the</p><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Indoor temperature and optimal heaters’ output</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/75238x48.png"/></fig><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Simulation results for the ASHP and TES</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/75238x49.png"/></fig><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Energy costs for each scenario and savings based on the reference</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Scenario</th><th align="center" valign="middle" >Electricity Costin 10 days (Euro)</th><th align="center" valign="middle" >Savings %</th></tr></thead><tr><td align="center" valign="middle" >Reference Case (PID)</td><td align="center" valign="middle" >28.00 [<xref ref-type="bibr" rid="scirp.75238-ref11">11</xref>]</td><td align="center" valign="middle" >-</td></tr><tr><td align="center" valign="middle" >Case 1: EMPC with Electric Heaters</td><td align="center" valign="middle" >26.79 [<xref ref-type="bibr" rid="scirp.75238-ref11">11</xref>]</td><td align="center" valign="middle" >4.3%</td></tr><tr><td align="center" valign="middle" >Case 2: EMPC with HP</td><td align="center" valign="middle" >7.39</td><td align="center" valign="middle" >73.6%</td></tr><tr><td align="center" valign="middle" >Case 3: EMPC with HP and TES</td><td align="center" valign="middle" >5.19</td><td align="center" valign="middle" >81.5%</td></tr></tbody></table></table-wrap><p>energy savings is small relative to the complexity of the system and the cost of the equipment, so the payback period for this system could be very long for the small buildings with small heat loads; for the future work, a business case could be developed by assuming an agent between the consumer level and the grid operators; this agent can aggregate all the distributed energy sources’ flexibility and provide ancillary services to the grid operators.</p></sec><sec id="s6"><title>Cite this paper</title><p>Awadelrahman, M.A.A., Zong, Y., Li, H.W. and Agert, C. (2017) Economic Model Predictive Control for Hot Water Based Heating Systems in Smart Buildings. Energy and Power Engineering, 9, 112-119. https://doi.org/10.4236/epe.2017.94B014</p></sec></body><back><ref-list><title>References</title><ref id="scirp.75238-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">(2012) Energy Policy in Denmark. DanishEnergy Agency.</mixed-citation></ref><ref id="scirp.75238-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">E. Union. (1995-2016). European Commission Website, Energy Efficiency, Buildings. https://ec.europa.eu/energy/en/topics/energy-efficiency/buildings</mixed-citation></ref><ref id="scirp.75238-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Chen,C., Wang,J.,Heo,Y. and Kishore,S.(2013) MPC-Based Appliance Scheduling for Residential Building Energy Management Controller.IEEE Transactions on Smart Grid, 4.https://doi.org/10.1109/tsg.2013.2265239</mixed-citation></ref><ref id="scirp.75238-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Zong,Y.,B&amp;ouml;ning,G. M., Santos,R. M., You, S.,Hu, J. and Han,X. (2016) Challenges of Implementing Economic Model Predictive Control Strategy for Buildings Interacting with Smart Energy Systems. Applied Thermal Engineering.</mixed-citation></ref><ref id="scirp.75238-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Arteconi,A., Hewitt, N. J., and Polonara,F.(2012) State of the Art of Thermal Storage for Demand-Side Management. Applied Energy,Alsevier.</mixed-citation></ref><ref id="scirp.75238-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Arteconi,A., Hewitt,N. J., and Polonara,F.(2012) Domestic Demand-Side Management (DSM): Role of Heat Pumps and Thermal Energy Storage (TES) Systems. Applied Thermal Engineering, Elsevier.</mixed-citation></ref><ref id="scirp.75238-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">V. ViessmannWerkeGmbH&amp;Co KG, (2012) Heating System and Domestic Ventilation System with Heat Pump Control Unit, Vitotronic 200, type WO1C, Oprating Instructions Vitronic 200.</mixed-citation></ref><ref id="scirp.75238-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Lavan,Z. and Thompson,J.(1977) Experimental Study of Thermally Stratified Hot Water Storage Tanks. Solar Energy, 19, 519-524. 
https://doi.org/10.1016/0038-092X(77)90108-6</mixed-citation></ref><ref id="scirp.75238-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Eicker,U.(2003) Solar Technologies for Buildings. John Wiley &amp; Sons Ltd, Chichester.</mixed-citation></ref><ref id="scirp.75238-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">BS EN 442-1. (2014) Specification for Radiators and Convectors.</mixed-citation></ref><ref id="scirp.75238-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">B&amp;ouml;ning,G. M.(2016) Model Predictive Control (MPC)-based Energy Management of Smart Buildings. MSc, Regenerative Energiesysteme, Fakult&amp;auml;t III - Prozesswissenschaften, TechnischeUniversit&amp;auml;t Berlin.</mixed-citation></ref></ref-list></back></article>