• Energy Research

The aim of the paper is to present an energy system optimization method, based on the Fuel Impact Formula, able to overcome typical limitations of previous formulations. In particular, the methodology allows time dependent production levels to be considered and on-off operation and presence-absence of any component to be modelled by means of binary decision variables and inequality constraints. These often happen if the synthesis and operation of Combined Heat and Power (CHP) systems are considered at the same time. Moreover, the effect of replacing actual not linear thermodynamic inputs-output relations of each component with linear, or even proportional ones, is highlighted in the paper. Finally, an example of application is shown, dealing with a multi-component DH and CHP system in the tertiary sector.

AbstractIn the context of district heating networks we consider a model for the distribution of energy through water (for heating or cooling) from a central power station to the consumers. We prove the well posedness of the system, by using the Banach Fixed Point Theorem together with stability estimates for reduced systems. Eventually we consider optimal control problems motivated by applications and we provide the existence of optimal controls in special situations.

Decentralized surplus feed-in of solar heat into a District Heating Network (DHN) is here addressed. The heat collected from solar panels located on rooftops of DHN connected buildings may either be used locally for domestic hot water and space heating or fed into the DHN. Two-way substations able to transfer heat from and into the network seem then to be required utilities. The present paper presents the specifications (60kW capacity, return-to-supply connection) and promising architectures of such two-way substation based on a previous analysis. A first-of-a-kind Modelica-based dynamic model of the substation together with the consumer and the solar field connected to it is then detailed. Two-day simulations considering real operating conditions of DHN were then performed. The results highlighted i) the good match between the periods of solar heat reinjection with the periods of low supply temperature and differential pressure and ii) the decisive benefit of the reinjection to increase the part of useful solar energy.

In this study, a holistic energy, economic and environmental assessment was performed on a prosumer-based district heating system, including scenarios with varying district size, retrofitting stages and system configurations. A modeling framework was built which comprises a thermal network design and simulation model; a building energy demand model for districts; and supply and storage technology models that allow assessing system solar fraction, equivalent annual cost and greenhouse gas emissions of district heating systems (DHS). Furthermore, the approach allows comparing the performance of a DHS with individual heating systems (IHS) for the district with the same set of technology options (rooftop-mounted solar thermal collectors, gas boilers and thermal storage tanks). The framework was applied to a Swiss case study. The results of the case study show that DHS often outperform IHS; however, parameters such as the district size, the retrofitting stage of buildings and the system configuration have impacts on the performance of the DHS. The most important parameter lies in the adequate selection of the storage volume over solar collector area ratio, which indicates that DHS solutions are only advantageous if they are properly sized. Smaller districts and districts with retrofitted buildings especially benefit from DHS solutions in terms of energy, economic and environmental performance. Maximum solar fractions of 50% (non-retrofitted case) and 63% (retrofitted case) were reached with the DHS solutions. Energies, 14 (4) ISSN:1996-1073

A data-set containing the structure, input parameters and optimization results of the heating networks used in benchmark of the paper "Non-linear Topology Optimization of District Heating Networks: A benchmark of Mixed-Integer and Adjoint Approaches". This data can be used to replicate the results of the publication.

Actuellement, les installations CAD sont souvent surdimensionnées par manque de connaissance pour garantir les besoins thermiques de tous les clients. Ce dimensionnement se fait pour un point de fonctionnement à la puissance maximale, puissance à fournir seulement de 1 à 5 jours par an. Une optimisation énergétique grâce à un logiciel pourrait permettre une économie des matières premières des conduites, des déperditions thermiques et de temps lié au dimensionnement. Le but de ce projet est de développer un algorithme itératif en C++, afin de calculer les pertes de charge en régime laminaire et turbulent selon le besoin de chaleur des clients dans un réseau ramifié. L’algorithme considère le réseau comme un arbre composé des points, selon un fichier d’entrée XML. Afin de simplifier la structure, seule la conduite d’aller est considérée. En plus seuls les points des divisons ou des consommateurs sont considérés comme importants dans une première version du Solver. Ces derniers sont nommés feuilles. Les entrées sont principalement l’énergie consommée par chaque sous-station, les longueurs des conduites, les propriétés des conduites et les températures d’aller et de retour du réseau. Les sorties du Solver se trouvent sous forme de fichier texte : débit, vitesse de l’écoulement, perte de charge, différence de température pour chaque branche, température et pression à chaque feuille. L’outil développé très simplifié permet d’imprimer les résultats pour chaque feuille principale dans la conduite d’aller et la visualiser dans QGIS. L’algorithme a été testé sur un cas d’étude réelle. Le résultat du Solver par rapport au tableur Excel actuellement employé par un bureau d’ingénieur a pu montrer son potentiel de précision et sa vitesse de résolution. Tous les résultats visuels sont cohérents par rapport aux résultats calculés.

Dans la continuité du projet de semestre 5 « Simulation et optimisation SIG du réseau CAD de Broc », il a été mis en évidence la nécessité d’optimiser la prise de mesures aux sous-stations ainsi qu’à la centrale de chauffe de Broc. Malheureusement, après avoir analysé les relevés de mesures aux sous-stations en conditions réelles, il s’est avéré que les valeurs ne peuvent pas être exploitées en vue de valider le Solver Eguzki. C’est pourquoi des mesures ont été prélevées directement sur les écrans des automates à une des sous-stations du réseau de chaleur à distance. Il faut savoir que les sous-stations composant le réseau CAD sont dépourvues d’instruments permettant de relever des mesures de pression. Il s’est avéré judicieux d’installer des sondes de pression à la sous-station précédemment citée en vue de valider les pertes de charge au niveau de l’échangeur et de la vanne 2 voies dépendant notamment du débit et du coefficient kvs. Après avoir valider le Solver Eguzki en effectuant diverses simulations se basant notamment sur les mesures susmentionnées, il a été possible de définir que l’algorithme pouvait atteindre une précision de ±12.8%. Ensuite, il a été question de chiffrer les gains économiques et énergétiques qui pouvaient être réalisés en exploitant le Solver Eguzki. En se basant sur des puissances horaires simulées et obtenues grâce au logiciel CitySim, diverses simulations ont été effectuées avec Eguzki. A l’aide de celles-ci, il a été possible de déterminer la nouvelle courbe MCR des pompes et de la comparer à celle actuellement exploitée à la centrale de chauffe de Broc. Les différents points de fonctionnement obtenus par simulation et l’équivalence pour le modèle de régulation actuel ont été insérés sur le site du fournisseur des pompes afin de déterminer la puissance électrique nette demandée par les pompes. Cette première optimisation a permis de définir que la courbe de régulation actuelle n’est pas adaptée et que les coûts électriques annuels peuvent être diminuer d’au moins 77.0% en utilisant la courbe MCR des pompes déterminée à l’aide du Solver Eguzki. Dans une deuxième phase, une optimisation dimensionnelle a été effectuée en se basant sur le jour type le plus critique déterminé sur la base des puissances horaires simulées avec CitySim. A partir de simulations effectuées avec Eguzki, les diamètres des conduites constituant le réseau CAD ont pu être modifiés en vue d’atteindre un équilibre optimisé entre les déperditions thermiques dans le sol et les pertes de charge que doivent remonter les pompes à la centrale de chauffe.