Application of a Numerical Simulation to the Estimation of Wind Loads on Photovoltaic Panels Installed Parallel to Sloped Roofs of Residential Houses

Many residential houses with sloped roofs are equipped with photovoltaic (PV) systems. In Japan, PV systems are generally designed based on JIS C 8955, which specifies wind force coefficients for designing PV panels. However, no specification is provided to the PV panels installed near the roof edges where high suctions are induced. When installing PV panels in such high-suction zones, we need to evaluate the wind loads on the PV panels appropriately, usually by performing a wind tunnel experiment. However, it is difficult to make wind tunnel models of PV panels with the same geometric scale as that for the building, e.g., 1/100, because the thickness of PV panels and the distance between PV panels and a roof are both several centimeters. Therefore, in the present paper a numerical simulation is applied to the estimation of pressures in the space between the lower surface of PV panels and the roof surface, called “layer pressures”, using the unsteady Bernoulli equation and the time histories of external pressure coefficients obtained from a wind tunnel experiment. An assumption of the weak compressibility of the air and an adiabatic condition is made for predicting the layer pressures from the flow speed through the gaps. The simulation method is validated by a wind tunnel experiment using a model of square-roof building.


Introduction
The objective of the present paper is to investigate the application of a numerical simulation based on the unsteady Bernoulli equation to the estimation of wind loads on PV panels installed parallel to sloped roofs of residential houses. Special attention is paid to a square roof. In the simulation, the time histories of external pressure coefficients on the roof, obtained from a wind tunnel experiment, were employed.
Power generation by photovoltaic (PV) systems has become popular worldwide due to the consideration for environmental conservation. In Japan, it is effective to install PV panels at a tilt angle of 20 • to 40 • with respect to the horizontal plane considering the power generation efficiency. Sloped roofs with a pitch of 20 • to 40 • are often used for residential houses in Japan. Therefore, PV panels installed parallel to the sloped roofs are very effective. However, being vulnerable to dynamic wind actions, they are often subjected to damage or collapse due to strong winds, such as those by typhoons and downbursts [1].
For reducing wind damage to PV systems, it is necessary to estimate the wind loads on PV panels accurately and to evaluate the wind resistant performance of PV systems appropriately. Many researchers have experimentally and/or numerically investigated the wind loads on tilted PV panels installed on the ground [2][3][4][5][6][7][8][9] or on flat roofs of large commercial buildings [10][11][12][13][14][15][16][17][18][19][20][21]. By comparison, only a few investigations have been made of the wind loads on PV panels installed parallel to sloped roofs of residential houses [22][23][24][25][26]. Geurts and Blackmore [22] performed full-scale and wind tunnel experiments to investigate the net uplift loads acting on a stand-off PV system installed parallel to the roof of an existing house with a roof pitch of β = 42 • . In the full-scale measurements, two dummy PV panels (1.6 m long, 0.8 m wide and 18 mm thick) were mounted on the roof. The distance H panel between the panel's lower surface and the roof surface was about 150 mm. The geometric scale of the wind-tunnel model was 1/100. The value of H panel was changed from 0.25 mm to 3 mm (from 25 mm to 300 mm at full scale). The results indicated that the net wind loads on the PV panels were substantially lower than the external loads on the roof surface. Furthermore, it was found that the effect of H panel on the net wind loads was relatively small. This feature may be related to such a high roof pitch as β = 42 • . Aly and Bitsuamlak [23] measured the wind forces on PV panels installed on gable roofs with two distinct pitches-i.e., 14 • and 22.6 • -in a wind tunnel. The geometric scale of the models was 1/15. The value of H panel was fixed to 0.1524 m at full scale. They tested three kinds of panel dimensions (small: 0.9144 m × 1.524 m, medium: 1.524 m × 2.4384 m, and big: 1.524 m × 2.7432 m) and four kinds of panel arrangements; the PV panels were installed almost all over the roof. Focus was on the effects of panel dimension and arrangement on the net wind loads of PV panels. Results indicated that the cladding loads on individual panel were smaller or larger in magnitude than those on the corresponding area of a bare roof (roof without PV panels), depending on the location and dimension of panels, as well as on the roof pitch β. Panels located close to the roof corners and edges were generally subjected to lower net pressures than the external pressures on the bare roof. Aly and Bitsuamlak recommended to avoid mounting PV panels in these zones, because the PV panels would be subjected to high suctions. Stenabough et al. [24] made a similar wind tunnel experiment using a gable roof model with β = 30 • . The geometric scale of the model was 1/20. The PV panels were modeled by flat panels with equivalent full-scale dimensions of 50 cm width, 145.5 cm length and 6 cm thickness. Focus was on the effects of the horizontal gap G between PV panels and H panel on the area-averaged wind loads of PV panels. It was found that larger G values and smaller H panel values yielded lower net wind loads due to pressure equalization, which resulted in the magnitude of net wind loads typically being lower than those for the bare roof surface. Naeiji et al. [25] investigated the net wind loads of PV panels installed on flat, gable and hip roofs using large models in a large wind tunnel, the cross-section of which is 6 m wide and 4 m high. To the authors' best knowledge, only this study dealt with a hip roof. The geometric scale of the models was 1/6. The panel was 2 m long, 1 m wide and 0.15 m thick at full scale. An investigation was conducted on the effects of building height H, panel clearance distance H panel and panel tilt angle β panel on the area-averaged wind loads of PV panels. Note that β panel was not necessarily equal to the roof pitch β. The value of H panel was either 0.3 m or 0.45 m at full scale, which is much larger than that of practical PV systems generally used in Japan. The results indicated that the critical wind direction generating the worst maximum or minimum peak wind force coefficient depended on the roof shape as well as on β panel .
Takamori et al. [26] conducted a wind tunnel experiment on the wind loads of PV panels installed parallel to gable roofs of low-rise buildings. The geometric scale of the models was 1/30. There was no horizontal gap between PV panels, i.e., G = 0 m. The PV panels were 1.06 m long and 0.86 m wide at full scale. The thickness of panel models was 3 mm (90 mm at full scale). The panels were not installed in the edge zones up to 0.3 m at full scale from the roof edges (eaves, gable and ridge). The values of β and H panel were changed from 10 • to 40 • and from 30 mm to 150 mm at full scale, respectively. Based on the results, they proposed positive and negative wind force coefficients for designing PV panels. It was noted that the magnitude of negative wind force coefficients might become larger than that of the proposed values when H panel > 100 mm.
Hip roofs are widely used for residential houses in Japan. Recently, PV panels are often installed parallel to these roofs, as shown in Figure 1. Such PV systems are generally designed based on JIS (Japanese Industrial Standard) C 8955 [27], which is based on the experimental results of Takamori et al. [26]. This standard provides wind force coefficients for designing PV panels installed on building's roofs. However, no specification is provided to PV panels located near the roof edges, up to 0.3 m from the edge. It is not recommended to install PV panels in such zones, because large up-lift forces are generated on the PV panels by flow separation at the roof edges [28,29]. When installing PV panels in such high-suction zones, we need to estimate the wind loads on the PV panels accurately. For estimating wind loads on structures we usually make a wind tunnel experiment. However, it is quite difficult to make wind tunnel models of PV panels with the same geometric scale as that of the building, such as 1/100, for example. The thickness t PV and the distance D cavity between the lower surface of PV panels and the roof surface are both as small as several centimeters. When the geometric scale of the wind tunnel model is 1/100, the values of t PV and D cavity should be less than 1 mm, such as 0.6 mm, for example. Furthermore, we have to measure the pressures on both the upper and lower surfaces of PV panels simultaneously for evaluating the net wind pressures (wind forces) on PV panels, which also makes it difficult to make wind tunnel models of PV panels. Indeed, some deformed models of PV panels are often used in wind tunnel experiments [26].
Wind 2022, 2, FOR PEER REVIEW designed based on JIS (Japanese Industrial Standard) C 8955 [27], which is based on the experimental results of Takamori et al. [26]. This standard provides wind force coefficient for designing PV panels installed on building's roofs. However, no specification is pro vided to PV panels located near the roof edges, up to 0.3 m from the edge. It is not recom mended to install PV panels in such zones, because large up-lift forces are generated on the PV panels by flow separation at the roof edges [28,29]. When installing PV panels in such high-suction zones, we need to estimate the wind loads on the PV panels accurately For estimating wind loads on structures we usually make a wind tunnel experiment However, it is quite difficult to make wind tunnel models of PV panels with the same geometric scale as that of the building, such as 1/100, for example. The thickness tPV and the distance Dcavity between the lower surface of PV panels and the roof surface are both as small as several centimeters. When the geometric scale of the wind tunnel model i 1/100, the values of tPV and Dcavity should be less than 1 mm, such as 0.6 mm, for example Furthermore, we have to measure the pressures on both the upper and lower surfaces o PV panels simultaneously for evaluating the net wind pressures (wind forces) on PV pan els, which also makes it difficult to make wind tunnel models of PV panels. Indeed, som deformed models of PV panels are often used in wind tunnel experiments [26]. In the present paper, we apply a numerical simulation based on the unsteady Ber noulli equation to the estimation of wind pressures on the lower surface of PV panel installed parallel to sloped roofs of residential houses. The pressure in the space between the lower surface of PV panels and the roof surface is called "layer pressure" in the presen paper. The layer pressure acts on the lower surface of PV panel as well as on the roo surface. In the numerical simulation we use the time histories of external pressure coeffi cients simultaneously measured at many points on the roof in a turbulent boundary laye (wind tunnel experiment). Among many kinds of sloped roofs, focus is on a square roo with no projection of eaves, because it is the simplest configuration of sloped roof. Two types of panel arrangements are tested. That is, PV panels are arranged with no gap be tween them (G = 0 m) in one case and with a small gap between them along the short side in the other case. The simulation method for evaluating the layer pressures is validated by a wind tunnel experiment, in which some modifications are made on the model of PV panels, as was done in previous studies [26]. Note that the main purpose of this paper i not to provide the wind force coefficients for designing PV systems but to investigate the application of the numerical simulation to the estimation of layer pressures. Similar sim ulations have been applied to the wind-load estimations of air-permeable double-laye roof systems, loose-laid roof-insulation systems, roofing tiles, permeable unit flooring decks and others [30][31][32][33][34][35]. To the authors' best knowledge, no previous studies have ap plied such a numerical simulation to the wind load estimation of PV panels installed par allel to sloped roofs of low-rise buildings.
The present paper consists of four sections. Section 2 describes the wind tunnel ex periment for measuring the external pressures on the roof and PV panels. Two kinds o building models are used; one has no PV panels and the other has PV panels installed In the present paper, we apply a numerical simulation based on the unsteady Bernoulli equation to the estimation of wind pressures on the lower surface of PV panels installed parallel to sloped roofs of residential houses. The pressure in the space between the lower surface of PV panels and the roof surface is called "layer pressure" in the present paper. The layer pressure acts on the lower surface of PV panel as well as on the roof surface. In the numerical simulation we use the time histories of external pressure coefficients simultaneously measured at many points on the roof in a turbulent boundary layer (wind tunnel experiment). Among many kinds of sloped roofs, focus is on a square roof with no projection of eaves, because it is the simplest configuration of sloped roof. Two types of panel arrangements are tested. That is, PV panels are arranged with no gap between them (G = 0 m) in one case and with a small gap between them along the short sides in the other case. The simulation method for evaluating the layer pressures is validated by a wind tunnel experiment, in which some modifications are made on the model of PV panels, as was done in previous studies [26]. Note that the main purpose of this paper is not to provide the wind force coefficients for designing PV systems but to investigate the application of the numerical simulation to the estimation of layer pressures. Similar simulations have been applied to the wind-load estimations of air-permeable double-layer roof systems, loose-laid roof-insulation systems, roofing tiles, permeable unit flooring decks and others [30][31][32][33][34][35]. To the authors' best knowledge, no previous studies have applied such a numerical simulation to the wind load estimation of PV panels installed parallel to sloped roofs of low-rise buildings.
The present paper consists of four sections. Section 2 describes the wind tunnel experiment for measuring the external pressures on the roof and PV panels. Two kinds of building models are used; one has no PV panels and the other has PV panels installed parallel to the roof. The characteristics of pressure distributions on the roof and PV panels are described. The wind pressure distribution on the PV panels is compared with that on the bare roof in order to investigate the effect of PV panels on the flow field around the roof. It is found that the pressure distribution on the PV panels is similar to that on the bare roof. Therefore, the time histories of wind pressure coefficients on the bare roof are used for the numerical simulation of the layer pressures in the following section. This is also because it is difficult to make models of PV panels with the same geometric scale as that for the building, as mentioned above. Brief discussion is made of the wind forces on the PV panels as well as on the roof. Then, Section 3 investigates the application of a numerical simulation based on the unsteady Bernoulli equation to the evaluation of layer pressures. First, the simulation method is briefly explained. Then, the simulation results for the wind loads on PV panels are compared with the experimental ones in order to validate the simulation method. The effect of PV panels on the wind pressures on the roof is also discussed. Finally, Section 4 summarizes the conclusions obtained from the present study

Wind Tunnel Model
The present study focuses on the PV panels installed parallel to the roof of a squareroof building, as shown in Figure 2. The width, depth and eaves height of the building are 11 m, 11 m and 8.2 m, respectively. The roof pitch is 25 • , which is often used for residential houses in Japan. Consequently, the roof top is 10.7 m high. The roof has no overhang. Three kinds of wind tunnel models, as shown in Figure 3, are used in the experiment. The models are made with a geometric scale of λ L = 1/75 using acrylic plates of 2 mm thickness. These models are named A, B0 and B1, respectively. Model A is not equipped with PV panels on the roof. Model B0 and B1 are equipped with PV panels almost all over a roof surface; the difference between Models B0 and B1 will be described later. The PV panels are also modeled by 2 mm thick acrylic plate. Therefore, the model is somewhat thicker than the practical ones considering that the geometric scale of the model is λ L = 1/75. However, this is a limitation in making models of PV panels. Similar modified models of PV panels were used in previous studies [26].
parallel to the roof. The characteristics of pressure distributions on the roof and PV panel are described. The wind pressure distribution on the PV panels is compared with that o the bare roof in order to investigate the effect of PV panels on the flow field around th roof. It is found that the pressure distribution on the PV panels is similar to that on th bare roof. Therefore, the time histories of wind pressure coefficients on the bare roof ar used for the numerical simulation of the layer pressures in the following section. This i also because it is difficult to make models of PV panels with the same geometric scale a that for the building, as mentioned above. Brief discussion is made of the wind forces o the PV panels as well as on the roof. Then, Section 3 investigates the application of a nu merical simulation based on the unsteady Bernoulli equation to the evaluation of laye pressures. First, the simulation method is briefly explained. Then, the simulation result for the wind loads on PV panels are compared with the experimental ones in order t validate the simulation method. The effect of PV panels on the wind pressures on the roo is also discussed. Finally, Section 4 summarizes the conclusions obtained from the presen study

Wind Tunnel Model
The present study focuses on the PV panels installed parallel to the roof of a square roof building, as shown in Figure 2. The width, depth and eaves height of the building ar 11 m, 11 m and 8.2 m, respectively. The roof pitch is 25°, which is often used for residentia houses in Japan. Consequently, the roof top is 10.7 m high. The roof has no overhang Three kinds of wind tunnel models, as shown in Figure 3, are used in the experiment. Th models are made with a geometric scale of λL = 1/75 using acrylic plates of 2 mm thickness These models are named A, B0 and B1, respectively. Model A is not equipped with PV panels on the roof. Model B0 and B1 are equipped with PV panels almost all over a roo surface; the difference between Models B0 and B1 will be described later. The PV panel are also modeled by 2 mm thick acrylic plate. Therefore, the model is somewhat thicke than the practical ones considering that the geometric scale of the model is λL = 1/75. How ever, this is a limitation in making models of PV panels. Similar modified models of PV panels were used in previous studies [26].  parallel to the roof. The characteristics of pressure distributions on the roof and PV panels are described. The wind pressure distribution on the PV panels is compared with that on the bare roof in order to investigate the effect of PV panels on the flow field around the roof. It is found that the pressure distribution on the PV panels is similar to that on the bare roof. Therefore, the time histories of wind pressure coefficients on the bare roof are used for the numerical simulation of the layer pressures in the following section. This is also because it is difficult to make models of PV panels with the same geometric scale as that for the building, as mentioned above. Brief discussion is made of the wind forces on the PV panels as well as on the roof. Then, Section 3 investigates the application of a numerical simulation based on the unsteady Bernoulli equation to the evaluation of layer pressures. First, the simulation method is briefly explained. Then, the simulation results for the wind loads on PV panels are compared with the experimental ones in order to validate the simulation method. The effect of PV panels on the wind pressures on the roof is also discussed. Finally, Section 4 summarizes the conclusions obtained from the present study

Wind Tunnel Model
The present study focuses on the PV panels installed parallel to the roof of a squareroof building, as shown in Figure 2. The width, depth and eaves height of the building are 11 m, 11 m and 8.2 m, respectively. The roof pitch is 25°, which is often used for residential houses in Japan. Consequently, the roof top is 10.7 m high. The roof has no overhang. Three kinds of wind tunnel models, as shown in Figure 3, are used in the experiment. The models are made with a geometric scale of λL = 1/75 using acrylic plates of 2 mm thickness. These models are named A, B0 and B1, respectively. Model A is not equipped with PV panels on the roof. Model B0 and B1 are equipped with PV panels almost all over a roof surface; the difference between Models B0 and B1 will be described later. The PV panels are also modeled by 2 mm thick acrylic plate. Therefore, the model is somewhat thicker than the practical ones considering that the geometric scale of the model is λL = 1/75. However, this is a limitation in making models of PV panels. Similar modified models of PV panels were used in previous studies [26].   Figure 4 shows the details of the model of PV panels for Model B1. The distance between the roof surface and the lower surface of PV panels is 1 mm. The PV panels are generally installed on horizontal rails mounted on pedestals, as shown in Figure 5. The horizontal rails are modeled by ribs 0.5 mm wide and 0.5 mm high. The pedestals are not reproduced in the model, because they may minutely affect the cavity flow under PV panels. It is assumed that the PV panels are arranged with gaps between them along the short sides. The gap width is tentatively assumed 10 mm at full scale in the present study. The net wind pressure (wind force) on PV panel is provided by the difference between the pressures on the upper and lower surfaces of the panel; the magnitude of net wind pressures may be reduced by pressure equalization due to the gaps [24]. Because the gap of 10 mm width cannot be reproduced as it is with a geometric scale of λ L = 1/75, it is replaced by a circular hole of 1 mm diameter, the area of which is almost equal to that of the gap. The circular hole is placed at the center of gap (see Figure 6). In Model B0, PV panels are arranged with no gap between them; that is, the PV panels are in contact with each other. In practice, Model B0 is obtained by covering the circular holes drilled on the PV panel model of Model B1 with thin adhesive tape.  Figure 4 shows the details of the model of PV panels for Model B1. The distance between the roof surface and the lower surface of PV panels is 1 mm. The PV panels are generally installed on horizontal rails mounted on pedestals, as shown in Figure 5. The horizontal rails are modeled by ribs 0.5 mm wide and 0.5 mm high. The pedestals are not reproduced in the model, because they may minutely affect the cavity flow under PV panels. It is assumed that the PV panels are arranged with gaps between them along the short sides. The gap width is tentatively assumed 10 mm at full scale in the present study. The net wind pressure (wind force) on PV panel is provided by the difference between the pressures on the upper and lower surfaces of the panel; the magnitude of net wind pressures may be reduced by pressure equalization due to the gaps [24]. Because the gap of 10 mm width cannot be reproduced as it is with a geometric scale of λL = 1/75, it is replaced by a circular hole of 1 mm diameter, the area of which is almost equal to that of the gap. The circular hole is placed at the center of gap (see Figure 6). In Model B0, PV panels are arranged with no gap between them; that is, the PV panels are in contact with each other. In practice, Model B0 is obtained by covering the circular holes drilled on the PV panel model of Model B1 with thin adhesive tape.      Figure 4 shows the details of the model of PV panels for Model B1. The distance between the roof surface and the lower surface of PV panels is 1 mm. The PV panels are generally installed on horizontal rails mounted on pedestals, as shown in Figure 5. The horizontal rails are modeled by ribs 0.5 mm wide and 0.5 mm high. The pedestals are not reproduced in the model, because they may minutely affect the cavity flow under PV panels. It is assumed that the PV panels are arranged with gaps between them along the short sides. The gap width is tentatively assumed 10 mm at full scale in the present study. The net wind pressure (wind force) on PV panel is provided by the difference between the pressures on the upper and lower surfaces of the panel; the magnitude of net wind pressures may be reduced by pressure equalization due to the gaps [24]. Because the gap of 10 mm width cannot be reproduced as it is with a geometric scale of λL = 1/75, it is replaced by a circular hole of 1 mm diameter, the area of which is almost equal to that of the gap. The circular hole is placed at the center of gap (see Figure 6). In Model B0, PV panels are arranged with no gap between them; that is, the PV panels are in contact with each other. In practice, Model B0 is obtained by covering the circular holes drilled on the PV panel model of Model B1 with thin adhesive tape.      Figure 4 shows the details of the model of PV panels for Model B1. The distance between the roof surface and the lower surface of PV panels is 1 mm. The PV panels are generally installed on horizontal rails mounted on pedestals, as shown in Figure 5. The horizontal rails are modeled by ribs 0.5 mm wide and 0.5 mm high. The pedestals are not reproduced in the model, because they may minutely affect the cavity flow under PV panels. It is assumed that the PV panels are arranged with gaps between them along the short sides. The gap width is tentatively assumed 10 mm at full scale in the present study. The net wind pressure (wind force) on PV panel is provided by the difference between the pressures on the upper and lower surfaces of the panel; the magnitude of net wind pressures may be reduced by pressure equalization due to the gaps [24]. Because the gap of 10 mm width cannot be reproduced as it is with a geometric scale of λL = 1/75, it is replaced by a circular hole of 1 mm diameter, the area of which is almost equal to that of the gap. The circular hole is placed at the center of gap (see Figure 6). In Model B0, PV panels are arranged with no gap between them; that is, the PV panels are in contact with each other. In practice, Model B0 is obtained by covering the circular holes drilled on the PV panel model of Model B1 with thin adhesive tape.       Figure 7 shows the location of pressure taps (diameter ϕ = 0.6 mm) installed on the roof and on the upper surface of PV panels. Because high suctions are induced near the roof edges, many pressure taps are arranged in the edge zones. Pressure taps are not installed on the lower surface of PV panels. Because the distance between the lower surface of PV panels and the roof surface is small, the pressure on the lower surface of PV panel can be replaced by that on the roof at the same location [23]. The model of PV panels is supported by brass pipes of 1 mm outer diameter connecting to the pressure taps, as shown in Figure 3d. The effect of the brass tubes on the cavity flow under PV panels can be neglected, because the outer diameter of brass pipes is as small as 1 mm and the cavity flow is regarded as a laminar flow [35,36].
Wind 2022, 2, FOR PEER REVIEW 6 Figure 6. Modelling of gaps between PV panels. Figure 7 shows the location of pressure taps (diameter φ = 0.6 mm) installed on the roof and on the upper surface of PV panels. Because high suctions are induced near the roof edges, many pressure taps are arranged in the edge zones. Pressure taps are not installed on the lower surface of PV panels. Because the distance between the lower surface of PV panels and the roof surface is small, the pressure on the lower surface of PV panel can be replaced by that on the roof at the same location [23]. The model of PV panels is supported by brass pipes of 1 mm outer diameter connecting to the pressure taps, as shown in Figure 3d. The effect of the brass tubes on the cavity flow under PV panels can be neglected, because the outer diameter of brass pipes is as small as 1 mm and the cavity flow is regarded as a laminar flow [35,36].

Wind Tunnel Flow
The experiment is conducted in an Eifel type wind tunnel at the Department of Architecture and Building Science, Tohoku University, which has a working section 1.4 m wide, 1.0 m high and 6.5 m long. The wind tunnel flow is a turbulent boundary layer. The profiles of mean wind speed Uz and turbulence intensity Iz at the location of the model's center without model are plotted in Figure 8a. The power-law exponent α for the mean wind speed profile is about 0.27 and the turbulence intensity IH at the mean roof height H (=0.126 m) is about 0.17. Figure 8b shows the normalized power spectral density function,  Figure 8b). According to the AIJ Recommendations for Loads on Buildings [37], the values of α, IH and Lx for Terrain Category III (typical of suburban exposure) are specified as 0.20, 0.26 and 58 m at full scale, respectively. Comparing the above-mentioned values of the wind tunnel flow with these specified values, we can find that α is larger, while IH and Lx are smaller. In particular, the value of Lx is as small as about 1/4 of the target value, considering that the geometric scale is 1/75. In wind tunnel experiments for low-rise buildings, it is generally difficult to satisfy the similarity for all of these parameters, particularly for Lx. Tieleman et al. [38][39][40] investigated the effects of these parameters on the characteristics of wind pressures on the roofs of low-rise buildings. They found that the fluctuating wind pressures were minutely affected by α. By comparison, the effect of IH was significant. Furthermore, the fluctuating pressures were less sensitive to Lx. Regarding the effect of Lx on the roof 14

Wind Tunnel Flow
The experiment is conducted in an Eifel type wind tunnel at the Department of Architecture and Building Science, Tohoku University, which has a working section 1.4 m wide, 1.0 m high and 6.5 m long. The wind tunnel flow is a turbulent boundary layer. The profiles of mean wind speed U z and turbulence intensity I z at the location of the model's center without model are plotted in Figure 8a. The power-law exponent α for the mean wind speed profile is about 0.27 and the turbulence intensity I H at the mean roof height H (=0.126 m) is about 0.17. Figure 8b shows the normalized power spectral density function, f S u ( f )/σ u 2 , of fluctuating wind speed at a height of z = 100 mm, where S u (f ) = power spectral density function, f = frequency, σ u = standard deviation of fluctuating wind speed, and L x = integral length scale of turbulence. It is found that the general shape of S u (f ) agrees well with that of the Karman-type spectrum with L x = 0.2 m (solid line in Figure 8b). According to the AIJ Recommendations for Loads on Buildings [37], the values of α, I H and L x for Terrain Category III (typical of suburban exposure) are specified as 0.20, 0.26 and 58 m at full scale, respectively. Comparing the above-mentioned values of the wind tunnel flow with these specified values, we can find that α is larger, while I H and L x are smaller. In particular, the value of L x is as small as about 1/4 of the target value, considering that the geometric scale is 1/75. In wind tunnel experiments for low-rise buildings, it is generally difficult to satisfy the similarity for all of these parameters, particularly for L x . Tieleman et al. [38][39][40] investigated the effects of these parameters on the characteristics of wind pressures on the roofs of low-rise buildings. They found that the fluctuating wind pressures were minutely affected by α. By comparison, the effect of I H was significant. Furthermore, the fluctuating pressures were less sensitive to L x . Regarding the effect of L x on the roof pressures they mentioned that the wind tunnel experiment could reproduce the fluctuating wind pressures almost satisfactorily, provided that the L x value of the wind tunnel flow was larger than 0.2 times the target value and the maximum length of the wind tunnel model. The wind tunnel flow used in this study satisfies this criterion. That is, the value of L x (≈0.2 m) of the wind tunnel flow is larger than 0.2 times the target value (=58 m/75 = 0.77 m) and the largest size of the building (=10.7/75 = 0.14 m). Tieleman [39] also found that a small-scale turbulence parameter S defined by Equation (1) played an important role on the flow simulation in wind tunnel experiments for low-rise buildings: where U = mean wind speed, and L B = characteristic length of the building. The wind pressure fluctuations were sensitive to S up to about 300, but they were minutely affected by S when S > 300. Therefore, the value of S should be larger than about 300. The value of S of the wind tunnel flow used in the present study was about 470, which was larger than 300. Note that U H and H are used for U and L B in Equation (1). Although the I H value of the wind tunnel flow is somewhat smaller than the target value, such a disagreement can be acceptable, because the main purpose of the present study is not to estimate the design wind loads on PV panels precisely but to investigate the application of the numerical simulation to the wind load estimation of PV panels.
Wind 2022, 2, FOR PEER REVIEW 7 pressures they mentioned that the wind tunnel experiment could reproduce the fluctuating wind pressures almost satisfactorily, provided that the Lx value of the wind tunnel flow was larger than 0.2 times the target value and the maximum length of the wind tunnel model. The wind tunnel flow used in this study satisfies this criterion. That is, the value of Lx (≈0.2 m) of the wind tunnel flow is larger than 0.2 times the target value (=58 m/75 = 0.77 m) and the largest size of the building (=10.7/75 = 0.14 m). Tieleman [39] also found that a small-scale turbulence parameter S defined by Equation (1) played an important role on the flow simulation in wind tunnel experiments for low-rise buildings: where U = mean wind speed, and LB = characteristic length of the building. The wind pressure fluctuations were sensitive to S up to about 300, but they were minutely affected by S when S > 300. Therefore, the value of S should be larger than about 300. The value of S of the wind tunnel flow used in the present study was about 470, which was larger than 300. Note that UH and H are used for U and LB in Equation (1). Although the IH value of the wind tunnel flow is somewhat smaller than the target value, such a disagreement can be acceptable, because the main purpose of the present study is not to estimate the design wind loads on PV panels precisely but to investigate the application of the numerical simulation to the wind load estimation of PV panels.

Experimental Procedure
The design wind speed UH at the mean roof height H (=9.45 m) is calculated based on the AIJ Recommendations for Loads on Buildings [37]. The 'Basic wind speed' U0 is set to 35 m/s as a representative value for the Main Island of Japan. The terrain category is assumed to be III. In practice, UH is calculated as 27.8 m/s. The value of UH of the wind tunnel flow is set to 8 m/s. The Reynolds number Re, defined in terms of UH and H, is approximately 6.7 × 10 4 . The blockage ratio Br of the model with respect to the wind tunnel's cross-section is about 1.9% at a maximum. The values of Re and Br satisfy the experimental criteria recommended in Wind Tunnel Testing for Buildings and other Structures [41]; that is, Re > 1.1 × 10 4 and Br < 5%. The velocity scale λV is approximately 1/3.5. Considering that the geometric scale is λL = 1/75, the time scale λT (=λL/λV) is calculated as approximately 1/22.
The wind direction θ, defined as shown in Figure 9, is varied from 0° to 180° at an increment of 15°, considering the symmetry of the model. The pressure taps are connected to pressure transduces (Wind Engineering Institute, MAPS-02, Midview City, Singapore) The wind pressure coefficient C p is defined by where P = pressure on the structure; P s = static pressure in the wind tunnel; and q H = mean velocity pressure (= 1 2 ρU H 2 , with ρ being the air density) of approach flow at the mean roof height H. The wind force (net wind pressure) acting on PV panel is provided by the difference between the pressures, P t and P b , on the upper and lower surfaces. The lower surface pressure P b is equal to the layer pressure P l under the PV panel. Therefore, the wind force coefficient, or the net wind pressure coefficient, C f , is provided by the following equation: where C pt represents the wind pressure coefficient on the upper surface of PV panel; and C pl is called "layer pressure coefficient" in this paper. When evaluating the statistical values of wind pressure and force coefficients, we apply ensemble averaging to the results of the 10 measurements conducted under the same condition. The distortion of measured fluctuating pressures due to tubing is compensated by using the frequency response function of the measuring system in the frequency domain. No moving average is applied to the time histories of wind pressure and force coefficients.
where P = pressure on the structure; Ps = static pressure in the wind tunnel; and qH = mean velocity pressure (= where Cpt represents the wind pressure coefficient on the upper surface of PV panel; and Cpl is called "layer pressure coefficient" in this paper. When evaluating the statistical values of wind pressure and force coefficients, we apply ensemble averaging to the results of the 10 measurements conducted under the same condition. The distortion of measured fluctuating pressures due to tubing is compensated by using the frequency response function of the measuring system in the frequency domain. No moving average is applied to the time histories of wind pressure and force coefficients.

Wind Pressures on the Roof without PV Panels
First, the characteristics of wind pressures on the square roof are investigated based on the results for Model A. Figures 10 and 11 respectively show the distributions of the maximum and minimum peak pressure coefficients, and , on the roof when θ = 0° and 45°. Note that and refer to the maximum and minimum peak values of wind pressure coefficients Cpe during a period of 10 min at full scale, which can be obtained from the time history of Cpe. Large positive pressures occur near the windward eaves at both wind directions. However, the magnitude is not so large. On the other hand, large suctions occur near the windward eaves when θ = 0° and near the ridge when θ = 45°. These suctions are induced by flow separation at the windward eaves or ridge. Figure  12a

Wind Pressures on the Roof without PV Panels
First, the characteristics of wind pressures on the square roof are investigated based on the results for Model A. Figures 10 and 11 respectively show the distributions of the maximum and minimum peak pressure coefficients,Ĉ pe andČ pe , on the roof when θ = 0 • and 45 • . Note thatĈ pe andČ pe refer to the maximum and minimum peak values of wind pressure coefficients C pe during a period of 10 min at full scale, which can be obtained from the time history of C pe . Large positive pressures occur near the windward eaves at both wind directions. However, the magnitude is not so large. On the other hand, large suctions occur near the windward eaves when θ = 0 • and near the ridge when θ = 45 • . These suctions are induced by flow separation at the windward eaves or ridge. Figure 12a,b, respectively show the distributions of the most critical maximum and minimum peak pressure coefficients,Ĉ pe,cr andČ pe,cr irrespective of wind direction, i.e., the maximum value ofĈ pe and the minimum value ofČ pe among all wind directions. The maximum value ofĈ pe,cr among all pressure taps is 0.66 when θ = 45 • and the minimum value of C pe,cr among all pressure taps is −4.1 when θ = 0 • . These values occur at pressure taps marked by white circles in Figures 10b and 11a. The present results for the wind pressure distributions on the square roof are consistent with those of previous studies investigating the wind pressure distributions on hip roofs [42][43][44][45][46][47][48].

of
, among all pressure taps is −4.1 when θ = 0°. These values occur at pressure taps marked by white circles in Figures 10b and 11a. The present results for the wind pressure distributions on the square roof are consistent with those of previous studies investigating the wind pressure distributions on hip roofs [42][43][44][45][46][47][48].  Figure 13 shows the distributions of the minimum peak pressure coefficients, and , on the roof of Model A and on the PV panels (upper surface) of Model B0 when θ = 0°, 45° and 90°. It is found that the distributions of and are similar to each other. Figure 14 and Figure 15, respectively, show the maximum and minimum peak pressure coefficients among all pressure taps on the roof of Model A and on the PV panels (upper surface) of Model B0 at each wind direction θ. Regarding the maximum peak pressure of , among all pressure taps is −4.1 when θ = 0°. These values occur at pressure taps marked by white circles in Figures 10b and 11a. The present results for the wind pressure distributions on the square roof are consistent with those of previous studies investigating the wind pressure distributions on hip roofs [42][43][44][45][46][47][48].  Figure 13 shows the distributions of the minimum peak pressure coefficients, and , on the roof of Model A and on the PV panels (upper surface) of Model B0 when θ = 0°, 45° and 90°. It is found that the distributions of and are similar to each other. Figure 14 and Figure 15, respectively, show the maximum and minimum peak pressure coefficients among all pressure taps on the roof of Model A and on the PV panels (upper surface) of Model B0 at each wind direction θ. Regarding the maximum peak pressure marked by white circles in Figures 10b and 11a. The present results for the wind pressure distributions on the square roof are consistent with those of previous studies investigating the wind pressure distributions on hip roofs [42][43][44][45][46][47][48].  Figure 13 shows the distributions of the minimum peak pressure coefficients, and , on the roof of Model A and on the PV panels (upper surface) of Model B0 when θ = 0°, 45° and 90°. It is found that the distributions of and are similar to each other. Figure 14 and Figure 15, respectively, show the maximum and minimum peak pressure coefficients among all pressure taps on the roof of Model A and on the PV panels (upper surface) of Model B0 at each wind direction θ. Regarding the maximum peak pressure  Figure 13 shows the distributions of the minimum peak pressure coefficients,Č pe andČ pt , on the roof of Model A and on the PV panels (upper surface) of Model B0 when θ = 0 • , 45 • and 90 • . It is found that the distributions ofČ pe andČ pt are similar to each other. Figures 14 and 15, respectively, show the maximum and minimum peak pressure coefficients among all pressure taps on the roof of Model A and on the PV panels (upper surface) of Model B0 at each wind direction θ. Regarding the maximum peak pressure coefficients, both results agree well with each other. Regarding the minimum peak pressure coefficients, on the other hand, the results for Model B0 are somewhat smaller in magnitude than those for Model A in oblique winds. Similar behavior can be observed in the results for θ = 75-90 • . This feature may be due to a small change in the characteristics of flow separation at the ridge or windward eaves by installing PV panels on the roof. That is, a change in the roof configuration by installing PV panels affects the separated flow, resulting in a change in the pressure distribution on the roof or on the PV panels. However, the effect is not so significant; the maximum difference is about 0.4. Considering that practical PV panels are thinner than the present model at full scale, it is thought that PV panels affect the flow around the roof only slightly. In other words, we can predict the values of layer pressures under PV panels based on the wind pressure coefficients on the roof obtained from a wind tunnel experiment using a building model without PV panels. coefficients, both results agree well with each other. Regarding the minimum peak pressure coefficients, on the other hand, the results for Model B0 are somewhat smaller in magnitude than those for Model A in oblique winds. Similar behavior can be observed in the results for θ = 75°-90°. This feature may be due to a small change in the characteristics of flow separation at the ridge or windward eaves by installing PV panels on the roof. That is, a change in the roof configuration by installing PV panels affects the separated flow, resulting in a change in the pressure distribution on the roof or on the PV panels. However, the effect is not so significant; the maximum difference is about 0.4. Considering that practical PV panels are thinner than the present model at full scale, it is thought that PV panels affect the flow around the roof only slightly. In other words, we can predict the values of layer pressures under PV panels based on the wind pressure coefficients on the roof obtained from a wind tunnel experiment using a building model without PV panels.   coefficients, both results agree well with each other. Regarding the minimum peak pressure coefficients, on the other hand, the results for Model B0 are somewhat smaller in magnitude than those for Model A in oblique winds. Similar behavior can be observed in the results for θ = 75°-90°. This feature may be due to a small change in the characteristics of flow separation at the ridge or windward eaves by installing PV panels on the roof. That is, a change in the roof configuration by installing PV panels affects the separated flow, resulting in a change in the pressure distribution on the roof or on the PV panels. However, the effect is not so significant; the maximum difference is about 0.4. Considering that practical PV panels are thinner than the present model at full scale, it is thought that PV panels affect the flow around the roof only slightly. In other words, we can predict the values of layer pressures under PV panels based on the wind pressure coefficients on the roof obtained from a wind tunnel experiment using a building model without PV panels.   In the present model, 26 rectangular PV panels are installed on the roof; the definition of panel numbers is shown in Figure 16. The PV panels labelled by red figures are subjected to large up-lift forces, as will be described later (see Table 1). Area-averaged wind force coefficient Cf,area is calculated for each PV panel. Because the resolution of pressure taps on the wind tunnel model is relatively coarse, we apply the cubic spline function to the experimental data of Cf in order to obtain the Cf values at the lattice points of a fine grid, from which the value of Cf,area is computed for each panel. Then, the maximum and minimum peak area-averaged wind force coefficients, , and , , are obtained from the time history of Cf,area for each PV panel. Figure 17 shows the minimum value of , among all panels at each wind direction, in which the results for Models B0 and B1 are plotted in order to investigate the effect of gaps between panels on , . Table 1 shows the PV panel that provides the minimum value of , among all panels at each wind direction. Note that the same wind direction provides the minimum values for both Models B0 and B1. It can be seen that the results for Models B0 and B1 show a similar trend. The minimum values occurred on panels located near the eaves or ridge. In particular, panels located near the windward eaves exhibited large-magnitude negative values (upward) when θ = 0°-30°. The largest one was −3.7 for Model B0 and −3.3 for Model B1, which occurred on Panel 4 at θ = 0°. It is found that the values for Model B1 are generally smaller in magnitude than those for Model B0. The largest difference between these two models is about 0.4 at θ = 0°. This feature may be due to the pressure equalization caused by the gaps between PV panels, as Stenabough et al. [24] indicated.   In the present model, 26 rectangular PV panels are installed on the roof; the definition of panel numbers is shown in Figure 16. The PV panels labelled by red figures are subjected to large up-lift forces, as will be described later (see Table 1). Area-averaged wind force coefficient C f,area is calculated for each PV panel. Because the resolution of pressure taps on the wind tunnel model is relatively coarse, we apply the cubic spline function to the experimental data of C f in order to obtain the C f values at the lattice points of a fine grid, from which the value of C f,area is computed for each panel. Then, the maximum and minimum peak area-averaged wind force coefficients,Ĉ f ,area andČ f ,area , are obtained from the time history of C f,area for each PV panel. Figure 17 shows the minimum value ofČ f ,area among all panels at each wind direction, in which the results for Models B0 and B1 are plotted in order to investigate the effect of gaps between panels onČ f ,area . Table 1 shows the PV panel that provides the minimum value ofČ f ,area among all panels at each wind direction. Note that the same wind direction provides the minimum values for both Models B0 and B1. It can be seen that the results for Models B0 and B1 show a similar trend. The minimum values occurred on panels located near the eaves or ridge. In particular, panels located near the windward eaves exhibited large-magnitude negative values (upward) when θ = 0-30 • . The largest one was −3.7 for Model B0 and −3.3 for Model B1, which occurred on Panel 4 at θ = 0 • . It is found that the values for Model B1 are generally smaller in magnitude than those for Model B0. The largest difference between these two models is about 0.4 at θ = 0 • . This feature may be due to the pressure equalization caused by the gaps between PV panels, as Stenabough et al. [24] indicated. In the present model, 26 rectangular PV panels are installed on the roof; the definition of panel numbers is shown in Figure 16. The PV panels labelled by red figures are subjected to large up-lift forces, as will be described later (see Table 1). Area-averaged wind force coefficient Cf,area is calculated for each PV panel. Because the resolution of pressure taps on the wind tunnel model is relatively coarse, we apply the cubic spline function to the experimental data of Cf in order to obtain the Cf values at the lattice points of a fine grid, from which the value of Cf,area is computed for each panel. Then, the maximum and minimum peak area-averaged wind force coefficients, , and , , are obtained from the time history of Cf,area for each PV panel. Figure 17 shows the minimum value of , among all panels at each wind direction, in which the results for Models B0 and B1 are plotted in order to investigate the effect of gaps between panels on , . Table 1 shows the PV panel that provides the minimum value of , among all panels at each wind direction. Note that the same wind direction provides the minimum values for both Models B0 and B1. It can be seen that the results for Models B0 and B1 show a similar trend. The minimum values occurred on panels located near the eaves or ridge. In particular, panels located near the windward eaves exhibited large-magnitude negative values (upward) when θ = 0°-30°. The largest one was −3.7 for Model B0 and −3.3 for Model B1, which occurred on Panel 4 at θ = 0°. It is found that the values for Model B1 are generally smaller in magnitude than those for Model B0. The largest difference between these two models is about 0.4 at θ = 0°. This feature may be due to the pressure equalization caused by the gaps between PV panels, as Stenabough et al. [24] indicated.   In the present model, 26 rectangular PV panels are installed on the roof; the definition of panel numbers is shown in Figure 16. The PV panels labelled by red figures are subjected to large up-lift forces, as will be described later (see Table 1). Area-averaged wind force coefficient Cf,area is calculated for each PV panel. Because the resolution of pressure taps on the wind tunnel model is relatively coarse, we apply the cubic spline function to the experimental data of Cf in order to obtain the Cf values at the lattice points of a fine grid, from which the value of Cf,area is computed for each panel. Then, the maximum and minimum peak area-averaged wind force coefficients, , and , , are obtained from the time history of Cf,area for each PV panel. Figure 17 shows the minimum value of , among all panels at each wind direction, in which the results for Models B0 and B1 are plotted in order to investigate the effect of gaps between panels on , . Table 1 shows the PV panel that provides the minimum value of , among all panels at each wind direction. Note that the same wind direction provides the minimum values for both Models B0 and B1. It can be seen that the results for Models B0 and B1 show a similar trend. The minimum values occurred on panels located near the eaves or ridge. In particular, panels located near the windward eaves exhibited large-magnitude negative values (upward) when θ = 0°-30°. The largest one was −3.7 for Model B0 and −3.3 for Model B1, which occurred on Panel 4 at θ = 0°. It is found that the values for Model B1 are generally smaller in magnitude than those for Model B0. The largest difference between these two models is about 0.4 at θ = 0°. This feature may be due to the pressure equalization caused by the gaps between PV panels, as Stenabough et al. [24] indicated.     0  15  30  45  60  75  90  105  120  135  150  165  180   Panel  4  3  2  3  3  23  9  23  19  19  9  19  18 The most critical maximum and minimum peak area-averaged wind force coefficients irrespective of wind direction,Ĉ f ,area,cr andČ f ,area,cr , in the edge and internal zones of Model B0 are compared with the specified values by JIS C 8955 [27] in Table 2. It should be mentioned that the Standard does not specify the wind force coefficients for PV panels located in the edge zone up to 30 cm from the edge; therefore, a direct comparison cannot be made for the edge zone in a strict sense. Furthermore, the Standard specifies the positive and negative wind force coefficients, C f ,positive and C f ,negative , and the gust effect factor G f . The values of G f × C f ,positive and G f × C f ,negative , respectively, correspond toĈ f ,area,cr anď C f ,area,cr , obtained from the present wind tunnel experiment. The values ofĈ f ,area,cr anď C f ,area,cr predicted from the C pe distributions on the roof of Model A are also presented in Table 2 for a comparative purpose. In the computation we assume that the net wind force on a PV panel is equal to that on the corresponding area of Model A (bare roof). The results indicate that the Standard provides conservative values of wind force coefficient for the panels in the internal zone. As might be expected, the value ofČ f ,area,cr for the panels in the edge zone is larger in magnitude than that in the internal zone. The value ofČ f ,area,cr for the panels in the edge zone predicted from the C pe distribution on the roof of Model A is somewhat larger in magnitude than those obtained from the direct measurements. Considering that the negative (up-lift) wind forces dictate the structural design of PV panels in practice, we can conclude that the design wind loads on PV panels can be estimated from the wind pressure coefficient distribution on the bare roof on the safer side. Table 2. The most critical positive and negative peak wind force coefficients irrespective of wind direction for panels in the edge and internal zones, compared with the specification of JIS C 9855 [27].

Positive Value
Negative Value

Edge Zone Internal Zone Edge Zone
Internal Zone  Figure 18 shows the distributions of the minimum peak pressure coefficients on the roof of Model B0 for some typical wind directions (θ = 0 • , 45 • and 90 • ) and for all wind directions. As might be expected, the magnitude is small in the internal zone, while it is large near the ridges. It can be found that the distributions are quite different from those for Model A (see Figures 12b and 13). Furthermore, the magnitude of the minimum peak pressure coefficients is found to be reduced significantly by installing PV panels. Figure 19 shows the minimum peak pressure coefficient among all pressure taps on the roof of Model A at each wind direction, compared with that of Model B1. We can see a significant reduction of the magnitude of suctions on the roof by installing PV panels on the roof when θ = 0 • , 15 • , 105 • and 120 • . In particular, the effect is very large when θ = 0 • .

Method of Simulation
The method of simulation used in the present paper is the same as that we employ in our previous paper [49], in which we applied this method to the evaluation of wi loads on PV panels horizontally installed on flat roofs. The original idea of this meth came from Amano et al. [30], Okada et al. [31] and Oh and Kopp [32]. In the simulati we estimate the layer pressures using the unsteady Bernoulli equation and the time his ries of the external wind pressure coefficients at the gap location, which were obtain from the above-mentioned wind tunnel experiment. The method is briefly explained he regarding the details, see Uematsu et al. [49].
The space under PV panels is divided into many sub-spaces called "Rooms" [33], shown in Figure 20. Note that the subject is similar to the problem of building inter pressures, when the space under PV panels is regarded as a building [30]. The pressure each Room is a function of time but assumed spatial uniformity. The unsteady Bernou equation is applied to the cavity flows between Rooms, as well as to the gap flows betwe the external space and Rooms. The driving forces of the gap and cavity flows are the pr sure differences between Rooms or between the external space and Room. The governi equations for the gap and cavity flows in the x, y and z directions (regarding the coor nate, see Figure 4) may be given by the following equations (see [49]):

Method of Simulation
The method of simulation used in the present paper is the same as that we emplo in our previous paper [49], in which we applied this method to the evaluation of w loads on PV panels horizontally installed on flat roofs. The original idea of this met came from Amano et al. [30], Okada et al. [31] and Oh and Kopp [32]. In the simulat we estimate the layer pressures using the unsteady Bernoulli equation and the time hi ries of the external wind pressure coefficients at the gap location, which were obtai from the above-mentioned wind tunnel experiment. The method is briefly explained h regarding the details, see Uematsu et al. [49].
The space under PV panels is divided into many sub-spaces called "Rooms" [33 shown in Figure 20. Note that the subject is similar to the problem of building inte pressures, when the space under PV panels is regarded as a building [30]. The pressur each Room is a function of time but assumed spatial uniformity. The unsteady Berno equation is applied to the cavity flows between Rooms, as well as to the gap flows betw the external space and Rooms. The driving forces of the gap and cavity flows are the p sure differences between Rooms or between the external space and Room. The govern equations for the gap and cavity flows in the x, y and z directions (regarding the coo nate, see Figure 4) may be given by the following equations (see [49]):  Figure 19. Minimum peak pressure coefficients among all pressure taps on the roof at each wind direction.

Method of Simulation
The method of simulation used in the present paper is the same as that we employed in our previous paper [49], in which we applied this method to the evaluation of wind loads on PV panels horizontally installed on flat roofs. The original idea of this method came from Amano et al. [30], Okada et al. [31] and Oh and Kopp [32]. In the simulation, we estimate the layer pressures using the unsteady Bernoulli equation and the time histories of the external wind pressure coefficients at the gap location, which were obtained from the above-mentioned wind tunnel experiment. The method is briefly explained here; regarding the details, see Uematsu et al. [49].
The space under PV panels is divided into many sub-spaces called "Rooms" [33], as shown in Figure 20. Note that the subject is similar to the problem of building internal pressures, when the space under PV panels is regarded as a building [30]. The pressure in each Room is a function of time but assumed spatial uniformity. The unsteady Bernoulli equation is applied to the cavity flows between Rooms, as well as to the gap flows between the external space and Rooms. The driving forces of the gap and cavity flows are the pressure differences between Rooms or between the external space and Room. The governing equations for the gap and cavity flows in the x, y and z directions (regarding the coordinate, see Figure 4) may be given by the following equations (see [49]): where U = flow speed (m/s); the subscripts of U represents the room location in a matrix form, for example, i,j U i,j+1 represents the speed of flow from Room (i, j) to Room (i, j + 1); subscript 'e' represents the external space; q H = reference velocity pressure (N/m 2 ); l = distance from the center of a Room to that of the next Room (or to the edge of the Room), parallel to the flow (m); e C = external pressure coefficient at the gap location; C = layer pressure coefficient; the subscripts of e C and C represent the Room location; C L = shape resistance coefficient for the cavity flow in the x or y direction; C Le = pressure loss coefficient of the gap (shape resistance coefficient) in the z direction, depending on the gap configuration; and ∆p = pressure loss due to friction (N/m 2 ). When the boundary of a Room corresponds to the periphery of PV panel, the subscript should be replaced by 'e' in Equations (4) and (5). Note that Equations (4) and (5) are applied to the cavity flows under PV panels or to the gap flows at the periphery of the PV system, while Equation (6) is applied to the flows through the gaps between PV panels. Although the values of C L and C Le depend on the gap configuration in a strict sense, they are assumed to be 1.0 for simplicity. This assumption can be accepted, because the width of gaps and the depth of cavity are relatively large. However, a detailed examination, such as a test with full-scale specimens [50], is required for validation.
ρ ρ  where U = flow speed (m/s); the subscripts of U represents the room location in a ma form, for example, represents the speed of flow from Room (i, j) to Room (i, j subscript 'e' represents the external space; qH = reference velocity pressure (N/m 2 ); l = tance from the center of a Room to that of the next Room (or to the edge of the Roo parallel to the flow (m); eC = external pressure coefficient at the gap location; C = l pressure coefficient; the subscripts of eC and C represent the Room location; CL = sh resistance coefficient for the cavity flow in the x or y direction; CLe = pressure loss co cient of the gap (shape resistance coefficient) in the z direction, depending on the gap figuration; and Δp = pressure loss due to friction (N/m 2 ). When the boundary of a R corresponds to the periphery of PV panel, the subscript should be replaced by 'e' in Eq tions (4) and (5). Note that Equations (4) and (5) are applied to the cavity flows unde panels or to the gap flows at the periphery of the PV system, while Equation (6) is app to the flows through the gaps between PV panels. Although the values of CL and CL pend on the gap configuration in a strict sense, they are assumed to be 1.0 for simpli This assumption can be accepted, because the width of gaps and the depth of cavity relatively large. However, a detailed examination, such as a test with full-scale specim [50], is required for validation. The cavity flow under PV panels is assumed laminar [35,36]. The external pres coefficients eC at the location of gaps in Equations (4)-(6) are obtained from the wind nel experiment mentioned in the previous section. Because the location of gaps does coincide with that of pressure taps on the wind tunnel model, a spatial interpolation u the cubic spline function is applied to the experimental data; the value at the cente each gap is used as a representative value for the gap. Considering the calculation l the above equations are solved by the 4th order Runge-Kutta method. Assuming the w compressibility of the air and an adiabatic condition, the internal pressure P in each R (layer pressure) can be obtained by the following equation (see [49]): The cavity flow under PV panels is assumed laminar [35,36]. The external pressure coefficients e C at the location of gaps in Equations (4)-(6) are obtained from the wind tunnel experiment mentioned in the previous section. Because the location of gaps does not coincide with that of pressure taps on the wind tunnel model, a spatial interpolation using the cubic spline function is applied to the experimental data; the value at the center of each gap is used as a representative value for the gap. Considering the calculation load, the above equations are solved by the 4th order Runge-Kutta method. Assuming the weak compressibility of the air and an adiabatic condition, the internal pressure P in each Room (layer pressure) can be obtained by the following equation (see [49]): where γ = heat capacity ratio of the air; P 0 = atmospheric pressure (N/m 2 ); V 0 = volume of Room (m 3 ); Q m = flow rate (m 3 /s) at gap m; M = total number of gaps; and t = time (s). The layer pressure at the next time step is calculated by the Euler method with a time increment ∆t of 1/8000 s. The sampling interval employed in the pressure measurements (Section 3) was approximately 0.036 s (=1/800 s × 29) at full scale, which is longer than ∆t. Therefore, a temporal interpolation using the cubic spline function is applied to the time histories of wind pressure coefficients obtained from the wind tunnel experiment. The wind force coefficient C f of PV panel is provided by Equation (3) mentioned above. Here, the values of C pt are obtained from the above-mentioned wind tunnel experiment, while those of C pl are provided by the numerical simulation. Figure 21 shows comparisons between experiment and simulation for the distributions of mean layer pressure coefficients C pl for Model B1 when θ = 0 • , 45 • and 90 • . Note that the distribution on a quarter of the roof (shaded area of Figure 9) is depicted in the figure. Figure 22 compares the simulated results with the experimental ones for C pl at each pressure tap when θ = 0 • , 45 • and 90 • ; the tap numbers are defined as shown in Figure 23. The agreement is generally good. When θ = 0 • , the simulation results at Taps 2-6 located near the windward eaves are somewhat larger in magnitude than experimental ones; the largest difference is about 0.2. However, the absolute values of C pl is relatively small at these pressure taps. A similar comparison is made for the minimum peak layer pressure coefficientsČ pl in Figure 24. When θ = 0 • , the simulation results are somewhat larger in magnitude than the experimental ones; the largest difference is about 0.6. The difference is relatively large at pressure taps near the windward eaves, which may be due to a difference in the geometry of the roof's windward edge where the flow separation occurs. In the wind tunnel experiment, PV panels are installed on the roof, while the simulation uses the wind pressure data on the bare roof.

Comparison with the Experimental Results
The simulation tends to overestimate the values of C pl and Č pl to some degree. When θ = 45 • , the agreement is generally good except for Tap 20 located near the roof top, where the local pressure may be affected by the ridge configuration significantly. When θ = 90 • , the simulation underestimates the magnitude ofČ pl at several pressure taps located near the eaves, such as Taps 3-5. This may be due to the effect of a difference in the eaves configuration. The result for Tap 20 is similar to that observed when θ = 45 • . However, the simulation captures the largest wind force coefficient at Tap 1. The data plotted in Figure 24 are the results obtained by applying ensemble averaging to the results of 10 times of simulations. Figure 25 shows the standard deviation (σČ pl ) ofČ pl obtained from these ten simulation data. Although the value of σČ pl is generally small, we can see some cases where it is relatively large, for example, σČ pl > 0.1. In these cases the absolute value of C pl is also large. Therefore, the coefficient of variation, defined by σČ pl / Ĉ pl , is smaller than 0.1 in most cases.

Concluding Remarks
The present paper investigates the wind loads on PV panels installed parallel to sloped roofs of residential houses, including the edge zones where large suctions are generated by flow separation. The main objective of the present paper is to investigate the application of a numerical simulation based on the unsteady Bernoulli equation to the

Concluding Remarks
The present paper investigates the wind loads on PV panels installed parallel to sloped roofs of residential houses, including the edge zones where large suctions are generated by flow separation. The main objective of the present paper is to investigate the Considering that the thickness of PV panels and the distance between the lower surface of PV panels and the roof surface are both smaller in practical PV systems than in the wind tunnel model, the effect of roof configuration at the edges, whether PV panels are installed or not, on the layer pressures is expected to be smaller than those observed here. Therefore, it can be concluded that the numerical simulation employed in the present paper is used for estimating the layer pressures appropriately.

Concluding Remarks
The present paper investigates the wind loads on PV panels installed parallel to sloped roofs of residential houses, including the edge zones where large suctions are generated by flow separation. The main objective of the present paper is to investigate the application of a numerical simulation based on the unsteady Bernoulli equation to the estimation of pressures in the space under PV panels (called 'layer pressures', here), in which the time histories of external pressure coefficients on the bare roof (roof with no PV panel) is used. Because it is difficult to make a wind tunnel model of PV panels with the same geometric scale as that of the building, we focus on such a numerical simulation.
First, we made a wind tunnel experiment using three models of a square-roof building, named Models 'A', 'B0' and 'B1'. Model A had no PV panels installed on the roof. PV panels were installed on the roof with no gap between them for Model B0 and with a small gap between them for Model B1. In any case, PV panels located near the roof edges were subjected to large up-lift forces. The wind pressure distribution on the upper surface of PV panels was similar to that on the roof of Model A (bare roof). This means that installing PV panels on the roof does not affect the flow around the roof significantly. This also implies that we can estimate the layer pressures for Models B0 and B1 from the external pressure coefficients on Model A. This finding is very important from a practical viewpoint. Indeed, various arrangements of PV panels are used. Once we have measured the wind pressures on the bare roof, we can apply the same data to the simulation of layer pressures for any arrangement of PV panels. A comparison for the wind force coefficients on PV panels between Models B0 and B1 indicated that the magnitude of wind force coefficients for Model B1 was somewhat smaller in magnitude than that for Model B0. This feature may be due to the effect of pressure equalization caused by gaps between PV panels, which were represented by small circular holes in the present wind tunnel model. The negative wind forces (up-lift forces), which generally dictate the structural design of PV panels, can be estimated on the safer side from the wind pressures on the roof of Model A (bare roof), assuming that the wind pressures on the upper surface of PV panels are the same as those on the bare roof and the wind pressure on the lower surface of PV panel is zero.
Then, we applied a numerical simulation based on the unsteady Bernoulli equation to the estimation of layer pressures. The layer pressure coefficients are estimated from the time histories of external pressures on Model A at the location of gaps. The simulation results for the layer pressure coefficients on Model B1 were compared with the experimental ones. A relatively good agreement between them was observed. The maximum difference in the minimum peak layer pressure coefficient was about 0.5. This result implies that the numerical simulation employed in the present study is useful for estimating the wind loads on PV panels installed parallel to sloped roofs. The present simulation model is based on the unsteady Bernoulli equation and the time histories of external pressure coefficients measured on the bare roof in a wind tunnel. Therefore, it cannot be applied to a case where the assumptions employed in the present model are not satisfied, for example, when the gap between PV panels and the distance between the lower surface of PV panels and the roof are so large that the flow speeds in the gaps and cavity are high. Furthermore, it cannot be applied to a case where the installation of PV panels on the roof affects the flow around the roof significantly, for example, when the height of PV panels is so high that the wind blows into the space between PV panels and roof. In such cases, a CFD simulation using ELS (Large Eddy Simulation) may be useful; however, it requires very fine mesh and high-spec computer for obtaining accurate Wind 2022, 2 147 solutions. It was found that the gap between PV panels might reduce the net wind forces on PV panels significantly due to the effect of pressure equalization. The optimum size from the viewpoint of wind load reduction can be obtained from a series of simulations, in which the gap size is changed over a wide range. This is the subject of our future study; a preliminary study has been carried out by Yambe et al. [51].