Optimising Embodied Energy and Thermal Performance of Thermal Insulation in Building Envelopes via an Automated Building Information Modelling (BIM) Tool

: Insulation systems for the ﬂoor, roof, and external walls play a prominent role in providing a thermal barrier for the building envelope. Design decisions made for the insulation material type and thickness can alleviate potential impacts on the embodied energy and improve the building thermal performance. This design problem is often addressed using a building information modelling (BIM)-integrated optimisation approach. However, one major weakness that lies in the current studies is that BIM is merely used as the source for design parameters input. This study proposes a BIM-based envelope insulation optimisation design tool using a common software Revit and its extension Dynamo to ﬁnd the trade-o ﬀ between the total embodied energy of the insulation system and the thermal performance of the envelope by considering the material type and thickness. In addition, the tool also permits data visualisation in a BIM environment, and automates subsequent material library mapping and instantiates the optimal insulation designs. The framework is tested on a case study based in Sydney, Australia. By analysing sample designs from the Pareto front, it is found that slight improvement in the thermal performance (1.3399 to 1.2112 GJ / m 2 ) would cause the embodied energy to increase by more than 50 times.


Introduction and Literature Review
The construction industry plays a prominent role in addressing energy use and emissions. In 2018, it was responsible for 36% of the total energy consumed and 39% of process-related carbon dioxide emissions [1]. Energy ensued from buildings can be categorised into two types; energy consumed at the operational stage and the energy capital of all building materials, referred to as embodied energy [2]. However, "energy" and "operational energy" are usually used interchangeably, while the embodied energy of materials is often ignored [3]. The consequence is that energy savings in the operational stage are achieved at the price of the accumulation of embodied energy [4]. Embodied energy accounts for a significant part in the total energy requirement [5]. In a domestic context, Treloar [6] indicated that embodied energy is 20 to 50 times the annual operational energy for most buildings in Australia. One holistic approach for minimising environmental impacts posed by the built environment is by considering materials' impact on the building performance and their embodied energy [7].

A Building Information Modelling (BIM)-Integrated Optimisation Approach to Support Decision-Making
In respect of the envelope insulation design problem, design parameters are largely dependent on the envelope assemblies, which are inevitably associated with tedious construction documentation in the early design phase. Design decisions made at the early design stage have significant impacts on the sustainability performance of final deliveries [23]. As a response to this, a Building Information Modelling (BIM) integrated approach is regarded as most suitable due to the comprehensive design information that BIM can store [24,25]. Krygiel et al. [26] also argued that BIM can effectively leverage material selections to achieve building performance and environmental goals.
Within the realm of BIM-integrated studies, some utilise BIM as a pure information platform to derive material data, such as insulation type and thickness, in order to calculate a series of performance or environmental impact indicators [27], or give a score to the design based on pre-determined weightings [28]. To account for more dynamic operational conditions, some studies link energy simulation processes to the BIM model. Cho et al. [29] attempted to do so by transferring the BIM model information in the format of gbXML schema (Green Building XML) to an energy simulation tool, Ecotect [30], and conducted a parametric analysis on the insulation type. Similarly, Jalaei and Jrade [31] developed a Revit plugin that exports the design to gbXML format for energy simulation in Ecotect. After Ecotect was discontinued in 2015 [32], Kim et al. [33] presented a workflow that maps an external material database to the material information in a standard Industry Foundation Class format (IFC) file and creates a new file input for subsequent energy simulation. Nonetheless, the premise of these approaches is that designers can always make high quality decisions effectively based on their past experience.
Some studies focus on the area of combining the simulation with optimisation process in a BIM environment to investigate material selections. Bank et al. [34] designed a decision-making framework by integrating a BIM model with a System Dynamic (SD) decision-making software, where a simulation-based optimisation on BIM design parameters are carried out. SD is a computer modelling method that simulates the feedback structure in a system, namely the components and the resulting behaviours [35], and it is believed to be suitable for sustainability design problems with highly dynamic parameters that are time or space dependent [34]. In this framework, the linkage between the BIM model, the SD tool AnyLogic [36], and energy simulation tools requires proficiency in the use of Visual Basic or C# and Java. In addition, the use of SD is very limited in the construction industry [37], showing a lack of feasibility for further adoption of the proposed framework in practice. Similarly, a BIM-based thermal performance multidisciplinary design optimisation framework was proposed by Welle et al. [38], which delivered an automated workflow to address technical barriers such as interoperability issues between modelling and analytical software and poorly coordinated import-export procedures. This workflow also involves building up customised connections with the external energy simulation programme EnergyPlus [39] and Radiance [40].
More recently, common BIM software such as Revit [8] has incorporated energy simulation and Lifecycle Analysis (LCA) as its integrated plugins. An optimisation package is also available in Revit via its Visual Programming Language (VPL) platform, Dynamo [41]. However, at the early design stage, as suggested by Shadram et al. [42], LCA is more often carried out at a highly detailed level considering the scope of the analysis and its poor interoperability with other BIM-driven design process. As for the insulation design problem, Shadram and Mukkavaara [43] developed an optimisation framework utilising BIM for placing insulations focusing on the embodied energy and the operational energy, which is obtained from the simulation results, featuring modelling software Revit and Grasshopper [44], a VPL platform of Rhino [45]. Even so, this framework is still highly relying upon an external database as the material inventory and the core of the framework for simulation and optimisation is carried out by the non-BIM software [46], Grasshopper, whereas the intrinsic characteristic of BIM as a data platform to store model and project information is underappreciated.
BIM, as the building design information database [25,47], should be used more than a platform where material information is simply read as optimisation input. The BIM model needs to stay in sync with the optimal design model in the way that any parameter and project-specific information changes are properly documented. Furthermore, by visualising the data in the VPL platform [48] after the optimisation, it is believed that users can make more informed decisions [49]. Given that Revit is a BIM platform that is adopted by a majority of construction engineers and architects around the world [46], there is still a need for an automated optimisation-based BIM framework as a design decision support system for selecting the most suitable insulation material and thicknesses in building projects. The aim of this research is therefore to develop a BIM-based optimisation framework for the building envelope insulation system, emphasising on the embodied energy and thermal performance. A BIM-based insulation design framework using Revit [50] is proposed, to enhance the decision-making process for designing envelope insulation system accounting for its embodied energy and thermal performance. Embodied energy used in this study is defined as the total primary energy ensued from the material extraction, manufacturing, and transportation within a cradle-to-gate boundary [18]. The objective is to deliver an automatic workflow that (i) finds the trade-off between the embodied energy and thermal performance while considering the insulation type and thickness; and (ii) places the insulations in the model; and (iii) maps insulation material information to the BIM platform.
The remainder of this paper is organised as follows. Section two first presents the proposed optimisation framework, and then gives the mathematical description of the optimisation problem; this section also provides a detailed illustration of framework components together with tasks carried out. Following that, the proposed framework is applied on a realistic case study to showcase the feasibility of it. Conclusions are drawn at the end.

Materials and Methods
The proposed framework contributes to the research area of BIM-integrated optimisation by providing a decision-making system for building envelope insulation design, complete with automated material documentation procedures in the BIM model. In this section, the detailed framework components are presented.

The Proposed Framework
An optimised design-based building framework for insulation selection and dimensions that places emphasis on the embodied energy and the thermal energy performance of the envelope insulation system has been developed. The proposed framework integrates the multi-objective optimisation insulation design problem into a BIM modelling environment, enabling a seamless workflow for optimisation setup, data visualisation, model modification, and material documentations. As shown in Figure 1, the overview of the framework, which is further broken down into specific details to be outlined later, is comprised of four main components: (i) BIM Software; (ii) Envelop Insulation Optimisation; (iii) Visualisation, and (iv) Nomination of Insulation Design. The BIM integration in this study is mainly two-fold. First, the intimacy between Dynamo and Revit allows the seamless transfer of optimisation parameters, such as existing envelope assemblies, quantity take-offs of each assembly layer, and the properties of assembly materials. This is achieved via the first module, BIM Software, where no manual work is required. However, it is expected that each envelope component is designed according to the construction specifications with correct  The BIM integration in this study is mainly two-fold. First, the intimacy between Dynamo and Revit allows the seamless transfer of optimisation parameters, such as existing envelope assemblies, quantity take-offs of each assembly layer, and the properties of assembly materials. This is achieved via the first module, BIM Software, where no manual work is required. However, it is expected that each envelope component is designed according to the construction specifications with correct assemblies and the material thermal information properly embedded. This allows Dynamo to extract the detailed assemblies of the building envelope and the quantities of the material used, as well as the corresponding material properties from Revit as the design parameters. Meanwhile, design variables, namely available insulation type and thickness to be tested, are exported to Dynamo from Comma-separated Value (CSV) files to initiate the optimisation. After the optimisation, Dynamo platform allows visualising all feasible design objective function results by interacting with the local directories specified by users where designated CSV files are stored and updating the envelope elements in Revit automatically via accessing Revit Application Programming Interface (API). Detailed explanations of these two integration processes are offered in Sections 2.2 and 2.4 respectively.
In this framework, major data transfer paths between the building model and the VPL platform; and the VPL platform and local directories are highlighted in Figure 1. Key features of the data interaction paths design lies in the characteristics of the envelope insulation problem that the design parameters are uninsulated layers that are held as constants in the Revit model while the design variables also require corresponding visual presentation and information documentation in the BIM model. Different types of tasks performed within the other components of Figure 1 are discussed below.

Envelope Insulation Optimisation
The second component, the Optimisation Module, is illustrated in Figure 2. Data Transfer Node 1 is responsible for bridging Dynamo with other platforms to fetch information for the envelope optimization problem. It involves preparing two classes of optimisation inputs, the first of which aims to read available insulation type and thickness from CSV files. The second class is the material information of envelope assemblies from the BIM model. As suggested in Figure 2, these design parameters can be retrieved from three Revit dialogs. Figure 3 is an illustrative example of how the information for a wood finish timber floor is displayed in Revit.     Figure 3A,B are the BIM dialogs of the Element Assembly and Material Thermal Assets as shown in Data Transfer Node 1. For each layer of the envelope component, its function, material applied, and thickness nominated need to be defined. Thermal assets that are relevant to this study are thermal  Figure 3A,B are the BIM dialogs of the Element Assembly and Material Thermal Assets as shown in Data Transfer Node 1. For each layer of the envelope component, its function, material applied, and thickness nominated need to be defined. Thermal assets that are relevant to this study are thermal conductivity and density, as shown in Figure 3B, and they are associated with the material at each layer. The material assembly and properties hosted can be obtained using the Dynamo commands group A and B illustrated in Figure 4A,B. Command C generates the floor material take-off schedule shown in Figure 3. Note that it is important to insert an insulation layer at the model building stage. This allows the nominated insulation area shown in the schedule, which will provide more accurate results in the subsequent optimisation.
Buildings 2020, 10, x FOR PEER REVIEW 7 of 23 conductivity and density, as shown in Figure 3B, and they are associated with the material at each layer. The material assembly and properties hosted can be obtained using the Dynamo commands group A and B illustrated in Figure 4A,B. Command C generates the floor material take-off schedule shown in Figure 3. Note that it is important to insert an insulation layer at the model building stage. This allows the nominated insulation area shown in the schedule, which will provide more accurate results in the subsequent optimisation. Task One in Figure 2 is where problem-specific parameters are specified, including the file path for exporting datapoints during the optimization process, optimisation setup, and the Heating Degree Days (HDD) and Cooling Degree Days (CDD) based on local climate data. These last two parameters are used to estimate the energy requirement when using a specific envelope system and Task One in Figure 2 is where problem-specific parameters are specified, including the file path for exporting datapoints during the optimization process, optimisation setup, and the Heating Degree Days (HDD) and Cooling Degree Days (CDD) based on local climate data. These last two parameters are used to estimate the energy requirement when using a specific envelope system and are further discussed in the following section. Optimisation Problem Description refers to inputs of parameters that specify distinct optimisation scenario including population size, number of iterations, and variable ranges. Variable ranges correspond to the discrete list input of insulation type and thickness. Population size and iteration times refer to the total number of individuals included at each iteration and when the optimisation will stop. Any tailoring to objective functions is also included in Task One.
The core of this component, Task Two, employs Optimo [47], a Non-dominated Sorting Genetic Algorithm (NSGA-II) optimisation package of Dynamo to conduct the optimisation. NSGA-II is a multi-objective evolutionary algorithm in essence and utilises nondominated sorting techniques [51]. With necessary information supplied through the data interaction and Task One, the optimisation calculation process is carried out. The detailed mathematical equations formulated for the envelope insulation design problem are introduced in the following section.
Optimisation Mathematical Model Table 1 summarises the notation set adopted to create the optimisation model proposed. In order to evaluate the primary design objectives, namely the embodied energy and the thermal performance, two sets of equations are formulated. The total embodied energy of the insulation system is calculated by summing up the embodied energy of the nominated insulation materials. The thermal performance is assessed by evaluating the energy requirement under static conditions using degree-days method, an approach commonly seen in insulation optimisation problems [20,52]. This method translates dynamic weather conditions into two local constant indicators, Heating Degree Days (HDD) and Cooling Degree Days (CDD), reflecting how much the daily average temperature deviates from a base comfortable thermal level over a certain period [53]. Table 1. Notation Set.

Sets c ∈ BE
Component c belongs to the building envelope system BE, including the roof, floors and external walls.
m ∈ M Set that contains available insulation materials.
t c ∈ thk c Set that includes a range of thickness to be tested for envelope component c.

A c m
The insulated area of component c using material m.
∆T Temperature difference.

HDD
Heating Degree days.

CDD
Cooling Degree days.
η Mechanical heating and cooling system efficiency.

R c in
The thermal resistance of the inner air film of component c.   The volume of an insulation layer using material m for component c. λ m Thermal conductivity of material m.

EE
Embodied Energy of the insulation system.

HT
Heat Transfer via the insulation system.
Q Heat transfer rate.
Within literature, one common approach to study the thermal resistance (R value) of wall assemblies, roof structure, and floor construction is by simple addition of the nominal thermal resistance value of each layer [8,[54][55][56][57].
However, this method is established on the assumption that heat flow in every cross section of the assembly is homogeneous [58]. In some occasions such as when framing elements are used, thermal bridges will form, providing a faster heat transfer path [59]. In this study, isothermal plane method is employed to estimate the total R value of the building envelope. This approach adds up the effective thermal resistance of each layer of a building assembly, when more than one heat flow path is presented in one layer, which is common when there exists an opaque-cavity structure, the thermal resistance of each heat path will be calculated in proportion to the area [60].
The first objective function Embodied Energy (EE) aims to minimise the embodied energy per unit envelope area [21], which is defined in Equation (1): where ρ m ·V m c calculates the mass of the insulation material m for envelope component c. EE m is the embodied energy coefficient given in MJ/Kg. The denominator denotes the total insulated area of the building envelope. V m c indicates the total volume of insulation material m used in component c, which is calculated as shown in Equation (2): where A m c is the area of insulation material m used in component c, measured in square metres, and t c is the corresponding insulation thickness. Considering different envelop components might have different requirements on insulation thickness, thk c specifies a range of thickness input for each component.
Equation (3) calculates the heat transfer rate where U is the overall thermal transmittance and ∆T is the temperature difference: Buildings 2020, 10, 218 9 of 23 Drawing upon this equation, Equation (4) utilises the degree-days approach to estimate the energy requirement [61,62]: In Equation (4), the energy requirement is expressed in GJ/m 2 . HDD and CDD are the heating and cooling degree-days for a specified region respectively, and η indicates the system efficiency. U is the overall heat transfer coefficient and is defined in Equation (5): U c is the heat transfer coefficient for a building envelope component. According to the Isothermal Plane method, it is defined as Equation (6) [60]: In this equation, R c in and R c out are constants for the each envelop component c, denoting the R-value for the inner and outer air film. ε c is a binary parameter which equals 0 when component c refers to the slab-on-ground floor system, and 1 otherwise [63]. The R c cont represents the total R-value of continuous layers, where only one heat flow path is identified, and therefore the resistance values are additive. One design constraint considered is related to the U c of the proposed insulation material. This is often regulated by local building codes [65][66][67]. This study adopts the maximum U-value allowed developed by the Australian national construction regulatory framework, the National Construction Code (NCC) [68], for different element assemblies on the building envelope. As a result, 0.24 W/m 2 ·K, 1 W/m 2 ·K and 0.36 W/m 2 ·K are set as the upper bound for roof, floor and wall construction, respectively [67].

Visualisation
The optimal solutions produced by NSGA-II can be visualised in the form of a Pareto front when no further improvements can be made for any variable without harming other variables [69]. By visualising the Pareto-optimal design solutions, designers are able to select among a series of equally good design options more effectively [70]. The visualisation component in the proposed framework aims to provide visual presentation of the distribution of all feasible objective function values documented during the optimisation, such that whether the process has achieved its pareto optimality is made obvious to decision makers. Feasible designs are those that do not violate the design constraint. Figure 5 demonstrates the Dynamo setup for visualising datapoints. All feasible results for the objective function Embodied Energy will be plotted on the X-axis and the Energy Requirement result is visualised on the Y-axis. Buildings 2020, 10, x FOR PEER REVIEW 10 of 23

Nomination of Insulation Design
The aim of the fourth Module of the framework is to update the insulation design of the building envelope based on the decision made with reference to the visualised datapoints from the optimisation process and the Pareto-optimal solutions. Meanwhile, the material library will be maintained through a process that maps the optimal insulation materials. This is achieved mainly via three main steps; (i) selection of the optimal insulation design from the Pareto-optimal list, (ii) mapping the corresponding thermal properties into the material library, and (iii) instantiation of the solution in the BIM model. Figure 6 depicts the proposed workflow. Most building materials displayed in BIM models are not merely a graphical representation [71]. BIM software has a collection of material visuals and thermal and physical properties as its native material library, which can be easily modified and appended on a project basis. Data Transfer Node 2 first iterates materials from the final solution list determined by the designer and if the material is already present in the library, its thermal properties will be modified according to the optimisation result, otherwise a new insulation type will be created with the corresponding properties, namely density and thermal conductivity embedded. Python scripts are used to interact with the Revit material library using the Revit API. The scripts can be found in Appendix A. Another set of commands duplicates the external walls, roof, and floor in the current BIM document, and then

Nomination of Insulation Design
The aim of the fourth Module of the framework is to update the insulation design of the building envelope based on the decision made with reference to the visualised datapoints from the optimisation process and the Pareto-optimal solutions. Meanwhile, the material library will be maintained through a process that maps the optimal insulation materials. This is achieved mainly via three main steps; (i) selection of the optimal insulation design from the Pareto-optimal list, (ii) mapping the corresponding thermal properties into the material library, and (iii) instantiation of the solution in the BIM model. Figure 6 depicts the proposed workflow.

Nomination of Insulation Design
The aim of the fourth Module of the framework is to update the insulation design of the building envelope based on the decision made with reference to the visualised datapoints from the optimisation process and the Pareto-optimal solutions. Meanwhile, the material library will be maintained through a process that maps the optimal insulation materials. This is achieved mainly via three main steps; (i) selection of the optimal insulation design from the Pareto-optimal list, (ii) mapping the corresponding thermal properties into the material library, and (iii) instantiation of the solution in the BIM model. Figure 6 depicts the proposed workflow. Most building materials displayed in BIM models are not merely a graphical representation [71]. BIM software has a collection of material visuals and thermal and physical properties as its native material library, which can be easily modified and appended on a project basis. Data Transfer Node 2 first iterates materials from the final solution list determined by the designer and if the material is already present in the library, its thermal properties will be modified according to the optimisation result, otherwise a new insulation type will be created with the corresponding properties, namely density and thermal conductivity embedded. Python scripts are used to interact with the Revit material library using the Revit API. The scripts can be found in Appendix A. Another set of commands duplicates the external walls, roof, and floor in the current BIM document, and then Most building materials displayed in BIM models are not merely a graphical representation [71]. BIM software has a collection of material visuals and thermal and physical properties as its native material library, which can be easily modified and appended on a project basis. Data Transfer Node 2 first iterates materials from the final solution list determined by the designer and if the material is already present in the library, its thermal properties will be modified according to the optimisation result, otherwise a new insulation type will be created with the corresponding properties, namely density and thermal conductivity embedded. Python scripts are used to interact with the Revit material library using the Revit API. The scripts can be found in Appendix A. Another set of commands duplicates the external walls, roof, and floor in the current BIM document, and then assigns them with a new name that indicates the total thickness and the type of each envelope component, such as "Roof_110", see Figure 7A, which also serves as the optimisation input. Figure 7B indicates the optimal insulation thickness determined by the designer.
Buildings 2020, 10, x FOR PEER REVIEW 11 of 23 assigns them with a new name that indicates the total thickness and the type of each envelope component, such as "Roof_110," see Figure 7A, which also serves as the optimisation input. Figure  7B indicates the optimal insulation thickness determined by the designer. Finally, the optimal insulation material and thickness will be assigned to the renamed duplicates using two nodes from the Dynamo package Orchid [72]. Figure 8 gives an example of how a new type insulation material named 'test1′ is populated in a drop-down list after mapping it to the material library via the Revit API. Accordingly, the thickness input should correspond to the optimal solution for the material 'test1.′ In addition, the element type and index can be fed by scripts shown in the previous graphs, as commands in Figure 7A have extracted the element type for each envelope component, and the index, representing where the insulation layer is inserted, has been specified in Figure 2A.

Case Study
The developed framework was utilised on the envelope insulation design of a house located in Sydney, Australia. The case study has a total floor area of 59.4 m 2 and is built on a slab-on-ground. The house has a pitched roof with ceiling lining under timber rafters and adopts the brick veneer external wall system, which is the most common wall construction method used in Australian dwellings [67]. Figures 9 and 10   Finally, the optimal insulation material and thickness will be assigned to the renamed duplicates using two nodes from the Dynamo package Orchid [72]. Figure 8 gives an example of how a new type insulation material named 'test1' is populated in a drop-down list after mapping it to the material library via the Revit API. Accordingly, the thickness input should correspond to the optimal solution for the material 'test1'. In addition, the element type and index can be fed by scripts shown in the previous graphs, as commands in Figure 7A have extracted the element type for each envelope component, and the index, representing where the insulation layer is inserted, has been specified in Figure 2A.
Buildings 2020, 10, x FOR PEER REVIEW 11 of 23 assigns them with a new name that indicates the total thickness and the type of each envelope component, such as "Roof_110," see Figure 7A, which also serves as the optimisation input. Figure  7B indicates the optimal insulation thickness determined by the designer. Finally, the optimal insulation material and thickness will be assigned to the renamed duplicates using two nodes from the Dynamo package Orchid [72]. Figure 8 gives an example of how a new type insulation material named 'test1′ is populated in a drop-down list after mapping it to the material library via the Revit API. Accordingly, the thickness input should correspond to the optimal solution for the material 'test1.′ In addition, the element type and index can be fed by scripts shown in the previous graphs, as commands in Figure 7A have extracted the element type for each envelope component, and the index, representing where the insulation layer is inserted, has been specified in Figure 2A.

Case Study
The developed framework was utilised on the envelope insulation design of a house located in Sydney, Australia. The case study has a total floor area of 59.4 m 2 and is built on a slab-on-ground. The house has a pitched roof with ceiling lining under timber rafters and adopts the brick veneer external wall system, which is the most common wall construction method used in Australian dwellings [67]. Figures 9 and 10 are the BIM model and the floor plan for the case study.

Case Study
The developed framework was utilised on the envelope insulation design of a house located in Sydney, Australia. The case study has a total floor area of 59.4 m 2 and is built on a slab-on-ground. The house has a pitched roof with ceiling lining under timber rafters and adopts the brick veneer external wall system, which is the most common wall construction method used in Australian dwellings [67].  The building envelope system is shown in Table 2. The assemblies are presented in the following order: exterior side → interior side. Where no construction details are available in the drawings, sustainable design specifications offered in the Australian government website YourHome [73] were referred to as a general guideline.   The building envelope system is shown in Table 2. The assemblies are presented in the following order: exterior side → interior side. Where no construction details are available in the drawings, sustainable design specifications offered in the Australian government website YourHome [73] were referred to as a general guideline.  The building envelope system is shown in Table 2. The assemblies are presented in the following order: exterior side → interior side. Where no construction details are available in the drawings, sustainable design specifications offered in the Australian government website YourHome [73] were referred to as a general guideline.
The HDD and CDD in Sydney adopted in this study are 743 and 556, measured with a base temperature of 18 degrees [75]. The system efficiency was assumed to be 80% [76]. In addition, the area ratio of the framing component for the brick veneer wall system and roof is assumed to be 20% [77,78]. Therefore, f c o and f c v shown in Equation (7) are assumed to be 20% and 80% respectively for both external wall and roof. Indoor Air Film 0.11 [74] A total of ten types of insulation were selected as optimisation input, including conventional insulation materials such as cellulose, mineral wool, a collective of insulation made of fibres [79], expanded purposes polystyrene (EPS), rock wool, one special type of mineral wool [79], polyurethane foam (PUR) [8], fibreglass batt [80] and emerging insulation like flax [81], recycled wool [82], and wood wool [8]. Their corresponding embodied energy coefficients following a cradle-to-ate approach, density and thermal conductivity were obtained from the Inventory of Carbon and Energy (ICE) database [18]. These design variables are presented in Table 3. Insulation thickness tested are some common dimensions seen in the market and this information was gathered from local suppliers, see Table 4. The indices in Tables 3 and 4 were used to identify the variables in the case study; for instance, ref. [1,10] was entered as the acceptable range for the variable insulation type.

Results and Discussion
The optimisation programme has a mutation probability of 0.01 and a crossover possibility of 0.9. The crossover and mutation distribution indices are 20 [24]. The Pareto front was obtained at the 150th generation with a population size of 200. The calculation time took no more than 10 min on a personal computer with a 2.5 GHz processor and 8 GB of RAM. Figure 11a gives a clear indication of the distribution of all feasible results in the Dynamo canvas that the black points from the final generation have formed the Pareto front. As a reminder, the horizontal axis represents the Embodied Energy function result, and on the vertical direction the value reflects the Energy Requirement function result. Note that in order to visualise this graph, geometry scaling was set as small in Dynamo, and a scale factor of 30 and 1/30 were applied to the Energy Requirement and Embodied Energy results respectively. For discussion purpose in this paper, the same group of results were also plotted on a 2-D scatter chart without geometry scaling to validate the Pareto optimality. This is illustrated in Figure 11b. In this graph, similarly, the lower the Energy Requirement value on the y-axis, the better the thermal performance of the insulation system. A total of 183 feasible solutions were obtained in the Pareto front. In addition to insulation design A and H, another six solutions have been identified in Figure 11b. The corresponding variables for solution A to H are listed Table 5. ER and EE are acronyms for the objective function Energy Requirement and Embodied Energy respectively. By contrasting designs A and H, it is noted that cellulose is opted for when considering the design objective Embodied Energy only as cellulose has the lowest embodied energy coefficient A total of 183 feasible solutions were obtained in the Pareto front. In addition to insulation design A and H, another six solutions have been identified in Figure 11b. The corresponding variables for solution A to H are listed Table 5. ER and EE are acronyms for the objective function Energy Requirement and Embodied Energy respectively. By contrasting designs A and H, it is noted that cellulose is opted for when considering the design objective Embodied Energy only as cellulose has the lowest embodied energy coefficient among all materials tested. In the meantime, design A is only associated with thinner insulation panels. However, the insulation dimension is bound by the design constraint on the system U value. From design A to H, there is a staggering difference in the Embodied Energy-this value of the design H is nearly 50 times of A. Insulation design B adopts the same material as A, while there is an increase in the thickness for floor insulation. Consequently, this leads to a 4.07% drop in Energy Requirement and a 16.56% increase in the embodied energy.
Moving towards designs with lower Energy Requirement, cellulose is replaced by fibreglass in design C, which has a slightly lower thermal conductivity. Thicker insulation panels were also proposed in this design. When comparing design C and D, it is noted that Embodied Energy doubled when rock wool is used in the floor assembly, even though it has better insulating properties and contains less embodied energy for a given mass. This is because rock wool's density is tenfold that of fiberglass. Similarly, solution E proposes PUR for the floor insulation, and achieved a 1.44% reduction on the estimated energy demand, while the energy capital for the insulation system soared from 33.1118 MJ/m 2 to 119.7994 MJ/m 2 (i.e., the difference is around 250%). PUR is known for its effective insulation property, but it also has the highest embodied energy coefficient value among all materials input in the case study.
It is also observed in the table that design H improved the Energy Requirement by around 2.66% compared to design C at the expense of storing 10 times more embodied energy. Meanwhile, in the case of designs A and B, more energy reduction is only associated with a 16.56% increase in embodied energy. This can be explained by the change in insulation dimension from solution A to B.
Another thing to note is that for all designs displayed in the table, between any two solutions, the change observed in Embodied Energy is always more vigorous than the decrease or increase in Embodied Energy. For instance, the Energy Requirement falls around 3.3% from design B to C, while the Embodied Energy has increased nearly four times. Biswas et al. [10] reported that the environmental impact associated with the embodied energy of the insulation material is secondary compared to the effect on reducing operation energy. The previous study focuses on a lifecycle assessment of the insulation material, which is not within the scope of this study. In fact, it was highlighted by Hammond et al. [18] that most insulation materials can pay back their embodied energy during their lifetime, in spite of the fact that common materials in the market have very poor energy capital. The Energy Requirement and Embodied Energy comparison results in this case study are in part in line with this argument and reaffirm that further improvements on the embodied energy for insulation materials should still be made.
Flax insulation, woodwool, and recycled wool are not chosen in any of the design. This is because the input data, namely the embodied energy coefficient, thermal conductivity, and density selected did not demonstrate significant advantages over other traditional materials since each material in reality has a wide range of variations. For instance, the density of the flax insulation mat could vary from 5 to 100 Kg/m 3 , and its thermal conductivity was reported to be in the range of 0.035 to 0.075 W/m·K [17]. Therefore, there is no one-size-fits-all insulation design strategy and the envelope insulation design framework developed in this study permits dealing with two critical insulation design considerations, namely the material type and thickness, while considering the environmental impacts.

Conclusions
A BIM-based envelope insulation design framework was proposed in this study. The proposed framework integrates a common BIM platform, Revit, with an optimisation model that is built within Dynamo to deal with two insulation design aspects: the insulation type and thickness for the roof, external walls, and floor, while addressing the embodied energy of the insulation system and the thermal performance of the building envelope. A Dynamo package that employs NSGA-II algorithm was utilised to solve the problem. Two objective functions that represent the embodied energy per insulated area (MJ/m 2 ) and the thermal energy requirement per envelope surface area (GJ/m 2 ) were formulated. Some other key features of the BIM framework include:

•
Plotting feasible objective function results on the Dynamo working canvas.

•
Mapping insulation material properties into the Revit material library.

•
Instantiating the optimal insulation design in the BIM model.
A case study was conducted using the developed tool and the results demonstrate that an improvement of the objective function Energy Requirement from 1.3399 to 1.2112 GJ/m 2 led to a drastic increase in the embodied energy contained by the insulation materials-the value has expanded 50 times more from 6.7115 to 360.0903 MJ/m 2 . Some limitations of the results generated lie in the insulation materials selected for the case study. The emerging insulation materials, such as flax, woodwool, and recycled wool, seem to be less preferred according to the optimisation results, as the data input in this case study are only representative of the properties for each insulation material, but cannot generalise all possible material variations, which in fact result in a wide range of acceptable values when discussing the properties of a given insulation material. Fire safety of insulation materials was also not accounted for. This is a significant design requirement as the building envelope system can effectively prevent the spread of fire [87]. According to the optimisation results, PUR seems to be an ideal option, while during a fire this material would pose potential health threats [88]. Additionally, the Energy Requirement adopts degree-days calculation approach based on static Sydney weather conditions. This leads to another limitation of this study that the objective function results are restricted to Sydney, Australia only, and would vary significantly depending on the location of the property tested. It is therefore critical to reiterate that the purpose of this study is not to make recommendations on the suitable insulation design but to showcase the applicability of the proposed optimisation framework.
One major limitation identified in previous BIM-based simulation or optimisation framework is that BIM is often only used as the data input source. This work contributes to the research area of BIM-integrated studies by providing an automated envelope insulation optimisation workflow with material mapping, design instantiation, and data visualisation features.
The weakness of this study is that it does not consider the whole lifecycle stage of the insulation materials, where the cost factor would also play a prominent role in the insulation design [89]. Furthermore, some future work could investigate the possibility of optimising the overall wall, roof, and floor assemblies in the BIM environment. This would lead to additional complexities as some assembly layers become optional and there is also a need to account for the structural performance [90].

Conflicts of Interest:
The authors declare no conflict of interest. Figure A1 demonstrates using the Revit API method IsNameUnique to evaluate whether the optimal insulation material is already documented in the Revit material library. [90]. Figure A1 demonstrates using the Revit API method IsNameUnique to evaluate whether the optimal insulation material is already documented in the Revit material library. The script in Figure A2 shows the procedure of updating the insulation material thermal conductivity and density in the Material Library based on the optimal design. The script in Figure A2 shows the procedure of updating the insulation material thermal conductivity and density in the Material Library based on the optimal design.