Start Submission

Reading: District cooling network optimization with redundancy constraints in Singapore


A- A+
Alt. Display

Technical Articles

District cooling network optimization with redundancy constraints in Singapore


Johannes Dorfner ,

Technical University of Munich, Arcisstraße 21, 80333 München, Germany, DE
X close

Patrick Krystallas,

SWM Infrastruktur, Emmy-Noether-Straße 2, 80992 München, Germany, DE
X close

Magdalena Durst,

Stadtwerke München, Emmy-Noether-Straße 2, 80992 München, Germany, DE
X close

Tobias Massier

TUM CREATE, 1 CREATE Way, #10-02 CREATE Tower, 138602 Singapore, Singapore, SG
X close


This study presents a mathematical optimization model for planning topology and capacity of a district cooling network. The model relies on a mixed-integer linear programming formulation to find the most economic network layout while satisfying redundancy criteria against unavailability of cooling stations. The simplicity of the formulation makes it easy to embed in other models or to extend it to other redundancy cases.

The model is applied to a case study in the central business district of Singapore. Results show that district cooling is a profitable option for Singapore, especially due to its constant high cooling demand that is currently satisfied mainly through decentralized cooling units. This result generalizes to tropical cities world-wide with high cooling demand density.
How to Cite: Dorfner, J., Krystallas, P., Durst, M. and Massier, T., 2017. District cooling network optimization with redundancy constraints in Singapore. Future Cities and Environment, 3, p.1. DOI:
  Published on 03 Jan 2017
 Accepted on 7 Dec 2016            Submitted on 2 Sep 2016


District cooling has been established in many countries. Especially in bigger cities, it offers an energy efficient alternative to individual generation of cooling power at the site of customers. Even in a city like Munich, Germany, which is located in a temperate climate zone, several district cooling systems have already been installed successfully and maintained.

In contrast to Munich, in Singapore, warm and humid weather prevails all year long, which results in much higher demand for cooling and dehumidification both in residential and non-residential buildings. However, in most buildings, cooling power is generated by using small chillers that use air for re-cooling which is very inefficient. As a result, much electricity is used to operate the chillers. In this case study, we model a district cooling network for the central business district (CBD) of Singapore. We present several cases and a comparison to conventional individual cooling.

Singapore’s CBD comprises mainly commercial buildings, e.g. shopping malls, hotels and office blocks, often with more than 50 floors. Thus, it stands in contrast to other areas of Singapore which comprise more residential or industrial buildings. Due to the different type of use of commercial buildings with the air conditioning running all day long and amounting to approximately 50% of each building’s energy demand, the introduction of district cooling could lead to tremendous energy savings.

Singapore already has a large underground district cooling system in the Marina Bay. The system contains two plants with a maximum of 330MW of cooling power. [1] According to the Singapore District Cooling Pte Ltd (SDC), a subsidiary of Singapore Power, energy savings amount to 40% which corresponds to the energy demand of 24 000 three-room flats operated by Singapore’s Housing Development Board (HDB).

State of projects

Over the last decades, district cooling networks have been developed in several cities throughout Europe. Researches result in data shown in Table 1. Paris for example already started to supply in 1991 [2, 3]. Vienna has set the aim of installing 200MW district cooling power until the year 2020 ([4], p. 2). Stockholm has probably the biggest district cooling network today with a connected load of already 330MW, ranging from small (8kW) to large (7000kW) customers with an average load of about 500kW each [5]. Munich’s district cooling network [6] is constantly growing and they have installed different kinds of cooling in order to supply the best technology for each purpose.

Table 1

State of projects in European cities

City Capacity Energy
[Sources] (MW) (GWh/a)
Berlin [27] 44
Helsinki [28, 29] 135 80
Munich downtown [6] 15
Munich groundwater [6] 11.4
Munich M-Campus (mixed form) [6] 4.39 4.4
Paris [2] 290 420
Stockholm [5] 330 460
Vienna [30] 113.3

State of research

This is a short summary of existing publications on district cooling in Singapore, district cooling in general, as well as similar or prominent approaches for district heating or more general multi-commodity district energy systems.

Cooling in Singapore

Apart from the aforementioned existing district cooling system [1], a recent study on Singapore [7] focused on comparing surface and subsurface cooling in tropical countries using an EnergyPlus model of a representative building. It concludes that an open-loop groundwater cooling system might be a better option for tropical cities like Singapore compared to conventional systems. The viability of district cooling in the CBD area of Singapore has already been assessed in a working paper [8] by Singapore Power Ltd.

District cooling

A thermal ice storage system and its influence on the cost-optimal design of district cooling distribution network is investigated [9] using non-linear optimization in two exemplary networks. It finds that ice storage can help reduce the required pipe diameters. A TRNSYS model is built to investigate a similar question [10] for a case study in Hong Kong; storage is not found to be an economically attractive option there.

A genetic algorithm is used [11] to assess the optimum shares of five building types to form a cooling load time series that is most suitable for district cooling. This approach was later embedded in an optimization scheme for a development project in South East Kowloon, Hong Kong. Another genetic algorithm was used by Feng and Long [12] to determine and compare piping layouts.

Khir [13] presents a mixed integer linear programming (MILP) model for designing district cooling system. It explicitly models flow rates, temperatures and pressures in a graph of arcs and nodes. Linearization of the physical constraints leads to a large amount of auxiliary variables. Both plant design and operation and network design can be optimized. However, to achieve a tractable problem size, those two subproblems are solved separately. However, this separation requires fixing the supply task. The paper also highlights the relationship with Steiner tree problems.

Södermann [14] presents a similar MILP model for designing district cooling systems in urban areas. It uses a linearized cost function for sizing equipment (plants, pipes, storages). A small set of representative load periods (8 in the example study) approximate the annual cooling load duration curve.

District heating

Due to the basic principle, district cooling and heating are comparable types of supply. In both types, with few central generation plants heat or cold is generated which is mostly transferred to water. This tempered water is transported through pipes to customers buildings. Having transferred the heat or cold to the buildings, the cooled or heated water returns to the generation plants where the process starts from the beginning. The biggest difference between both are the technologies used in generation plants, which is not the focus of this article.

Jamsek [15] presents a linear programming (LP) approach for determining cost-minimal pipe topology and capacity configuration for a given structure and cost data. A more conceptual, but worthwhile discussion of district heating and cooling systems is given by Rezaie [16].

Udomsri [17] presents and compares the different options to combine centralized and decentralized heating and cooling. It points towards solutions that use decentralized, thermally driven chillers to produce cooling on-site. This configuration can have advantages in regions with both heating and cooling demands throughout the year.

Multi-commodity district energy system

A detailed model on combined energy system planning for electric, heating and cooling demand is presented by Chinese from 2008 [18]. It presents a MILP model. A case study in Udine, Italy, is performed with 7 nodes and 41 individually weighted time steps to represent the seasonal demand characteristics of a full year.

A recent paper [19] investigates the case of combined district energy system for electricity, heating and cooling using a mixed integer linear optimization problem. Both spatial (7 nodes) and temporal (2 days in 48 time steps) resolution are low in the presented case study, but the non-linear performance of gas turbines is represented by a set of linearized constraints.

A broader review of computer models for renewable energy system analysis and optimization is provided by Connolly et al. [20]. Weber presents in her thesis [21] a comprehensive methodology to plan cost-optimal polygeneration energy systems.

This work

This paper presents the generalized version of a previously published model developed for planning of district heating distribution networks [22]. That model has then been adapted for optimization of district cooling systems and applied in a small case study for the city center of Munich [23].


This section describes a mathematical optimization model to represent the planning task for a minimum cost district cooling network. Its main input is a connected graph of street segments that represent the discrete possible network parts. Each segment – called edge in the following – has the attributes length and a cooling demand that may (but does not have to) be satisfied by the district cooling network. A subset of the graph’s vertices may provide cooling at fixed specific costs, representing (possible) cooling stations.

The model’s main contribution lies in a simple representation of redundancy constraints by means of artificial time steps (or periods) with pre-determined reduced availability of cooling stations. This technique can easily be embedded into other optimization models and can increase the robustness of obtained solutions against equipment downtime due to outages or maintenance.


A brief conceptual overview on the model’s inputs (parameters) and outputs (results) is given in Fig. 1. The presented optimization model, based on estimated (or known) cooling demand and a graph of possible segments to construct a cooling network, derives the cost-optimal size and topology of such a network. The objective function minimizes total costs (generation costs minus revenues), thus yielding the most profitable system design.

Fig. 1 

Optimization model overview. Inputs and outputs of the presented mathematical optimization model

To improve performance for large study regions, the model employs a linearized cost function for sizing the thermal pipe capacity. Depending on the parameterization, one can either cautiously overestimate the real estimated costs or try to fit as closely the observed cost function. In the latter case, one typically slightly underestimates the cost of medium capacities, while overestimating the cost of small and large capacities.

Mathematical description

The presented model minimizes the total costs ζ of generating and distributing cooling through a district cooling network. The key result is the size and topology of the network, represented by the value of network capacity (vector of all individual arc capacities ) in the optimal case:



A district or city is represented as a graph of vertices and arcs. This graph is derived from the street network. It should be derived in such a way to include all considerable locations for network pipes. The spatial resolution can be chosen as fine as required. Here, it is on the level of building blocks, i.e. a single street segment between two intersections.

Let V be the set of vertices vi, corresponding to connecting or terminating points of the graph. Set A of arcs then comprises ordered tuples of vertices aij=(vi,vj) with ij. A is symmetric, that means either both or none of the pair aij and aji are elements of A. For readability, the subscript ∘i is used to denote any parameter or variable that is defined over vertices viV, while ∘ij is used to denote a quantity defined over arcs aijA.

The set V0V defines so-called source vertices, which represent locations of possible cooling sources. Usually, the number of source vertices is low (≤10) compared to the size of the graph.

The neighbor sets Ni of a vertex vi are defined by the indices of all vertices that are connected to it by an arc: Ni={ kakiA }.

Time is represented by a set T of discrete time steps t, which represent a small number of representative operational situations. The minimum viable number of time steps is two: one for peak demand with a duration of one to several hours, a second for the average annual load with a duration of the remaining year. A more comprehensive choice are three to four time steps: The third can refer to the all-year base load, while the fourth can be used to represent a common intermediate load level. Additional time steps need to be introduced if redundancy requirements are to be stated. Refer to the discussion of the parameter availability below for more details.


The numerical given facts for this model are shown in Table 2. They are grouped by their defining domains. First are vertex parameters, then arc parameters, then constant or global parameters, and finally time-dependent parameters.

Table 2

Optimization model parameters

Name Unit Description
cfix S$/m Fixed investment costs
cvar S$/(kW m) Variable investment costs
com S$/(m a) Operation & maintenance costs
crev S$/kWh Revenue for cooling
wfix kW/m Fixed thermal losses
wvar kW/(kW m) Variable thermal losses
b Concurrence effect
q Connect quota
u 1/a Annuity factor for investment costs
S$/kWh Cooling costs at source vertices
kW Source vertex capacity
lij m Arc length
dij kW Arc peak demand
gij Existence of a pipe (1=yes, 0=no)
kW Maximum pipe capacity
st 1 scaling factor
wt h weight/duration
yit availability (1=yes, 0=no)

Vertices have a single parameter: their maximum power capacity . It is the thermal output power given in kW for that location. All non-source vertices have this parameter set to 0.

Arcs have three main attributes: their length lij (m), the thermal peak demand dij (kW) of their adjacent buildings and gij (binary), which states whether a pipe already exists in this arc.

As secondary parameter, indicates the maximum thermal power capacity (kW) in an arc, which can be derived from the maximum available pipe diameter. For arcs with existing pipes (gij=1), this value should be set according to the diameter of that existing pipe. As both arcs aij and aji refer to the same street segment, parameter values for both members are always identical.

Global parameters are technical and economic parameters that refer to the district cooling system as a whole.

Economic parameters are all costs and revenues. The investment costs are split into a fixed part cfix and a variable part cvar. The fixed part, given in S$/m, captures all costs that are not dependent on the capacity of the pipe to be built, mainly earth works. The variable part, given in S$/(kW m), refers to the capacity-dependent cost component of building a pipe. Both values must be tailored to the study area to compensate for differences in cost structure and availability of different pipe sizes. The same is true for the operation & maintenance (OM) cost parameter com in S$/m. In contrast to the fixed investment cost term, it is also to be paid for existing pipes. The cost for providing cooling here is dependent on the vertex to allow for representing cheap surface water cooling in contrast to more expensive re-cooling options.

The letter c and a subscript denote economic parameters. These are investment costs for building the pipe network cfix and cvar, maintenance com, costs of providing cooling, and revenue for delivering crev cooling to consumers. While cfix contains the size-independent costs (mainly earth works), cvar contains costs that are dependent on the thermal capacity (diameter) of the pipe.

Time-dependent parameters st and wt represent scaling factor (dimensionless) and weight or duration (h) of a time step t. A value of st=1 refers to peak demand, while smaller values correspond to moments of partial load. Together, these two parameters encode a discretized annual load curve.

Redundancy requirements can be stated in the model by setting the binary availability parameter yit to value 0 in certain time steps for 1,2,… source vertices. If one time step for each foreseen failure configuration is introduced, a full n−1,n−2,… failure safety (against cooling source outages) can be enforced in the produced solution. Clarification: in the conventional time steps discussed above in paragraph time step set, the value yit=1 is to be set for all source vertices viV0.


The main optimization task is to find values for the binary decision variable ξij. If its value is one, a pipe is built in the street segment corresponding to the arc aij. For each time step, the actual use of a given pipe is decided by setting the binary pipe usage decision variable ψijt. If its value is one, the pipe in arc aij is used in direction from vertex vi to vertex vj. Consequently, there must be a power flow into the pipe. Simultaneously, a value ξij=1 requires that the demand dij of this arc has to be satisfied at all times. The power flow variable at the other end of the pipe is called . Its value is determined by the ingoing power flow minus losses and demand.

The non-negative variable ρit represents thermal power output from a source vertex viV0 at time t in kW. It is limited in value by the source vertex capacity parameter and possibly by the availability parameter yit. Table 3 summarizes all model variables. The following line shows their domains, on which they are defined:

Table 3

Optimization model variables

Name Unit Description
ζ S$ Total system cost (Inv, O&M, Rev, Gen)
ξij Binary decision variable: 1 = build pipe
ψijt Binary decision variable: 1 = use pipe
kW Thermal power flow capacity into arc aij
kW Thermal power flow from vi into arc aij
kW Thermal power flow out of arc aij into vj
ρit kW Cooling in source vertex vi


Equations fall into two categories: The first is the cost function whose value is to be minimized. The second are constraints that codify all the previously discussed rules in mathematical form. By that, they define the region of feasible solutions, under which the solver selects a (close to) cost optimal solution.

The cost function value ζ is the sum of costs for providing cooling, minus revenue for that cooling. Network costs occur for building the network (annualized investment) and for operation & maintenance. The following derived parameters are introduced to shorten the definition of the cost function whose definition is given below: (S$/a) and (S$/(kW a)) cover network costs (investment, O&M). (S$/kWh) represents cooling generation costs, while (S$/h) represents revenue for provided cooling. With these, the cost function is equal to


The first summand forms a piece-wise linear function that is depicted in Fig. 2. The plot also depicts that the linear approximation possibly over- and underestimates the cost of either small or large capacities. Depending on the parameterization, one can either force an optimistic (by deliberate underestimation) or cautious (by deliberate overestimation) planning decision. The linear approximation exhibits unsteady jump from zero to the fixed cost ucfix+com for any positive value of for any non-existing (gij=0) pipe. The link between the binary decision variable ξij and the continuous pipe capacity variable is enforced through Eqs. (8), (12), and (14) below.

Fig. 2 

MILP pipe investment cost function. Investment and operation & maintenance costs for a single arc aij over the pipe capacity (thick, blue). The curved line (thin, green) indicates a representative possible real non-linear cost function

These are the definitions for the four derived parameters used in the cost equation above. The division by value 2, common to all three arc parameters, compensates for the double-representation of one street segment {vi,vj} as a pair two directed arcs aij, aji. More about this issue can be found at the explanation of the symmetry constraints (15) and (16) below. The cooling cost parameter is divided by the concurrence effect parameter b to remove the load reduction effect that is introduced in Eq. 10. In other words, the power flow that needs to be satisfied from a source vertex is lower by b than the energy flow that needs to be delivered. The division by b removes that difference in terms of costs.





The constraints formalize all physical laws and technical rules that need to be satisfied by a feasible solution. The first constraint concerns energy conservation in vertices vi, with respect to its neighbors Ni:


This inequality is thus a relaxed version of the law of energy conservation. Relaxed, because it allows for power to vanish at any vertex. As power does not come for free, any solution returned by the solver will satisfy this constraint to equality. For all non-source vertices (i.e. ρit=0), the difference between outgoing and incoming power must be smaller than or equal to zero. In source vertices, a positive difference may remain when it is met by an equal amount of input power ρit into the network at that vertex.

Demand satisfaction is the main arc constraint. It creates the logical connection between the discrete pipe usage decision ψijt and the demand satisfaction. Like the cost equation, it relies on two derived parameters xij and yijt. This constraint is also graphically explained in Fig. 3.

Fig. 3 

Sankey diagram of Eq. (8). Ingoing cooling power flow from vertex vi is reduced by variable losses (1−xij), fixed losses and cooling demand (yijt) before leaving to vertex vj as power flow

Parameter xij refers to losses wvar that are proportional to the amount of power flow into the arc aij. Parameter yijt refers to fixed thermal losses wfix that only depend on the pipe length, not on the power flow, and the time-dependent demand dij·st that is further reduced by the concurrence effect b and the connect quota parameter q. In notation:



Pipe capacity is a technical constraint that limits the power flow through a pipe by the built capacity. This is accomplished by first limiting the ingoing power flow into a pipe by the built pipe capacity . The second constraint is required to set the pipe usage decision variable ψijt. Its value is forced to 1 if a flow in direction from vi to vj at time step t is required. The reason for this seemingly artificial constraint is explained in the following constraint.



Unidirectionality of the power flow is required, as otherwise the solver would happily use a pipe’s capacity in both directions simultaneously, which would not be physically possible. Therefore, the pipe usage decision variable ψijt may only have the value 1 in one direction at any given time step t. This way, the direction of flow may change between time steps, but any given pipe may only be used in one direction at a given time.


Build capacity limits the pipe capacity variables to the maximum available diameter, whose capacity in kW is represented by parameter . At the same time, this constraint sets the value of the building decision variable ξij, if a non-zero pipe capacity is needed.


Symmetry of building decision are two constraints that enforce that the mathematically independent arcs aij and aji must have identical variable values for their pipe capacity and building decision ξij. This way, a pipe can then be used in both directions at different time steps.



Use if built is the last piece in the puzzle to link the building decision ξij to the pipe usage decision ψijt. Up to this point, nothing in the formulation requires that the cooling demand of a customer is satisfied at all times. This constraint enforces exactly that, by requiring that the sum of ψijt+ψjit is greater or equal to 1, if a pipe is built along that arc. Otherwise, the condition should not be enforced. This is done by calculating the expression (ξij+ξji)/2. It is equal to 1 if ξij=1, as Eq. (15) requires ξij=ξji. Otherwise, its value is equal to 0. In other words: if an arc is has a pipe, its demand must always be satisfied by a power flow from any of its two sides.


Source vertices is the last constraint. It limits the source vertex power flow ρit by its maximum allowed capacity . The availability parameter yitusually has value 1, except for special time steps in which the source vertex vi is made unavailable by setting yit=0:


Equations 2 to (18) together define a linear mixed-integer program, whose optimal solution is the desired vector of pipe capacities , accompanied by optimal cooling power flows.

Case study

The presented model is applied on a potential supply area in the central business district of Singapore, as shown in Fig. 4. The first part describes how the input parameters were derived. Then the assumptions for a sensitivity analysis against conventional per-building cooling are described.

Fig. 4 

Study region within Singapore. The green outline shows the location of the CBD study region within the extent of Singapore. It covers the Marina Bay and follows Orchard Road to the north west

Load profile estimation

Due to lack of access to a load profile from Singapore, a known cooling load from Munich was extrapolated based on weather data in Singapore.

The district cooling system of Munich has mostly non-residential customers like shopping malls, offices, congress centers and smaller retail trades. Its usage type composition is therefore rather similar to Singapore’s CBD. However, the buildings are much smaller and thus have a different surface to volume ratio and surface materials than the much higher buildings in Singapore’s CBD.

To record and to utilize the optimization potential of this expanding system, in a detail data analysis a question amongst others was examined on which factors the customers cooling demand is depending. Furthermore the collected data serve as basis for the evaluation of the cooling demand in other cities like Singapore. Based on a period of two years up to 2015 the influence of temperature, humidity, insolation and the customer behavior itself was analyzed.

The strongest influence on the cooling demand is the temperature. The humidity plays a subordinated role in Munich (but not in Singapore), cause of the temperate climate zone and the previous low dehumidification in the comfort zones. To count in this factor to the further customer and grid extension, the enthalpy of the outside air, as an index for the heat input to the buildings, is used as leading control parameter. The insolation has no considerable influence on the customers’ requirements because of the good insulation of the rather new buildings and the shading situation in the densely built city center.

The second strongest influence factor on the cooling demand of the Munich costumers are the business hours. Due to the simultaneous public traffic in the shopping areas and offices as well as the highly volatile temperature profile of a day, the cooling demand is five times higher during the day than at night or on a holiday. This characteristic day-night cycle are also recorded in other Central European cities ([24], pp. 8, 80).

Using the high cooling demand during the business hours and the enthalpy dependence, a trend was generated for the customers’ requirements. This trend line is combined with the climate records of Singapore to create an approximated annual load duration curve, shown in Fig. 5.

Fig. 5 

Annual load duration curve. Estimated duration curve (thin, green) and its discretization to 8 steps (thick, blue) as used by the optimization model. Step lengths correspond to time step weight wt and the relative load value to st

To reduce the temporal resolution for the optimization model to manageable sizes, this curve is discretized to 8 individually weighted time steps, also shown. One 24 h long time step represents peak demand, while the other 7 steps have durations ranging from 390 h to 2006 h.

Annual cooling demand

Cooling demand is estimated based on previous work conducted by Böhme et al. [25]. There, a geographic dataset of buildings with gross floor area and usage type was prepared. The total floor area by building type (e.g. residential, commercial, …) is summarized in Table 4. It also states the assumed design cooling loads, based on a mean design load of 112.5kW/m2, estimated from PengKee [8]. Building or area types with small total area were ignored. With those assumptions, a total summed peak cooling demand of 1241MW exists in the case study area.

Table 4

Total floor area and cooling peak load by building type

Building type Peak load Floor area
(W/m2) (103 m2)
Civic & community institution 112.5 529
Commercial 112.5 6685
Commercial & residential 100 1021
Educational institution 112.5 175
Health & medical care 125 40
Hotel 112.5 1605
Open space 0 22
Park 0 58
Place of worship 0 137
Reserve site 0 157
Residential 112.5 717
Residential/commercial (1st storey) 112.5 357
Sports & recreation 75 15
Transport facilities 0 29
Utility 0 43

Input data

The street graph is derived from OpenStreetMap data [26], processed using previously described [22] pre-processing steps using polygon skeletonization for simplifying the dense street network.

Cost data (given in S$) and technical parameters for this case study are summarized in Table 5. Pipe investment costs of 7000 S$/m are the main cost driver on the supply side. This value is scaled to one year with an annuity factor


Table 5

Parameter values in case study

Name Value Unit
cfix 7000 S$/m
cvar 8e-4 S$/(m kW)
com 80 S$/(m a)
crev 0.14 S$/kWh
wfix 0.01 kW/m
wvar 1e-8 kW/(kW m)
b 0.9
q 1.0
u 0.091 1/a
0.03 or 0.07 S$/kWh

with interest rate i=6.5% and depreciation duration n=20a, yielding a value of 0.091. Operation & maintenance are relatively cheap with annual 80 S$/m. Cooling cost depends on the location of the cooling station. The large triangles in Fig. 6 show stations with access to surface water, which allows cheap re-cooling with an assumed cost of 0.03 S$/kWh. The other six cooling stations have higher costs of 0.07 S$/kWh.

Fig. 6 

Demand and cooling stations. Location of cooling demand (color and line thickness of edges) and cooling stations (triangles). Large triangles (at vertices 57, 110, 154) correspond to cheap surface water cooling, while small triangles assume more expensive air re-cooling


Base case

The base case is designed to show today’s demand situation. In this scenario, the size and layout of a profitable district cooling network is to be determined.

In the base case, a revenue of 0.14 S$/kWh is assumed. A sensitivity analysis with reduced revenues in steps of 0.02 S$/kWh is also conducted to determine how sensitive the optimal network size is on the price for cooling.

Growing demand

In order to assess how the total cooling costs (generation and distribution) are affected by a possible demand growth in the study area, additional load is introduced in the edge (54, 114) in the south-eastern corner. It’s value is changed from 0MW to 200MW with steps of 50 MW. As this step creates more load than the existing cooling stations could satisfy while satisfying the redundancy constraint, all cooling stations’ capacity is increased by 50%.


Base case

The resulting network for the base case is shown in Fig. 7. It shows a district cooling network that satisfies almost all (1240MW of 1241MW) cooling demand. In other words: under the presented cost assumptions and a cooling revenue of 0.14 S$/kWh, district cooling is an economically viable strategy. A sensitivity check with a reduced cooling revenue of 0.10 S$/kWh still leads to an almost full saturation of 1237 MW supplied demand. The unsatisfied demands are mainly short, isolated edges with low demand (≤3MW) that cannot offset the network investment cost. Only when revenue is further reduced, between 0.08 and 0.06 S$/kWh, the optimal network drops from satisfying 1200 MW to 0 MW. This value range thus can be taken as a lower bound for profitability of a district cooling system (under the given cost assumptions).

Fig. 7 

Cost-optimal network layout. Line thickness corresponds to thermal capacity (kW) of each pipe

Figures 8 and 9 compare the cooling power flow in two of the 18 time steps considered for designing the network capacities in Fig. 7. The first situation shows the cost-optimal cooling power flow when all cooling stations are operational. The three cheap cooling stations (57, 110, 154) are delivering their full capacity of 150 MW each, while several of the other stations only run in partial load. The second situation in Fig. 9 shows how the redundancy requirement leads to the consideration of very different flow situations. Here, cooling station 13 now runs at full capacity and satisfies the middle of the study area.

Fig. 8 

Cost-optimal power flow during peak load. All three surface water cooling station operate at maximum capacity (150 MW each)

Fig. 9 

Power flow during outage of cooling station at vertex 57. In this situation, the large pipe capacities south of vertex 13 are needed to transmit backup cooling power from the western three cooling stations to the central bay area

Growing demand

The resulting cooling cost for gradually increasing the demand in edge (54, 114) for each case is shown in Fig. 10. The cooling cost here is calculated by dividing total costs for generation and network by the amount of provided cooling (in kWh). Cost for cooling grows slightly – along with demand – from approximately 0.045 S$/kWh to 0.048 S$/kWh in the case of additional 200MW peak cooling demand. This increase is caused by a higher utilization of the more expensive 0.07 S$/kWh cooling stations. The specific network costs actually exhibit a minor reduction as demand grows, but this effect is negligibly small here.

Fig. 10 

Specific cooling costs for growing demand with original 0.03 and 0.07 S$/kWh cooling stations. Cost increase is caused by higher utilization of expensive cooling stations

Discussion and future research

District cooling can be an attractive option for cities with constant high cooling load, compared to less efficient distributed cooling. The high demand density of Singapore’s CBD makes it a very well suited candidate. Direct access to surface and seawater provides enough cooling potential to satisfy a significant fraction of the present cooling loads.

As Singapore is planning to establish further dense commercial areas similar to the CBD with many office blocks, e.g. in Jurong, more regions will become suitable candidates for district cooling. Moreover, the potential of district cooling could also be evaluated for densely populated residential areas, e.g. Bukit Panjang or Punggol.

The redundancy requirements on cooling stations could need refinement, possibly by relaxing the full (n-1) capability for the whole system to only certain areas. An extension could introduce the same requirement not only to cooling station, but also to crucial parts of the network.


We express our gratitude to the Singapore Land Authority for supporting us with geospatial data.

This work was financially supported by the Singapore National Research Foundation under its Campus for Research Excellence And Technological Enterprise (CREATE) programme. This work was supported by the German Research Foundation (DFG) and the Technical University of Munich (TUM) in the framework of the Open Access Publishing Program.

Availability of data and materials

The optimization model’s implementation is available under the GPL 3.0 license at The model input data is available under the CC-BY license via doi:10.6084/m9.figshare.3582456.

Authors’ contributions

JD developed and applied the optimization model to the case study. PK lead his expertise in planning district cooling networks to provide input data and scenario definitions. MD contributed the method for deriving the annual load duration curve. TM provided data for the case study region and vetted model results against local experience. All authors collaborated for writing the article. All authors have read and approved the final submitted manuscript.

Competing interests

The authors declare that they have no competing interests.



    1. Singapore District Cooling PTE Ltd (2010) District Cooling as a Utility Service At Marina Bay Business District. Accessed 28 Jun 2016. 

    2. Terouanne D (2008) Climespace Cooling Paris, Naturally. Accessed 07 Sep 2015. 

    3. Inpal Energie (2010) 70 Km Fernkältenetz in Paris. Accessed 26 Jun 2016. 

    4. Kontrollamt der Stadt Wien (2012) Fernwärme Wien Gesellschaft m.b.H., Prüfung des Geschäftsfeldes Fernkälte. Tätigkeitsbericht KA IV - GU 45-4/12, Vienna Accessed 08 Sep 2015. 

    5. Hill A, Keilaniemi F (2012) District Heating an Cooling in Stockholm: Site Visit to the Ropsten Plant. Acessed 08 Sep 2015. 

    6. Stadtwerke München (2015) M-Fernkälte. Accessed 13 Mar 2016. 

    7. Liu, Y, Qin, XS and Chiew, YM Investigation on potential applicability of subsurface cooling in SingaporeAppl Energy2013103197206 

    8. Kee TP (2010) District cooling as an energy and economically efficient urban utility – its implementation at marina bay business district in singapore White paper, Singapore District Cooling Pte Ltd.Singapore. 

    9. Augusto, GL, Culaba, AB and Maglaya, AB Identification of design criteria for district cooling distribution network with ice thermal energy storage systemEnergy Procedia201579233238 

    10. Chan, ALS, Chow, TT, Fong, SKF and Lin, JZ Performance evaluation of district cooling plant with ice storageEnergy2006311427502762 

    11. Chow, TT, Chan, ALS and Song, CL Building-mix optimization in district cooling system implementationAppl Energy2004771113 

    12. Feng X, Long W (2008) doi:10.1109/icnc.2008.196 In: Fourth International Conference on Natural Computation.. Institute of Electrical & Electronics Engineers (IEEE), Jinan. 

    13. Khir, R and Haouari, M Optimization models for a single-plant district cooling systemEur J Oper Res201524726486583375850 

    14. Söderman, J Optimisation of structure and operation of district cooling networks in urban regionsAppl Therm Eng2007271626652676 

    15. Jamsek, M, Dobersek, D, Goricanec, D and Krope, J Determination of optimal district heating pipe network configurationWSEAS Trans Fluid Mech20105165174 

    16. Rezaie, B and Rosen, MA District heating and cooling: Review of technology and potential enhancementsAppl Energy201293210 

    17. Udomsri, S, Bales, C, Martin, AR and Martin, V Decentralized cooling in district heating network: System simulation and parametric studyAppl Energy201292175184 

    18. Chinese, D Optimal size and layout planning for district heating and cooling networks with distributed generation optionsInt J Energy Sector Manag200823385419 

    19. Ameri, M and Besharati, Z Optimal design and operation of district heating and cooling networks with cchp systems in a residential complexEnergy Buildings2016110135148 

    20. Connolly, D, Lund, H, Mathiesen, BV and Leahy, M A review of computer tools for analysing the integration of renewable energy into various energy systemsAppl Energy201087410591082 

    21. Weber CI (2008) Multi-objective design and optimization of district energy systems including polygeneration energy conversion technologies. PhD thesis, École Polytechnique Fédérale de Lausanne, Bucharest (Romania). 

    22. Dorfner, J and Hamacher, T Large-scale district heating network optimizationIEEE Trans Smart Grid20145418841891 

    23. Krystallas P, Dorfner J (2014) Expansion planning for district cooling networks using mathematical optimization: 38–42. EuroHeat Power Int 3:. 

    24. Urbaneck, T Kältespeicher: Grundlagen, Technik, Anwendung2012MunichDe Gruyter Oldenbourg 

    25. Boehme, P, Berger, M and Massier, T Estimating the building based energy consumption as an anthropogenic contribution to urban heat islandsSustain Cities Soc201519373384 

    26. OpenStreetMapLicensed Unter the Open Database License. Accessed 3 Aug 2016. 

    27. Vattenfall GmbH (2015) Kälte Für Berlins Mitte. Accessed 04 Sep 2015. 

    28. Riipinen M (2013) IEA Internationale Energy Agency Joint Workshop. Helsinging Energia. Accessed 03 Sep 2015. 

    29. Wallisch A (2012) Fernkälte in Europa. Wien Energie GmbH. Accessed 02 Sep 2015. 

    30. Energie W (2014) Wien Energie Energieversorgung Fernkälte. Accessed 31 Jul 2016.