1. Introduction
More than half of the world’s population is concentrated in cities that occupy 2% of the Earth’s surface [
1]. Buildings in cities represent 36% of the global final energy consumption [
2]. This energy demand is always increasing, and the lack of clean energy resources to cover it is a major challenge in reducing the urban carbon footprint. Moreover, the energy consumption of the residential sector is more than three times that of the non-residential sector [
3]. As people’s demand for environmental comfort grows, the proportion of “air-conditioned” buildings continues to rise [
4]. Reducing the energy consumption of residential buildings and improving the comfort of the indoor environment are important parts of solving the urban energy crisis and improving the quality of the human settlement environment. In addition, the urban-rural gap and globalization have significantly increased food mileage [
5]; however, it is still hugely challenging to fully obtain healthy and affordable food in many cities and regions [
6]. At present, vertical farming has been adopted to address the limitations of traditional agriculture, the shrinking area of arable land, and increasing food demand [
7,
8]. However, there is still a lack of research on the potential crop yields of building facades [
9,
10].
To address environmental and energy problems, a variety of facade design methods and concepts have been proposed in the field of architecture aimed at improving indoor environmental quality [
11], significantly reducing building energy consumption [
12], and generating a positive impact on local renewable resources, buildings, and human resource consumption [
13]. In recent years, the concept of productive facade systems has been introduced by Tablada et al. (2018), in which PV shading devices and vertical farming systems are integrated into building facades. Taking into account the balance between BIPV and BIA, PV shading is usually located in the upper part of the window, while vertical farming is located in the lower part, in order to adapt to the respective system requirements and ensure the accessibility of the indoor view [
14]. In addition, the integrated PV module design and modular planting slots are used to ensure the aesthetic effect of the building façade [
15].
Simulations and experimental results have shown that this integration can help achieve sustainability goals at urban and building scales by alleviating urban energy shortages and improving urban food security while improving indoor environmental quality [
9,
16].
However, predicting the yields of PV shading and vertical farming systems is complex. Previous studies were based on empirical formulas or simulation tools; however, more recent studies have shown that the calculated results of these methods and tools often have non-negligible errors compared to the actual results [
17]. This is because these two systems are affected by multiple environmental factors, such as solar radiation, ambient temperature, shadow, soil fertility, precipitation, and wind speed.
Theoretically, the solar radiation received by a PV module gets converted into electricity at a certain conversion rate; however, under actual conditions, this yield is reduced by both the temperature and the shadow of the PV cell. Regarding prediction of module temperature, Zhen and Nobre et al. (2013) reported that, in Singapore, the k value (Ross coefficient) of PV modules can reach as much as twice that of the same type of PV module but installed in a different environment [
18]. PV fabricators provide the PV temperature coefficient of power, which governs how strongly the conversion efficiency depends on temperature. However, after a module reaches a certain temperature, the coefficient has a large error compared to the actual situation, which causes the PV temperature coefficient of power to not meet the design accuracy requirements. Regarding shadows, they not only cause electrical mismatch losses in different parts of a PV module but also overheat the panel [
19].
In terms of vegetable yield, it is difficult to obtain accurate yield data using formulas or simulation methods due to several climatic variables and planting methods. Previous studies have often estimated vegetable yields based on the planting area and incident radiation, but this rough calculation method ignores other climatic variables and external factors, such as DLI, composts, and crop rotation [
20].
As problems in the building field become more complicated, the above-mentioned common methods gradually lose their advantages in terms of time, cost, and reliability. Simultaneously, machine learning methods represented by ANNs have gradually become crucial for seeking architectural solutions [
21]. ANNs can improve the accuracy of prediction results by using more complex functional relationships for modeling. In recent years, it has been used in building performance prediction [
22,
23] and has proven to be highly efficient for rapid assessment of building performance [
24,
25]. In particular, PV power yield prediction models have been proposed based on ANNs [
26,
27]. A multilayer perceptron (MLP)-based ANN architecture has been used to predict a PV plant power yield using three input elements [
28]. An ANN for PV panels has been developed, which uses the root mean square error (RMSE) method to evaluate the prediction quality of the ANN [
29]. Several commonly used input elements, such as global horizontal irradiance (GHI), azimuth and altitude angles, air temperature, module tilt angle and surface temperature, and date [
30,
31,
32], are mainly aimed at predicting the yield under environmental conditions such as PV power stations and rooftop PV. However, there is insufficient research on the prediction of PV yield under building facade installation conditions, especially with regard to solar PV shading. This problem becomes more apparent with the promotion of integrated photovoltaics.
Regarding the prediction of farming yield, Abrougui et al. (2019) used an ANN to predict organic potato yield [
33]. Their ANN had greater potential to estimate yield compared to a multiple linear regression model. Saad et al. (2009) proposed an ANN to predict rice yield [
34], and Pantazi et al. (2016) applied an ANN to predict wheat yield using soil and remote sensing vegetation indices as input parameters [
35]. The input elements typically include tillage systems and soil properties [
33]. Although ANNs are widely used for predicting farming yield, research on vertical farming is still lacking.
Several tools, such as MATLAB and Grasshopper, have a toolbox for ANNs. For Grasshopper, four known plugins adopt ANN algorithms: Dodo, Crow, Owl, and Lunchbox [
36,
37]. The 3D modeling software Rhinoceros 6.0 and its Grasshopper plug-in have been used to control geometric parameters and analyze indoor environmental performance [
38].
This study first compared the predictive modeling methods using ANN with different activation functions, proposed a relatively fast and reliable method for assisting productive facade design, and then explored the application of ANN predictions in practice to evaluate productive facade yield performance. The overall prediction process, shadow analysis, and solar radiation analysis were performed using the Grasshopper plug-in for the Rhinoceros platform. In addition, the indoor daylight and thermal environment of a residential unit were compared using this platform [
38,
39].
2. Materials and Methods
To obtain reliable PV electricity yield and vertical farming vegetable yield data, this study combined an ANN trained and tested by experimental data (given in
Appendix A) from the NUS-CDL Tropical Technologies Laboratory (T
2 Lab) of the National University of Singapore (NUS) [
27,
40]. The experimental data included the PV shading electricity yield (monocrystalline silicon, type: assembled by the Solar Energy Research Institute of Singapore/CIGS thin film, type: Solar Frontier SFL85-D) and vertical farming vegetable yield data (lettuce planting). Details of the T
2 Lab PV shading system are shown in
Table 1 and
Figure 1, and the performance indicators of the T
2 Lab envelope are shown in
Table 2.
The proposed ANN was trained using experimental data to obtain a predictive model that correlated the environmental and yield elements.
Figure 2 depicts two predictive models: Predictive model 1 (environmental and PV shading electricity yield elements) and Predictive model 2 (environmental and vertical farming vegetable yield elements). Through predictive modeling, this paper presents the yield on the southern productive facade of a typical public housing building (point block) in Singapore and analyzes indoor thermal and visual comfort.
The experimental data were from April 2019 to May 2019 (for PV shading,
Figure 3) and from December 2018 to June 2019 (for vertical farming,
Figure 4). The recorded data included: ① GHI (type: SPN1-A1815, W/m
2); ② ambient temperature (type: RM Young 92000, °C); ③ wind speed (type: RM Young 05103, m/s); ④ PV module’s temperature (type: RTD Pt 100, °C); ⑤ incident solar radiation (type: IMT Si-420TC, W/m
2); ⑥ output electricity recorded by a maximum power point tracking controller (type: BlueSolar charge controller MPPT, Wh/m
2); ⑦ vertical farming data including the photosynthetic photon flux density (PPFD, type: Quantum PAR Meter); ⑧ average dry weight of each plant in each planting cycle. The data collection frequency is once per minute (①–⑦).
In this study, two solutions were adopted to reduce the unreliable training data caused by the impact of environmental shadows on the PV module, which causes the hotspot effect. First, the uppermost PV module was selected to avoid shading between the PV modules. For example, there were nine monocrystalline silicon PV modules (assembled by the Solar Energy Research Institute of Singapore) facing east and west in the laboratory, and only three of the uppermost modules (module tilt angle: 40°) were selected. Second, we prevented the laboratory from blocking the PV modules. The CIGS thin-film PV modules (Type: Solar Frontier SFL85-D) in the north and south directions were all single-layer (module tilt angle: 20°), but the sun was located north of Singapore during the collection of the experimental data; therefore, only the north PV module was selected. Furthermore, the east, west, and north modules intercepted the experimental data from 6 a.m.–12 a.m. (6 h), 12 a.m.–6 p.m. (6 h), and 6 a.m.–6 p.m. (12 h). In this study, the experimental output power data of each PV module was converted into output electricity per unit area (Wh/m2). The predictive ANN was only applied to single-layer monocrystalline silicon and single-layer thin-film PV modules, owing to the aforementioned deliberate selection of experimental data. In addition, the electricity yield of PV shading was predicted for a variety of module angles because the predictive ANN was trained with substantial data that covered a significantly broad range.
The predictive model aimed to connect the environmental elements (independent variables) with the yield elements (dependent variables). This type of problem is mainly addressed in machine learning through supervised learning in training methods, specifically through an ANN in regression tasks. In this study, the back-propagation algorithm was used to train an ANN with an MLP structure [
30]. Two commonly used activation functions (the Bipolar Sigmoid and SoftPlus functions) in the back-propagation algorithm were compared to select the optimal function, and the multiple linear regression model of machine learning was chosen as a reference. In addition, ANN evaluations usually include training and testing errors. Training error expresses an ANNs quality of fit to the training data. If the training error is too large, the ANN has insufficient knowledge of the training data characteristics. Otherwise, the ANN has overlearned the characteristics of the training data if the training error is too small. Therefore, maintaining appropriate training errors is crucial for evaluating the quality of an ANN. In contrast to the training error, the test error characterizes the generalization ability of an ANN. In practical applications, the test error must be as small as possible.
The predictive models were calculated using the Grasshopper plug-in Dodo [
23,
41] and Rhinoceros 3D software [
42]. Dodo was created by Lorenzo Greco and features an extensive set of components for ANNs [
23]. Rhinoceros and Grasshopper provide a familiar operating platform for designers to establish predictive models and analyze the environmental performance of residential units.
2.1. Predictive Model 1: Environmental and PV Shading Electricity Yield Elements
Based on an ANN combined with the shadow loss coefficient, a predictive model was developed for environmental and PV shading electricity yield elements. The model contains two steps: the first establishes the ANN, and the second couples shadow influences. The second was training the ANN using the experimental data from the T2 Lab from April to May 2019. The input variables of the ANN were initially set to five environmental elements (independent variables: GHI on the roof, ambient temperature, wind speed, PV module temperature, and PV module incident solar radiation), and the output variable was only the PV shading electricity yield element (dependent variable: PV module output electricity per unit area (Wh/m2). In the second step, according to the mechanism of the PV module being affected by the shadow and the common methods to reduce shadow influence, the shadow influence was converted into the shadow loss coefficient, which was then coupled with the ANN model from the first step.
2.1.1. Correlation Analysis
Using a monocrystalline silicon PV module as an example, the Pearson correlation coefficient was used to analyze the correlation between the five environmental elements (independent variables) and one PV electricity yield element (dependent variable). We analyzed the experimental data through SPSS 22.0 from April to May 2019 containing 65,880 sets (3 modules × 61 d × 6 h × 60 min) from three monocrystalline silicon PV modules placed in the east and west orientations of the T
2 Lab. The results are listed in
Table 3. The correlation coefficients were all greater than 0.4 [
43] except “wind speed”, indicating a moderate or strong positive correlation between the independent and dependent variables. However, the PV module temperature and wind speed were excluded from the variable types after considering their time cost, feasibility, and relatively low correlation coefficients for the five independent variables. Finally, each set of experimental data contained three types of environmental elements (independent variables: GHI on the roof, ambient temperature, and PV module incident solar radiation) and one type of PV shading electricity yield element (dependent variable: PV module output electricity per unit area).
2.1.2. First Step: Building the ANN
ANN Training
The Latin cube sampling method was used to divide the actual experimental data into training (70%) and validation (30%) data [
44]. Each set of training data and validation data contained three independent variables: GHI on the roof, ambient temperature, PV module incident solar radiation, and one dependent variable: PV module output electricity per unit area. The categories and contents of each set of data are listed in
Appendix A Table A3 and
Table A4.
The root mean square error (RMSE) and mean absolute error (MAE) were used to analyze the training and test errors of different ANNs [
28,
45]. The other ANN parameter settings are listed in
Table 4.
① The experimental data of the monocrystalline silicon PV module were recorded by three PV modules in the east and west directions of the T
2 Lab of NUS from April to May 2019 for a total of 61 days. There were 65,880 sets of data (3 modules × 61 d × 6 h × 60 min), of which 46,116 and 19,764 sets were used as training and test data, respectively. (Due to space limitations, only part of the data are shown in
Appendix A Table A3).
② The source of the experimental data for thin-film PV modules was the same as that for monocrystalline silicon PV modules. The difference was that one thin-film cell was facing north, and there were 43,920 sets of data (1 module × 61 d × 12 h × 60 min), of which 30,744 and 13,176 sets were used as training and test data, respectively. (Due to space limitations, only part of the data are shown in
Appendix A Table A4).
ANN Evaluation
The RMSE is the square root of the average of the sum of squares of each data deviation from the true value. Its expression is
where
is the actual experimental value, and
is the predicted value of the ANN.
The MAE is the average value of the absolute error, which can better reflect the true situation of the predicted value error. It is expressed as
- (1)
Evaluation of the ANN for monocrystalline silicon PV modules.
As shown in
Figure 5 and
Table 5, the training and test errors of the three models (two ANNs and one multiple linear regression model) were relatively suitable. Therefore, there were no overfitting or underfitting situations in the model. Further analysis showed that the overall prediction quality was higher for the ANN than that for multiple linear regression. In terms of the activation function for the ANN, SoftPlus was more effective than Bipolar Sigmoid when considering RMSE; however, for MAE, the former was less effective than the latter.
Compared to MAEs, RMSEs give more weight to the abnormal points in the data. Therefore, the abnormal points in the data cause the RMSE to shift towards the abnormal points at the expense of other data points, thereby reducing the overall prediction quality of an ANN. Specifically, for MAE, the training and test errors of the SoftPlus activation function were 5.8% and 4.4% higher than those of the Bipolar Sigmoid function, respectively. For RMSE, the training and test errors of SoftPlus were 6.1% and 6.0% lower than those of the Bipolar Sigmoid function, respectively. Therefore, the prediction results of the ANN using the SoftPlus activation function had fewer abnormal points than those using the Bipolar Sigmoid, which is crucial. We also showed that the ANN using the SoftPlus activation function had a better prediction quality.
Figure 6 shows the comparison between the predicted and experimental results for evaluation. Considering fitting, the fit for the ANN model using the SoftPlus activation function was highest (R
2 = 0.66 for Bipolar Sigmoid; R
2 = 0.69 for SoftPlus; and R
2 = 0.60 for multiple linear regression) [
46]. Finally, the ANN with the SoftPlus function was selected as the ANN module for the monocrystalline silicon PV.
- (2)
Evaluation of the ANN for thin-film PV module
Similar to the evaluation of the monocrystalline silicon PV module,
Figure 7 and
Table 6 show that the training and test errors of the three models (two ANNs and one multiple linear regression model) were relatively suitable for the thin-film PV module. There was also no overfitting or underfitting of the model. Further analysis showed that the ANN had smaller errors and clear advantages over multiple linear regression. In addition, the MAE and RMSE values of the ANN using the Bipolar Sigmoid activation function were smaller than those for the SoftPlus activation function. Moreover,
Figure 8 shows that the ANN prediction results with the Bipolar Sigmoid activation function fit the experimental data better than the other two (R
2 = 0.88 for Bipolar Sigmoid; R
2 = 0.78 for SoftPlus; and R
2 = 0.52 for multiple linear regression) [
46]. Therefore, the ANN with the Bipolar Sigmoid activation function was selected as the ANN for the thin-film PV modules.
Establishment of the ANN
According to the model evaluation, the ANN with the SoftPlus activation and Bipolar Sigmoid activation functions were selected as the prediction models for the yield of the monocrystalline silicon and thin-film PV modules, respectively. The specific information about the models is listed in
Table 7.
2.1.3. The Second Step: Integrated Shadow Loss Coefficient
Considering the different characteristics of the two types of PV modules (monocrystalline silicon and thin-film cells), as well as the calculation complexity and time cost, the shadow loss coefficient is defined as
where
is the shadow loss coefficient,
(based on Rhino 3D shadow analysis) is the area of the shaded part of the PV module, and
is the area of the PV module.
2.2. Predictive Model 2—Environmental and Vertical Farming Vegetable Yield Elements
The experimental data were collected from six consecutive rounds of lettuce planting data (21 groups in total) during a half-year period from the T
2 Lab of NUS [
47] (
Appendix A Table A1 and
Table A2). Three factors were selected as the environmental elements (independent variables: average DLI of the planting cycle; time proportion of suitable illumination, which means the proportion of the total time that the PPFD received by the lettuce was in a suitable interval during the planting cycle; and soil fertility). Additionally, we selected the vertical farming vegetable yield element (dependent variable: average dry weight per plant). The experimental data were used to train the ANN and establish Predictive Model 2 (
Appendix A Table A5).
2.2.1. Correlation Analysis
The Pearson correlation coefficient method was used to analyze the three environmental elements (independent variables) and one vertical farming vegetable yield element (dependent variable).
Table 8 illustrates that the correlation coefficients between the environmental and vertical farming vegetable yield elements were all greater than 0.4 [
43]; therefore, there is a moderate correlation, which can be used as training data to train the ANN.
2.2.2. Building the ANN
There are three steps in establishing the prediction model: training, evaluation, and establishing ANN.
ANN Training
The experimental data processing and model evaluation methods were identical to those used for the PV section. There were 21 sets of training data, and each set of data contained three types of environmental elements: ① average DLI during the planting cycle (mol/m
2/d), ② time proportion of suitable illumination, and ③ soil fertility. In addition, there was one yield element for the vertical farming system, the average dry weight per plant (g). The mean ANN parameter settings are listed in
Table 9.
ANN Evaluation
Figure 9 and
Table 10 show that the training and test errors of the three models were relatively suitable. The model training results were similar to those of the PV section. Further analysis showed that the overall prediction quality of the ANN was better than that of the multiple linear regression. Specifically, considering the MAE, the training error of the Bipolar Sigmoid activation function was 2.6% higher than that of SoftPlus, but the test error was 44.8% lower. For the RMSE, the training error of the Bipolar Sigmoid was 19% higher than that of SoftPlus, but the test error was 55.1% lower. To further evaluate the quality of the prediction models, the prediction results were compared to the experimental results (
Figure 10), and it can be seen from the figure that the ANN prediction results using the Bipolar Sigmoid activation function fit the experimental data better than the other two (R
2 = 0.63 for Bipolar Sigmoid; R
2 = 0.61 for SoftPlus; and R
2 = 0.47 for multiple linear regression [
46] (The larger the R
2-value is, the closer the curve fits to the diagonal, the better the model’s prediction results).
Establishing ANN
According to the results of the model evaluation, the ANN with the Bipolar Sigmoid activation function was selected as the predictive model for the yield of vertical farming. The specific information about the model is listed in
Table 11.
2.3. Application of Predictive Models
As illustrated in
Figure 11, the research objects were selected from the Arcadia community in Punggol, Singapore. Two residential units were selected in the HDB building 289D (18 floors): one on the high zone (15th floor as the representative floor) and the other on the middle zone (9th floor as the representative floor). The building is a typical contemporary block high-rise residential building [
48]. Considering the layout of the surrounding area [
49], it is one of seven high-rise residential buildings rotated 30° from the north. The HDB 289D building is 58 m wide and 24 m deep. The spacing between buildings is 25–30 m in the front-to-back direction and 14 m in the left-to-right direction. In this study, we selected the two representative floors and discussed the yield performance of the productive facade on the south-southwest (SSW) facade in a real city context.
According to previous research, the size of the residential unit was set at 3.2 m wide, 4.5 m deep, and 2.8 m high [
50,
51] (
Figure 12). The prototype of the productive facade reflected the layout of the T
2 Lab [
52]. The upper part was a PV shading system (one-layer monocrystalline silicon or thin-film PV shading modules) with a tilt angle of 20°, considering the balance between the PV electricity yield and indoor light environment [
9]. The area of the PV module was 2.56 m
2 (0.8 m × 3.2 m), with the 0.8 m width of the PV module determined by the arrangement of four rows of 0.2 m PV cells to match the mainstream cell size development trend. The length of 3.2 m was the same as the width of the residential unit to maximize the potential of solar energy and improve the indoor environmental quality. Since the ANN prediction unit for the PV shading electricity yield was Wh/m
2, the final yield used the ANN prediction result multiplied by the area of the PV module. The lower part was a vertical farming system with three rows of planting units that allow 36 vegetables to grow, and the interval between each row of planting units was 0.2 m to meet the growing space requirements of vegetables. The parapet behind the vertical farming system was replaced by a customized opening to allow access to the crops [
16]. By installing a 3-layer planting unit on a rail that can slide up and down, the planting units can be maintained more conveniently. However, to simplify the model, an opaque wall was considered, as in the case of the T
2 lab testbed. The parameter settings of the residential unit were set according to the GreenMark [
52], and the layout parameters of the building planning are detailed in
Figure 12 and
Figure 13.
3. Results
First of all, using experimental data training and testing ANN to generate prediction results, which, when compared to the experimental results, testify that the prediction models using the Bipolar Sigmoid activation function fit the experimental data better (R2 = 0.88 for electricity yield and R2 = 0.63 for vegetable yields) than SoftPlus (R2 = 0.78 for electricity yield and R2 = 0.61 for vegetable yields) and MLR (R2 = 0.52 for electricity yield and R2 = 0.47 for vegetable yields). Therefore, the ANN-used Bipolar Sigmoid activation function is selected as the productive façade yield prediction model.
Then, this prediction model was applied in a case to evaluate the electricity and vegetable yields of the productive facade in the high and middle zones of a building facade facing the SSW. Using the Rhinoceros with Grasshopper plug-in platform, the prediction model obtained the annual PV shading electricity yield and annual vertical farming vegetable yield in the high and middle zones of the residential model (
Table 12). In general, the yields in the high zone were higher than those in the middle zone. The yield of the monocrystalline silicon PV module was higher than that of the thin-film PV module.
For the PV electricity yield, the annual electricity yield of the monocrystalline silicon and thin-film PV modules was 274.3 and 237.8 kWh in the high zone and 251.4 and 218.3 kWh in the middle zone, respectively. According to the Singapore Bureau of Statistics, the average monthly household electricity consumption of public housing in 2020 was 395.6 kWh (for four-bedroom residences) [
53]. According to the above data, the output of the productive facade can meet 23.1% (using a monocrystalline silicon PV module) or 20.0% (using a thin-film PV module) of the annual electricity demand from residential units composed of four-room residential units in the high zone. The middle zone can meet 21.2% (using a monocrystalline silicon PV module) or 18.4% (using a thin-film PV module) of the annual electricity demand.
In terms of vertical farming vegetable yield, the annual yield of a single residential unit in the high zone was 19.8 kg dry weight (calculated based on the average value of 73–95% of the moisture in leafy vegetables; 19.8 kg dry weight is about 99 kg of fresh vegetable weight [
54]), and a single residential unit’s annual yield in the middle zone was 17.1 kg dry weight (approximately 85.5 kg fresh vegetable weight). In 2014, the per capita consumption of leafy vegetables in Singapore was 96 kg [
55]. For the common four-room residential units in Singapore, the annual vegetable self-sufficiency rate of a four-person household can reach 103% (high zone) or 89.1% (middle zone) after the application of a productive facade. This is crucial to Singapore, where the current vegetable self-sufficiency rate is only 13%, and almost 90% of all vegetables are imported [
56].
4. Discussion and Conclusions
Based on the experimental data from the T2 Lab at NUS in Singapore, this paper developed a performance prediction method applied to ANN and presents the application of this method and the results of predicting the electricity and vegetable yield of a productive facade in a typical residential building in Singapore. First, based on different activation functions, an ANN was trained using the experimental data, and then the prediction models involving environmental and yield elements were built. Finally, due to the better performance of the other two activation functions, the yield prediction method of the productive facade was completed using ANN with the Bipolar Sigmoid activation function, which integrated the PV shading and vertical farming systems.
PV shading and vertical farming are increasingly being applied to architectural facade design and practice, but predicting PV electric energy and vertical farming vegetable yield is complex owing to the multitude of influencing parameters.
It should be noted that the integration of PV shading and vertical farming into building facades is a challenging task because of the need to consider compatibility between systems, urban shading effects, optimal design, architectural aesthetics, and so on. In the context of multi-objective decision-making, the optimal solution is often limited by the inability to obtain data on PV electricity yield and building facade vegetable output because these two are difficult to use with traditional methods such as simulation or formula derivation quickly and accurately. With PV shading and vertical farming being increasingly applied to building facade design and practice, the novel prediction method proposed in this research is helpful to solve the performance evaluation problem in the context of complex influencing factors. The novel and predictive ANN, based on the Rhinoceros and Grasshopper plug-in platforms, proposed in this paper allows designers to predict the yield performance of productive facades with similar layouts of PV shading and vertical farming systems. A relatively reliable method for predicting the output of productive facade resources was obtained in different tropical urban contexts. In addition, researchers can obtain satisfactory design prototypes by comparing the performance of different productive facade prototypes to aid in the design decision-making process.
Future research will focus on different weather conditions, building types, materials, and crops, and improving the accuracy of predictions. Specifically, the prediction process and results of this study require further analysis and verification. As the training data, namely the experimental data, were obtained from Singapore, the prediction model is currently only applicable in regions close to the equator. Moreover, considering the prediction process, different PV installation environments affect the temperature of the PV modules differently, which in turn affects the PV electricity yield. This paper discusses the electricity yield of PV shading in the facade installation environment of residential buildings only. In this research, we found that the prediction results of machine learning represented by the ANN might have “jumping points”, i.e., a small number of prediction results had large errors. Therefore, in this study, for the same situation, the prediction model was used to conduct multiple rounds of prediction, and the data with large errors were excluded by comparing the prediction results, but the time cost of the prediction process was also increased to a certain extent. Therefore, it is necessary to find ways to improve the predictive model from the perspective of tool iteration or algorithm updates.
CIGS thin-film was selected for this study, and more thin-film PV module types will be conducted. Finally, considering crop yield prediction, the experiment focused on shade-loving crops (lettuce) that make full use of low-light conditions. Therefore, the yield prediction of other crop types is another avenue for future research.
The next step is to examine how predictive models can be packaged and integrated into the current building facade design process. Through the application of measurement and satisfaction surveys, the prediction model is constantly improved, and a shared platform is built among designers, users, and manufacturers to achieve a balance between performance, user satisfaction, cost, and architectural aesthetics.
Energy and food security have always been bottlenecks in sustainable urban development. Through the rational use of solar energy resources, the productive facade, which integrates PV shading and vertical farming systems, provides feasible solutions for improving the quality of urban human settlements and promoting the development of urban green energy and food security.