A Multi-Parametric Mathematical Approach on the Selection of Optimum Insulation Thicknesses in Buildings

Detailed simulations have indicated optimum insulation thicknesses of walls’ insulation for a variety of cases. Simplified analytical relations have also been proposed to the same aim, allowing the extraction of more general results, with limited accuracy however, as imposed by mathematical simplicity requirements. In this sense, a variety of important parameters are ignored, such as: the existence of any glazing at the wall, the absorptance of the wall, the base temperature of the heated space which the referred to wall belongs to and its variation with insulation, the thermal characteristics (thermal capacitance, total heat losses coefficient) and the heat and solar gains of the heated space. An alternative analytical approach is consequently developed here, incorporating all above parameters and in this context accessing the wall as part of the whole heated space, instead of considering it solely as an isolated fabric element. The approach consists of a set of two implicit equations which are easily solved, and enables the investigation of the effects of all principal and secondary parameters on the optimum thickness. The ignorance even of the secondary of these parameters has proved to lead to significant errors.


Introduction
Thermal insulation of walls is regarded as a medium to highly cost-effective energy upgrade measure in buildings retrofits.In Greece, for instance, which lies at the southern part of Europe, insulation of external walls may cut the heating energy consumption down by 50% [1] and achieve pay-back in a range of ten years [2,3].In this framework, a tremendous potential of energy savings is formulated, considering that a great part of buildings were built before the two successive oil price crises in the 1970s and hence are insufficiently insulated.For instance, in Greece, only half of the dwellings have some kind of envelope thermal insulation [4].Some issues arise when retrofitting for insulation, with the principal being: the exact placement of the insulation (outside, inside, in cavity), the selection of insulation material and the required thickness for the insulating layer, with the last two interfering with each other.Insulation at buildings' envelopes is regarded considering cost-effectiveness, as we observe a critical thickness after which the achievable energy cost savings do not pay back the further increment in insulation cost.In this sense, an optimum (and at the same time, economically maximum) insulation thickness (OIT) is introduced.
Optimum insulation thickness can be the outcome of detailed simulation [5].Simplified analytical relations have also been proposed to the same aim, allowing the extraction of more general results, with limited accuracy however, as imposed by mathematical simplicity requirements.A great number Buildings 2017, 7, 15 2 of 26 of works on the determination of OIT have been elaborated in the past.An extended review on this topic is given in [6] and is not repeated here for economy of space.The problem of finding OIT is coped with the minimization of life cycle heating (and/or cooling) costs needed to compensate for the heat losses (and/or gains) through the considered fabric element.Both energy (fuel) and capital (insulation cost) expenses are accounted for and the issue renders to a single-objective optimization problem (minimization of the net present value of global costs): Minimize{(Annual heating/cooling costs)•PWF(N, r) + (Insulation cost)} (1) where PWF(N, r) = [(1 + r) N − 1]/[(1 + r) N •r] is the present worth factor for an N years' constant annual cash flow and a discount rate of r.Based on the above criterion, the researchers follow one of two distinguished alternative methods to find OIT, the first being a detailed simulation and the second an analytical relation based on simplification of the degree-days method.The application of detailed simulation leads to more accurate results by considering the effect of many factors.Nevertheless, the arising conclusions are site and case specific and, for this reason, of limited value to the designers/engineers.On the other hand, the proposed analytical relation ignores important factors, such as the effect of the insulation on the base temperature of the heated space and in this context it is set here under questionable accuracy.
Several researchers applied degree-days [7][8][9][10] and concluded to the same simple and explicit relation for the OIT H (the subscript H denotes the assumption of heating operation only): where DD H are the heating degree-days, C F the cost of fuel, H U its calorific value, n HS the efficiency of the heating system, k INS the thermal conductivity of the insulation material, C INS the cost of insulation and R W the initial thermal resistance of the wall (coefficient 293.94 stands for conversion of Watt-days to Joule).
On the other hand, Kayfeci et al. [11] found that the optimum insulation thickness may also depend on the cooling degree-days DD C and the coefficient of performance of the cooling equipment COP.In the same context, Bolatturk [12] proposed two alternative values for OIT one for the heating and one for the cooling operation, while Kurekci [13] proposed furthermore an OIT Y considering an all year around operation.In the same framework, a general relation to estimate OIT Y , through considering both the heating and cooling energy needs, was proposed [13][14][15]: where C E is the cost of electricity.
Although very simple and fast in estimating OIT, the degree-days concept and the subsequent simplified relations, Equations ( 2) and (3), have been criticized for inconvenience, by ignoring the effect of solar radiation and thermal mass [15,16].Actually, there are a few other parameters that have also been confirmed as affecting OIT significantly, but are not included yet in the above relations.For instance, wall orientation was found to affect OIT by 0.5 cm for cold climates like the one in Elazig, Turkey (2653 heating degree-days) [17] up to 1.6 cm in warmer climates such as in Tunis (550 heating degree-days) [15].Yuan et al. [18] investigated the optimum combination of external surface reflectivity and insulation thickness for several cities all over Japan thus demonstrating the effect of reflectivity on OIT.Other researchers examined alternatively the effect of external surface absorptance [17] and color [18] on OIT.In [16] it was examined the effect of shading on OIT, and demonstrated that full shading may decrease OIT by 3-3.5 cm (for a solar absorptance of 0.6).The significant effect of indoor temperature on OIT was also demonstrated [15], with the latter increasing remarkably with the winter base temperature.Other researchers investigated the effect of glazing area on OIT [19]; although they found that this does primarily affect the thermal load and consequently the achievable energy savings, it has a smaller effect on OIT which however can still reach up to 10% for the warmer areas of Turkey (region 1 of the Country's climatic zones).For any of the above reasons, the proposed simplified analytical relations may prove inaccurate, when comparing their results to more accurately calculated OIT by using, for instance, Complex Finite Fourier Transform [15].In this context, the need for amending the above relations, by incorporating additional parameters that proved to be of similar importance, arises.
In the previous publications on OIT, the optimization problem is regularly considered as an isolated one, focusing on the fabric element, and thus ignoring other data of the heated space which this element belongs to, such as the thermal transmittance of its other elements, the heat and the solar gains of the space, the thermal capacitance of the space etc.Only relatively recently, researchers took into consideration other data of the heated space [19].As a consequence they found that insulation is more effective (energy savings are larger and payback period is shorter, for the same insulation thickness) in buildings with smaller window areas.In their calculations they assumed a monthly time step and a utilization factor for the gains.More specifically, they investigated the effects of an alteration of the glazing percentage from 10% to 50% of the wall area and indicated that optimum insulation thickness varies between 0.036 and 0.087 m (for extruded polystyrene foam and use of natural gas).In this way they indicated that OIT should not be examined by considering the fabric element alone, but through considering the integration of the latter to the heated space as well, without however concluding to any mathematical expression.
Although having indicated the dependence of OIT on characteristics of the heated space, the above researchers have not concluded to any relations that could be used as a guide to this aim.A comprehensive parametric analysis has been elaborated to this aim [20] but it was entirely based on Equation (3) ignoring in this context other parameters of the heated space that also affect OIT.In the present work we exactly assume thermal insulation of a fabric element as a more generalized issue, by considering it as an interfering with the heated space element.As a consequence, we investigate the effects of factors not included within the previously suggested relations, like the total heat losses coefficient of the space, the heat gains and all relevant parameters affecting these quantities.To keep our approach simple and broadly applicable, we also attempt to solve the problem analytically.

Presentation of the Problem
The suggested analytical expressions Equations ( 2) and (3) include the heating degree-days which are usually calculated on a typical basis, e.g., for heating at 18 • C or 15.5 • C (used in the United Kingdom), or 65 • F (used in the United States) etc. [21].The base temperature of a building may vary from month to month, but it is satisfactorily accurate to neglect this variation and assume a mean base temperature for the whole winter period (this is indicated later in the text, in Section 5.1) performing the calculations on the basis of the heating season (as also foreseen, among others, in EN ISO 13790).Emphasis is given on climates where heating degree-days exceed cooling degree-days and heating is the major concern, as it happened, for instance, with all European Countries.
By adding insulation to a fabric element, the total heat losses coefficient of the heated space do also change and consequently the same follows suit with the base temperature of the space.The greater the gains are, the more significant the effect is, as is indicatively shown in Figure 1.The figure is based on detailed simulation results through using TRNSYS, regarding a 100 m 2 apartment in Athens, Greece.From the figure it becomes apparent that insulation does not play any role in the balance temperature only when there are no heat gains.In contrast, the existence of gains results in a lower base temperature (the performance curve is shifted downwards).Moreover, the addition of insulation does also reduce further the base temperature (provided there are heat gains).The degree-days method offers a mathematical explanation for this behavior [22].Anyway, the effect of the gains on the base temperature should be accounted for, when determining the optimum insulation thickness.

General Remarks
An analytical approach is followed here, in an attempt to conclude on more general results.To this aim, the degree-days method is also applied, and analytical relations are used (although approximate) e.g., for the heat gains utilization factor and for the heating degree-days.Although introducing some errors, these approximations have a minor impact on the determination of OIT due to the discrete thicknesses at which the insulation slabs are available.For instance, in Greece the respective thicknesses are 3, 5, 6, 7, 8, 10 and 12 cm; hence, an error smaller than ±5% seems acceptable, since within one standard error (68% probability in normal distribution [23]) do not fall more than one insulation thicknesses.
There are parameters which seem to have a minor effect on OIT.A sensitivity analysis will demonstrate whether these secondary parameters should also be taken into consideration in the optimization or not.We apply design of experiments [24], factorial simulation and analysis of variance (ANOVA), as it was similarly applied for envelope optimization in Malawi [25].
We assume the insulation of an opaque wall as a fabric element of a heated space where the latter generally constitutes part of a broader heated zone.Provided that the temperature distribution within the heat zone is uniform, no heat exchange takes place between the considered heated space and the rest of the zone.This is anyway desirable, to avoid overheating and thermal discomfort within the same heat zone, and can be easily attained (e.g., when heating is applied with radiators, by installing a thermostatic radiator valve (TRV) at each one of the heating elements and so on).

Consideration Regarding Insulation Costs
Fixed and variable costs are distinguished, when applying insulation slabs in retrofits.Fixed costs include the facilities for doing the job and the cost of placing the insulating panels on the wall, which are practically independent of the slab thickness.Insulation is placed preferably at the external side of the wall, hence the need to erect scaffolding; in addition, some materials are needed to secure installation and a decorative finishing touch (coat mesh, primer coat, render coat etc.).The space that is finally occupied by external insulation is not included in the built surface, a fact due to which there are no side costs (e.g., there is no loss of saleable floor area [26]).Labor costs and material costs to fix the wall insulation should be considered additional to the insulating slabs costs.All this however is not proportional to the volume (or equivalent to the thickness) of the insulation material, as is indicatively depicted for EPS in Figure 2. From this figure it becomes apparent that insulation material constitutes only a small part of the total cost, which is distinguished to a fixed

General Remarks
An analytical approach is followed here, in an attempt to conclude on more general results.To this aim, the degree-days method is also applied, and analytical relations are used (although approximate) e.g., for the heat gains utilization factor and for the heating degree-days.Although introducing some errors, these approximations have a minor impact on the determination of OIT due to the discrete thicknesses at which the insulation slabs are available.For instance, in Greece the respective thicknesses are 3, 5, 6, 7, 8, 10 and 12 cm; hence, an error smaller than ±5% seems acceptable, since within one standard error (68% probability in normal distribution [23]) do not fall more than one insulation thicknesses.
There are parameters which seem to have a minor effect on OIT.A sensitivity analysis will demonstrate whether these secondary parameters should also be taken into consideration in the optimization or not.We apply design of experiments [24], factorial simulation and analysis of variance (ANOVA), as it was similarly applied for envelope optimization in Malawi [25].
We assume the insulation of an opaque wall as a fabric element of a heated space where the latter generally constitutes part of a broader heated zone.Provided that the temperature distribution within the heat zone is uniform, no heat exchange takes place between the considered heated space and the rest of the zone.This is anyway desirable, to avoid overheating and thermal discomfort within the same heat zone, and can be easily attained (e.g., when heating is applied with radiators, by installing a thermostatic radiator valve (TRV) at each one of the heating elements and so on).

Consideration Regarding Insulation Costs
Fixed and variable costs are distinguished, when applying insulation slabs in retrofits.Fixed costs include the facilities for doing the job and the cost of placing the insulating panels on the wall, which are practically independent of the slab thickness.Insulation is placed preferably at the external side of the wall, hence the need to erect scaffolding; in addition, some materials are needed to secure installation and a decorative finishing touch (coat mesh, primer coat, render coat etc.).The space that is finally occupied by external insulation is not included in the built surface, a fact due to which there are no side costs (e.g., there is no loss of saleable floor area [26]).Labor costs and material costs to fix the wall insulation should be considered additional to the insulating slabs costs.All this however is not proportional to the volume (or equivalent to the thickness) of the insulation material, as is indicatively depicted for EPS in Figure 2. From this figure it becomes apparent that insulation material constitutes only a small part of the total cost, which is distinguished to a fixed and variable cost.Insulation material falls into the variable costs, together with the mechanical fixings, while all other components constitute the fixed costs.The total insulation costs (C INS,TOT in €/m 2 of insulated surface) get consequently the form: where

The Optimization Criterion
The usually applied criterion is the minimization of life cycle costs (see Equation ( 1)), and the optimization is applied based on the Life Cycle Cost Analysis LCCA [6].The optimization can alternatively be based on the marginal quantities of the insulation cost and the induced energy cost savings.The maximum insulation thickness which is economical to reach, is when the marginal benefits (energy costs savings) cannot counterbalance any longer the marginal cost imposed by the addition of more insulation.In this sense, optimum insulation thickness arises when the derivative of the accumulated benefits get equalized to the derivative of the insulation cost.
Actually, in retrofits where fixed and variable insulation costs are distinguished, a minimum insulation thickness that can pay back the fixed expenses and its own cost comes first (Figure 3).Well after that comes a maximum insulation thickness which is the economically optimum (e.g., greater IRR) and at the same time the thickness beyond which any addition of insulation is not paid back from the respective achievable savings.

The Optimization Criterion
The usually applied criterion is the minimization of life cycle costs (see Equation ( 1)), and the optimization is applied based on the Life Cycle Cost Analysis LCCA [6].The optimization can alternatively be based on the marginal quantities of the insulation cost and the induced energy cost savings.The maximum insulation thickness which is economical to reach, is when the marginal benefits (energy costs savings) cannot counterbalance any longer the marginal cost imposed by the addition of more insulation.In this sense, optimum insulation thickness arises when the derivative of the accumulated benefits get equalized to the derivative of the insulation cost.
Actually, in retrofits where fixed and variable insulation costs are distinguished, a minimum insulation thickness that can pay back the fixed expenses and its own cost comes first (Figure 3).Well after that comes a maximum insulation thickness which is the economically optimum (e.g., greater IRR) and at the same time the thickness beyond which any addition of insulation is not paid back from the respective achievable savings.

Commentary on the Most Important Factors that Affect Optimum Insulation Thickness
Optimum insulation thickness is a function of the building type, shape, orientation, construction materials, climatic conditions, insulation material and cost, energy type and cost, and the type and efficiency of air-conditioning system.Some of these factors are shortly commented below.
Heating degree-days: This is obviously the most important factor.High value of heating degree-days leads to greater achievable energy cost savings by applying insulation, and in this sense it is worth adding thicker layers of insulation.As a consequence, higher value of heating degree-days is reasonably expected to result in a greater optimum insulation thickness.
Cooling degree-days: The DDC also affect optimum insulation thickness.In colder climates the OITH for heating does not significantly differ from the value OITY regarding the heating and cooling operation, as it is depicted in Figure 4, which is based on published data [13].A linear relation between these quantities: (    ⁄ ) = 1.0 + 0.6 • (  ⁄ ) is apparently valid The coefficient in front of the degree-days ratio (value 0.6) results from the division of costs for cooling and heating, correspondingly:  = (  •   •   ) ( •   ) ⁄ .Indeed, by introducing HU = 34,485 kJ/m 3 , n = 0.90, CF = 0.385 $/m 3 from [12], CE = 0.07 $/kWh [8] and COP = 2.5 [12,27] we estimate indeed b = 0.6.As a consequence, a more general relation may arise:

Commentary on the Most Important Factors that Affect Optimum Insulation Thickness
Optimum insulation thickness is a function of the building type, shape, orientation, construction materials, climatic conditions, insulation material and cost, energy type and cost, and the type and efficiency of air-conditioning system.Some of these factors are shortly commented below.
Heating degree-days: This is obviously the most important factor.High value of heating degree-days leads to greater achievable energy cost savings by applying insulation, and in this sense it is worth adding thicker layers of insulation.As a consequence, higher value of heating degree-days is reasonably expected to result in a greater optimum insulation thickness.
Cooling degree-days: The DD C also affect optimum insulation thickness.In colder climates the OIT H for heating does not significantly differ from the value OIT Y regarding the heating and cooling operation, as it is depicted in Figure 4, which is based on published data [13].A linear relation between these quantities: (OIT Y /OIT H ) = 1.0 + 0.6•(CDD/HDD) is apparently valid The coefficient in front of the degree-days ratio (value 0.6) results from the division of costs for cooling and heating, correspondingly: b = (C E •H U •n HS )/(COP•C F ). Indeed, by introducing H U = 34,485 kJ/m 3 , n = 0.90, C F = 0.385 $/m 3 from [12], C E = 0.07 $/kWh [8] and COP = 2.5 [12,27] we estimate indeed b = 0.6.As a consequence, a more general relation may arise: Buildings 2017, 7, 15 6 of 26

Commentary on the Most Important Factors that Affect Optimum Insulation Thickness
Optimum insulation thickness is a function of the building type, shape, orientation, construction materials, climatic conditions, insulation material and cost, energy type and cost, and the type and efficiency of air-conditioning system.Some of these factors are shortly commented below.
Heating degree-days: This is obviously the most important factor.High value of heating degree-days leads to greater achievable energy cost savings by applying insulation, and in this sense it is worth adding thicker layers of insulation.As a consequence, higher value of heating degree-days is reasonably expected to result in a greater optimum insulation thickness.
Cooling degree-days: The DDC also affect optimum insulation thickness.In colder climates the OITH for heating does not significantly differ from the value OITY regarding the heating and cooling operation, as it is depicted in Figure 4, which is based on published data [13].A linear relation between these quantities: (    ⁄ ) = 1.0 + 0.6 • (  ⁄ ) is apparently valid The coefficient in front of the degree-days ratio (value 0.6) results from the division of costs for cooling and heating, correspondingly:  = (  •   •   ) ( •   ) ⁄ .Indeed, by introducing HU = 34,485 kJ/m 3 , n = 0.90, CF = 0.385 $/m 3 from [12], CE = 0.07 $/kWh [8] and COP = 2.5 [12,27] we estimate indeed b = 0.6.As a consequence, a more general relation may arise:  Effect of cooling degree-days to optimum insulation thickness (use of XPS insulating material and as a fuel the natural gas, data from [12]).
An amended value of DD C can be introduced in Equation ( 5), when adopting a different scenario of operation during the cooling period (e.g., unoccupied period due to summer holidays, intermittent operation etc.).For instance, in Greece it is assumed in the energy classification of buildings (based on ISO 13790) that the actual cooling energy needs are only 50% of the total cooling energy demand, hence the half value of DD C should be respectively inserted in the above equation etc. thus causing a limited effect on OIT Y .By taking into consideration that: (a) cooling degree-days affect OIT at a lesser degree than heating degree-days (b is generally much lower than unit) (b) for the weather conditions prevailing in the European Countries and more general in Countries with continental climates, heating accounts for a significant part of the total energy consumption (c) cooling load is only partly due to transmission gains through opaque surfaces (d) according to regulations, cooling may be considered covering part of the demand (e) optimum insulation thickness for year around operation can be easily extracted, by applying Equation ( 5) from optimum thickness for heating, we further consider in this work the determination of optimum insulation for heating operation only.
Existing insulation: The higher the existing insulation is, the lesser the room for further savings.Hence, preexistence of any insulation on the wall leads to a lower optimum insulation thickness.It is easily proved that the total insulation thickness should be anyway the same between pre-insulated and non-insulated walls, when pre-existing insulation (e.g., of thickness X PI ) is also accounted for in the total insulation thickness (provided that the same insulation material is considered).This is mathematically expressed by the relation: Such an issue is expected to arise when passing from already satisfactorily performing buildings to nearly zero energy buildings.Indeed, at an already insulated building, with X PI , the X OPT (X PI ) is expected to be much smaller and the question arising will be if it is worth proceeding with any additional insulation.In that case (but even more generally as well) the application of optimizing insulation thickness will not be enough but someone has to examine at first the profitability of the measure.If so, the fixed expenses should also be considered, and the following condition arises: NPV = (Annualy avoided energy costs)•PWF(N, r) − (Insulation cost) > 0 (7) Heat and Solar gains: Heat and solar gains result in a lower base temperature, then in reduced heating degree-days, in lower energy needs and in the end in a smaller OIT H value.The heat gains are more or less constant all through the year, while the solar gains vary.The application of insulation may reduce the solar gains through the opaque surface of the insulated wall, which constitutes however a minor part of solar gains in cases where the wall has some glazing surface on it.
Intermittency of operation: Intermittent heating is probably the rule in space heating.A night set-back temperature, and application of heating during occupancy hours only, lead to a lower mean temperature in the heated space and consequently to a reduced heating energy demand.The last fact has an impact on the economy of insulating, and reduces the OIT H , since the achievable energy cost savings in intermittent heating operation (for the same insulation thicknesses) will be less.At the same time, however, the mean internal temperature is also reduced, in a way probably affecting the thermal comfort of the occupants.A compromise is finally desirable between the insulation cost and the thermal comfort, but such an approach is beyond the scope of the present work, where continuous heating is the main consideration.
Absorptance: The solar radiation incident on the external surface of the fabric element affects heat exchange with the environment.This is expressed via the sol-air temperature, which is strongly related to the absorptance and emittance of the surface: where T A is the outdoor air temperature ( • C), F is the shading reduction factor for external obstacles on the surface, α ES is the solar absorptance of the external surface, I T is total solar irradiance incident on the surface (W/m 2 ), h ES is the combined convective and radiative heat transfer coefficient on the external surface (W/m 2 K), ε ES is the hemispherical emittance of the surface and ∆R is the difference between the long-wave radiation incident on the surface from the sky and surroundings and the radiation emitted by a blackbody at outdoor air temperature (W/m 2 ).The value of ∆R can be considered zero for vertical surfaces, because the radiant heat loss to the sky compensates the heat.Hence, solar gains through the opaque surface A* are equivalently expressed by the relation: where U* is the thermal transmittance of the considered fabric element.The higher the absorptance is, the higher the solar gains get, leading to respectively lower heating energy needs and in turn to lower values for OIT H . Economic factors: We distinguish cost factors (cost of insulation, cost of fuel) and financial factors (discount rate, life-time).All quantities are expressed in present values.The application of any loan can also be easily accounted for by converting all future installments on the debt to present values and adding them to the insulation cost.
The fact that the issue of specifying a suitable insulation thickness is multi-parametric arises as a conclusion of the above analysis of the factors affecting OIT.In addition, the wall should not be considered as isolated but as a constituent of the whole heated space.

Mathematical Background and Analytical Determination of the Optimum Insulation
Thickness OIT H

Formulation of the Problem
According to the criterion set in Section 2.4, OIT H will result as the solution of the equalization between marginal benefits and marginal costs.Based on Equation (4), the annuity of the induced cost by implementing a differential thickness of insulation dX INS is expressed, in present values, as: where N is the expected life time of insulation and r the discount rate.The annual benefits are respectively given by the relation: The estimation of the respective decrement of energy heating demand is not so straightforward, and is calculated as follows.According to the degree-days method, the annual energy demand of the heated space is: where DD H (T B ) are estimated at the base (or balance) temperature T B .The total heat loss coefficient TLC of the heated space is given by the relation: where V is the volume of the heated space and N AIR the infiltration rate through the space, expressed in air changes per hour; U i is the thermal transmittance and A i the respective surface area of any other external fabric element i of the considered space.
The fabric element for insulation is separately considered (not included in the summation term), unlike all other opaque surfaces of the heated space and openings which are included in the summation term.In the infiltration losses, only aeration through the openings between the space and the environment is considered.Aeration between the space under consideration and the other spaces of the zone do not contribute to the heating demand as long as the same internal temperature is sustained.The base temperature T B is given by the relation: where T SP is the set-point temperature, Q G the mean heat gains during the heating period and n G the gains utilization factor.The addition of insulation decreases the thermal transmittance of the considered element and, according to Equation ( 9), the same takes place with the respective solar gains through this opaque surface.As a consequence, the heat gains are distinguished regarding their variation with insulation to (i) the variable ones like the solar gains through the opaque surface which is considered for insulation Q GS * and (ii) the fixed ones Q GF which consist of the internal gains Q GI and the solar gains through any glazing and any other opaque surfaces of the space Q GS : The differential of the annual energy demand is consequently: where DD H (T B ) is the derivative of heating degree-days: By applying differentials in Equation ( 14) arises: For continuous heating, the heat gains utilization factor is given by the relation: 15, where τ is the time constant (in hours).A regression relation that gives the best fit to the utilization factor is [28]: with k G and D the regression parameters.According to a more recent study [29], parameter D ranges between (−0.02) and (−0.04) and in this sense is negligible within the accuracy limits of the present work.Parameter k G is given as a function of the thermal capacity per floor area practically ranging between 1.0 and 1.2 for typical constructions with time constants 10 and 20 h, but also for heavy constructions (see Figure 5) as long as ratio γ becomes large (this is expected with the addition of insulation).Setting D ≈ 0, Equation (20) gets the form: It is also valid: where with index o the values of quantities U* and TLC before the application of insulation at the considered fabric element are denoted.
Buildings 2017, 7, 15 10 of 26 The addition of insulation decreases the thermal transmittance of the element, according to the relation: The differential of U* becomes: Heating degree-days can be quite accurately written as a quadratic function of the base temperature [30]: with By using the above relations and applying simple algebra, the following relations arise: The addition of insulation decreases the thermal transmittance of the element, according to the relation: The differential of U* becomes: Buildings 2017, 7, 15 11 of 26 Heating degree-days can be quite accurately written as a quadratic function of the base temperature [30]: with and Here, T REF is the reference temperature used as the basis for the calculation of heating degree-days DD H (T REF ) and T MEAN the mean annual temperature of the site.
By introducing the above second degree approximation in heating degree-days for the estimation of dE H,D in Equation ( 28), and by equalizing marginal costs and marginal benefits (dC TOT = dB TOT ), the following implicit equation of U * finally arises: Coefficient 0.024 is introduced to convert W-days to kWh, assuming that the fuel cost is expressed per kWh of fuel calorific value.Here, n G is also a function of U * : The last set of equations, Equations (34) and (35), constitutes an implicit system of equations for the estimation of the optimum thermal transmittance U* that corresponds to the OIT H .

Application of a Trial and Error Process
The accurate solution of the implicit set of Equations ( 34) and ( 35) is of low practical value for rigid insulating slabs, as they are available at discrete thicknesses.Equation (34) has resulted by subtracting insulation costs from energy cost savings, and for this reason a positive value of F(U * ) does precisely mean that the benefits by insulating exceed the respective costs.In this way, the various available insulation slab thicknesses can be tried, and the greater one that still results in a positive value for F(U*) can finally be selected.

Successive Substitution
Equation (34) can be solved for U* to get the form: By setting initial value of U* = U 0 *, and introducing it consecutively to Equation (35) and to the right side of Equation (36), a new U* arises to the left side of Equation (36); this is also inserted in Equation ( 35) and to the right side of Equation (36) and so on.By these consecutive substitutions, convergence to the optimum value of U* is rapidly achieved (practically, two to three successive substitutions are sufficient).The whole approach is indicatively depicted in Figure 6.

Use of a Solver
The implicit set of Equations ( 34) and (35) can easily be solved by any solver like Mathematica or Maple (the latter is also integrated within MATLAB).Alternatively, the Excel software package (and its add-in solver) can also be used.

Approximate Graphical and Analytical Solutions
Graphical Solution Solar gains, through the opaque surface of the considered fabric element are expected to constitute a minor part of total gains.For this reason, their variation with the thermal transmittance of the element can be ignored, without introducing a significant error.Then, the gains Q GS * are included in the fixed gains: By removing Q GSO * from Equations ( 34) and (35), they get the following simpler forms: Setting: the system of equations ( 38) and (39) renders to a single equation of z: or equivalently to: By applying the Lambert function W to both parts of Equation (44) (the Lambert function is defined as the inverse relation of the function F(z) = z•e z ) we get: At last, by transforming z: the Equation (45) gets the final form: Equation ( 47) can be solved graphically in a diagram of Lambert function, or by using the discrete values of Lambert function given in Table 1.After finding the root z', we calculate z, then U* and finally the required OIT H .

Analytical Solution
The set of Equations ( 38) and (39) can be solved analytically only when the variation of the base temperature with the addition of insulation is ignored.Then, the set is solvable for U*: which, if solved for X INS , becomes practically identical to the previously suggested Equation (3) where the base temperature is now corrected to the actual balance temperature of the heated space (just before the application of insulation).

Analytical Solution
The set of Equations ( 38) and (39) can be solved analytically only when the variation of the base temperature with the addition of insulation is ignored.Then, the set is solvable for U*: which, if solved for XINS, becomes practically identical to the previously suggested Equation ( 3) where the base temperature is now corrected to the actual balance temperature of the heated space (just before the application of insulation).

Demonstrative Example
The present analytical approach is demonstrated through the following example.A heated space with dimensions of 5 × 5 × 3 m 3 (W × L × H), having an external wall facing to the west which additionally includes a double-glazed window, is in thermal equilibrium with adjacent, above and below heated spaces (see Figure 7).The wall is assumed to be consisted of three successive layers, an internal plaster layer of 0.02 m, a brick layer of 0.18 m and an external plaster layer of 0.02 m.Assuming internal and external thermal resistances of 0.13 m 2 -K/W and 0.04 m 2 -K/W, respectively, the initial thermal transmittance of the wall is estimated at 1.613 W/m 2 -K.The following are assumed for the calculation of internal gains.Population density 5 persons/100 m 2 , installed power for lighting 3.6 W/m 2 , use of lighting 6 h per day, sensible heat released 80 W/person, presence factor 0.75, installed capacity for equipment 4 W/m 2 , intermittency factor 0.5.According to them, the 24 h mean internal gains are 135 W.
In Table 2 the incident solar radiation values for various wall orientations in Athens are presented.In Table 3 the monthly temperatures and calculated mean solar irradiance on the west facing wall, the total gains (including internal gains, solar gains through glazing and through the opaque surface), the gains utilization factor and finally the base temperature for each month, are additionally depicted.In this sense, the months in which heating is necessary finally arise (the base temperature exceeds the monthly temperature).The mean value of monthly base temperatures of the heating period (from October to April, here) is calculated at 15.8 • C. Alternatively, through using the mean quantities for heat losses and gains for the whole heating period (presented in the last line of the same table), the base temperature would have been estimated at 15.5 • C, which reveals that the use of mean data for the whole heating period introduces some error.This small error however is negligible as equally concerning both the previously suggested analytical relation and the analytical approach proposed here, since neither applies variable-base degree-days.The necessary data for the calculations are further presented in Table 4.
In Table 5 other input necessary for the application of the relations proposed here, which are also calculated from the data of Tables 3 and 4 are presented.The application of the presented alternative methods to solve this problem is indicated below.

A. Solution by checking available thicknesses of insulating slabs
For each one of the available insulation thicknesses we calculate U* by using Equation ( 29), then n G by using Equation (35) and last F(U*) by using Equation (34), as shown in Table 6.According to these trials, it seems that there is no benefit to reach a thickness of 0.10 m, hence a slab of 0.08 m seems a safe selection instead.Noticeably, F(U*) is slightly negative at 0.10 m, thus if slabs were available at continuous thicknesses, the optimum thickness would be a value closer to 0.10 than to 0.08 m.Starting from the initial value of U O * = 1.613W/m 2 -K (for X INS = 0) we rapidly conclude to the final value of OIT H = 9.78 cm just after two consecutive substitutions, as indicated in Table 7.It is noticed here that we could alternatively start from the thermal transmittance set by the regulations (e.g., 0.5 W/m 2 -K) to reach even faster the final solution.

C. By applying a solver
By applying the solver of MATLAB we estimate OIT H = 9.7796676 cm.Hence, the 0.10 m slab is indeed marginally acceptable.

D. By applying a graphical solution
In this case, we include the solar gains through the opaque surface (8 W) to the fixed gains, and so total fixed gains reach the value of Q GF = 268 W. We estimate: Which can be graphically solved (see Figure 8) to give the root z' = 0.0278.Indeed, based on the figure we realize that the solution is between the z' values 0.02 and 0.03.We apply linear interpolation to the data of Table 3, for z' values 0.02 and 0.03: W(z') = 0.95z' + 0.006, and then by setting this equal to W(z') = 3.7380z' − 0.0769 we get z' = (0.0769 + 0.0006)/(3.7380− 0.95) = 0.0277977.
and the optimum insulation thickness: This graphical solution has an error of only 0.5% in comparison with the accurate solution attained through using a solver (0.0978 m).thus overestimating by +22% the optimum insulation thickness.

Effect of Orientation
All three alternative methods suggested (trial and error, successive substitution and graphical solution) lead quickly and satisfactorily to the optimum thermal transmittance calculated accurately with a solver (0.33 W/m 2 -K) and to the respective optimum insulation thickness (0.098 m).
The analytical relations instead, introduce a significant error (+22%) when a typical base temperature is used (e.g., 18 °C) which is reduced to 7.4% in this specific example, when the actual base temperature is applied.
and the optimum insulation thickness: This graphical solution has an error of only 0.5% in comparison with the accurate solution attained through using solver (0.0978 m).

E. By applying the analytical relation
By applying Equation (48) we calculate: which leads to the value OIT H = 0.1050 m which has an error of 7.4% (comparably to the accurate value of 0.0978 m).If a base temperature of 18 • C was used instead, then by using Equation ( 2) and introducing C F /H U = 0.08 €/kWh = 2.222•10 −8 €/J, we would calculate: thus overestimating by +22% the optimum insulation thickness.

Effect of Orientation
All three alternative methods suggested (trial and error, successive substitution and graphical solution) lead quickly and satisfactorily to the optimum thermal transmittance calculated accurately with a solver (0.33 W/m 2 -K) and to the respective optimum insulation thickness (0.098 m).
The analytical relations instead, introduce a significant error (+22%) when a typical base temperature is used (e.g., 18 • C) which is reduced to 7.4% in this specific example, when the actual base temperature is applied.
The calculated OIT for Athens (9.8 cm) is in agreement with relevant published results of detailed simulation (e.g., an OIT of 8 cm was calculated for Iraklion-it is at the Southern part of Greece, heating degree-days 1056 K-days-in the frame of energy efficient nZEB design [31]).The respective optimum thermal transmittance U* = 0.33 W/m 2 -K is almost 35% lower than the corresponding maximum permitted value set by the Greek regulation for the assumed climatic zone (U MAX = 0.50 W/m 2 -K).
The above assumed relations can be taken into consideration in relation to several parameters, such as the orientation which has been proved to be critical in OIT Y determination [32] but not regarded in the analytical expressions Equations ( 2) and (3).Indeed, let's assume that the wall has a due south orientation.The fixed gains and the variable solar gains are now 327.6 W and 27.6 W, respectively.Solving again as above, we find OIT H = 9.093 cm.
Last but not least, by assuming the north orientation, we now estimate the same heating period (seven months) but the fixed gains and the variable solar gains arise at present to 191.1 W and 8.0 W, respectively.Solving again as above, we find OIT H = 10.541 m.In this case, the installation of a 0.10 m insulating slab can be proposed.
The above estimations are depicted in Figure 9 together with the respective base temperatures, before and after the application of optimum insulation.It is firstly apparent that the assumption of a typical base temperature, such as 18 • C, is misleading as base temperatures are even 4 • C lower.Secondly, the effect of insulation on the base temperature which, for the assumed data, ranges from 0.7 to 1.0 • C also becomes apparent documenting once again the feasibility to take into consideration the base temperature variation when searching for the optimum insulation thickness.
Buildings 2017, 7, 15 19 of 26 The calculated OIT for Athens (9.8 cm) is in agreement with relevant published results of detailed simulation (e.g., an OIT of 8 cm was calculated for Iraklion-it is at the Southern part of Greece, heating degree-days 1056 K-days-in the frame of energy efficient nZEB design [31]).The respective optimum thermal transmittance U* = 0.33 W/m 2 -K is almost 35% lower than the corresponding maximum permitted value set by the Greek regulation for the assumed climatic zone (UMAX = 0.50 W/m 2 -K).
The above assumed relations can be taken into consideration in relation to several parameters, such as the orientation which has been proved to be critical in OITY determination [32] but not regarded in the analytical expressions Equations ( 2) and (3).Indeed, let's assume that the wall has a due south orientation.The fixed gains and the variable solar gains are now 327.6 W and 27.6 W, respectively.Solving again as above, we find OITH = 9.093 cm.
Last but not least, by assuming the north orientation, we now estimate the same heating period (seven months) but the fixed gains and the variable solar gains arise at present to 191.1 W and 8.0 W, respectively.Solving again as above, we find OITH = 10.541 m.In this case, the installation of a 0.10 m insulating slab can be proposed.
The above estimations are depicted in Figure 9 together with the respective base temperatures, before and after the application of optimum insulation.It is firstly apparent that the assumption of a typical base temperature, such as 18 °C, is misleading as base temperatures are even 4 °C lower.Secondly, the effect of insulation on the base temperature which, for the assumed data, ranges from 0.7 to 1.0 °C also becomes apparent documenting once again the feasibility to take into consideration the base temperature variation when searching for the optimum insulation thickness.

Revealing the Main Factors that Affect Optimum Insulation Thickness
Sensitivity analysis is further applied to the data and the results of the example, in discovering which are the most critical parameters affecting optimum insulation thickness in an area.To this aim, we vary (e.g., by ±20%) the main parameters of Equations (34) and (35).It is noticed that some of these quantities appear only in the form of a common factor, as is the case with the parameters under the square root of the equation, the effect of which is reasonably evaluated through the respective composite factor:

Revealing the Main Factors that Affect Optimum Insulation Thickness
Sensitivity analysis is further applied to the data and the results of the example, in discovering which are the most critical parameters affecting optimum insulation thickness in an area.To this aim, we vary (e.g., by ±20%) the main parameters of Equations (34) and (35).It is noticed that some of these quantities appear only in the form of a common factor, as is the case with the parameters under the square root of the equation, the effect of which is reasonably evaluated through the respective composite factor: In this context, we vary the considered parameters within the ranges shown in Table 8 and plot their effect on the optimum thermal transmittance in a normal probability diagram (Figure 10).In this context, we vary the considered parameters within the ranges shown in Table 8 and plot their effect on the optimum thermal transmittance in a normal probability diagram (Figure 10).From Figure 10, the composite quantity S arises as the most important parameter.This is quite reasonable for an economic optimization issue, as it incorporates all financial parameters.The strong volatility of energy prices, as is indicatively shown for natural gas in Figure 11, and in turn the variability of factor S, should be noticed here however.In this sense, any proposals for OIT should be thoroughly considered, as long as life cycle savings are accounted for.For instance, the drop in the price of natural gas from 0.08 to 0.06 €/kWh-GCV results in a 33% increase in factor-S leading to 15% higher U * value and finally to 16% smaller OITH (0.0819 m instead of 0.0978 m calculated earlier).From Figure 10, the composite quantity S arises as the most important parameter.This is quite reasonable for an economic optimization issue, as it incorporates all financial parameters.The strong volatility of energy prices, as is indicatively shown for natural gas in Figure 11, and in turn the variability of factor S, should be noticed here however.In this sense, any proposals for OIT should be thoroughly considered, as long as life cycle savings are accounted for.For instance, the drop in the price of natural gas from 0.08 to 0.06 €/kWh-GCV results in a 33% increase in factor-S leading to 15% higher U * value and finally to 16% smaller OIT H (0.0819 m instead of 0.0978 m calculated earlier).In addition, the heat gains, the incident solar radiation and the solar gain heat coefficient, have significant positive effects (on the optimum thermal transmittance).An increase of these factors leads to higher thermal transmittance values, or equivalently to lower values of economical insulation.On the other hand, the thermal transmittance of the window in the heated space and the infiltration rate through the space, have negative effects, as well.Both last parameters affect the total heat loss coefficient and in this sense their increase leads to lower thermal transmittance values.The glazing area seems to be effected neutrally.This is due to the fact that increased glazing leads to higher heat losses but at the same time to greater solar gains, which seem to counterbalance each other in this specific example.
The absorptance of the external surface of the wall does not seem to have a significant effect.The same is valid, regarding the initial value of the thermal transmittance of the wall.This last fact is very reasonable, meaning that the optimum total insulation thickness is independent of the initially placed insulation (see Equation (6) and relevant comments in Section 3).
The above results are also presented in the form of ANOVA, in Table 9 (printout of Design Expert ver.10 code).The assumed model (including the parameters B = U-window, D = Heat gains, E = Solar radiation, F = SHGC, H = Infiltration rate, J = Factor-S) scores an F-value of 4279.26 which implies that the model is significant.Values of "Prob > F" less than 0.0500 indicate that model terms are significant.In this case, B, D, E, F, H, J are indeed significant model terms.In addition, the heat gains, the incident solar radiation and the solar gain heat coefficient, have significant positive effects (on the optimum thermal transmittance).An increase of these factors leads to higher thermal transmittance values, or equivalently to lower values of economical insulation.On the other hand, the thermal transmittance of the window in the heated space and the infiltration rate through the space, have negative effects, as well.Both last parameters affect the total heat loss coefficient and in this sense their increase leads to lower thermal transmittance values.The glazing area seems to be effected neutrally.This is due to the fact that increased glazing leads to higher heat losses but at the same time to greater solar gains, which seem to counterbalance each other in this specific example.
The absorptance of the external surface of the wall does not seem to have a significant effect.The same is valid, regarding the initial value of the thermal transmittance of the wall.This last fact is very reasonable, meaning that the optimum total insulation thickness is independent of the initially placed insulation (see Equation (6) and relevant comments in Section 3).
The above results are also presented in the form of ANOVA, in Table 9 (printout of Design Expert ver.10 code).The assumed model (including the parameters B = U-window, D = Heat gains, E = Solar radiation, F = SHGC, H = Infiltration rate, J = Factor-S) scores an F-value of 4279.26 which implies that the model is significant.Values of "Prob > F" less than 0.0500 indicate that model terms are significant.In this case, B, D, E, F, H, J are indeed significant model terms.The resulting linear regression equation (with the parameters set in their coded values, ranging from −1 to +1), is: From this equation it becomes apparent that factor S has almost five times greater an effect than each one of the others.On the other hand, all other parameters together (heat gains, incident solar radiation/orientation), solar heat gain coefficient, thermal transmittance of the window and infiltration rate) may cause an equally significant effect similar to the one of factor S. Indeed, the summation of the absolute values of the respective regression coefficients of parameters B, D, E. F and H equals 29.36 × 10 −3 which is comparable to the coefficient of parameter J 31.70 × 10 −3 .For this reason, the effect of all these secondary parameters should not be neglected, when searching for the optimum insulation thickness.

Evaluating Divergences from the Previously Proposed Analytical Equation for OIT H
In continuation to the above sensitivity analysis, we further calculate for all arising combinations of parameters the respective optimum insulation thickness (OIT H ).These values are now compared to those estimated by the analytical expressions of Equation ( 2), where the actual base temperature is correctly introduced (however as already explained, by applying a typical base temperature like 18 • C the results are misleading).The comparison is depicted in Figure 12.From the figure it becomes apparent that the previously suggested analytical expression (Equation ( 2)) always overestimates the required insulation by almost +12%, as a mean value.There are cases however where this overestimation reaches even 20% leading to a greater thickness by +2 cm and becoming in this way very significant, when the discretization of the available thicknesses for the insulating slabs step 1-2 cm) is taken into consideration.
From this equation it becomes apparent that factor S has almost five times greater an effect than each one of the others.On the other hand, all other parameters together (heat gains, incident solar radiation/orientation), solar heat gain coefficient, thermal transmittance of the window and infiltration rate) may cause an equally significant effect similar to the one of factor S. Indeed, the summation of the absolute values of the respective regression coefficients of parameters B, D, E. F and H equals 29.36 × 10 −3 which is comparable to the coefficient of parameter J 31.70 × 10 −3 .For this reason, the effect of all these secondary parameters should not be neglected, when searching for the optimum insulation thickness.

Evaluating Divergences from the Previously Proposed Analytical Equation for OITH
In continuation to the above sensitivity analysis, we further calculate for all arising combinations of parameters the respective optimum insulation thickness (OITH).These values are now compared to those estimated by the analytical expressions of Equation ( 2), where the actual base temperature is correctly introduced (however as already explained, by applying a typical base temperature like 18 °C the results are misleading).The comparison is depicted in Figure 12.From the figure it becomes apparent that the previously suggested analytical expression (Equation ( 2)) always overestimates the required insulation by almost +12%, as a mean value.There are cases however where this overestimation reaches even 20% leading to a greater thickness by +2 cm and becoming in this way very significant, when the discretization of the available thicknesses for the insulating slabs (thickness step 1-2 cm) is taken into consideration.The ANOVA for the divergences in OITH, is presented in Table 10.The respective regression equation is: The ANOVA for the divergences in OIT H , is presented in Table 10.The respective regression equation is: The initial thermal transmittance is now revealed, as expected, because the higher the initial U-value the greater the required insulation to reach the optimum U* and consequently the greater the divergence between the two alternative estimations.Factor S is again the most important element but now glazing is revealed as an equally important reason for divergences between the two alternative estimations, as also expected (it affects the solar gains).Last, interaction terms are also noticed, but mainly including the most important parameters A and C. The high divergences, reaching up to +20% justify once again the need to consider all secondary parameters when estimating the optimum insulation thickness.Various insulation materials and their effect on OIT H are compared in the frame of the same example.Published data regarding cost and thermal conductivity of insulation materials are used to this aim [33].This data is presented in Table 11 together with the respective results.It becomes apparent that OIT H varies remarkably with the insulation materials, and the same is valid for the corresponding optimum thermal transmittance.The most effective of these alternative materials is extruded polystyrene, by succeeding the lower thermal transmittance (0.37 W/m 2 -K, according to the assumed data).The extraction of a rough guideline on the selection of the most cost-effective insulation material is furthermore attempted.To this aim, the resulted optimum thermal transmittances are correlated with the product C INS •k INS (the specific cost multiplied by the thermal conductivity of the insulation material), as it appears in factor S, Equation (58).Indeed, a strong correlation between optimum U* and C INS •k INS arises, as depicted in Figure 13.Conclusively, the insulation material with the lower product of its specific cost by its thermal conductivity is expected to lead to the lower optimum thermal transmittance, thus making it able to economically decrease the most the heating energy demand, and in this sense it becomes the most cost-effective one.
The extraction of a rough guideline on the selection of the most cost-effective insulation material is furthermore attempted.To this aim, the resulted optimum thermal transmittances are correlated with the product CINS•k INS (the specific cost multiplied by the thermal conductivity of the insulation material), as it appears in factor S, Equation (58).Indeed, a strong correlation between optimum U* and CINS•k INS arises, as depicted in Figure 13.Conclusively, the insulation material with the lower product of its specific cost by its thermal conductivity is expected to lead to the lower optimum thermal transmittance, thus making it able to economically decrease the most the heating energy demand, and in this sense it becomes the most cost-effective one.

Conclusions
The accuracy in using the degree-days method to estimate optimum insulation thickness OIT of an opaque fabric element can be further improved, by taking into consideration the base temperature variation with insulation.In this framework, all factors affecting the above base temperature are also accounted for, such as heat generation within the heated space (including the heat gains utilization factor and the thermal capacity of the space), solar gains through the transparent and opaque surfaces of the space (including the solar absorptance of the considered fabric element) and the total heat losses coefficient of the heated space (including thermal transmittances and areas of the rest of the elements, as well as infiltration).In this sense, optimization does not focus solely on the fabric element but regards it as an interfering part of the whole heated space to which it belongs.In this way a system of two implicit equations arises, which is easily solvable however, with or without using a solver (e.g., either graphically or by applying successive substitutions), to get the optimum thermal transmittance of the considered fabric element and the respective OIT.
The optimum insulation thickness for heating and cooling operation OITY is linearly correlated with the optimum thickness for heating only OITH.The regression coefficient depends on the ratio between cooling and heating degree-days and the relative cost of heating and cooling, respectively.In this context, the optimum insulation thickness for heating application can be firstly estimated followed by the optimum insulation thickness for the whole year operation being concluded, taking additionally into consideration any limitations set by the designer or the standards (e.g., consideration of summer holidays, partial application of cooling etc.).
The optimum insulation thickness (OITH) is estimated by equalizing the marginal insulation cost with the marginal induced benefits.In this sense, OITH is practically the insulation thickness, which leads economically to the maximum achievable energy savings.

Conclusions
The accuracy in using the degree-days method to estimate optimum insulation thickness OIT of an opaque fabric element can be further improved, by taking into consideration the base temperature variation with insulation.In this framework, all factors affecting the above base temperature are also accounted for, such as heat generation within the heated space (including the heat gains utilization factor and the thermal capacity of the space), solar gains through the transparent and opaque surfaces of the space (including the solar absorptance of the considered fabric element) and the total heat losses coefficient of the heated space (including thermal transmittances and areas of the rest of the elements, as well as infiltration).In this sense, optimization does not focus solely on the fabric element but regards it as an interfering part of the whole heated space to which it belongs.In this way a system of two implicit equations arises, which is easily solvable however, with or without using a solver (e.g., either graphically or by applying successive substitutions), to get the optimum thermal transmittance of the considered fabric element and the respective OIT.
The optimum insulation thickness for heating and cooling operation OIT Y is linearly correlated with the optimum thickness for heating only OIT H .The regression coefficient depends on the ratio between cooling and heating degree-days and the relative cost of heating and cooling, respectively.In this context, the optimum insulation thickness for heating application can be firstly estimated followed by the optimum insulation thickness for the whole year operation being concluded, taking additionally into consideration any limitations set by the designer or the standards (e.g., consideration of summer holidays, partial application of cooling etc.).
The optimum insulation thickness (OIT H ) is estimated by equalizing the marginal insulation cost with the marginal induced benefits.In this sense, OIT H is practically the insulation thickness, which leads economically to the maximum achievable energy savings.
Optimum insulation thickness is strongly correlated (positively) to the initial thermal transmittance of the wall; this may vary significantly (e.g., from 0.5 to 2.0 W/m 2 -K), as addition

Figure 1 .
Figure 1.Correlation of consumption to the mean monthly ambient temperature, to demonstrate the effect of insulation and heat gains on the balance temperature.

Figure 2 .
Figure 2. External insulation cost in retrofits, as a function of thickness.

Figure 2 .
Figure 2. External insulation cost in retrofits, as a function of thickness.

Figure 3 .
Figure 3. Variation of insulation costs, energy cost savings and accumulated net benefits with insulation thickness.

Figure 4 .
Figure 4. Effect of cooling degree-days to optimum insulation thickness (use of XPS insulating material and as a fuel the natural gas, data from[12]).

Figure 3 .
Figure 3. Variation of insulation costs, energy cost savings and accumulated net benefits with insulation thickness.

Figure 3 .
Figure 3. Variation of insulation costs, energy cost savings and accumulated net benefits with insulation thickness.

Figure 4 .
Figure 4. Effect of cooling degree-days to optimum insulation thickness (use of XPS insulating material and as a fuel the natural gas, data from[12]).

Figure 4 .
Figure 4. Effect of cooling degree-days to optimum insulation thickness (use of XPS insulating material and as a fuel the natural gas, data from[12]).

Figure 5 .
Figure 5. Gains utilization factor for various time constants, and estimations according to Equation (21) with k G = 1.2).

Figure 7 .
Figure 7. Sketch of the heated space.

Figure 7 .
Figure 7. Sketch of the heated space.

Figure 9 .
Figure 9. Optimum insulation thicknesses for the various orientations, according to the data of the example, and respective base temperatures before and after the insulation.

Figure 9 .
Figure 9. Optimum insulation thicknesses for the various orientations, according to the data of the example, and respective base temperatures before and after the insulation.

Figure 10 .
Figure 10.Revealing the most significant parameters affecting optimum thermal transmittance.

Figure 10 .
Figure 10.Revealing the most significant parameters affecting optimum thermal transmittance.

Figure 11 .
Figure 11.Variation of natural gas prices in Greece, during last three years.

Figure 12 .
Figure 12.Optimum thermal transmittance values estimated in the sensitivity analysis, as compared to the respective values by applying the analytical expression.

Figure 12 .
Figure 12.Optimum thermal transmittance values estimated in the sensitivity analysis, as compared to the respective values by applying the analytical expression.

Figure 13 .
Figure 13.Correlation between the product CINS•k INS and the resulted optimum thermal transmittance U* for the alternative insulation materials.

Figure 13 .
Figure 13.Correlation between the product C INS •k INS and the resulted optimum thermal transmittance U* for the alternative insulation materials.

Table 1 .
Values of Lambert W function.

Table 1 .
Values of Lambert W function.

Table 2 .
Incident solar radiation (kWh/m 2 -mo) in Athens, on vertical surfaces with various orientations.

Table 3 .
Estimation of the duration of the heating period and the respective mean gains.
1referring to the heating period (October to April).

Table 4 .
Data used in the example.

Table 5 .
Calculated quantities in the demonstrative examples.

Table 6 .
Trial and error method to find optimum thermal transmittance U*.

Table 7 .
Application of successive substitutions to find optimum thermal transmittance U*.

Table 8 .
Parameters considered in the sensitivity analysis and respective ranges of variation.

Table 8 .
Parameters considered in the sensitivity analysis and respective ranges of variation.

Table 9 .
Analysis of variance (ANOVA) for selected model for optimum U*.Variation of natural gas prices in Greece, during last three years.

Table 9 .
Analysis of variance (ANOVA) for selected model for optimum U*.

Table 10 .
Analysis of variance (ANOVA) of selected model for Delta(OIT H ).

Table 11 .
Evaluation of various alternative insulation materials.