Abstract
The microclimate affects the quality and efficiency of outdoor spaces of campuses, especially in the cold regions of China. In this paper, we propose a multiobjective optimization method to improve the thermal comfort of the outdoor environment of university campuses in severe cold regions. We used morphology data from 41 universities in the cold region of China to create a layout prototype of a campus cluster. Multiobjective optimization was used, and the effects of sunlight, solar radiation, and wind on the outdoor thermal comfort in winter were considered. A parameterized platform was established for the multiobjective optimization of the microclimate of the simplified model of the campus. A multiobjective optimization based on an evolutionary algorithm was used to obtain 108 groups of nondominated solutions. The optimum outdoor microclimate of the campus was obtained at a building density of 0.21–0.23, a plot ratio of 1.51–1.88, and a road width of 11–14 m. We recommend that buildings are designed based on the wind direction in winter and that the space between buildings is increased.
1. Introduction
The microclimate affects the quality and efficiency of open spaces on university campuses [1] and can be significantly improved by appropriate planning [2–4]. In severe cold regions, the thermal comfort of outdoor spaces on campuses is affected by multiple climatic factors, such as sunlight, solar radiation, and wind speed. Faulty planning can reduce the thermal comfort of the users. Therefore, in this study, we consider the relationship between several climatic factors and the layout of the university campus and propose a design strategy for campuses to optimize the outdoor microclimate. The strategy improves the use efficiency of campus open spaces in cold regions and promotes the sustainable development of the campus environment.
Extensive research has been conducted on the microclimate of urban blocks and street canyons [5, 6]. The users’ behavior at universities is similar to people in urban areas; however, the layout of a campus differs from that of urban blocks. The campus layout is characterized by many open spaces and large building clusters. The spatial morphology of the university campus influences the microclimate of the campus [7, 8], and architectural clusters have the greatest influence on the duration of direct solar radiation and the mean radiation temperature [9]. Unlike outdoor spaces on campuses, enclosed areas have good performance for reducing the heating energy consumption and enhancing natural ventilation and cooling [10]. Moreover, the structures on campuses have a significant impact on outdoor solar radiation and wind [11]. A study of university campuses in cold regions has shown that when the outdoor thermal comfort and heat dissipation of buildings in summer is considered in the design, the time when the outdoor temperature is perceived as uncomfortable is reduced by 25%, and the heat dissipation is reduced by 5% [12]. The main factors affecting the outdoor microclimate include sunlight, temperature, humidity, wind speed, surface material, and the sky view factor [13–15].
The urban climate map (UCM) system is a method to analyze the relationship between urban climate and urban morphology [16]. The concept of the UCM was first proposed by German scholars [17] and has been subsequently used in urban planning and local climate improvement [18]. Smith et al. found that a higher density of buildings and population exacerbated urban heat stresses [19]. Ren et al. analyzed the urban microclimate of Arnhem city using spatial data and found that land use was the most significant factor affecting the thermal load [20]. The UCM method is well suited to estimate the thermal environment of different regions according to the geographical characteristics and building characteristics in urban areas.
The optimization of the outdoor microclimate is affected by various climatic factors. The optimization of a single climatic factor may have adverse effects on other climatic factors. Therefore, the optimization of multiple climate factors has practical significance to improve the outdoor environment. Thus, the multiobjective optimization theory based on a genetic algorithm was used in this study, and we reviewed the literature on this topic [21]. Pareto proposed the Pareto solution set to describe the optimal state of resource allocation based on multiple objectives; the solution set is commonly depicted as a point distribution on a curve or surface [22]. The Pareto frontier can be solved using the strength Pareto evolutionary algorithm (SPEA) proposed by Zitzler [23]. A method of finding the Pareto optimal solution set was proposed by Wortmann et al. for use in architectural design [24]. Cui et al. established a microclimate optimization framework for the optimal design of an urban outdoor space in the cold regions in China [25]. Du et al. proposed a multistage optimization process to optimize the wind speed distribution and the thermal environment in urban canyons [26]. Waibel et al. presented a cosimulation framework for the optimization of building energy demand and solar potential [27].
The purpose of this study is to improve the outdoor microclimate on a university campus using campus layout design, which is achieved as follows:(1)Morphological classification is used to describe the characteristics of university campuses in a cold region in China (Heilongjiang Province, Jilin Province, Liaoning Province) and develop a layout prototype of campus clusters(2)The morphological characteristics of the campus are simplified, and a multiobjective optimization method combined with a genetic algorithm is used to conduct parametric linkage analysis and optimize the morphology(3)According to the optimized solution set and the established quantitative equation, we suggest planning strategies for campus layouts to guide sustainable environmental planning on campuses
2. Methodology
2.1. Principle of Multiobjective Optimization
The influence on the outdoor microclimate differs for different spatial morphological parameters. For example, increasing the height (H) to width (W) ratio of streets can reduce the air temperature and solar radiation and improve the thermal environment in summer, but excessive wind speed will result in thermal discomfort [28]. In addition, a morphological layout that is suitable for the summer thermal environment may not have advantages in the winter [29]. In cases where different parameters have to be optimized simultaneously, Pareto optimization is commonly used. The Pareto optimal solution set, which is also known as the nondominated solution set, provides a balance of all optimization objectives. The distribution of nondominated solutions is called the Pareto frontier. The goal of multiobjective optimization is to determine the Pareto frontier. The mathematical description is as follows:
As shown in equation (1), the variable feasible region is S and the corresponding target feasible region is Z = f (S). For a viable solution ∈S, it should be , so is the optimal solution for multiobjective optimization. If x ∈ S does not exist which makes , is called an efficient solution, and is the Pareto optimal solution of multiobjective programming.
2.2. Simulation-Based Optimization Workflow
For the design process, which generally uses simulation for performance prediction and optimization, the separation of the optimization and performance prediction processes often results in repetition and an increase in time cost. Thus, the simulation-based optimization workflow was established [30]. During the optimization process, some performance factors are often sacrificed to improve the performance of others, and the direction and effect of optimization are uncertain. The multiobjective optimization workflow allows for determining the complex competition among various performance factors and provides user-defined optimization [31].
The multiobjective optimization of the microclimate of the university campus cluster includes four steps, which are parametric modeling of the campus cluster, determination of the optimization objectives, calculation of the decision variables, and multiobjective optimization. We used a multiobjective optimization platform based on Rhinoceros and Grasshopper; the four steps are shown in Figure 1.(1)Determine the morphological parameters to be optimized based on the results of the morphology of the university campuses in the cold regions. Establish a cluster prototype parametric model that allows for the real-time adjustment of the parameters.(2)Determine the optimization objectives according to the characteristics of the cold climate, such as the wind and thermal environment of the outdoor space.(3)Select the decision variables according to the optimization objectives. The sunlight and solar radiation values of the outdoor space in the simulation were obtained using the Ladybug toolset [32]. The wind speed simulation was performed by importing the OpenFoam software [33] by the Butterfly toolset.(4)Load the multiobjective evolutionary algorithm using the Octopus tool, link the optimization objective variables with the parameters of the campus cluster prototype, and generate a large number of optimized cluster prototype models for the multiobjective optimization.

2.3. The Morphological Characteristics of the Campus and Its Simplification
A survey of campuses in the cold region of China showed that the campuses consist of a teaching area, accommodation area, service area, and sports area, which are arranged in clusters. In this study, the clusters of 41 campuses were classified based on the layout and scale. It was found that the clusters consisted of several types of building types, including enclosed buildings, semienclosed buildings, and multistory buildings. A grid was used to analyze the pattern. The types of buildings and shape simplification are shown in Table 1. The survey data showed that the width of the streets on the campus was in the range of 5–20 m, the grid size ranged from 45 m to 100 m, and the area of the clusters ranged from 200 m × 200 m to 400 m × 400 m. The distance between the buildings in the East-West direction was 13–20 m, and that of the North-South direction ranged from 20 m to 30 m and depended on the height of the buildings. The building density of the campus cluster was 0.2–0.5, and the building plot ratio of the clusters was significantly different, ranging from 0.8 to 5. The number of stories of the buildings on the campus was in the range of 2–24, and most ranged from 4 to 12.
The types of buildings were simplified into nine types, with a size of 100 m × 100 m. The simplified model and corresponding parameters are shown in Table 2.
Based on the survey data and simplified building units, a cluster prototype of the campus was constructed; it consisted of nine 100 m × 100 m building units, and the street width ranged from 5 m to 20 m. The nine types of campus building units and open spaces were randomly distributed at the site, resulting in 910 sets of clusters. Three cluster prototypes and the representative cases are shown in Table 3. The rapid generation of the cluster prototypes can be achieved by parametric modeling.
2.4. Optimization Objectives
Since this study focuses on adaptive planning of the campus layout in a severe cold climate, 41 university campuses were chosen in Northeast China (118°–135.2° E longitude, 38.7°–53.33°N latitude). The area has a temperate continental monsoon climate with long winters, low temperatures, and snow. The urban climate conditions in the severe cold region are shown in Table 4.
The objectives of this study were to maximize the sunlight duration and radiant heat gain of the outdoor space and minimize the wind speed to achieve the comfortable outdoor environment of the university campus. Three parameters were considered, namely, the average sunlight duration (ASLD), average solar radiation (ASR), and average wind speed (AWS). The optimization objectives are as follows:
ASLD: to increase the average sunlight hours in winter. The ASLD is the average value of sunlight hours on a cold day (ℏ) at the height of 0.9 m of the campus cluster; the unit is meter, and the values are obtained by simulation. The Octopus toolset calculates the minimum value of each objective value by default. The objective function of this parameter in the multiobjective optimization is α, α = −ℏ.
ASR: to increase the radiant heat gain of the site during the heating period (10 : 20–16 : 20). The unit of the parameter is kW·H/m2, and the values are obtained by simulation. The objective function of this parameter in the optimization is γ, γ = –I.
AWS: according to the characteristics of the Harbin urban wind field, the optimization objective is to decrease the average wind speed. The data are imported into the OpenFoam computing platform using the Butterfly toolset to simulate the wind speed distribution in the outdoor space of the cluster prototype. The prevailing wind direction in Harbin is WNW in winter, the initial wind speed is 5 m/s, the calculation grid size is 10 m × 10 m, and the wind speed is modeled at the height of 1.5 m. The objective function in the optimization is V.
2.5. Optimization Decision Variables and Constraints
The multiobjective optimization of the outdoor space microclimate of campus clusters is affected by various decision variables; the building density and plot ratio have the largest influences on the solar radiation of the outdoor space. The building height affects the solar radiation range and wind in the surrounding space, and the percentage of pavement area (PAVE) affects the ventilation performance of the external space. To simplify the optimization model, we focused on building morphology variables and cluster layout variables (Table 5).
2.5.1. Building Morphology Variables
After simplifying the types of campus buildings, nine types of buildings were extracted that cover different campus building types on university campuses in cold regions. The nine types of building modules were represented by 1–9, and 0 represented open spaces. In this manner, the layout types of different building types could be changed in the parametric module, as shown in Table 2, to create the layout of the campus cluster in the next step.
2.5.2. Cluster Layout Variables
The morphology of the university campus cluster is affected by many variables, such as the plane form, scale, and distance between buildings. Each of the nine 100 m × 100 m building areas could be assigned to 1 of the 9 building types in the parametric model to create different cluster morphologies. Thus, 910 sets of campus clusters were obtained. The building density ranged from 0 to 0.52, the plot ratio ranged from 0 to 5.62, and the width ranged from 6 m to 15 m, as shown in Table 5.
2.6. Parametric Simulation Modeling and Algorithm Settings
Based on the parametric modeling platforms Rhinoceros and Grasshopper, the campus cluster prototype model is integrated with the microclimate simulation of the outdoor space to establish a complete data flow and automatically generate the simulation model. In addition, Ladybug and Butterfly plugins were used to integrate the Grasshopper EnergyPlus Weather files (EPW), Dynamo, and OpenFoam. The optimization software platform is shown in Figure 2.

The EPW and Dynamo were used to calculate the sunlight duration and solar radiation, which allows for changing the simulation parameters, such as the building type and the road width of the cluster. The climate data of Harbin (45.77°N, 126.68°E), which is a typical city in the cold region of China, was chosen as the reference; the data were obtained from the EPW [29]. The location information was used to calculate the sunlight duration and the annual direct and scattered radiation values. The solar radiation data were imported into the simulation engine for the solar radiation calculation. The simulation-based optimization model is presented in Figure 3.

The optimization algorithm settings are listed in Table 6. Octopus provides two convergence mechanisms and three mutation mechanisms. The improved SPEA (SPEA-2) and polynomial mutation were used in the experiment to achieve good convergence and uniformity. The main parameter settings in octopus were elite ratio of 0.5, crossover probability of 0.8, mutation probability of 0.1, mutation ratio of 0.5, population size of 50, and maximum number of iterations of 40; the sunlight requirements were set as random constraints [35, 36].
3. Results and Discussion
3.1. Multiobjective Optimization of Campus Morphology
The multiobjective optimization calculation was performed for over 49 hr, and after 40 iterations, 108 nondominated solutions were obtained. As shown in Figure 4, the three axes represent the ASLD on cold days, the ASR during the heating period in winter, and the AWS of the outdoor space, respectively. Each point represents a group of solutions that correspond to a cluster prototype. The closer the objective values are to the origin, the better the solution performance is. The red-colored boxes represent the best performing solutions regarding the three optimization objectives.

(a)

(b)

(c)
As Figure 4(a) shows, each polygonal line represents a set of solutions, and the three objective values of ASLD, ASR, and AWS are mutually restricted. Figure 4(b) shows the Pareto front surface consisting of the nondominated solutions, and Figure 4(c) shows the distribution of the nondominated solutions in the coordinate system.
The spatial distribution of the nondominated solutions (Figure 5) indicates that the ASLD on cold days is positively correlated with the ASR in winter, and the relationship between the two parameters and the AWS is mutually restricted. The nondominated solutions of the ASLD range from 3.3 h to 8.07 h on cold days, and the ASR during the heating period ranges from 130.32 kW H/m2 to 235.74 kW H/m2. The AWS ranges from 1.76 m/s to 2.75 m/s (Figure 5).

(a)

(b)

(c)

(a)

(b)

(c)
Five typical cluster prototypes with nondominated solutions were analyzed (Table 7). Solution A has the optimal AWS, solution B has the optimal ASLD and ASR, and solutions C, D, and E are relatively balanced regarding the three optimization objectives. The optimization of the winter microclimate of the cluster can be achieved by changing the morphology parameters and the combination of multiple building units. Therefore, in the design, the morphological parameters should be considered to achieve a better microclimate.
3.2. Parameters of the Nondominated Solutions
The morphological parameters of the nondominated solutions were selected as the optimal value threshold (Figure 7). The building density was in the range of 0.17–0.31, and the median was 0.22. In half of the nondominated solutions, the building density was in the range of 0.21–0.23, indicating that more multiobjective optimization prototypes were obtained when using this value range as the building density of the cluster.

. (a)

(b)
The plot ratio ranged from 0.99 to 2.51, and the median was 1.63. Most of the values were in the range of 1.51–1.88, which can be used as the reference range of the parameters for cluster morphology optimization.
The average number of stories of the buildings was in the range of 3–8.29, and the median was 5.34. Half of the values in the nondominated solutions were in the range of 4.63–6.14. The optimized average number of stories of the building was between 4 and 6.
The road width of the cluster prototype ranged from 6 m to 15 m, and the median was 14 m. Most of the values were in the range of 11 m–14 m; thus, this range is suggested as the optimum range of the road width.
When the parameters of the campus cluster are within these ranges, the cluster morphology is optimized regarding the ASLD, ASR, and AWS in winter, and the microclimate of the outdoor space is optimized.
The building density (Figure 8) and road width (Figure 9) of the nondominated solution were proportional to the ASLD and ASR and inversely proportional to the AWS. The plot ratio (Figure 10) and the average number of stories (Figure 11) of the nondominated solution were inversely proportional to the ASLD and ASR and proportional to the AWS (Figure 10).

(a)

(b)

(c)

(a)

(b)

(c)

(a)

(b)

(c)

(a)

(b)

(c)
3.3. Recommendations Based on the Parameters of the Nondominated Solutions
According to the parametric analysis of the nondominated solutions, the optimal range of the building density is 0.21–0.23, and the types of qualified buildings include the following:(1)Type 7: multistory building (building density: 0.38)(2)Type 8: high-rise building (building density: 0.36)(3)Type 9: high-rise tower (building density: 0.16)
The optimal range of the plot ratio is 1.51–1.88, and the types of qualified buildings include the following:(1)Type 3: multistory semienclosed (plot ratio: 2.3)(2)Type 5: centralized (1) (plot ratio: 1.35)(3)Type 6: centralized (2) (plot ratio: 1.35)(4)Type 7: multistory building (plot ratio: 2.3), which can be used as the preferred building form (Table 2)
The simulation results of the optimal solutions showed higher wind speeds near the high-rise buildings; thus, multistory buildings should represent the main component of the campus cluster, such as type 7. The enclosed and semienclosed layout resulted in lower wind speeds of the outdoor space; when the longitudinal direction of the building was the same as the wind direction, the wind speed of the enclosed space was significantly reduced.
In the climate conditions of Harbin, the open spaces in the southeast of the campus provided advantages regarding sunlight and solar radiation, and the sunlight duration could be increased by increasing the building spacing. It is suggested to optimize the microclimate of the outdoor space by maintaining the average number of stories at 3–8.3 and the road width at 11–14 m. Green vegetation can be used to change the roughness of the surface to optimize the outdoor wind speed in winter.
3.4. Comparison of the Conditions before and after Multiobjective Optimization
3.4.1. Before Optimization
A model of Jilin Agricultural University was established using the Rhinoceros and Grasshopper platforms, and a simulation was performed with three optimization objectives: ASLD, ASR, and AWS. The results are shown in Figure 12.

(a)

(b)

(c)

(d)
Figure 12(a) shows the key points in the cluster space that were quantitatively analyzed. The sunlight duration of the cluster model on cold days was calculated using the Ladybug tool set. The ASLD was 0 to 18 h, with a mean value of 5.56 h. Figure 12(b) shows that the ASLD was low; 80% of the open space had an ASLD of 0 h to 7.2 h, and at the locations ③ and ④, the ASLD was close to 0 h. At ①, ②, and ⑥, the ASLD was between 0 h and 1.8 h. The data shows that the outdoor environment of the cluster requires improvement in terms of sunlight duration.
The ASR in winter ranged from 0 kW H/m2 to 357.6 kW H/m2, with a mean value of 162.58 kW H/m2. Figure 12(c) shows that at points ④ and ⑤, the ASR was near the minimum area, which was between 0 kW H/m2 and 143.0 kW H/m2. There was a greater distance between the buildings in the North-South direction, and the maximum value of the ASR was 321.84 kW H/m2. The distribution of the solar radiation in the cluster space was uneven, resulting in differences in the heat gain of the cluster space; this aspect requires further improvement.
The Butterfly toolset was used to simulate the wind environment of the cluster space. The initial wind speed was 5 m/s, and the dominant wind direction in winter is WSW. Figure 12(d) shows that the AWS in the cluster space is between 0 m/s and 4 m/s, with a mean value of 1.99 m/s. At points ⑦ and ⑧, the AWS was in the range of 0.58 m/s −2.48 m/s, representing the area with the lowest wind speed in the cluster space. At ①, ③, ④, and ⑤, the wind speed was higher, especially near the edge of the buildings, and reached 4 m/s, which substantially decreased the thermal comfort of the outdoor space in winter. It is necessary to adjust the morphology of the cluster and optimize the local wind environment.
3.4.2. After Optimization
The proposed design recommendations for Jilin Agricultural University take into account that the sunlight duration is short, such as at ③, ④, and ⑤, where high wind speeds of above 4 m/s were observed. The layout of the three centralized buildings in the cluster was adjusted, as shown in Figure 13(a). First, the distance between the buildings was changed from 6 m to 12 m, which increased the building spacing and the sunlight duration, as shown in Figure 13(b). Points ③, ④, and ⑤ were well optimized, and the solar radiation was also well optimized, as shown in Figure 13(c). Second, we adjusted the building form to an L-shape according to the wind direction to block the cold wind in winter, as shown in Figure 13(d). Since the average number of stories of the original building cluster was 5, which is within the optimal value range, the building height was not changed. In addition, the optimal sunlight duration of the outdoor space after optimization was 5.59 h, the sunshine radiation was 166.94 kW H/m2, and the AWS was 1.94 m/s. The optimized cluster morphology resulted in an improved microclimate, demonstrating the excellent performance of the optimization model. Moreover, the passage space was changed to a square space that provides a positive public outdoor space on campus.

. (a)

(b)

(c)

(d)
4. Conclusions
In this study, 41 university campuses in three northeastern provinces of China were surveyed, and the morphology was simplified. A prototype of the campus cluster was created with 9 types of building units, and the parametric model was established using the Rhinoceros and Grasshopper modeling method. Different building types and road widths were modeled, resulting in different building densities, plot ratios, and the average number of stories.
Three optimization objectives were used to optimize the outdoor microclimate of the campus cluster in winter, namely, maximizing the ASLD on cold days and the ASR during the heating period and minimizing the AWS. A parameterized platform was established for the optimization using the morphological parameters of the cluster prototype. The SPEA-2 was used for the multiobjective optimization, and 108 sets of nondominated solutions were obtained.
The ASLD of the open space on cold days ranged from 3.3 h to 8.07 h, the ASR was in the range of 130.32 kW H/m2–235.7 kW H/m2, and the AWS was 1.76 m/s–2.75 m/s. After multiobjective optimization of the three parameters, the building density of the cluster was 0.21–0.23, the plot ratio was 1.51–1.88, the average number of stories of the buildings was 4–6, and the street width was 11–14 m to achieve the optimum microclimate of the prototype of the campus cluster. Our design recommendations based on the optimization of the Jilin Agricultural University campus are as follows: the building shape should match the wind direction in winter, the space between buildings should be increased, the average number of stories should be 3–8.3, and the road width should be 11–14 m.
In the optimization of the campus outdoor space environment of Cold Regions University, there are many factors that have not been discussed in this study, such as orientation, vegetation, and vision, which will continue to be paid attention to in the future research.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors wish to thank members of the Architecture Creation Research Institute for the support of this work. This research was financially supported by the Beijing Postdoctoral Research Foundation (Contract no. 2021-ZZ-107) and National Natural Science Foundation of China (no. 51778008).