Determination of the Thermal Insulation for the Model Building Approach and the Global Effects in Turkey

Abstract

One of the most important considerations to be considered in the design of energy efficient buildings is the thickness of the insulation to be applied to the building. In this study the existing building stock in Turkey has been investigated depending on parameters such as the height and the area. A model building has been created covering all of these buildings. Fuel emission reduction of combustion system was calculated in the case of insulation applied to this model building. Heat loss of the existing building stock and exhaust emissions and the contribution to the country's economy with the model building methodology are also determined. The results show that the optimum insulation thicknesses vary between 3.21 and 7.12 cm, the energy savings vary between 9.23 US$/m² and43.95 US$/m², and the payback periods vary between 1 and 8.8 years depending on the regions. As a result of the study when the optimum insulation thickness is applied in the model building, the total energy savings for the country are calculated to be 41.7 billion US$. And also total CO₂ emissions for the country are calculated to be 57.2 billion kg CO₂ per year after insulation.

1. Introduction

Regulations regarding the thermal insulation of buildings utilize the calculation method stated in the TS 825 (Turkish standard) compulsory. A balance must be determined between the insulation thickness, cost, and the savings obtained. It is also changing dramatically due to various data (areas) for all buildings. Thermal insulation in buildings reduces energy consumption and provides fuel saving as well as other benefits such as eliminating condensation and mold formation on interior surfaces of walls. Thermal conductivity and first cost are important parameters in selecting an insulation material. Insulation materials are not independent energy production or conservation systems but are one of the complex structural elements that form a building's shell. In this sense, they cannot be evaluated in the same way that energy producing systems such as solar thermal systems or photovoltaic systems can; they must be evaluated as an integral part of a building's design and construction. Thirty years after the introduction of compulsory thermal insulation in most European countries, insulation materials still comprise the major means for improving a building's energy behavior. The use of insulation materials has increased, both in terms of the number of buildings being insulated and in terms of the minimum amount of insulation required by the national regulations [1]. In this study, the geometries of all of the buildings in each city of Turkey were examined. These buildings were classified for each climatic region, and model buildings were created. For the model buildings, heat losses, insulation thicknesses, and optimum insulation thicknesses were determined. The energy savings that would be achieved from modifications to all of these buildings were calculated. Effects of different insulation thicknesses on fuel consumption and thereby on emissions of pollutants such as CO₂ are evaluated.

When the studies in the literature regarding this issue are analyzed, the fuel costs per unit area are calculated according to degree-day values. Hasan [2] used life cycle cost analysis to determine optimum insulation thicknesses. The results showed that savings up to 21 US$/m² of wall area are possible for rock wool and polystyrene insulation. He determined payback periods of 1–1.7 years for rock wool and 1.3–2.3 years for polystyrene insulation, depending on the type of wall structure. Papadopoulos et al. [3] sorted old buildings by age and investigated the insulation payback period. Çomaklı and Yüksel [4] studied the optimum insulation thickness for the coldest cities in Turkey, and when the optimum insulation thickness was applied, a considerable energy saving of 12.11 US$/m² of wall was obtained in Erzurum. Dombayci et al. [5] reported a study in 2006, in which five different energy sources and two insulation materials were compared when calculating the optimum insulation thickness of wall for Denizli in Turkey. Sisman et al. [6] investigated optimum insulation thicknesses for different degree-day regions of Turkey, namely, Izmir, Bursa, Eskisehir, and Erzurum for a lifetime of N years, by using the present worth value of annual energy savings for insulated external walls. Mahlia et al. [7] found a relationship between the thermal conductivity and optimum thickness of insulation material and gave a correlation for practical use. Özel [8] investigated the optimal location and distribution of insulation in a wall by considering the time lag and decrement factor for various wall orientations in both summer and winter conditions in Elazıg, Turkey. Results showed that the best thermal performance was obtained when three equal insulation layers were used: the first piece placed on the outdoor surface of wall, the second piece placed within the wall, and the third piece placed on the indoor surface of wall. Bolattürk [9] in 2006 selected 16 cities from 4 climate zones in Turkey to analyze the optimum insulation thicknesses, energy savings, and payback periods for various fuels and climate zones. Five different fuels (including coal, natural gas, fuel oil, liquefied petroleum gas, and electricity) and polystyrene insulation material were considered. Results show that optimum insulation thicknesses vary between 2 and 17 cm, energy savings vary between 22% and 79%, and payback periods vary between 1.3 and 4.5 years, depending on the city and the type of fuel. Another study was later reported by the same author [10] in 2008, in which the effects of different base temperatures on the optimum insulation thickness were analyzed in Turkey's warmest zone with respect to the cooling and heating degree hours. For cooling loads, the optimum insulation thicknesses vary between 3.2 and 3.8 cm, the energy savings vary between 8.47 and 12.19 US$/m², and the payback periods vary between 3.39 and 3.81 years, depending on the cooling degree hours. For heating loads, optimum insulation thicknesses vary between 1.6 and 2.7 cm, energy savings vary between 2.2 and 6.6 US$/m², and payback periods vary between 4.15 and 5.47 years. Masoso and Grobler [11] present the findings of a case study with building simulation using the Energy Plus dynamic thermal simulation software, in which wall insulation was varied together with cooling set-point temperature in a hot and dry climate in Botswana. Six sets of simulations (cooling set temperatures between 22°C and 26.5°C) were performed and the effect of the set point on the insulation was discussed. Simultaneously, Kaynakli [12] investigated the variation in annual energy requirement of buildings with various architectural design properties considering long term and current outdoor air temperature records in Bursa. The optimum insulation thicknesses for Bursa varied between 5.3 and 12.4 cm depending on the fuel used. In a subsequent study in 2009 [13], the optimum insulation thicknesses of the external walls, energy savings, and payback periods over a lifetime of ten years were calculated using the P₁-P₂ method for five different fuels, four different insulation materials, and various cities in Turkey. The results show that the optimum insulation thicknesses vary between 1.06 and 7.64 cm, energy savings vary between 19 US$/m² and 47 US$/m², and payback periods vary between 1.8 and 3.7 years, depending on the city and the type of fuel. In another study, Ucar and Balo [14] examined the energy savings and payback periods in the four climatic regions of Turkey for four different insulation materials applied to the external walls in the cases of five different fuel sources. The longest payback period was determined to be 2.25 years in the city of Mersin, where natural gas is used for heating. Daouas et al. [15] determined the optimum insulation thicknesses by developing and using a “Complex Finite Fourier Transform” analytical model because fixed insulation thicknesses were applied throughout Tunisia. The optimum insulation thickness was determined to be 5.7 cm, corresponding to a payback period of 3.11 years and energy savings of 58%. Özkan and Onan [16] determined the optimum insulation thicknesses according to Turkey's regions and used different insulation materials, such as XPS, polystyrene foam, and rock wool, and examined different fuels such as natural gas and fuel oil. Kaynakli [17] compiled the optimum thickness of the thermal insulation material applied to the buildings and the studies determining its impact on energy consumption. Ozel [18] analyzed a system with a 20-year life and showed that the optimum insulation thickness value varies between 5.4 and 19.2 cm, the energy savings vary between 86.26 and 146.05 US$/m², and the payback period varies between 3.56 and 8.85 years. Bektas Ekici et al. [19] calculated the optimum insulation thickness, energy loss, and payback periods for four cities selected from different climatic regions according to the TS 825 (Antalya in the 1st region, Istanbul in the 2nd region, Elazığ in the 3rd region, and Kayseri in the 4th region). The optimum insulation thickness varied between 0.2 cm and 18.6 cm, energy savings varied between 0.038 US$/m² and 250.415 US$/m², and the payback period varied between 0.714 and 9.104 years. Badran et al. [20] performed a comparison between several wall, roof, and floor constructions and several levels of insulation thickness and the heating load calculations for those constructions (steady state analysis and dynamic behavior) were made. Zhou and Zhao [21] studied the energy-saving effects of different building envelop energy-saving technologies which were calculated and were compared with each other in the five major climatic regions of China. They used energy simulation with economic analysis to evaluate the optimum EPS insulation thickness for office buildings. Shekarchian et al. [22] determined the optimum insulation thickness influenced by local electricity tariff rate and the capital insulation outlays in Malaysia. And also some of the commonly used insulators available in the Malaysian market were analyzed. Ozel [23] determined the optimum insulation thicknesses, energy savings, and payback periods by using life cycle cost analysis. Two different wall structures and two different insulation materials (XPS and EPS) are used in calculations. Finally, they calculated fuel consumption and emissions of CO₂ and SO₂ by taking into consideration wall orientations for different structure and insulation materials. It is seen that optimum insulation thicknesses for south, north, east, and west orientations are obtained as 9.2, 10.2, 9.8, and 9.8 cm, respectively. Fang et al. [24] presented a technical-economic model for optimizing the thermal insulation thickness using Lagrange's method of undetermined multipliers. The approach was applied to a residential building in Shanghai in hot summer and cold winter zones and showed that the insulation on the west wall provided the largest savings, whereas the insulator on the south wall provided the smallest. Huang and Ye [25] established a mathematical model to calculate the optimum insulation thickness of external walls to estimate the life cycle cost including energy consumption and investment of insulation materials. Yu et al. [26] calculated the optimum thicknesses of five insulation materials (expanded polystyrene, extruded polystyrene, foamed polyurethane, perlite, and foamed polyvinyl chloride) with a typical residential wall using solar-air cooling and heating degree-day analysis. Total life cycle costs, savings, and payback periods were then calculated based on the life cycle cost analysis.

Unlike other studies, this investigation is the first one performed using the determination of the buildings’ heat losses. Through the determination of the building inventory in the country, the importance and negative impact of buildings without insulation on energy costs are asserted. Calculation of optimum insulation thickness has been made for natural gas for the model buildings in all climatic regions. And then, the effect of the alteration of insulation thickness on the pollutant emissions such as CO₂ gets out from the chimneys for the different regions which are analyzed. The emissions of CO₂ are decreased by maximum 44.7% for the cases in which optimum insulation material (XPS) and natural gas are used. Furthermore, the importance of this study is emphasized by examining the energy efficiency through the comparison of the optimum insulation thickness determined by the P₁-P₂ method and the minimum allowed insulation thickness for the current building inventory in Turkey.

The calculations were made for four degree-day regions of Turkey shown in Figure 1 with XPS (extruded polystyrene foam) as the insulation material. The density of the insulation material was 25 kg/m³ and thermal conductivity (λ) was 0.030 W/mK. Calculations of the optimum insulation thickness were made for natural gas.

Figure 1:

Four degree-day regions of Turkey.

Figure 2:

Monthly average outdoor temperature values (°C) used in heat loss calculation for the different regions.

Figure 3:

Monthly average solar radiation intensity used in the calculations (W/m²).

2. Annual Heating Energy Requirement Calculation Method

In our country, insulation thickness used in buildings is limited by heating energy requirements in TS 825 standard. Its method is similar to DIN 4108 and the method defined in DIN is adapted to our country. Turkey is divided into four regions according to climate conditions and these regions are shown in Figure 1. Annual heating energy requirements and the subequations are shown in (1)–(12). Descriptions and the formulas are shown, respectively.

The annual heating energy requirements for the model buildings are calculated as given in [27]:

Q_{year} = \sum_{j = 1}^{12} [H (T_{in} - T_{out, j}) - η_{j} (ϕ_{in} + ϕ_{s, j})] \times t .

(1)

Monthly average solar energy gain is

ϕ_{s, j} = \sum_{k}^{} r_{j} g_{j} I_{j, k} A_{gl, k} .

(2)

g is the solar energy permeation factor:

g_{j} = F_{w} \times g_{⊥} .

(3)

Unitless monthly average usage of heat gain factor (the ratio of the heat losses to the heat gains) is

η = 1 - e^{H (T_{out} - T_{in}) / (ϕ_{in} + ϕ_{s})} .

(4)

Specific heat loss (H) of the building is calculated by adding the heat loss which occurred in consequence of conduction and convection (H_tr) to the heat loss which occurred in consequence of ventilation (H_ven):

H = H_{tr} + H_{ven} .

(5)

Heat loss occurred in consequence of ventilation:

H_{ven} = 0.264 \times n_{a} \times V_{gross} .

(6)

Heat loss occurred in consequence of conduction and convection:

H_{tr} = U_{ew} A_{ew} + U_{gl} A_{gl} + 0.5 U_{fl} A_{fl} + 0.8 U_{ce} A_{ce} .

(7)

Rearranged version of (5) perimeter of the floor or ceiling area is shown by P:

H = U_{ew} P_{ce} h_{avg} (1 - ε) + ε U_{gl} P_{ce} h_{avg} + \dots + 0.5 U_{fl} A_{fl} + 0.8 U_{ce} A_{ce} + 0.264 \times n_{a} \times V_{gross} .

(8)

Ratio of the glazing area to the exterior wall surface area is shown by ε:

ε = \frac{A_{gl}}{A_{gl} + A_{ew}} .

(9)

The wall conductance U for a typical wall that includes a layer of insulation is

U_{ew} = \frac{1}{R_{in} + R_{cw} + R_{ins} + R_{out}} .

(10)

R_ins is the thermal resistance of the insulation layer:

R_{ins} = \frac{L_{ins}}{λ_{ins}} .

(11)

Rearranged version of (10) is

U_{ew} = \frac{1}{0.13 + 0.02 / 1 + 0.19 / 0.45 + L_{ins} / λ_{ins} + 0.03 / 1.6 + 0.04} .

(12)

The assumptions used in (1)–(12) are also shown in Table 1.

Table 1:

Assumptions used in model calculations.

Equation no.	Assumptions and descriptions

(1)	ϕ_in is the average interior heat gain; 5 W/m² per unit usage area is considered to be maximum as the interior gain in residences, schools, and the buildings equipped normally (such as office buildings). Monthly average outdoor temperature values are stated in Figure 2 and monthly average interior temperature value is taken as 19°C.

(2)	r is the monthly average shading factor of the transparent surfaces and isconsidered to be 0.6 due to the fact that generally buildings are located around uninsulated buildings or there is shading caused by trees in Turkey. I values are taken from TS 825 according to direction (k) and month (j). Values of the monthly average solar radiation intensity used for calculation of the solar radiation coming on the windows are given in Figure 3.

(3)	F_w is the correction factor for windows and is considered to be 0.8. g_⊥ is the solar energy permeation factor measured under laboratory conditions for the rays striking the surface vertically and is considered to be 0.75 for colorless insulation glass [27].

(6)	Air changing ratio (n_a) is taken as constant “0.8” according to TS 825 which is similar to DIN 4108 for all regions (natural ventilation).

(7)	Total surface area of exterior wall, floor area or ceiling, and window areas are changing according to model buildings. Heat transfer coefficients recommended for ceilings, floors, and windows according to the regions are received from TS 825 (Turkish standard 825) and they are shown in Table 2 [27].

(8)	The height of each floor in the country is architecturally 2.8 m; this value was taken to be fixed. The floor and ceiling area are taken to be equal to each other in all buildings.

(9)	Glazing area ratio (ε) is taken as constant 0.20 for the model buildings.

(10)	R_in and R_out are the inside and outside air film thermal resistances, respectively; R_cw is total thermal resistance of the composite wall materials without the insulation. Recommended design values for air film thermal resistances on the inner and outer surfaces of a building are R_in = 0.13 m²K/W and R_out = 0.04 m²K/W. Exterior wall is considered to be 2 cm interior plaster, 19 cm hollow brick, insulation material, and 3 cm exterior plaster, respectively.

(11)	L and λ are the thickness and thermal conductivity of the insulation material, respectively. Thermal conductivity (λ) is taken to be 0.030 W/mK for XPS.

Table 2:

Values of the heat transfer coefficients (W/m²K).

Region	U_ce (W/m²K)	U_fl (W/m²K)	U_gl (W/m²K)

1st	0.45	0.70	2.40 (double glazing window covered with low-e)
2nd	0.40	0.60
3rd	0.30	0.45
4th	0.25	0.40

Minimum insulation thickness is calculated by (13) according to TS 825 for all regions:

L_{s}^{} = \frac{\sum_{month} A_{ew}^{} λ Δ T^{}}{Q_{s}^{} \times 0 . \bar{44} V_{gross} + \sum_{month} (η (ϕ_{in} + ϕ_{s}^{}) - M^{} Δ T^{})} - R_{ww} λ .

(13)

3. Determination of the Building Inventory and the Model Building

The building distribution in each city was determined by using the building inventory data taken from the State Institute of Statistics (Prime Ministry, Republic of Turkey) for 81 cities in Turkey, divided into four climatic regions. All of the cities specified in different colors in Figure 1 were grouped, and a model building for each region was created.

Weighted average is used and the average story numbers for cities are calculated by

s_{avg}^{} = \frac{[\sum_{1}^{s} {Bn}_{s} \times s]}{{Bn}_{total}} .

(14)

When the story number distribution in Turkey is examined, the highest story value (3.38) is observed in Istanbul, which is located in the second degree-day region. The lowest story value (1.23) is observed in Iğdır city, which is located in the third degree-day region.

When the distribution of the floor area in Turkey is examined, the highest floor area value (195.85 m²) is observed in Ankara city, which is located in the third degree-day region. The lowest floor area value (116.04 m²) is observed in Aydın city, which is located in the first degree-day region.

Whereas model buildings are created for every region, the average floor area was determined by means of (15). Here, m represents the number of buildings, which varies based on the range of the base area, and n represents the cities found in this region:

A_{f l_{avg}}^{} = \frac{\sum_{1}^{m} [A_{fl}^{m} \times \sum_{1}^{n} {Bn}_{n}^{m}]}{{Bn}_{total}} .

(15)

Likewise, for the floor height of the buildings, the Turkish city average was found using (16) for every region and is based on the cities and the distribution of the number of buildings. Building height is determined by multiplying story number by the height of the floor which is taken to be constant 2.8 as shown in Table 1:

h_{avg}^{} = \frac{\sum_{1}^{m} [h_{}^{m} \times \sum_{1}^{n} {Bn}_{n}^{m}]}{{Bn}_{total}} .

(16)

And the total building number used in (14)–(16) is calculated by

{Bn}_{total} = \sum_{1}^{m} \sum_{1}^{n} {Bn}_{n}^{m} .

(17)

The variation of the floor areas of the buildings by region is given in Figure 4 and the variation of the building story number is shown in Figure 5. Based on these data, the specifications of the model buildings created for the thermal calculation are given in Table 3.

Table 3:

Specifications of the model buildings and building inventory.

Region	Average floor area (m²)	Average story number	Total building number

1st	145.63	1.93	1630329
2nd	157.90	2.36	3504662
3rd	160.85	1.86	2019946
4th	159.79	1.59	683738

Figure 4:

Distribution of the number of buildings in four climatic regions based on the floor area.

Figure 5:

Distribution of the number of buildings in four climatic regions based on the story number.

4. Cost Analysis and Energy Saving, Payback Period Correlations

The P₁-P₂ method of calculating life cycle savings is particularly useful. Therefore, operation and investment costs of the building have been calculated according to the P₁-P₂ method [28]. Investment cost of the insulation is calculated by

C_{ins} = C_{I} L P_{fl} h_{avg} (1 - ε) .

(18)

Where C_I (US$/m³) is the unit cost of the insulation and equals 162.6 US$/m³ for XPS. Annual fuel cost of the building heating system is found by the use of

C_{A, f} = \frac{Q_{year}}{{LHV}_{f} η_{b}} C_{f},

(19)

where LHV_f is the low heating value of the fuel, equals 34485 kJ/Nm³, defines efficiency of the heating system, and is 90% and C_f is the unit cost of the fuel and equals 0.44 US$/Nm³. After determination of the insulation and fuel costs, total net energy saving of the building is calculated by the use of the following equation according to the P₁-P₂ method:

S = Δ E_{f} C_{f} P_{1} - C_{I} P_{2} L P_{fl} h_{avg} (1 - ε) .

(20)

The first cost of thermal insulation may be evaluated as capital investment. It is possible to calculate the present worth value of the amount of net energy savings via insulation using the P₁-P₂ method [28]. This is a practical, well-known method and can be used for optimizing the size of insulation of external walls [8, 10, 14]. P₁ is the present worth factor of a series of N future payments with interest rate (I) per period accounting for the time value of money with inflation rate (d). P₁ and P₂ are used for a building insulation application as a nonincoming process:

P_{1} = \frac{1}{d - i} [1 - {(\frac{1 + i}{1 + d})}^{N}],

(21)

where P₁ is the ratio of the cycle of fuel cost savings to the first year fuel cost savings and N is the lifetime of the investment and where 10 years is considered in the calculations. In this study, the average yearly inflation rate is considered to be 7.66% and the interest rate is considered to be 9.76%, according to the 2011 data of the Turkish Central Bank:

P_{2} = 1 + P_{1} M_{s} - R_{v} {(1 + d_{t})}^{- N},

(22)

where P₂ is the ratio of the life cycle expenditures incurred because of the additional capital investment to the initial investment, M_s is the ratio of the annual maintenance and operation cost to the original first cost, and R_v is the ratio of the resale value to the first cost. P₂ can be taken as 1 if the maintenance and operation cost is zero [29].

In the definition of the specific heat loss of the building (8), if the expression apart from the thermal resistance of the exterior wall is shown by M and the yearly heating energy requirement of the building is defined by (1), net energy saving of the model building is found by

\begin{matrix} S = - (\sum_{j = 1}^{12} [(M + A_{ew} \frac{1}{R_{ww} + L_{ins} / λ_{ins}}) Δ T \\ + (e^{- (M + A_{ew} U_{ew}) Δ T / (ϕ_{in} + ϕ_{s})} - 1) (ϕ_{in} + ϕ_{s})] \times \frac{86400 \times 30}{1000}) \times {({LHV}_{f} η_{b})}^{- 1} \times C_{f} P_{1} + C_{A, f, uni} P_{1} - C_{I} P_{2} L A_{ew} . \end{matrix}

(23)

The payback period of insulation cost, N, can be calculated by setting the net savings function to zero:

\begin{matrix} N = (\ln (1 - C_{I} P_{2} L P_{fl} h (1 - ε) {LHV}_{f} η_{b} (d - i) \times (C_{A, f, uni} {LHV}_{f} η_{b} - C_{f} \sum_{j = 1}^{12} [(M + A_{ew} \frac{1}{R_{ww} + L_{ins} / λ_{ins}}) Δ T + (e^{- (M + A_{ew} U_{ew}) Δ T / (ϕ_{in} + ϕ_{s})} - 1) \times (ϕ_{in} + ϕ_{s})] \times {\frac{86400 \times 30}{1000})}^{- 1})) \times {(\ln \frac{1 + i}{1 + d})}^{- 1} . \end{matrix}

(24)

5. Optimum Insulation Thickness

Optimum insulation thickness is the value where the net energy saving is maximum. Therefore, if the derivative of the net energy saving defined by (23) is taken according to the insulation thickness (L) and then is made to be equal to zero, the following equation giving the optimum insulation thickness is found:

\begin{matrix} L_{opt} = (\sum_{j = 1}^{12} \frac{λ Δ T 2592 C_{f} P_{1}}{{LHV}_{f} η_{b} C_{I} P_{2}} {\times (1 - e^{(M + A_{ew} (1 / (R_{ww} + L_{ins} / λ_{ins}))) Δ T / (ϕ_{in} + ϕ_{s})}))}^{0.5} - R_{ww} λ . \end{matrix}

(25)

Temperature, solar radiation, and heat transfer coefficients vary according to the degree-day regions of Turkey, so that investment costs, operating costs, and optimum insulation thickness are determined using the MATLAB program in this study.

6. Relationship between Environment and Optimum Insulation Thickness

Turkey has 2.86 m³ reserve natural gas at the end of 2006 and the consumption in 2006 was 30.5 billion cubic meters so that our country imported 19.65 bcm from Russia [30]. Specifications of the Russian gas are taken into consideration and the general combustion formula for the complete combustion is given by

\begin{matrix} (0.85 C H_{4} + 0.07 C_{2} H_{6} + 0.03 C_{3} H_{8} + 0.02 C_{4} H_{10} + 0.01 C_{5} H_{12} + 0.018 N_{2} + 0.002 O_{2} + 0.002 C O_{2}) + 1.1 a (O_{2} + 3.76 N_{2}) ⟶ x C O_{2} + z H_{2} O + (4.136 a + t) N_{2} + w O_{2} . \end{matrix}

(26)

Excess air factor is taken to be 1.1 for the natural gas. Carbon dioxide which is discarded into the atmosphere and causes greenhouse effect is an accepted emission as a result of good combustion. The total emission of CO₂ (kg CO₂/year) is calculated by

\begin{matrix} M_{C O_{2}} = \frac{44 x}{M_{natural gas}} M_{f}, \\ M_{C O_{2}} = \frac{44 x}{M_{natural gas}} [\frac{Q_{year}}{LHV η_{b}}], \end{matrix}

(27)

where M_f is the fuel consumption in the boiler system per year and M is the weight of the natural gas.

7. Analysis Results and Discussion

The annual heating energy requirement of the model building for each region was calculated from the data taken from the State Institute of Statistics (Prime Ministry, Republic of Turkey). The energy savings, payback period, and optimum insulation thickness were also determined.

7.1. Model Buildings Results

The annual heating energy requirements for all model buildings according to the regions are shown in Figure 6 for the case using XPS insulation and natural gas. It can be seen that the yearly heating energy requirement decreases due to an increase in the insulation thickness in buildings and that the yearly heating energy requirement changes proportionally with region.

Figure 6:

Annual heating energy requirements for all model buildings.

Changes in the energy savings (US$) per model building for all regions are shown in Figure 7 according to the insulation thickness. The energy savings in the 1st region are calculated to be 2000 US$, whereas the savings would be 8000 US$ in the 4th region. The peak points of the curves show the optimum insulation values. Figure 8 shows the change in the payback period versus insulation thickness for all regions. With increasing insulation thickness, the payback period calculated by the P₁-P₂ method increases as well.

Figure 7:

Energy savings (US$) per model buildings.

Figure 8:

Change of payback period according to insulation thickness.

The energy savings (US$/m²) of a building heated by natural gas, the payback period of the insulation application, and the optimum insulation thickness values for the four regions are shown in Table 4. The optimum insulation thickness calculated using the P₁ − P₂ method with natural gas fuel varies between 3.21 cm and 7.12 cm, depending on the region; the energy savings per area vary between 9.23 US$/m² and 43.95 US$/m².

Table 4:

Analysis results.

Region	Optimum insulation thickness (cm)	Payback period (year)

1st	3.21	1.50
2nd	4.94	1.26
3rd	6.06	1.12
4th	7.12	1.00

7.2. Global Effects and Environmental Impacts

The total energy savings and CO₂ emissions per square meter determined from the total number of buildings were calculated for each model building, as shown in Table 5. For the case in which the optimum insulation thickness is applied in the model building, the total energy savings for each region are seen in Table 5 and the total energy savings for the country are calculated to be 41.7 billion US$.

Table 5:

Energy savings and CO₂ emissions.

Region	Total building number	Energy savings per area	Energy savings per model building	Total energy savings (billion US$)	Total CO₂ emission after insulation (million kg·CO₂/m²·year)	Total CO₂ emission before insulation (million kg·CO₂/m²·year)	Reduction in total CO₂ emission (million kg·CO₂/m²·year)

1st	1630329	9.2	1916	3.12	29.2	47.1	17.9
2nd	3504662	21.6	5630	19.73	105.4	186.2	80.8
3rd	2019946	32.2	6672	13.48	82.3	148.1	65.9
4th	683738	44.0	7838	5.36	36.4	65.9	29.4

Total	7838675	106.9	22057	41.7	253.3	447.2	194.0

The change of the CO₂ emission values according to insulation thickness applied to the model buildings in four climatic regions in Turkey is shown in Figure 9.

Figure 9:

Change of CO₂ emission per model building according to insulation thickness.

The reduction ratios in the CO₂ emission values are given in Table 6 when the optimum insulation thickness was applied in the model buildings. The highest reduction is observed in the 4th and 3rd regions and the minimum reduction is seen in the 1st region by 38%.

Table 6:

Natural gas CO₂ emission results for model buildings.

Region	Natural gas CO₂ emission (kg·CO₂/m²·year)
Region	Without insulation	After insulation	Percent %

1st	28.86	17.90	38.0
2nd	53.13	30.07	43.4
3rd	73.33	40.73	44.5
4th	96.33	53.29	44.7

The CO₂ emissions are shown in Table 7 after and before the insulation application for the current building inventory in Turkey. When the reduction in total CO₂ emissions is examined, the first region has the minimum value and the second region has the maximum value.

Table 7:

Natural gas CO₂ emission results for all buildings in Turkey.

Region	Total CO₂ emission after insulation (billion kg·CO₂/year)	Total CO₂ emission before insulation(billion kg·CO₂/year)	Reduction in total CO₂ emission(billion kg·CO₂/year)

1st	6.06	9.77	3.71
2nd	27.52	48.62	21.10
3rd	17.05	30.70	13.65
4th	6.50	11.75	5.25

7.3. Comparison of the Model Building Results

Model results are compared with the latest optimization studies and the results are given in Table 8. Parameters such as interest rate, inflation rate, wall resistance, cost of fuel, insulation material, and thermal conductivity are different from each other in the calculations, so that results show variation when making comparison. Bektas Ekici et al. [19] used different wall components thicker than those used in this study. Ozel [18] made the model calculations for both heating season and cooling season. And also Ozel [18] accepted lifetime as 20 years. The calculation method is similar with Özkan and Onan [16] but they used constant geometric parameters in their study. In our model, calculations are made for the model building created by the data taken from the State Institute of Statistics (Prime Ministry, Republic of Turkey). Ucar and Balo [13] determined their results according to sandwich wall with insulation material so that their results are lower than this study. Dombayci et al. [5] took insulation cost 29 $/m³ in the case; they used expanded polystyrene as an insulation. However, calculations in this study are done using 162 $/m³, the recent price of extrude polystyrene which is used as an insulation material, so a difference between energy-saving values occurs. Overall, the results shown in Table 8 are similar to this model results.

Table 8:

Comparison of the model building results data with the latest studies.

Description	Author	Year	1st region	2nd region	3rd region	4th region

Optimum insulation thickness (cm)	Bektas Ekici et al. [19]	2012	2.1	3.2	4.1	4.5
	Ozel [18]	2012	—	—	8.2	—
	Özkan and Onan [16]	2011	4.02	7.04	7.46	8.66
	Ucar and Balo [13]	2009	0.5	1.06	1.82	3.14
	Bolattürk [10]	2008	1.6–2.7	—	—	—
	Kaynakli [12]	2008	—	5.3	—	—
	Sisman et al. [6] (coal)	2007	3.3	4.7	6.1	8
	Bolattürk [9]	2006	1.9–3.3	3.6–4.7	4.8–5.9	5.8–7.9
	Dombayci et al. [5]	2006	—	—	7.6	—
	Present study		3.21	4.94	6.06	7.12

Payback period (year)	Bektas Ekici et al. [19]	2012	3.8	2.9	2.4	2.2
	Ozel [18]	2012	—	—	5.18	—
	Özkan and Onan [16]	2011	2.21	2.07	1.65	1.49
	Ucar and Balo [13]	2009	3.14	3.77	2.75	2.49
	Bolattürk [10]	2008	4.15–5.47	—	—	—
	Kaynakli [12]	2008	—	—	—	—
	Sisman et al. [6] (coal)	2007	2.82	2.28	1.89	1.54
	Bolattürk [9]	2006	2.68–4.51	2.10–2.52	1.85–2.08	1.61–1.86
	Dombayci et al. [5]	2006	—	—	1.55	—
	Present study		1.50	1.26	1.12	1

Energy saving ($/m²)	Bektas Ekici et al. [19]	2012	9.3	21.4	35.1	43.4
	Ozel [18]	2012	—	—	125.52	—
	Özkan and Onan [16]	2011	15.01	32.88	48.14	64.65
	Ucar and Balo [13]	2009	19.07	—	26.13	—
	Bolattürk [10]	2008	2.20–6.54	—	—	—
	Kaynakli [12]	2008	—	8.64	—	—
	Sisman et al. [6] (coal)	2007	12.8	22.1	35	56.7
	Bolattürk [9]	2006	—	—	—	—
	Dombayci et al. [5]	2006	—	—	—	—
	Present study		9.2	21.6	32.2	44

8. Conclusion

The optimum insulation thicknesses of model buildings, the energy savings over a lifetime of 10 years, and the payback periods are calculated for the four different climatic regions in Turkey. The results show that the optimum insulation thicknesses vary between 3.21 and 7.12 cm, the energy savings vary between 9.23 US$/m² and 43.95 US$/m², and the payback periods vary between 1 and 8.8 years depending on the regions. The results show that the total energy savings vary between 3.12 and 19.73 billion US$ in Turkey. New correlations are given to determine the optimum insulation thickness for all concept buildings or models. Heat gains, solar radiation, and economical data are affecting the optimum insulation thickness. When the building inventory is examined, the number of buildings in the 1st region is approximately the same as the number of buildings in the 3rd region. For this reason, the distribution of the buildings is a parameter affecting the energy savings. Another factor is the application of the insulation material at the optimum thickness. Energy savings will be reduced if the insulation application is not made at the determined thickness. The savings seen in the country within a short time of 10 years are approximately forty billion dollars and will increase further when the financial returns, energy planning, and investments of the country are considered. For this reason, energy efficiency in buildings is important. When the reduction in total CO₂ emissions is examined for all buildings in Turkey, the first region has the minimum value and the second region has the maximum value. One tree will prevent the release of 2.3 kg CO₂. The twenty-three percent of the total forest area in Turkey is recommended before the insulation applications. It should be noted that, after the insulation, it could be decreased to 13.5%. An increase in the insulation thickness will reduce the energy consumption in buildings. As a result of this, the emission rates of combustion products decrease such as CO₂. In the case in which natural gas is used as a fuel and optimum insulation material (XPS) is applied, the emissions of CO₂ are decreased by maximum 44.7% in the 4th region and minimum 38.0% in the 1st region.

Footnotes

Nomenclature

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

References

Papadopoulos

A. M.

, “State of the art in thermal insulation materials and aims for future developments,” Energy and Buildings, vol. 37, no. 1, pp. 77–86, 2005.

Hasan

, “Optimizing insulation thickness for buildings using life cycle cost,” Applied Energy, vol. 63, no. 2, pp. 115–124, 1999.

Papadopoulos

A. M.

Theodosiou

T. G.

, and Karatzas

K. D.

, “Feasibility of energy saving renovation measures in urban buildings—the impact of energy prices and the acceptable pay back time criterion,” Energy and Buildings, vol. 34, no. 5, pp. 455–466, 2002.

Çomaklı

and Yüksel

, “Optimum insulation thickness of external walls for energy saving,” Applied Thermal Engineering, vol. 23, pp. 473–479, 2003.

Dombayci

O. A.

Gölcü

, and Pancar

, “Optimization of insulation thickness for external walls using different energy-sources,” Applied Energy, vol. 83, pp. 921–928, 2006.

Sisman

Kahya

Aras

, and Aras

, “Determination of optimum insulation thicknesses of the external walls and roof (ceiling) for Turkey's different degree-day regions,” Energy Policy, vol. 35, no. 10, pp. 5151–5155, 2007.

Mahlia

T. M. I.

Taufiq

B. N.

, and Masjuki

H. H.

, “Correlation between thermal conductivity and the thickness of selected insulation materials for building wall,” Energy and Buildings, vol. 39, no. 2, pp. 182–187, 2007.

Özel

, “Dynamic approach and cost analysis for optimum insulation thicknesses of the building external walls,” Journal of the Faculty of Engineering and Architecture of Gazi University, vol. 23, no. 4, pp. 879–884, 2008.

Bolattürk

, “Determination of optimum insulation thickness for building walls with respect to various fuels and climate zones in Turkey,” Applied Thermal Engineering, vol. 26, no. 11–12, pp. 1301–1309, 2006.

10.

Bolattürk

, “Optimum insulation thicknesses for building walls with respect to cooling and heating degree-hours in the warmest zone of Turkey,” Building and Environment, vol. 43, no. 6, pp. 1055–1064, 2008.

11.

Masoso

O. T.

and Grobler

L. J.

, “A new and innovative look at anti-insulation behaviour in building energy consumption,” Energy and Buildings, vol. 40, no. 10, pp. 1889–1894, 2008.

12.

Kaynakli

, “A study on residential heating energy requirement and optimum insulation thickness,” Renewable Energy, vol. 33, no. 6, pp. 1164–1172, 2008.

13.

Ucar

and Balo

, “Effect of fuel type on the optimum thickness of selected insulation materials for the four different climatic regions of Turkey,” Applied Energy, vol. 86, no. 5, pp. 730–736, 2009.

14.

Ucar

and Balo

, “Determination of the energy savings and the optimum insulation thickness in the four different insulated exterior walls,” Renewable Energy, vol. 35, no. 1, pp. 88–94, 2010.

15.

Daouas

Hassen

, and Aissia

H. B.

, “Analytical periodic solution for the study of thermal performance and optimum insulation thickness of building walls in Tunisia,” Applied Thermal Engineering, vol. 30, no. 4, pp. 319–326, 2010.

16.

Özkan

D. B.

and Onan

, “Optimization of insulation thickness for different glazing areas in buildings for various climatic regions in Turkey,” Applied Energy, vol. 88, no. 4, pp. 1331–1342, 2011.

17.

Kaynakli

, “A review of the economical and optimum thermal insulation thickness for building applications,” Renewable and Sustainable Energy Reviews, vol. 16, no. 1, pp. 415–425, 2012.

18.

Ozel

, “Cost analysis for optimum thicknesses and environmental impacts of different insulation materials,” Energy and Buildings, vol. 49, pp. 552–559, 2012.

19.

Bektas Ekici

Aytac Gulten

, and Aksoy

U. T.

, “A study on the optimum insulation thicknesses of various types of external walls with respect to different materials, fuels and climate zones in Turkey,” Applied Energy, vol. 92, pp. 211–217, 2012.

20.

Badran

A. A.

Jaradat

A. W.

, and Bahbouh

M. N.

, “Comparative study of continuous versus intermittent heating for local residential building: Case studies in Jordan,” Energy Conversion and Management, vol. 65, pp. 709–714, 2013.

21.

Zhou

and Zhao

, “Optimum combinations of building envelop energy-saving technologies for office buildings in different climatic regions of China,” Energy and Buildings, vol. 57, pp. 103–109, 2013.

22.

Shekarchian

Moghavvemi

Rismanchi

Mahlia

T. M. I.

, and Olofsson

, “The cost benefit analysis and potential emission reduction evaluation of applying wall insulation for buildings in Malaysia,” Renewable and Sustainable Energy Reviews, vol. 16, no. 7, pp. 4708–4718, 2012.

23.

Ozel

, “Thermal, economical and environmental analysis of insulated building walls in a cold climate,” Energy Conversion and Management, vol. 76, pp. 674–684, 2013.

24.

Fang

, and Liu

, “Economic thermal insulation thickness analysis of airconditioned buildings,” journal of Tsinghua University, vol. 42, no. 5, pp. 625–628, 2002 (Chinese).

25.

Huang

and Ye

, “Economical optimum of insulation thickness on external wall of energy saving building,” Build Energy Environment, vol. 24, pp. 73–76, 2005 (Chinese).

26.

Yang

Tian

, and Liao

, “A study on optimum insulation thicknesses of external walls in hot summer and cold winter zone of China,” Applied Energy, vol. 86, no. 11, pp. 2520–2529, 2009.

27.

Turkish Standard 825 (TS 825), Official Gazette Number 27019, Ankara, 2008.

28.

Duffie

A. J.

and Bechman

W. A.

, Solar Engineering of Thermal Processes, John Wiley & Sons, New York, NY, USA, 2nd edition1980.

29.

Söylemez

M. S.

and Ünsal

, “Optimum insulation-thickness for refrigeration applications,” Energy Convers Manage, vol. 40, pp. 13–21, 1999.

30.

Ghosh

K. T.

and Prelas

M. A.

, Energy Resources and Systems: Volume 1: Fundamentals and Non-Renewable Resources, Springer, 2009.