Economic Analysis of Heat Distribution Concepts for a Small Solar District Heating System

: One challenge in today’s district heating systems is the relatively high distribution heat loss. Lowering distribution temperatures is one way to reduce operational costs resulting from high heat losses, while changing the distribution system from steel pipes to plastic pipes and changing the heat distribution concept can reduce investment costs. The result is that the overall life cycle cost of the district heating system is reduced, leading to the improved cost competitiveness of district heating versus individual heating options. The main aim of this study was to determine the most cost-efﬁcient distribution system for a theoretical solar district heating system, by comparing the marginal life cycle cost of two different distribution systems. A secondary aim was to determine the inﬂuence of the employed pipe type and insulation level on the marginal life cycle cost by comparing detailed economic calculations, including differences in pipe installation costs and construction costs, among others. A small solar-assisted district heating system has been modeled in TRNSYS based on a real system, and this “hybrid” model is used as a basis for a second model where a novel distribution system is employed and the heating network operating temperature is changed. Results indicate that a novel distribution concept with lower network temperatures and central domestic hot water preparation is most efﬁcient both from an energy and cost perspective. The total life cycle costs vary less than 2% for a given distribution concept when using different pipe types and insulation classes, indicating that the investment costs are more signiﬁcant than operational costs in reducing life cycle costs. The largest difference in life cycle cost is observed by changing the distribution concept, the novel concept having approximately 24% lower marginal life cycle cost than the “hybrid” system.


Introduction
This section sets the premise for the work presented in this article by treating the background of the topic and placing the paper in context by reviewing previous literature.The aims and scope are also outlined, together with a graphical presentation of the overall method used to answer the research questions.

Background
District heating (DH) was invented to recycle energy that would otherwise go to waste.With the onset of the industrial revolution and the employment of heat energy derived from coal, the first generation of DH networks naturally utilised high-temperature steam as a medium for heat transport.However, these networks led to high distribution heat losses and low system efficiency, a deficit of the distribution technology that endured until the third generation of DH that we have today.Presently, operating temperatures are significantly lower than in the past, allowing for a host of new industrial and renewable heat sources alike, while the pipe insulation has improved to reduce heat losses.Although Energies 2022, 15, 4737 2 of 29 the fourth generation of DH technology is still being studied by research in the field, the foreseen improvements include even lower operating temperatures, enabling a wide variety of waste-heat sources and a larger share of renewables [1].
Space heating (SH) and domestic hot water (DHW) are two low-temperature application areas often supplied by DH systems in the EU, particularly in the urban areas.The residential sector is responsible for about 50% of the final energy demand in the EU, out of which the vast majority is SH and DHW demand, which is supplied mainly by fossil fuels.The DH share of the EU28 heating demand was about 12% in 2015, out of which only a very small (<1%) share was supplied by solar [2].The temperature requirements of SH and DHW preparation correspond well to the operating temperature range of solar thermal (ST) energy, and most (>90%) of the installed ST systems are used in small and medium-sized hot water systems in the residential and public sector.However, by the end of 2019, less than 3% of all ST collectors in Europe were installed in large-scale systems (e.g., for DH) [3].Although as much as 22% of total energy use for heating and cooling in the EU28 was supplied by renewable energy by the end of 2019 [4], only a small (~2%) part of the thermal energy demand in European buildings was covered by ST [5].
Solar thermal (ST) energy has been used in several DH systems with success, the earliest examples of which were realised in the 1980s.In recent years, SDH has matured into a viable alternative for DH companies around Europe.Denmark is today the market leader, closely followed by Germany, Sweden and Austria.The economic feasibility of SDH in these countries has largely been determined by energy policy so that ST energy has achieved parity or even surpassed fossil fuel alternatives in terms of cost per unit of heat [3].However, the economic viability of DH has always been subject to line heat density, favouring urban areas and limiting its use in rural and suburban environments.Paradoxically, with the development of low-energy housing and political aims to reduce residential sector carbon footprints [6], the heat density of present and future residential areas is set to reduce significantly, which increases the fraction of distribution heat losses and reduces the DH system efficiency.This reduces the economic competitiveness of DH in comparison to other heating methods, and due to this, research efforts are necessary to identify the most efficient distribution options available for low-energy housing areas [1].Despite previous efforts to identify essential improvements in the future district energy system [7], case studies on implementation are still needed for proof of concept [8].
Today's solar-assisted district heating (SDH) plants typically employ low-temperature flat-plate solar collectors with gas or biomass-fired boilers to achieve a solar-fraction (SF) of 20-50% in the annual heat supply.Solar fractions in the upper half of this range are typically achieved by employing centralised seasonal storage, whereas the values in the lower half are usually achieved by using centralised or decentralised diurnal storage [9,10].The use of decentralised short-term storage may allow for more strategic use of the heating network and increase the efficiency of the DH system by enabling the use of ST energy locally and reducing distribution heat losses associated with circulation flows needed to uphold the target network operating temperatures.This may be particularly true in so-called hybrid networks, where 3rd and 4th generation DH technology is combined using both high-temperature steel pipes and low-temperature plastic pipes.Operational experiences from these types of systems as well as preliminary economic evaluations show that these networks have advantages over conventional 3rd generation DH networks, although it has been indicated that there may be economic advantages to the use of pure low-temperature plastic networks for heat distribution when compared to the hybrid networks [11,12].This is the premise of this study.

Previous Work
Early research into distribution networks for suburban housing areas indicated that plastic pipes could be a cost-efficient way to extend networks of higher heat demand densities.Studies in Sweden particularly focused on how to increase the economics of sparse district heating and concluded that plastic pipes were well suited to increase Energies 2022, 15, 4737 3 of 29 profitability in areas with low heat density, especially if used in a secondary network working as an extension to a conventional primary network [13].This conclusion was re-iterated recently when it was stated that cross-linked polyethylene (PEX) pipes encased in evacuated polystyrene (EPS) (abbreviated EPSPEX culvert) have the lowest operating costs for distribution systems in areas of low heat demand density [14].Another study linked the EPSPEX culvert to best practice in low-temperature district heating (LTDH) [10], while a recent study stated the use of PEX culverts as an essential improvement in 4th generation DH systems [7].
Until the present, the pipe type traditionally employed in DH networks has been made from steel and, although modern pipes employ polyurethane (PUR) foam for insulation, the developments in distribution technology have been limited to increasing the insulation level.One of the precursors for this may be the re-emergence of PEX pipes for distribution systems after their use was largely discontinued in the 80s [13].However, it has been indicated that the cost difference of these pipes compared to their steel counterparts has been reduced [15], which may have limited their employment.In recent technical developments, a new pipe type has emerged for low-temperature applications.The pipe is employed for LTDH, using a PUR insulated polyethylene (PE) pipe with some of the best features of steel and PEX pipes combined [16].Nonetheless, as the use of this pipe is part of an ongoing research project, the research results are still not known, but the intention is to provide the industry with a new solution for future DH systems.
However, as much as new technology can be a solution to an old problem, old technology used in new ways can be a solution too.The GRUDIS (Swedish for GRUppcentralDIstri-butionsSsystem) distribution concept from the 1980s was developed as a solution for sparse DH networks, where conventional steel pipes would prove too expensive.As such, it is an old solution to the relatively new problem of low heat density in urban environments.Due to its characteristics, it may be considered a 4th generation DH technology and should be suitable for low-temperature district heating (LTDH), according to desired characteristics specified in [7,17].The main characteristics are plastic (PEX) pipes and domestic hot water circulation (DHWC) with direct draw off DHW volume from the pipe without hydraulic separation by heat exchanger (HEX).The advantages of this concept are much the same as for the previously mentioned PE-pipe, in that it features simple installation without welding, which leads to faster and less costly installation.
A previous report [13] has recommended GRUDIS for use mainly in secondary networks, hydraulically separated from conventional (steel-pipe) primary networks with higher working temperatures and pressures.A case study has previously been made on such a system, which is the basis for the study presented here.The initial performance and inner workings of this system were studied in detail in [18], and the performance was later verified in [19].Subsequent studies have focused on potential energetic improvements in the use of the concept by moving towards a pure LTDH system and removing the hightemperature (HT) primary network altogether [20].Detailed continuations of these studies have solidified in increasing detail that, although an LTDH system approach has similar energetic performance as a combined network, it could be the economically preferred solution [11,12].The premise of both these studies includes focusing on the GRUDIS concept and the employment of PEX pipes.The advantages presented for GRUDIS, therefore, automatically include the use of plastic pipe distribution.
Although numerous studies exist on the detailed performance of steel DH pipes and optimal properties of these, few known studies treat detailed economical calculations on the use of PEX pipes versus steel pipes.According to one report [13], several research reports were made in the 1990s on this topic, although only one of these has been found [21].This study states that the costs for material, installation and groundwork, are lower for PEX pipes than for steel and that, due to diminishing cost savings for large pipe dimensions, these are better suited for low heat-density areas.One report published a decade later makes a simple comparison of EPSPEX versus steel pipes.However, this is based solely on the difference in catalogue values for specific pipe heat loss for "equivalent" pipe dimensions Energies 2022, 15, 4737 4 of 29 and an assumed cost for the heat energy loss.The EPSPEX culvert has a 20% lower annual cost than twin steel pipes, although more detailed calculations are recommended [22].The employed installation costs are also stated to be "rough" estimates, which reduces the generality of the results.
The essence of these few studies was that there is a cost-benefit when using plastic pipe systems in terms of lower heat losses and installation costs.However, despite this, the employment of PEX pipes is still rare in modern DH systems.Older studies have shown that the investment costs of the heating system are the driving factor of the system life cycle cost (LCC) [23].This fact may be well known by industry professionals, so the costs of distribution pipes and associated heat losses may have been neglected.This would explain the low adoption of PEX in the DH sector.Furthermore, presently, with the development of LTDH, the focus of the research has largely focused on operational and technological aspects of the LTDH system and, to a negligible degree, on the distribution network.Surprisingly, studies with a more holistic approach-those that evaluate heating systems and distribution networks together-have not been found.Therefore, studies that holistically highlight the impact of the chosen distribution concept, as well as the most appropriate pipe type, are necessary.
The novelty of this study lies in the presentation of detailed cost calculations both for an existing and novel system solution, comprising a combination of 3rd and 4th generation DH system technology, as well as for a proposed 4th generation system solution for LTDH.Furthermore, the study examines various pipe types suitable for these system solutions and presents detailed costs to quantify the size of these relative to the total system cost.Lastly, the influence of insulation level on the cost competitiveness of the two system solutions is evaluated to clarify the benefits of one over the other.None of this has been found in previous studies, which is why this study is considered a significant contribution to knowledge in the field of solar-assisted district heating.

Aim and Scope
The main research questions of this paper are: how large are the differences in marginal LCC when using two different distribution concepts for a hypothetical DH system, and how much do the pipe type and insulation level of the employed pipes influence this marginal LCC?
An existing solar DH (SDH) system is used as a basis for a theoretical SDH system, and a simulation model with similar technical specifications is developed.A hybrid system combines a HT steel pipe primary network and a LT PEX pipe secondary network, where the GRUDIS system with DHWC and central DHW preparation in intermediate substations is employed for secondary distribution.An alternative system model employs PEX pipes and GRUDIS distribution throughout the system.Simulations of the two distribution concepts are made using a reference insulation level for both to establish the impact of changing distribution systems on the energy performance of the system.A parametric study is made where the two systems are simulated, varying insulation levels for all distribution pipes and pipe types for PEX pipes, aiming to determine the most techno-economic solution.This is done by providing energetic considerations, followed by calculations of marginal LCC for the two distribution concepts using a reference insulation level and pipe type, as well as for several variants with different pipe types and/or insulation levels, to determine the most economic distribution concept and insulation level.Finally, because the assumed values for the financial boundary conditions have large uncertainty, a sensitivity analysis is made by a variation of financial parameters to determine the potential influence of these on the results from the parametric study.
The study is based on Swedish conditions, concerning both weather and costs.The DH system studied has a low line heat density (LD) as it supplies a load comprised of low-energy houses with a special form of heat distribution.This study does not include conventional DH distribution with pre-insulated steel pipes, higher operating temperatures, and a local heat exchanger for DHW preparation, as this has been covered in a previous Energies 2022, 15, 4737 5 of 29 study and, according to cost estimations, was much more costly than the distribution concepts covered in this study.One reason is that the house substation in conventional DH systems is more costly than the one needed in a GRUDIS system, as the GRUDIS system needs no local DHW preparation.
In lack of a commonly adopted terminology to describe hydraulic pipe networks, this paper uses specific terminology to differentiate between pipe types.The distribution network consists of "main" pipes that are connected in series to supply the (parallel) "branches" of the system.Generally, branch pipes supply any load > 1 house, whereas individual houses are supplied by "service" pipes.The pipe systems as a whole are termed "conventional" when referring to pre-insulated steel pipes, "GRUDIS" when referring to plastic pipes, and "hybrid" when referring to a combination of the two.

Overall Method
Figure 1 shows how the work in this study was executed in order of sequence.A case study is made of a solar-assisted DH system, from which two theoretical DH systems are defined using different distribution systems; one hybrid and one GRUDIS.The DH networks were dimensioned by using the technical specifications of the case study to determine the design heat load in the different areas of the DH network, from which flow rates were calculated using hydrodynamic theory.The design heat losses of the pipes were calibrated to design values by simulating the pipe performance under the same specific test conditions listed in manufacturer catalogues.Based on this, models are made in TRNSYS 17 (v.17.02.0005) [24], which includes houses, house substations, intermediate substations, primary and secondary distribution networks, as well as solar collectors and a boiler central (BC).These are denoted variants H1-H8 and G1-G4, where H and G stand for "hybrid" and "GRUDIS", respectively.
Energies 2022, 15, x FOR PEER REVIEW 5 of 30 conventional DH distribution with pre-insulated steel pipes, higher operating temperatures, and a local heat exchanger for DHW preparation, as this has been covered in a previous study and, according to cost estimations, was much more costly than the distribution concepts covered in this study.One reason is that the house substation in conventional DH systems is more costly than the one needed in a GRUDIS system, as the GRUDIS system needs no local DHW preparation.
In lack of a commonly adopted terminology to describe hydraulic pipe networks, this paper uses specific terminology to differentiate between pipe types.The distribution network consists of "main" pipes that are connected in series to supply the (parallel) "branches" of the system.Generally, branch pipes supply any load > 1 house, whereas individual houses are supplied by "service" pipes.The pipe systems as a whole are termed "conventional" when referring to pre-insulated steel pipes, "GRUDIS" when referring to plastic pipes, and "hybrid" when referring to a combination of the two.

Overall Method
Figure 1 shows how the work in this study was executed in order of sequence.A case study is made of a solar-assisted DH system, from which two theoretical DH systems are defined using different distribution systems; one hybrid and one GRUDIS.The DH networks were dimensioned by using the technical specifications of the case study to determine the design heat load in the different areas of the DH network, from which flow rates were calculated using hydrodynamic theory.The design heat losses of the pipes were calibrated to design values by simulating the pipe performance under the same specific test conditions listed in manufacturer catalogues.Based on this, models are made in TRNSYS 17 (v.17.02.0005) [24], which includes houses, house substations, intermediate substations, primary and secondary distribution networks, as well as solar collectors and a boiler central (BC).These are denoted variants H1-H8 and G1-G4, where H and G stand for "hybrid" and "GRUDIS", respectively.A parametric study is made where the insulation level is varied for the two distribution concepts, either by changing the insulation series and/or the pipe technology/manufacturer, and the energy balance (EB) for all system variants is calculated based on the simulation results.Furthermore, the marginal LCC of all system variants is calculated by including differences in initial capital costs (ICC) of construction due to differences in system configuration, based on cost estimates from entrepreneurial cost calculation software Wikells [25] and cost tenders from pipe manufacturers.Finally, a sensitivity analysis is conducted by changing the interest rate (IR) of the initial capital costs (ICC) for finance and by changing the annual fuel cost increase (FCI) for boiler fuel.
The simulation models have been simplified to make the study more general by scaling common subsystem components such as houses and intermediate substations to give the heat load of the entire subsystem.Furthermore, the distribution network has been modeled by using pipe elements to represent all pipes of the same size and type to give the same heat loss as the design value of the network.A parametric study is made where the insulation level is varied for the two distribution concepts, either by changing the insulation series and/or the pipe technology/manufacturer, and the energy balance (EB) for all system variants is calculated based on the simulation results.Furthermore, the marginal LCC of all system variants is calculated by including differences in initial capital costs (ICC) of construction due to differences in system configuration, based on cost estimates from entrepreneurial cost calculation software Wikells [25] and cost tenders from pipe manufacturers.Finally, a sensitivity analysis is conducted by changing the interest rate (IR) of the initial capital costs (ICC) for finance and by changing the annual fuel cost increase (FCI) for boiler fuel.
The simulation models have been simplified to make the study more general by scaling common subsystem components such as houses and intermediate substations to give the heat load of the entire subsystem.Furthermore, the distribution network has been modeled by using pipe elements to represent all pipes of the same size and type to give the same heat loss as the design value of the network.
This article is structured such that some content comprises a summary of information from detailed content contained in a data repository to improve readability.Therefore, every section where more details are available will have the data repository cited specifically [26].

Case Study
Two simulation models have been made in this study, both based on a hypothetical DH system.The systems supply 100 single-family houses, divided into four identical areas (A1-A4) of 25 houses each.This hypothetical system is based on the real system at Vallda Heberg, Sweden [10,[18][19][20]27,28], which has a hybrid distribution system.The two distribution systems configured for the hypothetical DH system are presented in this chapter.

Solar District Heating System Layout
Detailed information on pipe lengths and methods for calculating pipe sizes is available in an Appendix in a data repository related to [26].
Figure 2 [12] shows an overview of (a) the GRUDIS and (b) the hybrid distribution network.Evacuated tube collectors (ETC, red rectangle) are located on the boiler central (BC), which supplies the housing area directly in the GRUDIS system.In contrast, it supplies intermediate substations (SS1-SS4) through a conventional DH primary network in the hybrid system.These intermediate substations, in turn, supply each housing area (note that 1 house represents 5 houses) through a GRUDIS secondary network.Each housing area is named according to the substation it is connected to (only A1 is shown).In Figure 2, this conventional piping is marked in orange and GRUDIS piping in purple.A green solar culvert connects flat plate collectors (FPC, green rectangle) to the BC or intermediate substations.Note that the boiler central is larger in the GRUDIS system (stifled line) due to the absence of intermediate substations.
This article is structured such that some content comprises a summary of information from detailed content contained in a data repository to improve readability.Therefore, every section where more details are available will have the data repository cited specifically [26].

Case Study
Two simulation models have been made in this study, both based on a hypothetical DH system.The systems supply 100 single-family houses, divided into four identical areas (A1-A4) of 25 houses each.This hypothetical system is based on the real system at Vallda Heberg, Sweden [10,[18][19][20]27,28], which has a hybrid distribution system.The two distribution systems configured for the hypothetical DH system are presented in this chapter.

Solar District Heating System Layout
Detailed information on pipe lengths and methods for calculating pipe sizes is available in an Appendix in a data repository related to [26].
Figure 2 [12] shows an overview of (a) the GRUDIS and (b) the hybrid distribution network.Evacuated tube collectors (ETC, red rectangle) are located on the boiler central (BC), which supplies the housing area directly in the GRUDIS system.In contrast, it supplies intermediate substations (SS1-SS4) through a conventional DH primary network in the hybrid system.These intermediate substations, in turn, supply each housing area (note that 1 house represents 5 houses) through a GRUDIS secondary network.Each housing area is named according to the substation it is connected to (only A1 is shown).In Figure 2, this conventional piping is marked in orange and GRUDIS piping in purple.A green solar culvert connects flat plate collectors (FPC, green rectangle) to the BC or intermediate substations.Note that the boiler central is larger in the GRUDIS system (stifled line) due to the absence of intermediate substations.

Heat Supply
Table 1 shows an overview of the heat supply for the hybrid and GRUDIS distribution system simulated in this study.The DH is supplied with heat by a 300 kW woodpellet boiler and 108 m 2 evacuated tube collectors (ETC).Both boiler and solar collectors are connected to a 15 m 3 storage in the boiler central (BC).However, the location of flatplate collector (FPC) arrays and connected storage tanks varies according to the distribution concept.In the hybrid system-each intermediate substation (see Figure 2) features roughly 155 m 2 of FPC, divided into about 44 m 2 roof-mounted on the substation and about 111 m 2 ground-mounted (GM) in the vicinity of the substation, both arrays connected to a 15 m 3 storage.In the GRUDIS system, roughly 155 m 2 of GM collectors are located in the vicinity of each housing area.These large GM collector arrays can be seen as extensions of the GM arrays in the hybrid system, with the major difference being that heat from these is stored in 60 m 3 of storage tanks located in the boiler central.The storage employed consists of

Heat Supply
Table 1 shows an overview of the heat supply for the hybrid and GRUDIS distribution system simulated in this study.The DH is supplied with heat by a 300 kW wood-pellet boiler and 108 m 2 evacuated tube collectors (ETC).Both boiler and solar collectors are connected to a 15 m 3 storage in the boiler central (BC).However, the location of flat-plate collector (FPC) arrays and connected storage tanks varies according to the distribution concept.
Table 1.Heat supply-Overview of the heat supply configuration for the two distribution concepts investigated in this study.BC = boiler central, SS = Substation, GM = ground-mounted.

Parameters
Hybrid GRUDIS In the hybrid system-each intermediate substation (see Figure 2) features roughly 155 m 2 of FPC, divided into about 44 m 2 roof-mounted on the substation and about 111 m 2 ground-mounted (GM) in the vicinity of the substation, both arrays connected to a 15 m 3 storage.In the GRUDIS system, roughly 155 m 2 of GM collectors are located in the vicinity of each housing area.These large GM collector arrays can be seen as extensions of the GM arrays in the hybrid system, with the major difference being that heat from these is stored in 60 m 3 of storage tanks located in the boiler central.The storage employed consists of Energies 2022, 15, 4737 8 of 29 smaller, off-the-shelf, 5 m 3 accumulator tanks.This means that 15 m 3 storage consists of three storage tanks, while 60 m 3 consists of 12 storage tanks.

Distribution Systems
Two district heating systems supplying the same load type and size have been configured in this study, each using a different distribution concept.The energy system supply is the same for both systems, featuring the same installed solar collector area, solar heat storage volume, and boiler capacity.The distribution network route lengths are also the same, although the pipe type and technology vary.Table 2 shows the configuration of district heating subsystems in the two distribution concepts studied.A reference to the figure showing the subsystem schematic is written below the name of each distribution concept.Note that the design ground temperature is 10 • C. One advantage when using the GRUDIS concept is a fairly simple house substation, which lowers the investment costs.However, a drawback is that solar heat harvested in the distributed collector is stored centrally, which requires locating the solar buffer storage volume to the boiler central.This requires making the boiler central a bit larger, which induces some additional construction costs.Furthermore, the central storage location may lead to additional heat losses when solar heat is transported over a longer distance, which can lower the available solar energy and decrease system performance.

Hybrid Distribution Concept
Figure 4 [12] shows a schematic illustration of the hybrid heat distribution concept, with the boiler central to the right, an intermediate substation in the middle, and a passive single-family house to the left.The hybrid distribution system (see Figure 2b) consists of one boiler central and four intermediate substations.In contrast to the GRUDIS system The hybrid system comprises a combination of 3rd-and 4th-generation DH technology, using a novel system design to address the mismatch between the high operating temperatures seen in conventional DH systems and the lower operating temperatures wanted to increase the performance of ST collectors.The GRUDIS system is a suggestion for an alternative system configuration that better matches the optimal temperatures for increased performance of ST collectors and can be viewed as a 4th generation DH technology.The two systems are described more closely with a schematic in Sections 2.3.1 and 2.3.2.
(see Section 2.3.1), the solar buffer storage is located centrally and distributed in the substations.There are FPCs located on the substation roofs and in GM arrays nearby the substations, which supply heat to the substations.The boiler central contains a wood-pellet boiler and building integrated ETCs, which provide the intermediate substations with auxiliary heat to reach the target supply temperature of the building stock through a conventional steel pipe primary network.Heat is supplied to the buildings through a GRU-DIS secondary network, which was explained in Section 2.3.1.As the secondary network in this system employs the GRUDIS concept, one advantage of this configuration is a simple house substation with lower investment costs than in a conventional type of system.However, one disadvantage of the hybrid distribution concept is large heat losses in the distribution pipes and substations, which increases running costs.

GRUDIS
Figure 3 [12] shows a schematic illustration of the GRUDIS distribution concept, with the boiler central and central DHW preparation.In the GRUDIS concept, hot water is circulated in plastic pipes in a 2-pipe system [13], similar to how it is done in a DHWcirculation loop.Roof-mounted ETCs and ground-mounted FPCs supply a centrally located storage tank, which contains spiral heat exchangers used partially for DHW pre-heating or circulation pre-heating.A wood-pellet boiler supplies additional heat when required to reach the target supply temperature of the buildings.When DHW is tapped in the house substations, it is drawn directly from the circulation loop, and the draw-off volume is immediately replaced by make-up cold water in the boiler central.Space heat is supplied through floor heating in the bathroom of the houses and a ventilation air heat exchanger.
One advantage when using the GRUDIS concept is a fairly simple house substation, which lowers the investment costs.However, a drawback is that solar heat harvested in the distributed collector is stored centrally, which requires locating the solar buffer storage volume to the boiler central.This requires making the boiler central a bit larger, which induces some additional construction costs.Furthermore, the central storage location may lead to additional heat losses when solar heat is transported over a longer distance, which can lower the available solar energy and decrease system performance.

Hybrid Distribution Concept
Figure 4 [12] shows a schematic illustration of the hybrid heat distribution concept, with the boiler central to the right, an intermediate substation in the middle, and a passive single-family house to the left.The hybrid distribution system (see Figure 2b) consists of one boiler central and four intermediate substations.In contrast to the GRUDIS system (see Section 2.3.1), the solar buffer storage is located centrally and distributed in the substations.There are FPCs located on the substation roofs and in GM arrays nearby the substations, which supply heat to the substations.The boiler central contains a wood-pellet boiler and building integrated ETCs, which provide the intermediate substations with auxiliary heat to reach the target supply temperature of the building stock through a conventional steel pipe primary network.Heat is supplied to the buildings through a GRUDIS secondary network, which was explained in Section 2.3.1.
As the secondary network in this system employs the GRUDIS concept, one advantage of this configuration is a simple house substation with lower investment costs than in a conventional type of system.However, one disadvantage of the hybrid distribution concept is large heat losses in the distribution pipes and substations, which increases running costs.

Boundary Conditions and Heat Demand
This section contains some of the boundary conditions relevant to the simulations.More details are found in the data repository [26].
Figure 5 [12] shows a graph of cold water inlet supply temperature, monthly average ambient temperature, and global insolation for a typical meteorological year (TMY), together with the simulated monthly mean specific heat demand for a single-family house, in Kungsbacka, Sweden.The annual amount of global horizontal insolation is 984 kWh/m 2 .The simulated annual specific heat demand (SH + DHW) is ~38.5 kWh/m 2 , with values ranging from 5.6 kWh/m 2 (January) to 1.1 kWh/m 2 (July).The simulated total annual heat demand is about 542 MWh.The employed ground temperatures are not shown here but are calculated by the distribution pipe model during simulation.This is done using a sine function based on annual average surface temperature plus an amplitude representing the variation between highest and lowest daily average temperatures during the year, which was taken from the daily ambient temperature profile.

Boundary Conditions and Heat Demand
This section contains some of the boundary conditions relevant to the simulations.More details are found in the data repository [26].
Figure 5 [12] shows a graph of cold water inlet supply temperature, monthly average ambient temperature, and global insolation for a typical meteorological year (TMY), together with the simulated monthly mean specific heat demand for a single-family house, in Kungsbacka, Sweden.The annual amount of global horizontal insolation is 984 kWh/m 2 .The simulated annual specific heat demand (SH + DHW) is ~38.5 kWh/m 2 , with values ranging from 5.6 kWh/m 2 (January) to 1.1 kWh/m 2 (July).The simulated total annual heat demand is about 542 MWh.The employed ground temperatures are not shown here but are calculated by the distribution pipe model during simulation.This is done using a sine function based on annual average surface temperature plus an amplitude representing the variation between highest and lowest daily average temperatures during the year, which was taken from the daily ambient temperature profile.

DHW Profiles
The DHW profiles for the district heating system were generated with DHWcalc (v.2.02), using a 3-min random distribution [30].The draw-off profiles are based on a statistical profile of a multifamily building using 4 tap categories and model the seasonal variations by considering holidays.
Two DHW profiles are used-each one for 50 houses.The cold water inlet supply temperature is modeled using a sinus function for a temperature curve with an annual average of 12 °C and amplitude of 4 °C (January being the coldest-see Figure 5).

DHW Profiles
The DHW profiles for the district heating system were generated with DHWcalc (v.2.02), using a 3-min random distribution [30].The draw-off profiles are based on a statistical profile of a multifamily building using 4 tap categories and model the seasonal variations by considering holidays.
Two DHW profiles are used-each one for 50 houses.The cold water inlet supply temperature is modeled using a sinus function for a temperature curve with an annual average of 12 • C and amplitude of 4 • C (January being the coldest-see Figure 5).

Method
Section 3.1 of this chapter outlines the modeling approach.This has been provided more extensively in a previous study [12], and details of underlying calculations made to model the system can be found in the data repository [26].Section 3.2 provides an overview of the parametric study, while Section 3.3 provides the fundamental equations and inputs for calculating life cycle costs.Section 3.4 provides an overview of variables altered in the cost sensitivity analysis.Finally, Section 3.5 explains the employed key performance indicators (KPI) used to analyse the results of simulations and calculations.

Modelling Approach for Deriving Energetic KPIs
Figure 6 [12] shows the simplified model for the hybrid system, including the intermediate substation (SS) model and building model, together with applied pipe sizes.Due to software limitations in TRNSYS on the number of component outputs and the long simulation times related to models exceeding these, both simulation models used in this study are simplified in that they employ scaling of loads.The DH network modeled is symmetric (see Figure 2), so such a simplified approach can be assumed to affect the results insignificantly as this study focused mainly on comparing energy use between system variants.

Modelling Approach for Deriving Energetic KPIs
Figure 6 [12] shows the simplified model for the hybrid system, including the intermediate substation (SS) model and building model, together with applied pipe sizes.Due to software limitations in TRNSYS on the number of component outputs and the long simulation times related to models exceeding these, both simulation models used in this study are simplified in that they employ scaling of loads.The DH network modeled is symmetric (see Figure 2), so such a simplified approach can be assumed to affect the results insignificantly as this study focused mainly on comparing energy use between system variants.
The hybrid system model can be divided into three subsystem models: • Building and house substation model (SH load × 25 = one housing area).
The hybrid system was modeled by scaling up the house heating loads by a factor of 25 and connecting this load at the far end of a series of pipe elements that represents the secondary (GRUDIS) network.Each intermediate substation (SS) is scaled up by a factor of 2 to give the total load of the system, as seen from the BC.The simulated load was used to calculate the primary network flow rates by using the simulated supply and return temperatures at the substation, resulting in a realistic load in the primary network.[12].BC = Boiler central.SS = substation.

Common Components for All System Models
Table 3 lists the main TRNSYS components employed and parameter settings.shown in the substation models where these are used [12].BC = Boiler central.SS = substation.
The hybrid system model can be divided into three subsystem models:

•
Building and house substation model (SH load × 25 = one housing area).
The hybrid system was modeled by scaling up the house heating loads by a factor of 25 and connecting this load at the far end of a series of pipe elements that represents the secondary (GRUDIS) network.Each intermediate substation (SS) is scaled up by a factor of 2 to give the total load of the system, as seen from the BC.The simulated load was used to calculate the primary network flow rates by using the simulated supply and return temperatures at the substation, resulting in a realistic load in the primary network.

Common Components for All System Models
Table 3 lists the main TRNSYS components employed and parameter settings.[32].Data from [12].

Name Component Type Main Parameters Descriptions
Weather data Type 15 Kungsbacka, Sweden (57.5

Building and Space-Heating Model
A two-zone building model is used, and its geometry is based on that of a real house.The house model takes into account external gains from sunlight through the windows as well as internal gains from occupants and electricity use.The inputs to the building model are supplied in building information files in the data repository [26].
The house is heated by an air heating circuit, where heat is supplied from the distribution network to a mechanical ventilation system with heat recovery, in addition to year-round passive floor heating.

Distribution Pipe Model
Both twin and single pipes use a buried horizontal pipe model, although only the twin pipe includes the thermal resistance of the casing material (HDPE80).For PEX pipes, the EPS insulation and the casing are the same, so no casing resistance is modeled.
The heat losses were calibrated against catalogue data for all pipes to give all of the pipes a common theoretical reference point and allow for a more even comparison between pipes.

Reference Pipe Types and Combinations
There are two different pipe types used in the pipe networks in this study, pre-insulated steel pipes and PEX pipes.For steel pipes, one manufacturer is used, whereas, for PEX pipes, two manufacturers are used-both employing a different type of insulation.For the sake of objectivity, manufacturer names are omitted in the content of this article but are mentioned in the acknowledgments section together with references.Generic names have instead been adapted, such as pipe manufacturer 1 (M1), 2 (M2), and 3 (M3).Pipe specifications (pipe dimensions, insulation thermal conductivity, etc.) are found in the data repository [26] for all manufacturers.
For the steel pipes, two insulation classes are used (series 1 and series 2).The dimensions and material parameters used were provided by M1.For PEX pipes, two different technologies are used; EPSPEX and pre-insulated PEX (PI-PEX).The EPSPEX culvert is made by M3, and only the standard insulation class has been used, as additional insulation classes are manufactured only upon request.This culvert comprises a pair of PEX pipes enclosed in a square/rectangular EPS-casing and is available with twin pipes for all pipe sizes used in this study.The PI-PEX pipes used are made by M1 and M2, which both produce pipes in two insulation classes.However, PI-PEX pipes are predominantly available as twin pipes up to pipe size DN63, above which pairs of single pipes must be used.
A reference combination of pipe types and insulation levels has been chosen for the two distribution concepts in order to enable a comparison between various pipe types and combinations in a parametric study (see Section 3.2).The reference pipe type and insulation level for the two distribution concepts are shown in Table 4.There is no common terminology for steel and PEX pipes regarding insulation class, wherefore the used terminology varies between the two pipe types.Proper modeling of steel and PEX pipe heat losses was ensured by calibrating the specific heat losses [W/m] against values provided in manufacturer product catalogs for a specified set of boundary conditions.For copper pipes, the theoretical overall heat transfer coefficient was calculated, and the specific heat loss was derived from this for manufacturerspecified boundary conditions.This approach was chosen as the copper pipes were used in the solar thermal system and located above ground so that the catalogue values were invalid.Table 5 overviews the boundary conditions used for the pipe calibration process in TRNSYS.These values have been specified according to EN13941 or EN15632, and the listed heat transfer coefficients (λ-values) have been calculated by the pipe manufacturer according to formulae specified in the standard employed (depends on the manufacturer).
Table 5. Boundary conditions-Input values used for calibration of pipe heat losses in TRNSYS; catalogue heat transfer coefficients λ for pipe, soil insulation, and ground, as well as supply/return temperatures and soil cover.

Pipe Type
T supply The calibration process was performed in TRNSYS with a simulation period and time-step of 744 h and 0.05 h, respectively.The range of pipe sizes used in the distribution network and solar culvert were simulated iteratively under the specified boundary conditions.The insulation heat transfer coefficient (λ ins ) was adjusted for each pipe size until the simulated and catalogue or calculated (copper pipes) values for specific heat loss were equal.The calibrated values of the heat transfer coefficients, which were the values used in the simulation of the DH system, can be found in the calibration datasheet.This datasheet (MS Excel), which includes calculations for the copper pipes and the product catalogues containing information used in the calibration process, has been added to the data repository [26].

Parametric Study
Table 6 lists the pipe combinations used to simulate variants of the hybrid distribution concept, listing pipes according to their application in either the primary or secondary network.Table 7 shows an overview of the pipe combinations used to simulate variants of the GRUDIS distribution concept.Each distribution concept has one reference pipe combination/insulation class, introduced in Section 3.1.4.The EPS-casing has a square cross-section, which means that the insulation class is not easily comparable to the insulation standards used for other pipes.No definition of the insulation class has been found, so this is defined as "not applicable" (NA).The EPS-casing has a square cross-section, which means that the insulation class is not easily comparable to the insulation standards used for other pipes.No definition of the insulation class has been found, so this is defined as "not applicable" (NA).
For the GRUDIS distribution concept, the pipe combinations used are limited to the two different PEX pipe technologies, so there are only four variants.
For all pipe sizes used in this study, M1 and M2 only provide single pipes for the largest pipe sizes used, DN90 and DN110.However, M2 offers two levels of insulation for all sizes, whereas only M1 provides a lower insulation level for the sizes DN63 and DN110.Therefore, the insulation level can be mixed and consist of two different insulation levels when using PI-PEX pipes from M1.This has been shown in Tables 6 and 7 as a "mixed" insulation class.

Calculation of Life Cycle Cost
The hybrid system modeled in this study was based on a real SDH system, which has been treated extensively in previous publications.One IEA report used this system as a "best-case" example for solar-assisted DH.Economic information provided cost estimates of SEK 340 million (EUR ~34 million) for the entire residential area, including the heating system, out of which the solar thermal system cost was estimated at SEK 3.3 million (~EUR 0.33 million) [10].However, these costs are not certain and detailed costs on both the solar energy system and the district heating subsystems are not included.Due to the difficulty in obtaining accurate and specific investment costs for the different parts of the system, the calculation of per-unit costs of heat such as LCOH is considered unfeasible in this study.
Instead, the LCC is adopted for use in this analysis.The LCC is calculated for the entire project lifetime and is the total cost when summarizing investment and capital costs, operation and maintenance (O&M) costs, re-investment costs, and residual values.Re-investment costs apply to the malfunction or replacement of components with a shorter technical lifetime than the project lifetime, while residual values are related to any components with a remaining technical lifetime at the end of the project lifetime.
The life cycle costs (LCC) can be defined according to Equation (1) [34]: where ICC is the initial capital cost, A is the annual O&M cost, R is the re-investment cost, and Res is the residual value.The LCC is calculated as the net present value (NPV) over the project lifetime.
The primary interest of this study was to compare two distribution concepts and, in total, 12 variants of distribution systems in terms of LCC.For this purpose, it was assumed that the annual maintenance costs were the same for both distribution concepts, while-for simplicity-the re-investment cost and residual value were zero.These three equation variables have therefore been neglected when using Equation (1).Furthermore, when evaluating the ICC, only the differences in construction costs between distribution concepts are included.To do this, the costs of the houses and solar collectors, as well as the boiler central illustrated in the hybrid concept (see Figure 2), are used as a basis, and only the difference in construction costs is calculated.Other ICCs related to the pipe network and financial services are included in full, as the distribution system is the main focus of this article.This is considered appropriate, as the costs of construction elements common to all variants only contribute to the absolute (total) costs, but not to evaluate the difference in costs of one distribution variant relative to another.Therefore, Equation (1) can be reformulated as presented in Equation ( 2): where LCC mrg designates a marginal LCC.This is because of the omission of common construction costs, so the LCC mrg becomes lower than the total (or actual) LCC.However, as a means of comparison, LCC mrg is considered a sufficient measure.For reasons of brevity, the marginal LCC will be referred to as simply LCC for the remainder of this article.

Summary of Economic Boundary Conditions
Table 8 shows an overview of cost types and categories included in the LCC (Equation ( 2)) for the GRUDIS and hybrid distribution concepts.These costs will be explained in more detail in the following subsections.It should be noted that the installation costs of the solar thermal system are assumed to be the same in both concepts so that only the pipe network costs of the solar culvert are included.This also applies to the costs of solar storage tanks and associated plumbing works.Furthermore, future increases in boiler fuel costs are not accounted for due to the high uncertainty associated with forecasting this.It has been assumed that the boiler fuel costs increase at the same rate as the target inflation rate, which should be 2% p.a. [36].This means that the resulting NPV of boiler fuel costs may be somewhat conservative.Detailed cost calculation summaries for the costs included in Equation ( 2) and Table 2 can be found in the data repository [26].

Initial Capital Costs (ICC)
Supplementary information about calculations is available in Appendix A of the data repository [26].
The ICCs are comprised of the initial construction costs for buildings, HVAC and electric costs for heating system components, installation costs of the pipe network, and financial costs associated with raising initial capital through a loan.Note that the HVAC costs do not include the cost of solar storage tanks and plumbing work related to the connection of these.The costs for storage and associated plumbing work are assumed to be the same for both distribution concepts, as the number of tanks is identical for both concepts.The cost calculation software Wikells [25] was employed to properly calculate the investment costs related to the construction of additional buildings needed for the different distribution concepts.The program was used under instruction from an experienced calculator (Rickard Åkesson, personal communication) [40].
As shown in Table 8, different construction and pipe network costs apply according to the distribution concept.
The applicable construction and pipe network costs calculated are outlined for each distribution concept as follows:

GRUDIS concept-higher boiler central and pipe installation costs
The GRUDIS distribution concept features a large boiler central to house all solar storage tanks.As such, the costs for the boiler central are larger than in the hybrid system.However, the building envelope of the boiler central is assumed to comprise the same materials as those used in the intermediate substations of the hybrid system.
The pipe network in the GRUDIS system comprises only PEX pipes, and for the preinsulated pipe type from M1 and M2, the largest employed pipes are only available as single pipes.The purchased length of pipe, therefore, varies depending on the pipe type and is lower for variant G4 than the others.Furthermore, double pipe trenches increase pipe installation costs for variant G1-G3 particularly, but also generally for the GRUDIS concept, due to parallel installation of the solar culvert with the main pipes (see Figure 2), which increases excavation costs.

Hybrid concept-additional costs for Intermediate substation and pipe network
The hybrid distribution concept features intermediate substations that act as a hydraulic separation between the primary and secondary network, in addition to providing space for solar storage and roof area for solar FPCs.There are four of these substations in the hybrid network, and the calculated additional investment costs that apply for the construction of one intermediate substation have been divided into the categories of building, HVAC, and electric costs.The HVAC components include a heat exchanger and other necessary components such as temperature sensors for the control system, an expansion vessel, a circulation pump, and internal piping for connecting storage, heat exchangers, and solar collectors.The labour costs associated with the installation of these components are also included.Building costs include walls, roof, and base slab, while electric costs are for lighting and connecting various components demanding power, such as control systems, pumps, etc.The primary network of the hybrid system consists of steel main and branch pipes.The costs of these are generally higher than for PEX pipes, as both the pipe material cost and the installation cost is higher for pipes of similar size.

Annual Operation and Maintenance Costs (A)
Supplementary information is available in Appendix B of the data repository [26].Maintenance costs were discussed with a representative for HVAC consulting company Andersson & Hultmark (P.A.-Jessen, Personal communication) [38], responsible for the design of the Vallda Heberg (hybrid) system on which this study is based.
The maintenance costs assumed are those related to the boiler, which include inspection, cleaning, and repairs.It was assumed that three person-hours are needed weekly for a boiler technician, resulting in an annual cost of approximately 50,000 SEK.Maintenance costs for house substations and intermediate substations were not included, as the house substations were assumed to be the responsibility of the house owners, and maintenance of the intermediate substations was assumed to be unscheduled, i.e., only performed when needed.
The calculations of boiler fuel costs were based on numbers for 2019 in an official report [39] from the Swedish Energy Authority.The calculation of annual boiler fuel demand was done by considering the theoretical boiler efficiency and the simulated boiler energy delivered to the fluid stream per hour, according to Equation (3).
where Q fuel, boiler is the boiler fuel supply, Q boiler is the boiler energy to a fluid stream, and η boiler is the calculated boiler efficiency at a given time and load.The boiler efficiency as a function of the load factor was calculated using Equation (4) [41]: where η boiler is the calculated boiler efficiency, η boiler,100% is the boiler efficiency at full load, K (=0.140) is an empirical constant, and β is the boiler load factor in percent.The boiler load factor is the quotient of Q boiler and maximum boiler capacity (Q max = 300 kW).The boiler modulation was calculated based on the minimum turndown ratio (MTR) of the boiler used in the original DH system on which this study is based (see Section 2).For every hourly value where the boiler output was above zero but lower than the minimum capacity, the output value was set to the minimum, and boiler fuel use was calculated accordingly.The hourly values were averaged over a 3 h period to account for potential fluctuations in boiler capacity within one hourly time step as a result of boiler start/stop and modulation behavior.

Cost Sensitivity Analysis
Due to the uncertainty in variables like interest rates, boiler fuel price, and cost estimates, it was interesting to study how the results change when these variables change.A cost sensitivity analysis has been made using a variation of these economic boundary conditions:

•
Economic lifetime: 20 years reference and 30 years extended.
A boiler fuel price increase of 2% corresponds to the inflation target of many central banks and represents a minimum increase in fuel price.The conversion of the energy system to mitigate climate change may increase the demand for biomass fuels for heating purposes, so it is interesting to study how higher fuel prices affect the LCC.Regarding the pipe installation costs for EPSPEX, it should be noted that these are uncertain.

Key Performance Indicators (KPI)
The key figures used for evaluating results on system performance are shown in

Results
This section is divided into three sections; Section 4.1 presents the energy balance for the reference variant of GRUDIS and the hybrid concept, while Section 4.2 shows the LCC for the same.Section 4.3 presents how energy balance and LCC changes when insulation level changes, while Section 4.4 shows how LCC changes under different economic boundary conditions.A boiler fuel price increase of 2% corresponds to the inflation target of many central banks and represents a minimum increase in fuel price.The conversion of the energy system to mitigate climate change may increase the demand for biomass fuels for heating purposes, so it is interesting to study how higher fuel prices affect the LCC.Regarding the pipe installation costs for EPSPEX, it should be noted that these are uncertain.

Key Performance Indicators (KPI)
The key figures used for evaluating results on system performance are shown in Table 9

Results
This section is divided into three sections; Section 4.1 presents the energy balance for the reference variant of GRUDIS and the hybrid concept, while Section 4.2 shows the LCC for the same.Section 4.3 presents how energy balance and LCC changes when insulation level changes, while Section 4.4 shows how LCC changes under different economic boundary conditions.From an energy perspective, the two distribution concepts seem to be performing similarly, both with regard to the energy supply and the energy demand.The main differences are observed in the solar storage and boiler central (BC) losses, where a reduction in one is accompanied by an increase in the other and vice versa.The solar storage loss is higher in the hybrid concept, partly due to the distribution of the storage volume into smaller entities located in intermediate substations around the network.These substations contain more internal piping and other hydraulic components than the centralised solar storage solution used in the GRUDIS concept, which leads to higher associated losses.Furthermore, as the FPC solar energy is only used for heat losses in the secondary network, some of the stored solar energy could be lost due to low demand in periods of high supply.On the other hand, the boiler central loss is higher in the GRUDIS concept, mostly due to the lower operating temperatures employed, leaving more of the heat in the boiler and ETC accumulator tanks unused.However, the combined losses in solar storage and boiler central are still about 11 MWh lower in the GRUDIS concept, although this is mostly reflected in lower FPC yield and has a low effect on the boiler energy demand.It is, therefore, not clear from a purely energetic perspective which one of these distribution concepts is preferable, considering the reference pipe combination (hybrid) and insulation class.

Life Cycle Cost
Figure 8 shows an overview of the calculated LCC of the two reference variants for the two distribution concepts simulated in this study.It is clear that there is a significant LCC gap between the GRUDIS and Hybrid distribution concepts.The value for LCC is SEK 15.5 M for G1 and SEK 20.4 M for H1, i.e., the LCC of G1 is 24% lower than H1.Using only PEX pipes give lower distribution pipe costs, although higher solar pipe costs lead to similar total pipe costs (SEK ~5.7 M) for the two concepts.As the boiler fuel costs also vary little between the concepts, the combined effect of higher construction costs and added financial costs as a result of this becomes an economic disadvantage of using the hybrid concept.The construction costs for expanding the boiler central in the GRUDIS concept are roughly SEK 0.4 M and make out about one-ninth (~11%) of the SEK 3.9 M applicable for the construction of four intermediate substations in the hybrid system.Due to this, the financial costs in the GRUDIS system are about 1.4 M SEK lower than those in the hybrid system, a difference of 35%.Based on these results, it is clear that the GRUDIS distribution concept is more economical.
From an energy perspective, the two distribution concepts seem to be performing similarly, both with regard to the energy supply and the energy demand.The main differences are observed in the solar storage and boiler central (BC) losses, where a reduction in one is accompanied by an increase in the other and vice versa.The solar storage loss is higher in the hybrid concept, partly due to the distribution of the storage volume into smaller entities located in intermediate substations around the network.These substations contain more internal piping and other hydraulic components than the centralised solar storage solution used in the GRUDIS concept, which leads to higher associated losses.Furthermore, as the FPC solar energy is only used for heat losses in the secondary network, some of the stored solar energy could be lost due to low demand in periods of high supply.On the other hand, the boiler central loss is higher in the GRUDIS concept, mostly due to the lower operating temperatures employed, leaving more of the heat in the boiler and ETC accumulator tanks unused.However, the combined losses in solar storage and boiler central are still about 11 MWh lower in the GRUDIS concept, although this is mostly reflected in lower FPC yield and has a low effect on the boiler energy demand.It is, therefore, not clear from a purely energetic perspective which one of these distribution concepts is preferable, considering the reference pipe combination (hybrid) and insulation class.

Life Cycle Cost
Figure 8 shows an overview of the calculated LCC of the two reference variants for the two distribution concepts simulated in this study.It is clear that there is a significant LCC gap between the GRUDIS and Hybrid distribution concepts.The value for LCC is SEK 15.5 M for G1 and SEK 20.4 M for H1, i.e., the LCC of G1 is 24% lower than H1.Using only PEX pipes give lower distribution pipe costs, although higher solar pipe costs lead to similar total pipe costs (SEK ~5.7 M) for the two concepts.As the boiler fuel costs also vary little between the concepts, the combined effect of higher construction costs and added financial costs as a result of this becomes an economic disadvantage of using the hybrid concept.The construction costs for expanding the boiler central in the GRUDIS concept are roughly SEK 0.4 M and make out about one-ninth (~11%) of the SEK 3.9 M applicable for the construction of four intermediate substations in the hybrid system.Due to this, the financial costs in the GRUDIS system are about 1.4 M SEK lower than those in the hybrid system, a difference of 35%.Based on these results, it is clear that the GRUDIS distribution concept is more economical.

Parametric Study
The parametric study shows how the energy balance (Section 4.3.1)and LCC (Section 4.3.2) vary with the employed insulation level of the pipe network for the GRUDIS and hybrid distribution concept.
Energies 2022, 15, 4737 20 of 29 4.3.1.Energy vs. Insulation Level Figure 9 shows how the supplied and lost energy changes with insulation level for the GRUDIS concept when compared to reference variant G1 (not shown).Variant G2 shows an increase in insulation level from "standard" to "plus" for the pipes (M2) used in the reference variant, which results in a reduction of 18% in distribution losses and 5% in boiler energy.Variant G3 and G4 represent a decrease and increase in insulation level, respectively.Thus, the standard insulation level in pipes by M1 is lower than that of M2 or M3.This is indicated by lower distribution heat losses for G2 and G4 and higher losses for G3 compared to the reference variant.For variant G3, the increase in distribution losses is about 16% and boiler energy 4%, while the decrease in losses is 7% and boiler energy is 2% for variant G4.

Parametric Study
The parametric study shows how the energy balance (Section 4.3.1)and LCC (Section 4.3.2) vary with the employed insulation level of the pipe network for the GRUDIS and hybrid distribution concept.

Energy vs. Insulation Level
Figure 9 shows how the supplied and lost energy changes with insulation level for the GRUDIS concept when compared to reference variant G1 (not shown).Variant G2 shows an increase in insulation level from "standard" to "plus" for the pipes (M2) used in the reference variant, which results in a reduction of 18% in distribution losses and 5% in boiler energy.Variant G3 and G4 represent a decrease and increase in insulation level, respectively.Thus, the standard insulation level in pipes by M1 is lower than that of M2 or M3.This is indicated by lower distribution heat losses for G2 and G4 and higher losses for G3 compared to the reference variant.For variant G3, the increase in distribution losses is about 16% and boiler energy 4%, while the decrease in losses is 7% and boiler energy is 2% for variant G4.  7 for overview.
It is evident from Figure 9 that the change in boiler supplied energy is lower than the change in distribution heat loss, largely due to changes in the solar energy yield.The average boiler energy is 591 MWh, less than a percent more than for the reference variant G1.For variant G2, the lower distribution heat loss leads to lower yield during the solar heating season as there is a surplus of heat, which is seen in the increase in solar storage loss.For variant G3, the heat losses increase.This leads to an increase in solar yield and simultaneously lower losses from boiler central and solar storage, which indicates that the increased loss is covered by solar.For variant G4, the distribution losses decrease but not sufficiently to increase the boiler central and solar storage loss, implying that all of the solar heat is utilised as efficiently as in variant G1.Furthermore, the solar yield is affected insignificantly, which indicates that the solar energy harvested during the daytime is just enough to cover the load and heat loss, giving equal changes in boiler energy and distribution loss.
Figure 10 shows how the losses and supplied energy in the hybrid distribution concept change with insulation levels relative to the reference variant H1 (not shown).System variants H2-H4 show increasing insulation levels for the reference combination of steel  7 for overview.
It is evident from Figure 9 that the change in boiler supplied energy is lower than the change in distribution heat loss, largely due to changes in the solar energy yield.The average boiler energy is 591 MWh, less than a percent more than for the reference variant G1.For variant G2, the lower distribution heat loss leads to lower yield during the solar heating season as there is a surplus of heat, which is seen in the increase in solar storage loss.For variant G3, the heat losses increase.This leads to an increase in solar yield and simultaneously lower losses from boiler central and solar storage, which indicates that the increased loss is covered by solar.For variant G4, the distribution losses decrease but not sufficiently to increase the boiler central and solar storage loss, implying that all of the solar heat is utilised as efficiently as in variant G1.Furthermore, the solar yield is affected insignificantly, which indicates that the solar energy harvested during the daytime is just enough to cover the load and heat loss, giving equal changes in boiler energy and distribution loss.
Figure 10 shows how the losses and supplied energy in the hybrid distribution concept change with insulation levels relative to the reference variant H1 (not shown).System variants H2-H4 show increasing insulation levels for the reference combination of steel pipes (M1) and PEX pipes (M2), while variants H5-H8 employ either M1 or M3 PEX pipes but alternates the insulation level of the steel pipes.An overview of the different variants is found in Table 6.
pipes (M1) and PEX pipes (M2), while variants H5-H8 employ either M1 or M3 PEX pipes but alternates the insulation level of the steel pipes.An overview of the different variants is found in Table 6.6 for an overview.
From Figure 10, it is clear that increasing the insulation level for steel pipes (i.e., H1  H2) has a smaller effect on distribution losses than when increasing insulation in PEX pipes (i.e., H1  H3).Employing M1 PEX pipes as in variant H6 corresponds to decreasing insulation level relative to H1 as the losses increase.On the other hand, the use of M3 (EPS) PEX pipes as in variant H7 is seemingly equivalent to using M2 PEX pipes with standard insulation level, as in variant H1.Variant H4 leads to the largest reduction in distribution losses (−28%) and boiler energy demand (−7%).The average boiler energy for all variants is 580 MWh, which is one percent less than for the reference variant H1.The standard deviation is about 3%, which indicates that the insulation level is of minor importance for the boiler energy demand.
By comparing Figures 9 and 10, and recalling that the distribution loss and boiler supplied energy were similar for reference variants G1 and H1 (see Figure 7), it can be seen that variant H4 is the most energy-efficient of any variant, regardless of the distribution concept.Furthermore, it appears that steel pipes with a low insulation level in the primary network of the hybrid system, as in variant H3, have nearly the same heat loss as if these pipes were PEX pipes with a higher insulation level, as in variant G2 (see Table 10).However, the additional overall heat loss in G2 is due to the additional network length necessary for the stem pipe sizes DN90 and DN110, which are only available as single pipes.If not for this, variant G2 would have the lowest losses of all variants.Nevertheless, variant H4 shows the lowest overall heat loss, followed by H3 and G2.That is, two variants of the hybrid concept lead to lower network heat losses than the best GRUDIS network.However, despite this, the variant G2 uses less boiler energy than variant H3, although the difference is too insignificant to give a definitive answer which is energetically preferable when considering uncertainty.6 for an overview.
From Figure 10, it is clear that increasing the insulation level for steel pipes (i.e., H1 → H2) has a smaller effect on distribution losses than when increasing insulation in PEX pipes (i.e., H1 → H3).Employing M1 PEX pipes as in variant H6 corresponds to decreasing insulation level relative to H1 as the losses increase.On the other hand, the use of M3 (EPS) PEX pipes as in variant H7 is seemingly equivalent to using M2 PEX pipes with standard insulation level, as in variant H1.Variant H4 leads to the largest reduction in distribution losses (−28%) and boiler energy demand (−7%).The average boiler energy for all variants is 580 MWh, which is one percent less than for the reference variant H1.The standard deviation is about 3%, which indicates that the insulation level is of minor importance for the boiler energy demand.
By comparing Figures 9 and 10, and recalling that the distribution loss and boiler supplied energy were similar for reference variants G1 and H1 (see Figure 7), it can be seen that variant H4 is the most energy-efficient of any variant, regardless of the distribution concept.Furthermore, it appears that steel pipes with a low insulation level in the primary network of the hybrid system, as in variant H3, have nearly the same heat loss as if these pipes were PEX pipes with a higher insulation level, as in variant G2 (see Table 10).However, the additional overall heat loss in G2 is due to the additional network length necessary for the stem pipe sizes DN90 and DN110, which are only available as single pipes.If not for this, variant G2 would have the lowest losses of all variants.Nevertheless, variant H4 shows the lowest overall heat loss, followed by H3 and G2.That is, two variants of the hybrid concept lead to lower network heat losses than the best GRUDIS network.However, despite this, the variant G2 uses less boiler energy than variant H3, although the difference is too insignificant to give a definitive answer which is energetically preferable when considering uncertainty.11 gives an overview of how the absolute costs change when insulation level changes, compared to reference variant G1 and H1 for the GRUDIS and hybrid distribution concept, respectively.According to the results, variant G4 has the largest net cost reduction compared to G1, which should largely be due to the difference in network length between variants G1-G3 and variant G4, where the latter uses twin pipes in larger dimensions (DN90 and DN110) in contrast to the former that use single pipes.The resulting increase in pipe costs due to a longer network offsets the low heat losses and correspondingly lower boiler energy costs.For the hybrid system, variant H4 yields the largest net cost reduction compared to H1. heat loss , calculated according to network length L, for the different distribution variants simulated in this study.Figure 11 gives an overview of how the absolute costs change when insulation level changes, compared to reference variant G1 and H1 for the GRUDIS and hybrid distribution concept, respectively.According to the results, variant G4 has the largest net cost reduction compared to G1, which should largely be due to the difference in network length between variants G1-G3 and variant G4, where the latter uses twin pipes in larger dimensions (DN90 and DN110) in contrast to the former that use single pipes.The resulting increase in pipe costs due to a longer network offsets the low heat losses and correspondingly lower boiler energy costs.For the hybrid system, variant H4 yields the largest net cost reduction compared to H1.  6 for an overview.

Variant
However, because the scale of the axis is two orders of magnitude lower than the scale in Figure 8, the relative influence on the total costs is in the lower one-digit range.Because the largest difference between variants within each distribution concept is less than 2%, the difference between variants within each concept is too insignificant to conclude whether it is worthwhile to invest in one over the other, particularly when considering the uncertainties in simulations and economic boundary conditions.Therefore, the average (based on all variants simulated) discrepancy in LCC for the GRUDIS distribution concept and the Hybrid concept is equal to the discrepancy between G1 and H1-24%.According to this, it can be concluded that the GRUDIS concept is more economic than the hybrid concept, regardless of insulation level.6 for an overview.
However, because the scale of the axis is two orders of magnitude lower than the scale in Figure 8, the relative influence on the total costs is in the lower one-digit range.Because the largest difference between variants within each distribution concept is less than 2%, the difference between variants within each concept is too insignificant to conclude whether it is worthwhile to invest in one over the other, particularly when considering the uncertainties in simulations and economic boundary conditions.Therefore, the average (based on all variants simulated) discrepancy in LCC for the GRUDIS distribution concept and the Hybrid concept is equal to the discrepancy between G1 and H1-24%.According to this, it can be concluded that the GRUDIS concept is more economic than the hybrid concept, regardless of insulation level.

Sensitivity Analysis
The sensitivity analysis investigates the effect of variation in the boundary conditions of the economic calculations, as described in Section 3.4.
Figure 12 shows the difference in LCC between distribution concepts for reference variants G1 and H1, as well as how the LCC change (increases) when economic boundary conditions are varied for reference variants G1, H1, and the most energetically efficient variants G4 and H4.According to the results, the absolute difference in LCC between the reference concepts remains quite similar regardless of fuel cost increase (FCI), although it increases with increasing IRs.This can be explained by observing that the relative cost increase is different for the hybrid and GRUDIS concepts.Costs increase slightly faster for GRUDIS variants G1 and G4 than they do for hybrid variants H1 and H4-with increasing fuel cost-while on the other hand, the GRUDIS variants are less sensitive to increasing IRs.This results from the differences in relative shares made out by individual ICCs in the LCC (see Figure 8).
of the economic calculations, as described in Section 3.4.
Figure 12 shows the difference in LCC between distribution concepts for reference variants G1 and H1, as well as how the LCC change (increases) when economic boundary conditions are varied for reference variants G1, H1, and the most energetically efficient variants G4 and H4.According to the results, the absolute difference in LCC between the reference concepts remains quite similar regardless of fuel cost increase (FCI), although it increases with increasing IRs.This can be explained by observing that the relative cost increase is different for the hybrid and GRUDIS concepts.Costs increase slightly faster for GRUDIS variants G1 and G4 than they do for hybrid variants H1 and H4-with increasing fuel cost-while on the other hand, the GRUDIS variants are less sensitive to increasing IRs.This results from the differences in relative shares made out by individual ICCs in the LCC (see Figure 8).Because the absolute boiler fuel costs are very similar for the two distribution concepts, and these costs make up a significantly larger relative share of the LCC for the GRU-DIS concept than for the hybrid concept, a similar absolute increase in boiler fuel costs has a larger relative impact on the increase in LCC.In the GRUDIS concept, the boiler fuel costs make out about 38% of the LCC, compared to 28% for the hybrid concept.Conversely, because the financial ICC is significantly larger in the hybrid concept than in the GRUDIS concept, a similar increase in IR has a larger relative impact on the LCC.In the hybrid concept, construction and financial ICC make out about 66% of the LCC, compared to 54% in the GRUDIS concept.
These results indicate that regardless of the boundary conditions, the GRUDIS distribution concept remains economically preferable and that the absolute cost savings are consistent over the studied range of boundary conditions.

Discussion
The discussion section mainly focuses on parts of the method and their influence on the overall results to reveal uncertainties in the results.Because the absolute boiler fuel costs are very similar for the two distribution concepts, and these costs make up a significantly larger relative share of the LCC for the GRUDIS concept than for the hybrid concept, a similar absolute increase in boiler fuel costs has a larger relative impact on the increase in LCC.In the GRUDIS concept, the boiler fuel costs make out about 38% of the LCC, compared to 28% for the hybrid concept.Conversely, because the financial ICC is significantly larger in the hybrid concept than in the GRUDIS concept, a similar increase in IR has a larger relative impact on the LCC.In the hybrid concept, construction and financial ICC make out about 66% of the LCC, compared to 54% in the GRUDIS concept.
These results indicate that regardless of the boundary conditions, the GRUDIS distribution concept remains economically preferable and that the absolute cost savings are consistent over the studied range of boundary conditions.

Discussion
The discussion section mainly focuses on parts of the method and their influence on the overall results to reveal uncertainties in the results.

Modelling
The modeling approach can have a significant influence on the simulation results, and this is discussed here.

Model Simplifications
Several model simplifications have been brought up in Section 3.1:

•
One housing area to represent two housing areas.

•
One house model to represent all houses.

•
Load connected to the supply at the end of the network.
Each of these simplifications influences the result to a minor extent, which has been extensively covered in other studies relating to this system model [11,12].The main conclusion from these was that because the modeling approach is the same for both distribution concepts, the potential effects of the modeling approach affect the results similarly for both.Hence, although the absolute values for energy quantities presented may not be entirely accurate for each respective concept, the differences in values between concepts should be considered accurate.

Boiler Modulating Behaviour
The boiler modulation in this study was modeled by post-processing of hourly simulation output.This required simplifying the data processing method, such as adapting a 3 h moving average (MA) value to account for hourly fluctuations in boiler output.This approach is considered appropriate due to the presence of buffer storage that provides peak shaving so that the boiler rarely works at minimum capacity.The absence of storage would render proper modeling of modulation more significant.
According to information from a technician working on the system used as a case study for the hypothetical system modeled in this study, the cycling time for a pellet boiler of the employed size is 15-20 min.That is, the time from the boiler turns off below the minimum turndown ratio, is started again, and arrives at the supply target temperature.The control response time is usually in a matter of minutes, so it is possible with 3-4 boiler cycles within one hour, resulting in potentially high instantaneous loads to give a low average hourly load.In particular, a jump in the hourly supply power from a value of, e.g., 75 kW to 0 kW-or vice versa-is unrealistic, as the boiler would down-modulate as far as possible before turning off to avoid extensive cycling.In general, all fluctuations in the part-load behaviour should be of importance, such as those that move between minimum load (22.5 kW) and up to about 1/3 boiler capacity, where the boiler efficiency stabilises.In this load interval, the efficiency varies the most, and, therefore, it seems most important to model the modulation behaviour correctly here.
Nonetheless, using a 3 h MA instead of the hourly values seems to have little influence on the calculated boiler fuel use.For the GRUDIS concept, the 3 h MA leads to a roughly 2% higher annual fuel use than when using hourly values.For the Hybrid concept, the 3h MA results in 1-3% higher annual fuel use than when using the hourly values.Variations in the base value for the MA in the range 3-6 h give less than a 1%-point difference in these values, which limits the importance of modeling boiler modulation more precisely.

Cost Calculations
The economic boundary conditions have the largest influence on the overall results.These are therefore subject to extra scrutiny when evaluating the results.This section focuses on the main uncertainties underlying the results and their influence on these.

Calculations from Price List vs. from Supplier
The pipe network prices were collected by contact with the sales department of each respective manufacturer.For M2 and M3, a network schematic and list of pipe sizes was supplied to a sales representative, and a tender was received with a part specification list.For M1, a price list was handed out with information about applicable discounts for the type of pipe (steel or PEX).The tender of M2 was used as a basis for the choice of pipe parts, such as joints and reductions, as these are similar to the pipe type (PI-PEX) for both manufacturers.
Two issues that arise with the employed method of cost information gathering are: 1.
Tenders from suppliers are without vendor margins (below market price).

2.
Possibility of tenders being undervalued (potential market strategy).
At first glance, both of these issues seem to be two sides of the same issue.However, issue 1 refers to the challenge of knowing how the tender relates to the market prices for pipes.Issue 2, on the other hand, refers to the possibility that manufacturer sales representatives-knowing the intended use of the cost information-could deliver lower tenders than would be the case for a client.
Although it is impossible to establish whether this happened, the small differences in pipe costs between different manufacturers observed in the LCC (see Section 4.3.2) implies that this is not the case.The similarity in costs further implies that the difference in cost calculation approach for the different manufacturers (i.e., discounts in M1 prices vs. tenders) does not influence the prices significantly either.It can be speculated upon whether the discounts omitted by M1 were a part of a market strategy to be competitive in the market, a strategy in which the calculated pipe costs in this study would seem to confirm as effective.Nevertheless, aside from speculation, the small pipe cost differences seem to reflect a healthy market with sufficient competition between manufacturers.On account of this, it is deemed improbable that the above-mentioned issues are important to the results.

Cost Level in Calculation Software Wikells
The employed cost calculation software is designed by collecting costs for vast amounts of entrepreneurial costs from various projects and companies.The cost calculations that are delivered are mostly estimates.According to professional users, every company delivering tenders with the program has its own marginal cost added to the project based on the experience from previous projects.According to these users, the actual cost may differ as much as 20-30% for some projects, which is a significant overhead.The "standard cost" or average entrepreneurial cost is more of a guideline than a real estimate.
The construction costs included in this study are based on the design of real buildings and using real wall constructions, HVAC components, and electrical installations.Access to 2017 cost data from an entrepreneur involved in building one intermediate substation equal to the ones modeled in this study shows that the actual costs of one intermediate substation were 51% more expensive than the estimate from Wikells after adjusting for inflation.This is a significant difference, and it seems that the single most important deviation from the calculated costs is due to the HVAC installation material and installation costs.Adjusting the calculated HVAC costs from Wikells per hour to correspond to the same amount of work hours as in the actual substation increases the total substation costs by 44% if all other Wikells costs are kept.
Such a massive uncertainty in costs unquestionably leads to questionable results from a cost analysis.Nonetheless, the costs of an intermediate substation only apply to the Hybrid concept.As this has already been identified as the most costly concept, the potential additional cost does not change the results.On the other hand, the uncertainties poised by using Wikells do give some room for questioning the apparent difference in LCC for the two concepts, as this may be different from that which is presented.However, even with a 51% increase in the ICC of construction (seeFigure 8), the result on the LCC of the GRUDIS system is insignificant compared to the influence on the LCC of the Hybrid system, even when disregarding the financial cost increase.As such, the apparent difference in LCC between distribution concepts should do anything but decrease.

Fuel Costs
The fuel cost used as the basis for the LCC analysis was taken from costs statistics made available by the Swedish Energy Authority (Energimyndigheten) [39].The cost development of refined and unrefined biomass fuels for the period 2011-2020 shows a falling trend in the period 2011-2017 (except for 2016), before a rise from 2018 to 2020.The average fuel cost in this period is 282 SEK/MWh, ranging from 252-332 SEK/MWh.As such, the year-to-year cost of fuel shows high variability due to various factors affecting availability, which makes it difficult to determine a representative value.However, the transition towards increased use of renewable energy is expected to increase the demand for biomass at least until 2030, which could result in increased costs [42].The cost ultimately chosen was based on a single year (2019) as this was equal to that of 2011 (305 SEK/MWh) before the costs started declining, and it is roughly 8% higher than the average for the historical period.This, therefore, includes some margin for higher costs in the future as the sensitivity analysis investigates the potential impacts of various annual fuel cost increase rates.Notably, the influence of the fuel cost is insignificant for the difference between distribution concepts and, therefore, the results between these.Still, it has an observable impact on the LCC and the relative share of LCC cost components.Nevertheless, seeing that the boiler fuel consumption is rather similar between the concepts, the effect of this on the results is assumed insignificant.

•
A hypothetical district heating system has been modeled in TRNSYS using two distribution concepts; a hybrid comprised of 3rd and 4th generation DH technology and GRUDIS, which may be considered 4th generation DH technology.

•
A range of pipe insulation levels has been simulated for each concept, and the average boiler-supplied energy is similar for both the GRUDIS and Hybrid distribution concepts.

•
Increase/decrease in heat losses due to change in insulation levels primarily affects the boiler supplied energy and, to a minor extent, solar contribution.The reduction in boiler supplied energy is at best 5% for the GRUDIS concept and 7% for the Hybrid concept.

•
The lowest pipe network losses are found for a hybrid system variant using the highest available insulation level for both steel pipes and PEX pipes.The GRUDIS variant with the lowest pipe heat loss uses the PEX pipes with the highest insulation level.

•
The LCC is highest for the hybrid concept, being 20.4 M SEK on average for all variants simulated, compared to 15.5 M SEK for all variants of the GRUDIS concept-a difference of the average of 24%.

•
An annual increase in fuel cost of 2-6% has a negligible influence on differences in LCC between the hybrid and GRUDIS concept, whereas an interest rate of 2-6% has a significant impact on the differences due to the larger initial capital cost associated with the construction and consequently, financing of, the hybrid system.

•
Simplifications made in the modeling approach appear to have an insignificant influence on the results, as these are based on inter-comparison between systems built on the same assumptions.

•
The influence of the cost calculation method appears to have little to no effect on the difference in LCC between the distribution concepts, as different types of cost input for the calculation of pipe network costs give similar results.

Future Work
As this study has focused on the best distribution concept for a new DH system with a pre-determined heat demand and solar array size, a future study should look at determining the optimum solar thermal system integration for a pre-determined distribution concept and heat demand.For DH systems built anew, solar thermal integration should be a requirement.The GRUDIS concept could have other system configurations that would lead to higher solar fractions and hence, better economy.Furthermore, due to the market share of DH already being so high in many parts of the world, there is a large amount of existing 3rd generation DH systems currently running at high temperatures and possibly wholly or in part on fossil fuels.This means installing solar can significantly reduce carbon emissions and/or heating costs, should solar heat replace other fuel.For these systems, the most techno-economic integration of solar thermal needs to be investigated to provide updated guidelines on the present economic conditions for solar-assisted DH.One related aspect of this is the influence of increasing demand on the system as the DH network expands, either due to urbanisation or the integration of low heat-density residential areas into the network.Both how the network is expanded and optimal integration of the solar thermal into the expanded system is of interest, as low-temperature networks could yield higher solar fractions and lower costs.However, more studies are required to have conclusive information on this.

Figure 1 .
Figure 1.Overall process-flowchart showing how the work in this study was executed.

Figure 1 .
Figure 1.Overall process-flowchart showing how the work in this study was executed.

Figure 2 .
Figure 2. System layout-schematics of (a) GRUDIS network and (b) hybrid network.Evacuated tube collectors (ETC) are located on the boiler central (BC), and flat plate collectors (FPC) are located on either intermediate substation (SS) and/or ground [12].

Figure 3 .
Figure 3. GRUDIS distribution-Simple schematic of GRUDIS distribution system with boiler central (BC) and single-family house.Ground-mounted flat-plate collectors (FPC) and evacuated tube collectors (ETC) on the BC [12].

Figure 4 .
Figure 4. Hybrid distribution-Simple schematic of the hybrid heat distribution concept with boiler central (BC), intermediate substation, and single-family house.Roof and ground-mounted flat plate collectors (FPC), evacuated tube collectors (ETC) on boiler central [12].

Figure 4 .
Figure 4. Hybrid distribution-Simple schematic of the hybrid heat distribution concept with boiler central (BC), intermediate substation, and single-family house.Roof and ground-mounted flat plate collectors (FPC), evacuated tube collectors (ETC) on boiler central [12].

Figure 6 .
Figure 6.Simplified approach-Schematic showing the hybrid system simulation model and applied pipe sizes.Copper pipes in the solar thermal system are not shown.Scaling factors (2x and 25x) are shown in the substation models where these are used [12].BC = Boiler central.SS = substation.

Figure 6 .
Figure 6.Simplified approach-Schematic showing the hybrid system simulation model and applied pipe sizes.Copper pipes in the solar thermal system are not shown.Scaling factors (2x and 25x) are shown in the substation models where these are used [12].BC = Boiler central.SS = substation.

Figure 7
Figure 7 shows the simulated annual energy balance for the two distribution concepts modeled in this study with the reference pipe combinations introduced in Section 3.1.4.The legend used refers to the KPIs listed in Section 3.2.

Figure 7
Figure 7 shows the simulated annual energy balance for the two distribution concepts modeled in this study with the reference pipe combinations introduced in Section 3.1.4.The legend used refers to the KPIs listed in Section 3.2.

Figure 7 .
Figure 7. Energy balance; simulation results for the reference pipe combination/insulation class of the two distribution concepts modeled in this study, showing energy input is positive and output is negative.

Figure 7 .
Figure 7. Energy balance; simulation results for the reference pipe combination/insulation class of the two distribution concepts modeled in this study, showing energy input is positive and output is negative.

Figure 8 .
Figure 8. LCC analysis-Overview of the LCC of the two reference variants for the two distribution concepts simulated in this study with an economic lifetime of 20 years.

Figure 8 .
Figure 8. LCC analysis-Overview of the LCC of the two reference variants for the two distribution concepts simulated in this study with an economic lifetime of 20 years.

Figure 9 .
Figure 9. Energy balance difference-Relative changes in lost and supplied energy compared to the reference variant (G1) for the variants of the GRUDIS concept.The bold line indicates a change in pipe manufacturer; G1 and G2 employ M2, while G3 and G4 employ M1 and M3, respectively-see Table7for overview.

Figure 9 .
Figure 9. Energy balance difference-Relative changes in lost and supplied energy compared to the reference variant (G1) for the variants of the GRUDIS concept.The bold line indicates a change in pipe manufacturer; G1 and G2 employ M2, while G3 and G4 employ M1 and M3, respectively-see Table7for overview.

Figure 10 .
Figure 10.Energy balance difference-absolute changes in supplied and lost energy compared to the reference variant (H1) for the variants of the hybrid distribution concept.The bold line indicates a change in the manufacturer of PEX pipes; H1-H4 employs PEX M2, while H5-H8 employs M1 and M3-See Table6for an overview.

Figure 10 .
Figure 10.Energy balance difference-absolute changes in supplied and lost energy compared to the reference variant (H1) for the variants of the hybrid distribution concept.The bold line indicates a change in the manufacturer of PEX pipes; H1-H4 employs PEX M2, while H5-H8 employs M1 and M3-See Table6for an overview.

Figure 11 .
Figure 11.Cost difference-Absolute cost changes for the variants of the GRUDIS and hybrid distribution concept relative to the reference variant (G1 and H1, respectively).The bold line separates variants of different distribution concepts, while the stifled line separates reference pipe manufacturer(s); G2 and H2-H4 employ PEX M2, while G3 and G4 and H5-H8 employ M1 and M3-See Table6for an overview.

Figure 11 .
Figure 11.Cost difference-Absolute cost changes for the variants of the GRUDIS and hybrid distribution concept relative to the reference variant (G1 and H1, respectively).The bold line separates variants of different distribution concepts, while the stifled line separates reference pipe manufacturer(s); G2 and H2-H4 employ PEX M2, while G3 and G4 and H5-H8 employ M1 and M3-See Table6for an overview.

Figure 12 .
Figure 12.Boundary conditions-Cost difference [M SEK] between reference concepts G1 and H1, as well as relative increase [%] in total costs for variants G1, G4, H1, and H4, with changing fuel cost increase (FCI) and interest rate (IR).Relative values are calculated using FCI 2%/IR 2% as reference.

Figure 12 .
Figure 12.Boundary conditions-Cost difference [M SEK] between reference concepts G1 and H1, as well as relative increase [%] in total costs for variants G1, G4, H1, and H4, with changing fuel cost increase (FCI) and interest rate (IR).Relative values are calculated using FCI 2%/IR 2% as reference.

Table 1 .
Heat supply-Overview of the heat supply configuration for the two distribution concepts investigated in this study.BC = boiler central, SS = Substation, GM = ground-mounted.

Table 3 .
Main TRNSYS components and parameter settings [31].Non-standard components are referenced individually, TESS components are described in

Table 4 .
Reference pipe type and insulation class for the GRUDIS and hybrid distribution concept.

Table 6 .
Hybrid variants-Overview of different pipe combinations used to simulate variations (H1-H8) of the hybrid distribution concept.The pipe types are listed according to their application in either the primary or the secondary network.Abbreviations: ins = insulation and ref. = reference.

Table 7 .
GRUDIS variants-Overview of different pipe combinations used to simulate variations (G1-G4) of the GRUDIS distribution concept.Abbreviations: ins = insulation and ref. = reference.

Table 8 .
LCC types-Overview of cost types included in the LCC for the hybrid and GRUDIS distribution concepts.The LCC cost category is included in parenthesis under the cost type.

Table 9 Table 9 .
KPI-overview over key performance indicators used in this study.

Table 9 .
KPI -overview over key performance indicators used in this study.

Table 10 .
Heat loss-Overview of the total distribution of heat loss