Multi-Criteria Optimisation of an Experimental Complex of Single-Family Nearly Zero-Energy Buildings

: The Directive 2010 / 31 / EU on the energy performance of buildings has introduced the standard of “nearly zero-energy buildings” (NZEBs). European requirements place the obligation to reduce energy consumption on all European Union Member States, particularly in sectors with signiﬁcant energy consumption indicators. Construction is one such sector, as it is responsible for around 40% of overall energy consumption. Apart from a building’s mass and its material and installation solutions, its energy consumption is also a ﬀ ected by its placement relative to other buildings. A proper urban layout can also lead to a reduction in project development and occupancy costs. The goal of this article is to present a method of optimising single-family house complexes that takes elements such as direct construction costs, construction site organisation, urban layout and occupancy costs into consideration in the context of sustainability. Its authors have analysed di ﬀ erent proposals of the placement of 40 NZEBs relative to each other and have carried out a multi-criteria analysis of the complex, determining optimal solutions that are compliant with the precepts of sustainability. The results indicated that the layout composed of semi-detached houses scored the highest among the proposed layouts under the parameter weights set by the developer. This layout also scored the highest when parameter weights were uniformly distributed during a test simulation.


Introduction
Designing and building a complex of single-family houses, with a holistic focus on sustainable development, and construction and occupancy costs, is a problem that is very difficult to solve. It is a multi-aspect problem, in which decision variables include, among others: selecting construction materials, building services, the parameters of the shape of the buildings, their placement on the site, the sequence in which individual buildings are to be constructed and the completion deadlines for individual tasks. During the occupancy stage, one should add an analysis of the impact of the change in energy costs. When considering the entire life-cycle of the buildings, renovation and repair costs should also be included, in addition to any remodelling and demolition expenses. In recent years, the focus on the energy efficiency and pollution emission of buildings has grown exponentially and has had an immense impact on the real estate market and property values. The paper [1] presents an analysis of the profitability of investing in a selection of three systems based on renewable energy sources, using a semi-detached house as an example.

The Scale of the Problem
According to the Eurostat Agency [30], in 2017, as much as 33.6% of the population of the European Union resided in detached single-family houses, while 24.0% of the population lived in semi-detached single-family houses. According to data provided by the Polish General Construction Inspector's Office [31], in 2018, a total of 98,915 building permit applications were filed for residential buildings, out of which 93,714 building permits were issued for single-family houses (almost 94.8%), which constituted a 3% increase from the previous year [31]. The results of statistical studies presented by BPD Europe BV [32] indicate that 12% of respondents from The Netherlands, 16% of respondents from France and 17% of respondents from Germany who had been searching for a new place of residence pointed to a single-family house as a preferred option in a situation in which they would have to move to a different dwelling. This proves that single-family buildings remain popular in many developed countries, while in countries such as Poland they continue to enjoy very high popularity. The design of single-family building complexes can therefore be considered an important field, both to planners who design urban layouts for individual projects, as well as designers who work on projects composed of many such buildings.
The necessity of considering parameters associated with energy efficiency is dictated by, among other things, the current energy policy vector of the European Union, expressed by such documents as Directive 2010/31/EU on the energy performance of buildings [33]. This EU document introduced the following definition of nearly-zero energy buildings (NZEB): nearly zero-energy building' means a building that has a very high energy performance, as determined in accordance with Annex I. The nearly zero or very low amount of energy required should be covered to a very significant extent by energy from renewable sources, including energy from renewable sources produced on-site or Energies 2020, 13, 1541 4 of 30 nearby". However, every European Union Member State has adopted its own NZEB parameters, which are adapted to local conditions. The implementation of the Directive into Polish state law has introduced a new standard of near zero-energy buildings (NZEBs). The document that has been amended to include the standard, as stipulated in the Directive, is the Ordinance on the Technical Conditions that Must be Met by Buildings and their Placement [34]. In the Polish definition of NZEBs there are two criteria. The first criterion is ensuring proper building envelope thermal insulation. The second criterion for NZEB buildings in Poland applies to Primary Energy (PE), which cannot exceed a set value for different types of buildings. For single-family residential buildings without cooling installations, the maximum PE value is 70 kWh/(m 2 /annum). The buildings analysed in the article meet both criteria, thus being compliant with NZEB standards as stipulated in Polish regulations.
Selecting the appropriate orientation of buildings, their shape, construction materials, building services and shading system solutions leads to lower energy demand and costs and affects the comfort of use of the buildings themselves, reducing overheating. This subject has been discussed in, among others, [35][36][37][38][39][40][41][42]. It can therefore be argued that the application of research findings in the field of energy efficiency and thermal comfort optimisation in the design of buildings and complexes that enjoy significant popularity on the real estate market should be incorporated into current design methodology. It can also be argued that by increasing the energy efficiency of buildings forming market-preferred layouts one can contribute to increasing the popularity of such solutions in areas where detached or semi-detached houses are the most popular.

Cost Estimation Methods
Calculating the overall costs of construction is a difficult task. Typically, the overall costs of construction are divided into indirect and direct costs. However, this is a major oversimplification. Situations that can significantly affect the expenditures of the contractor can take place over the course of working on a construction project. Some of these are difficult to predict (e.g., natural disasters, economic and social changes), while others can be prepared for during a project's planning stage. The predicted additional costs can include penalties for failing to meet planned deadlines or penalties for failing to maintain the continuity of construction crew operation.
Furthermore, every project is subjected to certain technological and organisational constraints, such as the time during which construction crews or buildings are available, or technological and organisational requirements concerning ongoing work performed by a given crew on a specific building. Numerous cost estimation or total construction project cost optimisation methods have been developed. They usually utilise metaheuristic methods, such as simulated annealing, hybrid algorithms [43] and genetic algorithms [44]. Mathematical programming methods are also used, such as the Decision model for planning material supply channels in construction [45]. In [46], an innovative approach was proposed, in the form of a stochastic decision network. It enables the integration of random and decision-based aspects of planned projects. This approach allows the decision-maker to obtain a solution that is optimal in terms of the expected time and cost. They can do so by defining their preferences as to the results of the planned project and setting their risk aversion accordingly. In the case of multi-criteria analysis of time and cost, the decision-maker can also define weight values for the expected project time and cost. The methods listed above do not incorporate all of the essential costs of carrying out construction projects, such as direct costs, indirect costs, penalties for failing to meet contractual deadlines, as well as the costs of discontinuities in the work of construction crews. They also do not take into account technological and organisational constraints. Furthermore, metaheuristic methods involve lengthy calculations when compared to linear programming.

Previously Used Methods
Multi-criteria decision-making analysis methods that have been applied in urban design and spatial planning have been predominantly focused on large-scale planning. As Afshari and Vatanparast [6] Energies 2020, 13, 1541 5 of 30 have indicated, such proposals utilised the Preference Ranking Organisational method for Enrichment Evaluation (PROMETHEE and PROMETHEE-2) methods in combination with Hasse Diagram Technique (HDT), the AHP method, Fuzzy AHP, Weighted Linear Combination (WLC), a combination of the AHP and TOPSIS methods, the ELECTRE-2 method (ELimination Et Choix Traduisant la REalité), the EXPROM-2 method (an extended version of the PROMETHEE method), as well as evolutionary algorithms. Guarini, Battisti and Chiovitti [23] also pointed to the possible use of Multi-attribute Utility Theory (MAUT) and Measuring Attractiveness by a Categorical Based Evaluation methods (MACBETH).
The selection of the multi-criteria analysis method depends on the specificity of the decision problem at hand. For instance, in [47] the author referred to studies in which they used the ELECTRE III and AHP methods to analyse the possibility of using an alternative version of the BIPOLAR method to solve the problem of selecting a construction company for partnering cooperation during a construction project. The study included a calculation example and verified the results of previous studies by applying an alternative BIPOLAR method version. Ultimately, its use produced the same end result as when using two other methods.
Multi-criteria methods that incorporate co-dependencies between decision-making criteria and alternatives have seen considerable development in recent years. In [48] the authors analysed a multi-criteria problem of selecting a new function for historical buildings in light of the proposed selection criteria. Apart from the economic aspect and cultural heritage, these criteria also included sustainability. The Weighted Influence Non-linear Gauge System method (WINGS), which was extended to incorporate the uncertainty of expert opinions and their aggregation, was applied to model the structure of the problem and analyse the dependencies between alternatives and criteria. The method had to be extended because of the specificity of the problem discussed in the paper. The majority of previously used methods did not analyse interdependencies between decision-making criteria despite such interdependencies actually existing, as proven by the authors. Similarly, the authors of [49] proposed a multi-criteria hybrid model, using the Decision Making Trial and Evaluation Laboratory method (DEMATEL) and the Analytic Network Process (ANP) to select a utility function for the purpose of adapting the building of 'Stara Polana', located in Zakopane. Both in this publication and in [50], a group of experts and specialists used a questionnaire to assess dependencies between factors that affect the criterion of societal benefits and the benefits associated with preserving cultural heritage. These were some of the parameters selected for analysis and rating in association with the investigated alternative forms of adaptive reuse of the Stara Polana heritage site in Zakopane. The remaining parameters, associated with, among others, acoustics, vibrations or thermal comfort, were studied using specialist technical instruments. These interdisciplinary studies became a basis for a multi-criteria analysis performed using a proposed hybrid model of selecting a form of use for a heritage site adaptation project.

The Significance of the Method Proposed in the Article
The innovative method proposed in this paper is an important contribution to the professional toolkit of designers and can significantly aid real estate developers in tailoring their offering. The method presented by the authors includes elements and criteria that are important during many stages of the design process: the conceptual urban design and site plan stage, the conceptual and technical architectural stage, as well as the structural system design and building services design stages. It also pertains to project cost assessment in the context of building single-family residential NZEBs, which are an important part of the housing sector. This approach can be considered a holistic perspective of the design task and the project as a whole, which is both necessary and desirable from the point of view of sustainable design. The necessity of this approach has been confirmed by, among others, Garau and Pavan [51] in a broader context. The method developed by authors assumes a multi-criteria analysis that includes aspects of sustainability, as well as the costs of the construction and occupancy of a complex of single-family houses. Sustainability parameters, such as the energy demand of the buildings and greenhouse gas Energies 2020, 13, 1541 6 of 30 emissions (depending on the given building variants and their siting) have been determined through analyses using the EnergyPlus tool [52]. The energy analysis took the life-cycle stages of the buildings into consideration. The cost parameter was determined by including direct and indirect costs, the costs of construction (taking technological and organisational constraints into consideration), as well as the costs of building occupancy-the authors performed calculations for different building types and urban layout alternatives. Indirect and direct costs, as well as construction costs, were determined using a specifically-tailored optimisation model and the concept of priority scheduling. Occupancy costs were calculated using energy demand and the discount rate. Using previously selected parameters as a basis, the authors carried out a multi-criteria analysis using the Weighted Aggregated Sum Product Assessment method (WASPAS). This made it possible to determine the optimal solution for the problem in question.

Method Overview
The design optimisation task analysed by the authors is an actual design problem. The analysed complex of forty experimental buildings was, at the time of the writing of thepaper, planned to be constructed in the municipality of Mogilany near Krakow, Poland. The developer must select one of several urban layout alternatives, which is to be optimal in terms of cost and energy demand minimisation, in addition to being compliant with the precepts of sustainable development and real estate market preferences. The urban layouts under analysis can differ in terms of building siting, as well as the size or type of buildings (detached, semi-detached, terraced). Criteria values were determined for these layouts. The criteria adopted for the analysis included: overall construction cost, building energy demand, greenhouse gas emissions associated with building occupancy, building thermal comfort and an assessment of the urban layout by potential buyers. After setting parameter values for all alternatives, the decision-maker (developer) can set the weights of each parameter. After preparing the data in the listed manner, the authors conducted a multi-criteria optimisation using the WASPAS method. A block scheme of the proposed method will be presented in Figure 1. The method of setting the values of individual criteria will be presented in the following sections of the article. The method developed by authors assumes a multi-criteria analysis that includes aspects of sustainability, as well as the costs of the construction and occupancy of a complex of single-family houses. Sustainability parameters, such as the energy demand of the buildings and greenhouse gas emissions (depending on the given building variants and their siting) have been determined through analyses using the EnergyPlus tool [52]. The energy analysis took the life-cycle stages of the buildings into consideration. The cost parameter was determined by including direct and indirect costs, the costs of construction (taking technological and organisational constraints into consideration), as well as the costs of building occupancy-the authors performed calculations for different building types and urban layout alternatives. Indirect and direct costs, as well as construction costs, were determined using a specifically-tailored optimisation model and the concept of priority scheduling. Occupancy costs were calculated using energy demand and the discount rate. Using previously selected parameters as a basis, the authors carried out a multi-criteria analysis using the Weighted Aggregated Sum Product Assessment method (WASPAS). This made it possible to determine the optimal solution for the problem in question.

Method Overview
The design optimisation task analysed by the authors is an actual design problem. The analysed complex of forty experimental buildings was, at the time of the writing of thepaper, planned to be

Overall Construction Cost
Despite the contribution of the aforementioned methods, the approach prepared by the authors of this paper features numerous advantages. Its greatest advantage is that it describes all of the included cost dependencies (direct costs, indirect costs, penalties for failing to meet contractual deadlines, penalties for failing to ensure continuity in the work of construction crews) and constraints in the form of a linear programming model. This allows one to quickly find the optimal solution (e.g., using the Simplex method). If we allow the possibility of altering the sequence of building the structures, the rapid calculation of the goal function in the discrete optimisation algorithm can be considered a significant advantage. The proposed model permits the overlap of the work of construction crews on different plots, which is a good reflection of the working conditions of an actual construction site. Furthermore, the model allows for flexible adherence to technological and organisational constraints. This means that the model permits certain departures from the constraints defined by the designer. However, these departures are optimal (minimal). This is of particular significance in the case of projects with a large (even exceedingly large) number of constraints, as one can encounter a situation where it is impossible to meet all of them. The model's disadvantage is that it describes time-cost dependencies in linear form, which is a considerable simplification. Furthermore, the sequence of building the structures cannot be changed in the model. Both of these disadvantages will be studied further by the authors. The individual sections of the model have been discussed in [53,54].
The model has been defined as follows: Given: where: with decision variables: Goal function: where: K indir = Z n,m ·k indir (7) Energies 2020, 13, 1541 8 of 30 Constraints: In order for the presented model to be applied, the following data is required(1)-(3): • tgr i,j -limit completion time for task performed on building i by crew j; • kgr i,j -limit costs of task performed on building i by crew j; The parameters are additional given data in the model, enabling the establishment of technological and organisational constraints. They have been described in the article at length [52]. The model's decision variables (4) are as follows: • t i,j -task duration (i,j); • R i,j -task initiation time (i,j); • Z i,j . -task completion time (i,j); • p i -delay value relative to the planned deadline for building i; • in addition to variables that make it possible to establish technological and organisational constraints: CK OL i,j , CK OU i,j , CK BL i,j , CK BU i,j . The goal function (5) is the sum of all essential project costs. Direct cost K dir (6) is the sum of the costs of all works performed by all crews on all buildings, depending on the completion time of a given Energies 2020, 13, 1541 9 of 30 task. Direct cost includes the cost of labour, materials and equipment. It is assumed that a reduction in the task completion time, from a standard completion time to a limit time, is possible, but requires additional expenses. The dependency between the cost of completing a given task depends on the set standard and limit costs and is assumed to be linear. Therefore, the popular Critical Path Method Cost (CPM-COST) approach was incorporated into the model [55].The CPM-COST method enables the calculation of minimum direct costs, assuming the ability to shorten task completion time by incurring additional expenses when a fixed contractual deadline is imposed. Element (6) of the goal function and constraints (12)- (18) apply the assumptions of the CPM-COST method, extended to include flexible time couplings. Had the model included only (6) and (12)-(18), with parameters CT and CK removed, with the addition of the condition of meeting a contractual deadline, it would have been a classical CPM-COST model. The values a i,j and b i,j can be calculated using equations (25) and (26).
Indirect cost K indir (7) depends on the value of indirect per unit cost or the completion time of the entire multiple-structure project. Indirect per unit cost can be assessed using either a detailed or indicator-based method.
Costs associated with failing to meet planned deadlines K p (8) are set as the sum of the products of the per unit cost of failing to meet a deadline and the length of the delay for all buildings. the per unit cost of failing to meet a planned deadline can vary between buildings. The penalties for such are defined in the contract between the developer and the buyer and typically range between 0.05 and 1% of the contract value for every day of delay. For instance, in [56] the authors define the cost of failing to meet contractual deadlines as ca. 0.1% of contract value per day of delay. Additional conditions depend on the country and its legal environment.
Costs associated with pauses in the operation of crews K c (9) are calculated as the sum of the products of the per unit cost of pauses in crew operation and the length of the pause for all crews. The pause length is calculated on the basis of the earliest time the crew can complete its tasks on all buildings, the earliest time of initiation of work on the first building and the task completion times of the crew on all buildings.
In addition, the goal function includes the costs associated with failing to meet flexible couplings between tasks performed on successive buildings performed by the same crew K CO (10) as well as flexible time couplings between tasks performed by successive crews on each building K CB (11). The impact value of these couplings on the goal function for both types of couplings is calculated as the sum of the product of failing to meet the lower value of the flexible time coupling and the unit value for coupling failure, as well as the product of failing to meet the upper value of the flexible time coupling and the per unit value for coupling failure.
The following constraints have been included in the model. Constraints (12) and (13) constrain the task completion time to the limits between limit time and normal time. Constraint (14) enables the calculation of task completion times. Constraints (15)-(18) make it possible to determine the value of failing to meet appropriate flexible time couplings. Constraint (19) makes it possible to calculate the value of failing to meet the deadline of completing a given building. Constraints (20)-(23) implement the assumptions concerning the availability of every building and every crew. All decision variables used in the model are to take on non-negative values (24).
The goal function and all constraints are linear, which makes this a linear programming model. The Simplex method was used to solve the model [57].
The model described above is used to optimise a pre-set sequence of constructing buildings. In general, this sequence can be modifiable. Changing the building completion sequence affects the time and cost of completing the project. The authors plan to expand the model so that it will be able to Energies 2020, 13, 1541 10 of 30 optimise the sequence of constructing buildings. The proposed model was implemented using the Python programming language. The program can be obtained from the authors.

Energy Consumption
One of the criteria of assessing near-zero energy buildings (NZEBs) that are applicable in Poland is Primary Energy (PE), which cannot exceed the values stipulated in [34], and which differ depending on building type. For single-family residential buildings without cooling installations, the maximum PE value is 70 kWh/(m 2 /annum). The Primary Energy (PE) values for the analysed buildings are lower than the maximum values allowed by Polish law. Primary Energy values predominantly refer to the building's sources of heating and electric power. Final Energy (FE) is a more precise indicator of material and building services solutions, as it results from a building's overall energy performance balance. This is why this indicator was incorporated into the multi-criteria analysis presented in the article.
The energy consumption simulation for all of the analysed alternatives was performed using DesignBuilder, a program enabling dynamic simulations. It is based on an hourly calculation method. The calculation engine of the program is based on EnergyPlus. DesignBuilder uses EnergyPlus format hourly weather data to define external conditions during simulations. These hourly weather data sets are data derived from hourly observations at the chosen location-in this case, Krakow, Poland.

CO 2 Emissions
A building consumes various natural resources, including water, materials and energy, and releases many pollutants and emissions during its life-cycle, i.e., from the raw material extraction to the final treatment of waste from the building's demolition [58]. The authors assume that emissions generated by construction material production, the construction of the buildings or their demolition, will be highly similar for both building variants (detached, semi-detached). This is why only occupancy-related CO 2 emissions were incorporated into the analysis, particularly as it is concordant with building energy performance calculations stipulated in Polish law [34].
The CO 2 emission in EnergyPlus is calculated from Fuel Totals using a locally published CO 2 emission factor for the specific region-in this case, Krakow. The consumed fuel (electricity for lighting and gas for DHW and Heating) is multiplied by the carbon emission factor. Based on emissions factors entered for a specific location, EnergyPlus calculates the mass or volume of CO 2 (carbon dioxide) [52] The amount of CO 2 produced depends on the carbon content of the fuel, and the amount of heat produced depends on its carbon and hydrogen content. Because natural gas, which is mostly methane (CH 4 ), has a high hydrogen content, combustion of natural gas produces less CO 2 for the same amount of heat produced from burning other fossil fuels [59].

Thermal Comfort
The methodology for determining thermal comfort is based on PN-EN ISO 7730 [60]. The thermal insulation of clothing was determined based on the PN-EN ISO 9920:2009 standard [61].
According to the ANSI/ASHRAE Standard 55-2010, Thermal comfort analyses were processed for indoor winter clothing (I cl = 1 clo) for a typical winter week period, and summer clothing (I cl = 0.5 clo) for a typical summer period. The l clo value of 0.7 is assumed for transitional periods over the year. The insulation of clothing was determined as the value for the transitional season of clothing worn indoors I clo = 0.7 (clo). Determining the PMV (Predicted Mean Vote) indicator requires the calculation of the physical parameters of the analysed space, while considering energy expenditure and the thermal insulation properties of the clothes of the persons working inside. The calculated PMV value is compared with a seven-point psychophysical scale of thermal sensation which has been developed by Fanger [62] (Table 1). A comfortable environment in terms of microclimate (thermally neutral) can be considered to be within the 0.5 < PMV < +0.5 range [63]. Thermal comfort sensations on the Fanger scale have been presented in Table 1. Thermal comfort evaluation analysis was performed for two periods: • A typical summer week-a week identified by the weather data translator software as being typical of summer. • A typical winter week-a week identified by the weather data translator software as being typical of winter.
The authors assumed that these periods will be suitable for the purpose of this paper as they represent weather conditions that are repeatable and most likely to happen every year. The Polish thermal comfort standard [60] does not indicate specific periods of the year for measurements.
Design Builder uses three different mathematical models to calculate Thermal Comfort: the one by Fanger (the Fanger Comfort Model), the Pierce Foundation model (the Pierce Two-Node Model), and the model developed by researchers from Kansas State University (the KSU Two-Node Model). For the purpose of this paper, the Fanger model was chosen and implemented for thermal comfort analysis based on the equation from ISO 7730. All PMV values are taken from dynamic simulation outputs.
The results of the thermal comfort simulation analysis will be compared with in-situ measurements when similar conditions occur once all of the buildings have been completed. One of the hottest days of the year was selected for the comparative analysis, i.e., 12.07.2018. The simulation was performed for the living rooms of each building.

Urban Layout
The selection of the urban layout of a housing complex must take the preferences of future residents into consideration. In order to determine the optimal layout, a survey was conducted among the potential users of the housing complex's buildings. The respondents were presented with three urban layouts (Layout A, Layout B and Layout C) and asked to select the optimal urban layout in terms of housing comfort, their preferences and convenience. The survey was filled out by 31 middle-aged persons (30-50 years of age). The survey featured a three-point scale: a rating of 1; a rating of 2; a rating of 3. A rating 1 was given to the layout considered to be the most desirable by the respondent, a rating of 2 was given to the next-to-best layout and a rating of 3 was given to the least desirable of the layouts from the point of view the adopted assessment criteria.

MCDM Analysis-The WASPAS Method
The parameter sets selected by the authors (construction cost, energy demand, CO 2 emission, thermal comfort, urban layout assessment) are diverse and cannot be directly compared. This is why the authors decided to utilise a multi-criteria decision-making approach. There are numerous MCDM methods, such as AHP, ANP, ELECTRE, BIPOLAR, TOPSIS and many others [64,65]. For this specific problem, the authors selected the WASPAS method. It is a method that combines two other methods: the Weighted Sum Model (WSM) and the Weighted Product Model methods (WPM) [66][67][68]. The WASPAS method was selected because of its mathematical simplicity and capacity to provide highly precise results when compared to the WSM and WPM methods. The WASPAS method is currently widely accepted as an effective decision-making support tool [68]. It is continuously being developed (e.g., by integrating a fuzzy approach) [69,70]. The WASPAS method requires a decision matrix to be prepared: X = x i,j m×n where x i,j is the value of alternative i in relation to criterion j, m is the number of alternatives and n is the number of criteria. Afterwards, matrix X is to be standardised according to the following formulae: for criteria which are stimulants In order to calculate the importance of alternative i, the following formula is to be used: The optimal value of parameter λ can be calculated based on the formula below: can be calculated from the following formulae: The approximate variance value for the standardised criteria is calculated using the formula: The criteria adopted for the analysis using the WASPAS method have been presented in Figure 2. Of note is the weight selection process in the WASPAS method. The proposed method has been formulated with the developer in mind and it is they who make autonomous decisions concerning weight values, based on their own preferences. The method is a decision-making support tool for use by the developer and models their preferences. In specific cases, the developer can use a different method of setting weights, e.g., using experts' opinions. The WASPAS method was applied by the authors using a Microsoft Excel spreadsheet. It is a simple solution that can be helpful to the developer. All that needs to be done is to set the weights and individual numerical values for parameters. The spreadsheet will calculate a score for each layout and produce a ranking of said layouts. The Excel file can be obtained from the authors.  In specific cases, the developer can use a different method of setting weights, e.g., using experts' opinions. The WASPAS method was applied by the authors using a Microsoft Excel spreadsheet. It is a simple solution that can be helpful to the developer. All that needs to be done is to set the weights and individual numerical values for parameters. The spreadsheet will calculate a score for each layout and produce a ranking of said layouts. The Excel file can be obtained from the authors.

Case Study-Housing Complex in Mogilany near Krakow, Poland
The planned housing complex which is the subject of the case study in this article is planned to be built on a 3-hectare plot in Libertów, in the municipality of Mogilany near Krakow, in the Lesser Poland Voivodship (Poland). The site is adjacent to an expressway which leads in the direction of Zakopane, a popular tourist destination located in the mountains. There is a lot of traffic along this road, particularly during weekends and holidays. The location of the experimental housing complex has been shown in Figure 3. In specific cases, the developer can use a different method of setting weights, e.g., using experts' opinions. The WASPAS method was applied by the authors using a Microsoft Excel spreadsheet. It is a simple solution that can be helpful to the developer. All that needs to be done is to set the weights and individual numerical values for parameters. The spreadsheet will calculate a score for each layout and produce a ranking of said layouts. The Excel file can be obtained from the authors.

Case Study-Housing Complex in Mogilany near Krakow, Poland
The planned housing complex which is the subject of the case study in this article is planned to be built on a 3-hectare plot in Libertów, in the municipality of Mogilany near Krakow, in the Lesser Poland Voivodship (Poland). The site is adjacent to an expressway which leads in the direction of Zakopane, a popular tourist destination located in the mountains. There is a lot of traffic along this road, particularly during weekends and holidays. The location of the experimental housing complex has been shown in Figure 3. The site plan for the plot (MU1) has been shown in Figure 3b. Figure 3c shows an aerial view of the plot. The housing complex was assumed to be composed of 40 single-family residential buildings, with various building layouts (detached and/or semi-detached buildings).
Three proposals of the urban layout of the complex's buildings were analysed (Layout A, B and C). They will be presented in Figure 4. The expressway runs near the plot, to its south-east (SE). To the north of the plot is the direction towards Krakow, while from the south, to the holiday resort-Zakopane. The urban layouts that were analysed are a reflection of the preferences of persons who buy single-family houses in Poland in terms of their view of property desirability. The site plan for the plot (MU1) has been shown in Figure 3b. Figure 3c shows an aerial view of the plot. The housing complex was assumed to be composed of 40 single-family residential buildings, with various building layouts (detached and/or semi-detached buildings).
Three proposals of the urban layout of the complex's buildings were analysed (Layout A, B and C). They will be presented in Figure 4. The expressway runs near the plot, to its south-east (SE). To the north of the plot is the direction towards Krakow, while from the south, to the holiday resort-Zakopane. The urban layouts that were analysed are a reflection of the preferences of persons who buy single-family houses in Poland in terms of their view of property desirability. Their selection is the result of the actual housing market situation in Poland. The proposed urban layouts are options that the developer offers for sale to their clients on the open market. The authors are aware of the precepts of sustainable design concerning the planning of small-scale urban complexes, but in this case the intent is to provide an iteration of the method for a more market-oriented approach. The authors do not rule out preparing a modification of the method that would feature other sustainability metrics in the future.  The site plan for the plot (MU1) has been shown in Figure 3b. Figure 3c shows an aerial view of the plot. The housing complex was assumed to be composed of 40 single-family residential buildings, with various building layouts (detached and/or semi-detached buildings).
Three proposals of the urban layout of the complex's buildings were analysed (Layout A, B and C). They will be presented in Figure 4. The expressway runs near the plot, to its south-east (SE). To the north of the plot is the direction towards Krakow, while from the south, to the holiday resort-Zakopane. The urban layouts that were analysed are a reflection of the preferences of persons who buy single-family houses in Poland in terms of their view of property desirability. Their selection is the result of the actual housing market situation in Poland. The proposed urban layouts are options that the developer offers for sale to their clients on the open market. The authors are aware of the precepts of sustainable design concerning the planning of small-scale urban complexes, but in this case the intent is to provide an iteration of the method for a more market-oriented approach. The authors do not rule out preparing a modification of the method that would feature other sustainability metrics in the future.     The housing complex will be composed of buildings with different floor areas and layouts. The types of individual buildings will be presented in Figure 5.  Table 2 features a listing of the number of detached and semi-detached buildings for each layout alternative.
The housing complex will be composed of buildings with different floor areas and layouts. The types of individual buildings will be presented in Figure 5.   The housing complex will be composed of buildings with different floor areas and layouts. The types of individual buildings will be presented in Figure 5.  The buildings under analysis were assumed to be built using timber framing technology with thermal insulation composed of wood wool. The ground floor plan, assumed to have a usable floor area of 110 m 2 , in a detached building, is to be presented in Figure 5a, and in a semi-detached building-in Figure 5b. The geometric models of the buildings subjected to the analysis have been presented in Table 3. B1 type buildings are linked with a common wall that separates the bathrooms of both buildings.
The cost of individual buildings was calculated based on "Biuletyn cen obiektów budowlanych (BCO) część I a-obiekty kubaturowe" (in English: Construction price bulletin, Part Ia-buildings) [71]. It is a cost bulletin that lists the mean prices of buildings in Poland. The BCO can be used to estimate the cost of constructing buildings or their fragments for the purposes of cost assessment or preparing financial schedules for construction projects. Direct cost was assessed using the construction price bulletin and includes labour, construction materials and machinery. The completion time for individual works was estimated on the basis of available material costs catalogues. The values for completion times and costs were assessed for both types of buildings featured in the case under analysis. Table 4 includes the values for a detached building with a usable floor area of 110 m 2 . The buildings under analysis were assumed to be built using timber framing technology with thermal insulation composed of wood wool. The ground floor plan, assumed to have a usable floor area of 110 m 2 , in a detached building, is to be presented in Figure 5a, and in a semi-detached building-in Figure 5b. The geometric models of the buildings subjected to the analysis have been presented in Table 3. B1 type buildings are linked with a common wall that separates the bathrooms of both buildings. The cost of individual buildings was calculated based on "Biuletyn cen obiektów budowlanych (BCO) część I a-obiekty kubaturowe" (in English: Construction price bulletin, Part Ia-buildings) [71]. It is a cost bulletin that lists the mean prices of buildings in Poland. The BCO can be used to estimate the cost of constructing buildings or their fragments for the purposes of cost assessment or preparing financial schedules for construction projects. Direct cost was assessed using the construction price bulletin and includes labour, construction materials and machinery. The completion time for individual works was estimated on the basis of available material costs catalogues. The values for completion times and costs were assessed for both types of buildings featured in the case under analysis. Table 4 includes the values for a detached building with a usable floor area of 110 m 2 .
In order to determine the overall construction costs, the following parameter values were assumed: • indirect construction costs, estimated at 400 €/day; Indirect costs include the general costs of construction (wages for the site director, works managers, construction site infrastructure costs), as well as management costs (the cost of maintaining the project, company offices, costs of office work, etc.). Indirect costs were assessed as per the construction price bulletin. The bulletin lists the value of indirect costs at a level of ca. 15% of the total cost. The cost per day was calculated using an estimate of the project's duration. • costs of failing to meet planned deadlines, estimated at 150 €/per day of delay; Penalties for failing to meet planned deadlines are defined in contracts between the developer and the buyer In order to determine the overall construction costs, the following parameter values were assumed: • indirect construction costs, estimated at 400 €/day; Indirect costs include the general costs of construction (wages for the site director, works managers, construction site infrastructure costs), as well as management costs (the cost of maintaining the project, company offices, costs of office work, etc.). Indirect costs were assessed as per the construction price bulletin. The bulletin lists the value of indirect costs at a level of ca. 15% of the total cost. The cost per day was calculated using an estimate of the project's duration. • costs of failing to meet planned deadlines, estimated at 150 €/per day of delay; Penalties for failing to meet planned deadlines are defined in contracts between the developer and the buyer of the building. There is no official data on this subject and in the experience of the authors such penalties in Poland can range from 0.05% to even 1% of a contract's value per day of delay. It should be noted that penalties amounting to 0.5% of a contract's value per day are not considered unjustified by Polish courts of law. Furthermore, it should be noted that the penalty is not applied when the developer cannot be held responsible for the cause of the delay. It is impossible to factor in unforeseen events that can cause delays and for which the developer is not responsible during the construction planning stage, which is why this aspect was ignored in the paper. A penalty value of 0.2% per day per individual building was assumed in the paper. • it was assumed that agreements with some contractors specified penalties for failing to ensure the continuity of crew operation. For the third crew, a 150 € penalty for each day of pause in the work was assumed, with 200 € per day assumed for the fourth crew. The penalty value for failing to ensure crew work continuity is a direct result of contracts with specific contractors. If a contractor performs a highly specialised section of the work, these penalties will be higher. This is why the third crew, which was assumed to assemble the innovative roof structure, and the fourth crew, which was assumed to handle building services installation, were assigned such high penalties for work discontinuity. • it was determined that the first building is to be completed within 90 days. Every subsequent building was to be completed 14 days after the completion of the previous one. Planned deadlines are a direct consequence of contracts or agreements between potential buyers and the developer.
In addition, it was assumed that technological continuity was to be maintained between footing, foundation wall construction and insulation application, as well as the construction of load-bearing walls, decks, roof structure and roofing. Table 4 includes the standard and limit costs of construction work, as well as the standard and limit completion times for individual works.
The optimisation process produced schedules for each of the analysed layouts that were optimal in terms of cost. For layout A, an overall cost of 3034 thousand euro was obtained, along with an overall completion time of 629 days. For layout B, the costs amounted to 2928 thousand euro, while the overall completion time was 629 days. For layout C, the costs amounted to 3233 thousand euro and the overall completion time was 680 days. The completion time was taken into consideration when calculating the costs of failing to meet planned deadlines for each building.
The buildings under analysis were designed to have a timber frame structure. The external partitions of the analysed buildings were assumed to be built using timber framing technology. The technology selection was primarily dictated by the developer's guidelines, who had the technical infrastructure enabling them to manufacture prefabricated timber elements. Figure 6 is to show a cross-section of the external wall.
Energies 2020, 13, x FOR PEER REVIEW 18 of 30 infrastructure enabling them to manufacture prefabricated timber elements. Figure 6 is to show a cross-section of the external wall. U heat transfer coefficients [W/(m 2 K)] for elements of the buildings' envelope were calculated in accordance with the methodology outlined in [72,73]. The U heat transfer coefficients have been listed in Table 5, along with the standard requirements for NZEBs applicable in Poland [34].    [72,73]. The U heat transfer coefficients have been listed in Table 5, along with the standard requirements for NZEBs applicable in Poland [34]. The windows of the buildings meet Polish requirements for NZEBs. They were partially fitted with external blinds. The geometry of the buildings was also subjected to optimisation. Their massing is simple and thermal bridges were minimised. Their orientation relative to the cardinal directions is determined by local development guidelines. The structure of the buildings, material choices and their mass were determined by the developer. The air-tightness of the external envelope of the buildings was assumed to be n 50 = 1.5 [1/h]. In Polish regulations concerning building design [34] there are two requirements concerning building air-tightness. For buildings with gravitational ventilation it is a value of 3.5 air exchanges per hour [1/h], while for buildings with mechanical ventilation, the number of air exchanges per hour is 1.5 [1/h]. This is why the authors assumed a value of 1.5 [1/h]. The analysis was performed during the buildings' design stage. At present, the first buildings have already been built. Air-tightness inspections performed on site indicate that the precise construction of the buildings has ensured a much greater air-tightness than initially assumed during the design stage. The air-tightness achieved during construction has a value of n 50 = 0.6 [1/h], and is at a level expected of passive buildings [34]. The internal temperatures were assumed in accordance with the PN-EN 12831 standard [74].
Occupancy-related (metabolic) gains were defined in accordance with the ASHRAE 140/BESTEST standard [75], depending on the form of use of the spaces.
The buildings were assumed to be equipped with mechanical supply and exhaust ventilation with heat recovery. Supply air streams were assumed as in [76].
The Central Heating installation was designed in the form of gas-powered heating using a dual purpose gas-powered condensation boiler with a total efficiency rating of 0.95.
The indoor air temperature was assumed per PN-EN 12831:2006. Instalacje ogrzewcze w budynkach-Metoda obliczania projektowego obciążenia cieplnego-a Polish standard for heating installations governing methods of calculating heat loads during the design stage [74]. Per this standard, spaces for permanent occupancy by persons without outdoor clothing should be designed to have a temperature of +20 • C, while spaces meant for occupancy by persons potentially without clothing (such as bathrooms), should have a temperature of +24 • C. This temperature difference can be obtained by properly setting a building's thermostats or by using BMS installations. Concerning the schedule, it was assumed that such temperatures should be kept 24 h per day throughout the entire year. In residential buildings with building automatics systems it is much easier to assume periods of lowered temperature (e.g., at night) or when residents leave the building to work. However, most buildings in Poland do not feature elaborate BMS installations. The analyses assumed the buildings to have such systems, which is why no temperature drops were accounted for.
Design Builder Software uses Occupancy model data which defines the number of people in a space and occupancy times according to a set schedule. The Occupancy schedule setting is also used to control internal gains and/or HVAC systems. In DesignBuilder, schedules can be used to define:

•
Occupancy times • Equipment, lighting HVAC operation • Heating and Cooling temperature setpoints For the purpose of the simulation, the authors used a slightly modified Design Builder template with default data for residential buildings occupancy was used, to better reflect Polish living conditions, hours of work etc.
Lighting was designed in compliance with the PN-EN 12464-1:2012 standard [77]. Lighting was assumed to be energy-efficient, with a value of 3.3 [W/m 2 ] per 100 [lux]. The lighting intensity level was assumed as in [78].
The remaining parameters were assumed in accordance with the principles of construction, technical knowledge and applicable standards.

Energy Consumption of Central Heating and Domestic Hot Water Preparation
The yearly final energy consumption for the purposes of heating, ventilation, domestic hot water preparation and lighting is to be presented in Figure 7. The yearly final energy consumption for the analysed buildings has been presented in Table 6.
The Central Heating installation was designed in the form of gas-powered heating using a dual purpose gas-powered condensation boiler with a total efficiency rating of 0.95.
The indoor air temperature was assumed per PN-EN 12831:2006. Instalacje ogrzewcze w budynkach-Metoda obliczania projektowego obciążenia cieplnego-a Polish standard for heating installations governing methods of calculating heat loads during the design stage [74]. Per this standard, spaces for permanent occupancy by persons without outdoor clothing should be designed to have a temperature of +20 °C, while spaces meant for occupancy by persons potentially without clothing (such as bathrooms), should have a temperature of +24 °C. This temperature difference can be obtained by properly setting a building's thermostats or by using BMS installations. Concerning the schedule, it was assumed that such temperatures should be kept 24 h per day throughout the entire year. In residential buildings with building automatics systems it is much easier to assume periods of lowered temperature (e.g., at night) or when residents leave the building to work. However, most buildings in Poland do not feature elaborate BMS installations. The analyses assumed the buildings to have such systems, which is why no temperature drops were accounted for.
Design Builder Software uses Occupancy model data which defines the number of people in a space and occupancy times according to a set schedule. The Occupancy schedule setting is also used to control internal gains and/or HVAC systems. In DesignBuilder, schedules can be used to define: • Occupancy times • Equipment, lighting HVAC operation • Heating and Cooling temperature setpoints For the purpose of the simulation, the authors used a slightly modified Design Builder template with default data for residential buildings occupancy was used, to better reflect Polish living conditions, hours of work etc.
Lighting was designed in compliance with the PN-EN 12464-1:2012 standard [77]. Lighting was assumed to be energy-efficient, with a value of 3.3 [W/m 2 ] per 100 [lux]. The lighting intensity level was assumed as in [78].
The remaining parameters were assumed in accordance with the principles of construction, technical knowledge and applicable standards.  In the case of the 110 m 2 house (B1/110), its external wall touches its "twin", which results in an adiabatic zone and a heat transfer value of zero for these walls. This results in a lower final energy  Table 7 presents the yearly energy consumption for heating, ventilation and domestic hot water preparation for the entire complex. The high operational energy demand has its source in the Polish definition of NZEBs. The general definition of nearly-zero energy buildings is provided in the Energy performance of buildings directive 2010/31/UE as "'nearly zero-energy building' means a building that has a very high energy performance, as determined in accordance with Annex I. The nearly zero or very low amount of energy required should be covered to a very significant extent by energy from renewable sources, including energy from renewable sources produced on-site or nearby". However, every European Union Member State has adopted its own NZEB parameters, adapted to local conditions. In Poland, NZEBs are defined by the Ordinance on the Technical Conditions that Must be Met by Buildings and their Placement [34]. In the Polish definition of NZEBs there are two criteria. The first criterion is ensuring the proper building envelope thermal insulation. The heat transfer coefficients for external partitions in the analysed buildings are compliant with the values for NZEB buildings in Poland. The second criterion is for NZEBs in Poland pertains to Primary Energy value (PE), which cannot exceed a value set for different types of buildings. For single-family residential buildings without cooling, the maximum PE value is 70 kWh/(m 2 annum). In the article, the authors based their analyses on operational energy (OE), even though it is not listed in the requirements of the Polish NZEB definition. However, operational energy is more useful for the analyses performed by the authors than primary energy, which is a criterion for NZEBs in Poland. A definition of near-zero energy buildings for Poland has been included in the article. PE values were compared with maximum values for the Polish standard for NZEBs. The PE values for the analysed buildings are lower than the maximum values stipulated in Polish law, despite the fact that they are not fitted with renewable energy sources. This is why the authors define the buildings under analysis as NZEBs in their discussion. Table 8 shows the energy demand values for type A and type B buildings.

CO 2 Emission in the Analysed Building Types
Calculation of CO 2 emissions was based on [52,59]. The calculation of CO 2 emissions is shown in Table 9. The yearly CO 2 emission for the entire complex (40 buildings) has been presented in Table 10. CO 2 emissions are associated with the amount of fuel consumed for heating, as well as other types of fuel, such as electric power associated with, among other things, lighting. The low CO 2 emission of the buildings is also a product of their structural system. A reinforced concrete structure is associated with greater CO 2 emissions relative to a timber structure [78].

Thermal Comfort in the Analysed Building Layouts
The thermal comfort assessment analysis was carried out in a typical summer week, characteristic of climate data considered in this article (Southern Poland, around Krakow). The authors assumed that this period would be appropriate for the purposes of this paper, because it represents weather conditions that are repetitive and will most likely occur every year at a given location (this is also an option resulting from the software). The thermal comfort analysis includes building use schedules and HVAC system operation, which are discussed in the article. Graph 11 shows the daily average PMV values that are used for comparison between the buildings analyzed. Figure 8 presents the average values of thermal comfort for the analysed buildings. CO2 emissions are associated with the amount of fuel consumed for heating, as well as other types of fuel, such as electric power associated with, among other things, lighting. The low CO2 emission of the buildings is also a product of their structural system. A reinforced concrete structure is associated with greater CO2 emissions relative to a timber structure [78].

Thermal Comfort in the Analysed Building Layouts
The thermal comfort assessment analysis was carried out in a typical summer week, characteristic of climate data considered in this article (Southern Poland, around Krakow). The authors assumed that this period would be appropriate for the purposes of this paper, because it represents weather conditions that are repetitive and will most likely occur every year at a given location (this is also an option resulting from the software). The thermal comfort analysis includes building use schedules and HVAC system operation, which are discussed in the article. Graph 11 shows the daily average PMV values that are used for comparison between the buildings analyzed. Figure 8 presents the average values of thermal comfort for the analysed buildings. The buildings were designed as low energy and have a high degree of air tightness. In addition, they are not equipped with external solar screens, so they receive high solar heat gains thanks to large glazed surfaces. The PMV value of 0.7-1.27 during a sunny and warm summer week can be considered a predictable result [79].
Overheating in buildings, as simulations show, is a real problem in heavily insulated, airtight buildings with windows facing south. It was assumed that no shading elements, such as blinds, venetian blinds or brise-soleils would be installed on buildings. Such assumptions were intentional and are part of another research project, which aims to determine the extent to which solar barriers improve the comfort of using buildings with almost zero energy consumption. The buildings analysed in the article did not have solar barriers, which is why their interiors were overheated.
In order to improve the thermal comfort of the analysed buildings, additional shading mechanisms should be added (external louvers, blinds, etc.). The shape of greenery around the The buildings were designed as low energy and have a high degree of air tightness. In addition, they are not equipped with external solar screens, so they receive high solar heat gains thanks to large glazed surfaces. The PMV value of 0.7-1.27 during a sunny and warm summer week can be considered a predictable result [79].
Overheating in buildings, as simulations show, is a real problem in heavily insulated, airtight buildings with windows facing south. It was assumed that no shading elements, such as blinds, venetian blinds or brise-soleils would be installed on buildings. Such assumptions were intentional and are part of another research project, which aims to determine the extent to which solar barriers improve the comfort of using buildings with almost zero energy consumption. The buildings analysed in the article did not have solar barriers, which is why their interiors were overheated.
In order to improve the thermal comfort of the analysed buildings, additional shading mechanisms should be added (external louvers, blinds, etc.). The shape of greenery around the buildings should be adjusted to increase the amount of air supplied to the mechanical ventilation or air conditioning installations. The authors calculated the average value of thermal comfort for each system to compare different urban layouts. The thermal comfort value for system A is 1.075, 1.21 for system B and 0.94 for system C.

Assessment of the Analysed Urban Layouts
The results of the survey performed among the potential residents of the experimental housing complex have been presented in Figure 9. The number of persons listing a given layout as the most desirable one is indicated by the blue column (no. 1), the red column shows the number of persons who rated the given layout as the second-most desirable (no. 2), while the number of persons who rated the layout as the least desirable has been shown in green (no. 3).
air conditioning installations. The authors calculated the average value of thermal comfort for each system to compare different urban layouts. The thermal comfort value for system A is 1.075, 1.21 for system B and 0.94 for system C.

Assessment of the Analysed Urban Layouts
The results of the survey performed among the potential residents of the experimental housing complex have been presented in Figure 9. The number of persons listing a given layout as the most desirable one is indicated by the blue column (no. 1), the red column shows the number of persons who rated the given layout as the second-most desirable (no. 2), while the number of persons who rated the layout as the least desirable has been shown in green (no. 3). The mean result of the surveys and the ranking of each layout has been presented in Table 11. The respondents who selected Layout C as the most desirable reported that it was the layout that provided the most privacy (no adjacent neighbour). This layout was ranked as the most desirable by the greatest amount of respondents (16). Layout B was reported to be the least desirable by the largest amount of respondents (19). The respondents who rated this layout as the most desirable (rank 1) stressed their opinion as to lower construction and occupancy costs in semi-detached buildings when compared to other alternatives. Table 12 will present all of the results for each layout. As mentioned previously, the weights for each criterion, as included in the analysis, were set by the developer. They reflect the developer's preferences concerning each criterion. The developer decided that the most important criterion was total construction cost (a weight of 0.55). As the buildings will be sold to individual clients, their price was a key factor. The developer also took thermal comfort (a weight of 0.2) and energy consumption (0.15) into consideration, as this was intended to encourage potential clients and raise The mean result of the surveys and the ranking of each layout has been presented in Table 11. The respondents who selected Layout C as the most desirable reported that it was the layout that provided the most privacy (no adjacent neighbour). This layout was ranked as the most desirable by the greatest amount of respondents (16). Layout B was reported to be the least desirable by the largest amount of respondents (19). The respondents who rated this layout as the most desirable (rank 1) stressed their opinion as to lower construction and occupancy costs in semi-detached buildings when compared to other alternatives. Table 12 will present all of the results for each layout. As mentioned previously, the weights for each criterion, as included in the analysis, were set by the developer. They reflect the developer's preferences concerning each criterion. The developer decided that the most important criterion was total construction cost (a weight of 0.55). As the buildings will be sold to individual clients, their price was a key factor. The developer also took thermal comfort (a weight of 0.2) and energy consumption (0.15) into consideration, as this was intended to encourage potential clients and raise the price of buildings. The developer saw CO 2 emissions and urban layout assessment (each with a weight of 0.05) as the least important. Other weight values would produce different multi-criteria analysis results, which is why the authors also performed an analysis of the impact that weights can have on scores, as presented later in the paper. The scores for each urban layout produced by the analysis have been presented in Table 13. The scores have been presented on the scale of the entire complex (40 buildings). The authors used the WASPAS method to obtain the results. This made it possible to determine the optimal solution in light of the analysed factors (assessment criteria). We also performed a simulation for different criteria weight values using the WASPAS method. The randomly generated weights and their respective results are presented in Table 14.

Discussion and Conclusions
The results presented in Table 13 can aid the decision-maker in making decisions in selecting an architectural design alternative. The results were obtained for weights as defined by the developer and formulated to their individual preferences. An initial analysis of standardised values can provide valuable information about each alternative. Alternative B appears to clearly outrank the two other alternatives (highest value in four parameters). It also appears that alternative A is better than alternative C, as it has higher Q i scores in all parameters except thermal comfort. The WASPAS method confirms this analysis. Alternative B was selected as the best, followed by alternative A, with alternative C ranked the lowest. The differences in Q i scores were not significant (a difference of 1.5% and 2.8% relative to alternative B, for alternatives A and C, respectively).
An analysis of different weight values on the final results of the overall analysis was also performed (see Table 14). Generated with a uniform distribution, the weights and values from the Case Study were used to calculate Q i score values in the WASPAS method. Statistically, alternative B was indicated to be the best (it ranked first in 80% of cases). Of note is that alternative C (ranked first in 20% of cases) was shown to be markedly better than alternative A (no first ranks). Therefore, it can be concluded that alternative A was dominated by alternatives B and C, which indicates that it can be excluded from the decision-making process. Alternative C outranked alternative B when a sufficiently high value (ca. 0.3) was assigned to the weight associated with thermal comfort. It is interesting that alternative B received the lowest rank in 10% of cases. It was even outranked by alternative A when thermal comfort was assigned a weight with an appropriately high value (ca. 0.3) and urban layout assessment was assigned a low value (ca. 0.05). However, alternative C was shown to outrank the remaining alternatives each time when this was the case.
In the presented Case Study, only three alternatives of the design of the housing complex were formulated. The proposed model can be used to support decision-making problems with a greater number of alternatives. The method can also support the decision-maker in rejecting dominated alternatives, reducing the number of analysed alternatives.
It should be noted that analysing such complex cases, which incorporate parameters with different units, value ranges and significance, is impossible without using multi-criteria analysis methods. Despite the significant contribution of these methods to analysis, the decision-maker is irreplaceable. It is up to the individual to quantify their preferences and expectations and it is they who make the final decision. The method presented in the paper can only support them in this.
The work presents pioneering research associated with formulating a method of determining the optimal urban layout of a complex of single-family houses. The study took various cost factors into consideration. The method presented in the article is compliant with the concept of sustainable development. The article does not refer to the results of the work of other scholars as no direct counterparts were found in the literature. The discussion presented in the article covers the results of a multi-criteria assessment meant to aid in the selection of an optimal urban layout of timber single-family houses while taking different types of building layouts into account. The discussion concerning the results of the study took into consideration energy demand, the direct costs of construction and its organisation, occupancy costs, thermal comfort, energy consumption, greenhouse gas emissions and the results of a survey concerning preferences as to the proposed architectural and urban layouts conducted among potential home buyers. The authors assumed a goal function in the form of a sum of all essential costs, which, along with the constraints, is linear. The Simplex method was used to solve the linear model. Sustainability parameters were determined through analyses employing EnergyPlus code. Afterwards, using the procured data, a multi-criteria optimisation was performed, using the WASPAS method. The results of this optimisation made it possible to report a solution that can be considered optimal in light of the adopted assessment criteria.

Future Research Directions
Future analyses should automatically include the possibility of selecting construction materials, the shape of the buildings themselves and their placement. The user would define permissible alternatives and a computer program would, through discrete optimisation, set parameters that could prove useful in multi-criteria optimisation.