<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta>
      <journal-title-group>
        <journal-title>Journal of Solar Energy</journal-title>
      </journal-title-group>
    </journal-meta>
    <article-meta>
      <title-group>
        <article-title>Parametric Virtual Design-based Multi-Objective Optimization for Sustainable Building Design</article-title>
      </title-group>
      <contrib-group>
        <aff id="aff0">
          <label>0</label>
          <institution>Cunyi Cheng, Jelena Ninic, Walid Tizani University of Nottingham</institution>
          ,
          <country country="UK">United Kingdom</country>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2019</year>
      </pub-date>
      <volume>22</volume>
      <fpage>230</fpage>
      <lpage>245</lpage>
      <abstract>
        <p>With the development of Building Information Modelling aiming for automatic, the automating of sustainability analysis will be a certain requirement in the future. Unlike the current research stream, this paper investigates a novel approach of directly linking parametric architectural models to sustainability optimization through an automatic design-through analysis workflow with the support of parametric virtual design techniques. Data required for optimization was directly extracted from architectural models and then being evaluated. After the optimization, architectural models are automatically updated with optimal parameters through parametric virtual modelling steps. A case study was carried out with optimization of daylight and energy performance for a residential building. The result demonstrates the advantages of directly using architectural models instead of energy models and the possibility of further development.</p>
      </abstract>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction</title>
      <p>later design stages, both the working load of preparing analytical model, and the impact scale
of parameter change due to sustainability optimization, is significantly increased.
In recent years, many researchers have investigated different automatic approaches for
sustainability designs (Wang, et al., 2006; Tuhus-Dubrow &amp; Krarti, 2010; Evins, 2013; Hachem,
et al., 2013). However, the current state-of-the-art in this field has significant limitations for
practical design tasks. Most of the studies have been performed based on conceptual energy
models that cannot directly be used for architectural design in the BIM design process, or based
on architectural models with a limited set of parameters, which are too simple to provide critical
ideas for the building’s initial design. For example, Asl, et al. (2015) optimized the daylight
and thermal performances of a manually predesigned building in Revit with only changing
window size and window type as parameters. Backer (2017) optimized the Daylight Factor,
Solar Heat Gains and Thermal Energy Losses with varying room depth, window size and
building orientation for a highly simplified L-shape building. The research outcomes of these
papers are therefore only a step towards parametric optimization of sustainability, which is
however not detailed enough for informing an initial practical design.</p>
      <p>In conclusion, both the official manual optimization method provided by Autodesk® and the
automatic methods being explored in recent research are useful but not efficient or automated
enough considering the development of higher maturity level BIM in the future. In real-design
conditions, it is more reasonable to directly extract the building geometry and construction
information from architectural models and use them for analysis. After optimization, the model
can be automatically updated, with the optimal parameters acquired from calculation. However,
since BIM has grown out of the design tools for interactive object-based parametric design
(Eastman, et al., 2008), and recently developed tools enable the combination of the BIM
objectbased concept with parametric modelling techniques by also embedding some intelligence in
the relationship between objects (Boeykens, 2012). This could be used to bridge the limitations
in sustainable building optimization. Therefore, in this paper, we attempt to overcome the
barrier of current software and implementation by using parametric virtual design technique.
Instead of using conceptual masses to evaluate the building, by combining parametric virtual
design with optimization techniques, the feasibility and efficiency of linking multi-objective
sustainable analysis directly to an architectural model is investigated.</p>
    </sec>
    <sec id="sec-2">
      <title>2. Methodology</title>
      <p>In this paper, we propose an approach for the sustainable building design loop by extracting
information from an architectural model, running automated optimisation of building
sustainability in terms of daylight and energy performance, and updating the original model
using parametric virtual design technique. The performance of this approach is explored with a
predefined case study of a target building. The methodology includes parametric virtual design,
energy performance assessment, sensitivity analysis and optimisation. To cover the relevant
areas in the initial sustainability design, parameters of both building geometry and construction
materials are considered. Prior to optimization, a sensitivity analysis is performed to decrease
the range of parameters, which helps to accelerate the convergence. Then, during the
optimization, similar to Generative Design (Zarzycki, 2012), possible solutions are investigated
using an optimization algorithm within predefined ranges, and a range of optimal solutions is
obtained through an iterative process. Then, the results are automatically applied with a
parametric virtual design technique to generate an example of an optimal building model.</p>
    </sec>
    <sec id="sec-3">
      <title>2.1 Linking Architectural Model and Sustainability Analysis</title>
      <p>An obvious disadvantage of using architectural models in optimization tasks is that standard
design software currently spends too much time in automatically regenerating models with
large amounts of architectural elements. However, two alternative approaches can be explored
to avoid this shortcoming. The first approach is by applying parametric virtual design, where a
conceptual geometry model (e.g. surfaces created in Dynamo) is constructed following the
same algorithm as the architectural model. In this case, these two models could change
synchronously with the variation of parameters. The conceptual geometry model can be applied
in the optimization process without generating an actual design model. After optimization, the
same optimal architectural model can be built by directly applying the optimal parameters. This
method works quickly with simple calculations, but appears to be insufficient for the analysis
software, which require a large amount of specific data to finish the calculation. The second
approach is to use an architectural model during the optimization, but only for the critical part
instead of the whole model. This means that only instances described with the parameters
relevant for the analysis will be considered, while other excess elements which are not
influencing the result or not worth studying will be neglected. An example for selecting the
critical parts of the model is demonstrated in the case study. This approach can be adopted for
different types of buildings, such as residential and office buildings, where the most important
rooms for sustainability design are normally identical, with a large amount of repetition inside
the building. In this case, only simplified models appear in the optimization with limited
possible variations. In this research, the second approach was explored because it is more
suitable for practical applications and can integrate with 3rd-party analysis tools.</p>
    </sec>
    <sec id="sec-4">
      <title>2.2 Optimization Algorithm</title>
      <p>The optimization algorithm adopted in this research is a Genetic Algorithm, specifically the
Non-dominated Sorting Genetic Algorithm-II (NSGA-II) (Deb, et al., 2002). This is a powerful
multi-objective optimization tool for environmental design (Kheiri, 2018). The function of
Genetic Algorithm is an improved version of Evolutionary Algorithms. As shown in Figure 1,
an initial list of samples is randomly created between predefined lower and upper bounds, then
a Fitness Function is used to calculate the series of initial objectives corresponding to the
variables. The Fitness Function in this research is a script including steps of extracting geometry
and material properties, evaluating daylight and energy performance, and exporting results for
further data processing. Then, initial variables and objectives are grouped into a list and
regarded as the initial generation, which is then input into a while loop for optimization. Inside
the loop, new parameter values are created after crossover and mutation. The Fitness Function
is applied again for generating the new generation and continued until the iteration number
approaches a predefined boundary. A special advantage for NSGA-II is that the elitism samples
among each child generation are compared with the last elitism before the creation of the next
generation, so that the variables in the next generation will automatically be close to the value
of the last elitism (Kheiri, 2018). In this case, the convergence is accelerated as the impact of
useless samples among a generation is significantly decreased. The result of optimization is
plotted and the best trade-off between the thermal and daylight performance is found. This is
followed by a further weighted-sum optimization to optimize the Normalized Average
Objective (NAO) (Qu &amp; Suganthan, 2010) calculated from the two performances. After this,
the difference between each sample in the Pareto front is compared with the optimal NAO value
to select the final best result.
The daylight simulation in this paper is evaluated by estimating the in-door lux values (lx). For
the design purpose, this value needs to be high, to reduce the electricity costs spent by the
lighting system. Professional 3rd-party analysing software has been used for this step. In this
research, the climate environment data for location in London is used. Building properties
include glass-to-wall ratio, façade orientation, glass transmittance and opaque material surface
properties such as roughness. In order to be compared with energy performance, the objective
Daylight Factor (DF), which is calculated based on the average in-door lux value, is introduced
to the optimization system as given in the following equation:
Where:
n is the total amount of individual lx values.
2.4 Energy Analysis

=
0.01

∑ 
 =1

In this research, the Solar Heat Gains (SHG) through glassing areas in summer and
Transmission Heat Losses (THL) in winter are calculated. For the design purpose, the value of
these factors needs to be lower, in order to reduce the energy consumption. SHG gives the total
solar heat absorbed by the walls and transmitted by the glass (Baker, 2017).
α is the absorptivity of glass defined;
U is the overall heat transfer coefficient for glass extracted from Revit W/(m2K);
h0 is the external heat transfer coefficient set as constant in this case study;
The total radiation can be calculated from direct radiation, diffuse radiation and reflected
radiation.
2.3-1
2.4-1
  =   cos[180° − ( +  )] ∙ cos( − | −  |) +   +  
d is the declination on the selected dates;
Σ is the tilt angle of surface (90°for vertical walls);
ξ is the wall azimuth angle;
ρg is the wall reflectivity;</p>
    </sec>
    <sec id="sec-5">
      <title>2.5 Sensitivity Analysis</title>
      <p>A sensitivity analysis has been carried out before the optimization is carried out in order to find
the individual impact of each parameter on the final sustainable performance. During the
sensitivity analysis, the target parameter is varied, while other parameters are kept constant. By
doing this, the range of each parameter between predefined upper and lower bond can be
Thermal Losses (TL) is the thermal energy diffused out of the building in winter when the
environment temperature is lower than the room temperature, in which the Transmission Heat
Losses (THL) are the most direct way. THL gives the heat directly transferred outside the
building through the building envelope, including walls and glass:</p>
      <p>( (
= ∑
 =1
/
), ∙</p>
      <p>)( 
 
,
−   )(1 + 100)
Where:
U is the heat transfer coefficient extracted from wall properties in Revit W/(m2K);
Tin is the temperature in room taken as 20 °C;
Tout is the temperature outside the room taken as 0 °C in winter;
α is the proportion of thermal bridges %;
Then, a normalised Thermal Factor (TF) is calculated following the general structure presented
below to make a comparison with DF. After finding the pareto front between DF and TF, the
final NAO could be calculated as shown:
 =  1( 1) 1 ±  2( 2) 2 ±  3( 3) 3 ± ⋯ ±   (  ) 
 
1
reduced to a reasonable value based on the evaluated results. Since the overall range of
parameters is decreased, the required number of generations for NSGA-II to reach the same
optimal result is significantly reduced. As a result, the convergence of optimization is
accelerated. Moreover, by comparing the resultant range of different parameters, the most
influential one can be found. In conclusion, it is efficient to run a sensitivity analysis before the
optimization. However, it needs to be considered that the total time cost of sensitivity analysis
is significantly increased with an increasing number of parameters in the optimization task.</p>
    </sec>
    <sec id="sec-6">
      <title>2.6 Parametric Virtual Design</title>
      <p>After getting the optimal range of parameters from sustainability optimization, values are
chosen and input to the script for generating building models with a parametric virtual design
tool. The building model is generated based on the coordinates of several control points. An
algorithm for positioning these points was developed such that the building profile is changed
by adjusting predefined geometrical parameters, for example, rotation of several control points
according to the orientation of one façade. Construction properties such as the insulation type
and the glass type are controlled by the family type parameters of the architectural elements.
Hence, the whole building, including geometric and construction parameters, is generated in a
fully automatic way, by importing parametrized family instances onto a parametrically
generated layout.</p>
    </sec>
    <sec id="sec-7">
      <title>3. Implementation</title>
      <p>As shown in Figure 2, following the general methodology briefly described in the previous
section, the parametric virtual design and multi-objective optimization is carried out using
Autodesk Revit® and its visual-programming plug-in Dynamo (Dynamo, 2019). With the
support of open-source packages in Dynamo, the optimization task with the Genetic Algorithm
is processed with the built-in package Optimo (Asl, et al., 2015) containing NSGA-II and the
daylight performance is analyzed by the package Honeybee included in the Ladybug Tools
(Mackey &amp; Roudsari, 2017).</p>
    </sec>
    <sec id="sec-8">
      <title>4. Case Study</title>
      <p>A case study was carried out to investigate the feasibility and efficiency of applying sustainable
optimization directly with an architectural model. The target building is a standard residential
building with 12 floors as shown in Figure 3.
As shown in Figure 4, following the general parametric virtual design method described in
previous section, the specific architectural model in this case study is generated. The algorithm
performs the following steps: i) generation of initial control points and lines for one section; ii)
and iii) rotations and symmetry operations for other sections; iv) generation of the reference
lines and points for the architectural elements; v) generation of the initial architectural layer;
and, vi) extraction of critical sections and generation of whole building.
As shown in Figure 5, the profile of each floor is identical and includes three critical parts: west,
middle and east. They are three residential areas in this building that are also the focus areas in
this case study. Daylight and energy consumption in the remaining areas, such as stair cases
and corridors, are not considered because normally no heating/cooling would be provided, and
natural lighting is not normally required in these areas.
As mentioned in the methodology, the critical parts are then selected as the target model for
optimization. The parameters for optimization in the middle critical part are shown in Table 1.
The discrete variables are default settings for architectural elements and materials provided by
the software. The continuous ones, are independently defined values, which are required for the
sustainability performance evaluation but not provided by software. Here, the size of the
windows controls the glass-to-wall ratio. Moreover, the angle to project north is used to
represent the façade orientation. By doing this, we avoid the rotation of the model during the
optimization process.
Wall Heat Transfer Coefficient</p>
      <p>0.0876 0.1802 0.0876 0.1802</p>
      <p>Glass Heat Transfer Coefficient 1.9873 6.7018
With these initial range of parameters, a sensitivity analysis was performed for each of the
critical building parts. A filtered range of parameters has been obtained and applied to the final
optimization to get the final range as also shown in Table 1. It can be observed that the optimal
results for discrete variables such as glass transmittance and glass heat transfer coefficient are
reasonable as the balanced point of higher daylight quality and lower SHG and TL are achieved.
However, some of the continuous variables such as wall surface specularity directly converged
to the upper bound, which seems unreasonable. This is because the discrete variables are all
defined by the type of architectural elements, so they work like passive variables and cannot
change independently. In this case, for the glass transmittance reaching the upper bound, the
glass heat transfer coefficient reaching the lower bound would be the theoretically best solution.
However, the type of glass satisfying these two values does not exist in practice. This is one
strong advantage of directly extracting information from architectural models employing real
objects, so this unrealistic scenario can be avoided. However, the obtained results for
independent continuous variables were relatively valuable, since they are not limited by the
practical condition. The error mentioned above could happen if all the parameters are simply
used as semantic values to run an analysis without connecting them with real physical objects.
The drawback of this result actually proved the advantage of directly extracting information
from architectural models.
Finally, the values of objectives during convergence are plotted as DF versus TF and the optimal
trade-off between these two objectives has been found. After applying a further weighted-sum
optimization, figures of NAO for three critical parts are plotted. One example of the running
result for the middle part of the building is shown in Figure 6.
After getting the optimal range of parameters, one combination of parameters is selected and
input into the custom node for generating an optimal building.</p>
    </sec>
    <sec id="sec-9">
      <title>5. Conclusion</title>
      <p>In this paper, we proposed a design loop for multi-objective optimization of building
sustainability based on the parametric virtual design approach. The methodology of this concept
includes setting an automatic link between architectural models and sustainability assessment
tools, optimization of sustainable performance based on the genetic algorithm and finally an
automatic update of the architectural model with the determined optimal parameters. The case
study presented in this paper, demonstrated the potential of the proposed methods to perform a
fully automated design-through-analysis-optimization workflow. Furthermore, from the
presented case study, the behavior of passive and independent variables revealed the potential
risk obtaining unrealistic results when extracting the geometry and semantics from non-related
sources, if the input variables are not carefully constrained. Often these constraining steps
require professional knowledge and experience in environmental engineering. However, this
problem could be solved by developing a method or a tool for direct use of all the information
(geometry and semantics) from a single architectural project. By doing this, all the parameters,
which are passive variables, are defined by a real physical property of the construction material
and the element size. Even though there are limitations in the existing practical
implementations, this research area has a large potential for further improvement. Firstly, the
enrichment of the environmental database of the standard design tools considering the
requirements of different analysis tools, will allow for more user-friendly and efficient analysis.
Secondly, the programming environment associated with these tools (open-source and
developer-friendly) offers the advantage for exploration of different optimization strategies,
since they show different performance for different optimization problems and conditions
(Kheiri, 2018). Finally, the strong functionality of parametric virtual design techniques could
significantly decrease the workload for the manual update of the design models. However,
developing a general algorithm for most common geometrical transformations of building
forms that are compatible with the optimization tools, poses a significant challenge for future
research. In conclusion, the design approach proposed in this paper has demonstrated that an
automated design optimization loop based on a real design model is feasible and also offers a
huge space for further developments.
Asl, M. R., Stoupine, A. and Zarrinmehr, S. (2015). Optimo: A BIM-based Multi-Objective Optimization Tool
Utilizing Visual Programming for High-Performance Building Design. BIM - Applied, Volume 1, pp. 673-682.
Asl, M. R., Zarrinmehr, S., Bergin, M. and Yan, W. (2015). BPOpt: A framework for BIM-based performance
optimization. Energy and Buildings. Energy and Buildings, Volume 108, pp. 401-412.</p>
      <p>
        <xref ref-type="bibr" rid="ref1">AUTODESK. Autodesk Insight (2019</xref>
        ). [Online] Available at: https://insight360.autodesk.com/oneenergy
[Accessed 11 April 2019].
      </p>
      <p>
        <xref ref-type="bibr" rid="ref1">AUTODESK. Autodesk Revit (2019</xref>
        ). [Online] Available at: https://www.autodesk.co.uk/ [Accessed 11 April
2019].
      </p>
      <p>
        <xref ref-type="bibr" rid="ref1">AUTODESK. Autodesk Revit (2019</xref>
        ). Create the Energy Model. [Online]
https://knowledge.autodesk.com/support/revit-products/learnexplore/caas/CloudHelp/cloudhelp/2019/ENU/Revit-Analyze/files/GUID-A77EF18D-0EB8-47D3-A6536370D9F1374A-htm.html [Accessed 11 April 2019].
      </p>
      <p>Available</p>
      <p>Baker, N. (2017). Modelling and Analysis of Daylight, Solar Heat Gains and Thermal Losses to Inform the Early
Stage of the Architectural Process. Stockholm: KTH Royal Institue of Tcchnology.</p>
      <p>Boeykens, S. (2012). Bridging building information modeling and parametric design. In: eWork and eBusiness in
Architrcutral, Engineering and Construction: ECPPM . London: Taylor and Francis Group, pp. 453-458.
Deb, K., Pratap, A., Agarwal, S. and Meyarivan, T. (2002). A Fast and Elitist Multiobjective Genetic Algorithm:
NSGA-II. IEEE TRANSACTIONS ON EVOLUTIONARY COMPUTATION, 6(2), pp. 182-197.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          <string-name>
            <surname>AUTODESK. Green Building Studio</surname>
          </string-name>
          (
          <year>2019</year>
          ). [Online] Available at: https://gbs.autodesk.com/GBS/ [Accessed 11
          <string-name>
            <surname>April</surname>
          </string-name>
          <year>2019</year>
          ].
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>