3.1. Flow and Pressure Fields
Figure 3 shows vertical profiles of the horizontally averaged streamwise velocity
, TKE
, and vertical Reynolds stress
for six selected cases of Solid, A12C12, A22C22, A32C32, A12C32, and A32C12. Reference profile data of velocity and Reynolds stress are obtained from Coceal et al. [
28]. The vertical Reynolds stress is calculated based on gradient transport theory with the eddy viscosity. The values inside the block are excluded from spatial averaging to compare flows of outdoor spaces. The wind speed is normalized by the reference wind speed defined at
. TKE and Reynolds stress are scaled by the friction velocity defined as
, where
indicates the domain height, and
represents the air density. Profiles from Coceal et al. [
28] in the figures are determined by direct numerical simulation. As can be seen in the figure, these six case profiles of mean wind speed, TKE, and Reynolds stress are quite similar to one another, both within and above the canyon region. This result indicates that the effect of air flow from the outlet of the wall openings to the outdoor area is minimal in the spatially-averaged sense because of the large differences between the mean wind speed and TKE in the canyon region compared to the outlet flow from the wall openings. Though other cases are not shown in the figure, almost identical profiles are obtained for all statistics.
In general, the streamwise wind speed profiles over
show good agreement with those of Coceal et al. [
28]. On the other hand, the profiles inside the canyon show discrepancies: the profiles here increase almost linearly with height within the canopy layer, whereas those of Coceal et al. [
28] show a slight convex curve below
and a concave curve above
. Furthermore, the Reynolds stress profiles of the present study shows an upward convex curve below
; however, that value determined by Coceal et al. [
28] shows a large gradient near the bottom and top of the canopy. Though there is no reference data for TKE in the canopy layer for the square array to compare, it might be expected that the underestimations of TKE can happen in the canopy layer, as explained in Xie and Castro [
29], resulting in an almost linear slope of the velocity inside the canyon due to insufficient eddy diffusivity. Although the Launder Sharma
k-ε model cannot reproduce velocity profiles inside the canopy for streamwise wind speed with high accuracy in this calculation situation, comparable velocity fields both above and within the canopy layer were reproduced by RANS simulation.
Exterior flow fields are shown in
Figure 4 and
Figure 5.
Figure 4 shows sectional views of the flow fields at the spanwise centres of the buildings for three cases that are typical of the flow patterns among all 11 cases: solid, A12C12, and A32C32. In the figure, only the contours and vectors below
are shown for clarity. The solid case in
Figure 4a shows that recirculation flows exist between buildings. These flow patterns are well known as skimming flow [
20]. The formation of the skimming flow can also be conformed in
Figure 5, which shows the horizontal sectional flow fields of the solid case at
. The high velocity flows in the open spaces between
to
are introduced into the canyon region between two buildings. These flows cause high-speed regions along the right and left faces of the windward wall, resulting in the counter-rotating cortex pair reported by Coceal et al. [
28]. The sectional flow field presented in
Figure 4b for A12C12, where the openings are located in the upper row, shows a jet flow from the windward opening into the room. The jet is neither vertical nor perpendicular to the windward wall, but flows diagonally into the room because of the recirculation flow in the outside canyon region. Moreover, a flow is also observed from the leeward opening to the canyon region behind the building. The leeward jet merges with the recirculation flow outside the room. Although the jet flows can be observed at both windward and leeward openings, the jet flows nonetheless have minimal effects on the outside recirculation flow. The sectional flow field presented in
Figure 4c for A32C32, where the openings are located in the lower row, shows that, in this case, ventilation flow appears to be very weak. A weak jet perpendicular to the opening can be seen at the bottom edge of the windward opening. As in A12C12, the presence of the openings has very little effect on the outdoor flow fields.
Section 3.2 presents a more detailed discussion of the flows in the openings.
The wall pressures due to the flow fields near the buildings are important to estimate the ventilation efficiencies of the buildings; therefore,
Figure 6a,b show the pressure coefficient distributions of the Solid case’s windward and leeward walls, respectively. The rectangles on the figures represent the opening locations of the other cases. On the windward wall, the pressure values generally become larger at higher locations. On the upper half of the buildings, larger pressure values can be observed along the edges, whereas on the lower half, pressure is higher near the spanwise centre. This is because fast winds blow against the building’s upper parts, resulting in a large dynamic pressure contribution to the wall pressure. Conversely, in the building’s lower parts, wind flows toward the spanwise centre along with the windward wall due to counter vortices, which are attributed to the small wall pressure because of the large dynamic pressure of the flow itself.
Figure 6b shows the pressure distribution on the leeward wall. In contrast to the windward wall, the relative values of the pressure coefficients on this face are considerably smaller, and some negative values can be observed. In particular, the pressure values at higher positions are significantly smaller because of recirculation flow in the canyon, which generates flows away from the leeward wall.
Although these wall pressure distributions on cube faces in a block array are qualitatively similar to those obtained by LES [
10] and wind-tunnel experiments [
29], we have to state that RANS simulation may not capture even relative distributions of the pressure distribution. According to Xie and Castro [
30], the vertical gradient of laterally-integrated pressure coefficients near the block roof level of cubes become steeper by the standard
-
model than those obtained by LES. Zaki et al. [
29] reported pressure coefficients of cubes in block arrays for three arrangement patterns with five packing densities. By interpolating their data to 25% packing density in order to compare our simulation results, the pressure coefficient
, where
indicates the block-face averaged pressure differences between front and back faces, which is estimated as 0.12. Contrary, the coefficient becomes 0.05 in the present simulation for solid block, meaning that pressure values of the present simulation are underestimated by more than 60%.
According to these flow and pressure distributions, the RANS simulation reproduces qualitative outdoor flow and pressure fields for the purpose of generating skimming flow interacting with surrounding building arrays, although quantitative difference can been significant especially for turbulence statistics and block wall pressure.
These discrepancies indicate that further quantification of the ventilation rates based on wall pressure values and turbulence statistics might cause considerable differences as compared with those determined by LES or DNS. However, the detailed analyses of various combinations of opening positions can advance a qualitative understanding of the relationship between flow fields and ventilation efficiency. In addition, error estimation of TKE can show how the mean flow and turbulence can affect the ventilation rates. Therefore, results of the RANS-type simulation should be examined well in the next section.
3.2. Flows in Building Openings
As can be seen in the recirculation flow patterns presented in
Section 3.1, the presence of openings has very little effect on the flow fields of outdoor spaces. However, flows introduced through openings may differ considerably depending on the opening positions, with some positions promoting more effective indoor ventilation. Therefore, this section discusses flow distribution patterns in the openings of windward and leeward walls in order to categorize flow introduction types through openings.
Figure 7 shows vertical sectional flow fields in windward openings at the spanwise centres of the three typical cases.
Figure 7 only illustrates flow distributions within openings. The horizontal axis
represents the coordinates from the windward wall normalized by wall thickness
D. Though both of
xz and
xy cross-sectional views were confirmed, the
xz cross-sections are suitable for the following flow categorization. In A12C12, where the openings are in the upper row, wind blows diagonally downward through the openings in most of the regions with large wind speeds, which is why the flows are introduced into the room diagonally, as shown in
Figure 4b. In addition, a triangular reverse flow region can be observed near the upper corner of the interior side. In contrast, in A32C32, where the openings are in the lower row, the wind direction is almost downward as it approaches the exterior side after passing through the opening, and then the flow returns to the bottom edge of the opening. As a result, only air flows from the lower parts are introduced into the room with a wind direction that is normal to the opening. In most of the openings, a reverse flow can be observed though the opening of the windward wall. Similar reverse flows are seen in the case of A31C31, where openings are also located in the lower row, as shown in
Figure 7c. In this case, a vortex appears at
mm; air flows are introduced into the room below the vortex, whereas flows are ejected from the room above the vortex.
In the same manner, vertical sectional flow fields in the windward openings of all cases were visualized and classified into the following two types: the one-way type is characterized by air flows that are unidirectional, such as the downward diagonal of A12C12, and the two-way type is characterized by air flows that are introduced into the room at the bottom edges of the openings and ejected out from the room in the upper regions of openings, such as A32C32 and A31C31. Following this classification, the three cases of A32C32, A31C31, and A31C33 were two-way, and the other seven cases were one-way (schematic figures are summarized in Figure 9).
Flow patterns in the leeward openings were also categorized as either the one-way or two-way type.
Figure 8 shows the vertical sectional flow fields in the leeward openings for two typical cases. Horizontal axis
represents the coordinates from the leeward wall normalized by the wall thickness
D. In A12C12, flows with almost homogeneous wind speeds blow perpendicular to the opening. In contrast, the wind direction of A31C31 is almost parallel to the openings and upward. The upward flow in the opening volume is generated by flow introduced at the bottom edge of the opening from both sides of the room and from the leeward canyon. Accordingly, air flows are ejected into both the room and leeward canyon from the opening volume at the top edge. Eventually, all cases classified as either the one-way type or two-way type at the windward opening were also classified as the same type at the leeward opening.
3.3. Various Ventilation Rate Definitions
In order to quantify the ventilation efficiency differences due to opening locations, four ventilation rates are calculated in this section: net ventilation rate , gross ventilation rate , estimated cumulative and average instantaneous ventilation rate , and the conventional ventilation rate . These ventilation rates are defined in detail as follows:
The conventional ventilation rate
was calculated from wall pressure differences
of the Solid case as expressed in the following equations:
Here,
is the air density;
is the flow rate coefficient (or discharge coefficient) (=0.6) (White, 2010);
(=0.02
2 m
2) as the front and rear opening areas, respectively; and
is the effective opening area. Since the conventional method calculates the introduction flow speed by pressure differences between openings between based on the Bernoulli’s theorem, flow in indoor and outdoor spaces are assumed to diffuse perfectly, indicating openings should be small enough. In addition, correct flow rates need to be determined based on flow introduction angles. According to previous research, conventional methods estimate ventilation rates reasonably when wall porosity is less than 10% [
2]. Advantages of the conventional methods are easiness of determining the ventilation rates as first approximation only based on wall pressure distribution, though correct flow rates are needed. In the present estimation, pressure difference
is determined by the pressure distribution of solid conditions.
As shown in
Section 3.2, there are reverse flows in the openings in the several cases classified as two-way ventilation and, therefore, these reverse flows must be considered to accurately quantify the ventilation efficiency. The net ventilation rate
regards the reverse flow as negative ventilation. In other words,
is determined by the streamwise mean wind speed
times the opening area
as expressed as follows:
Here, and represent the streamwise mean wind speed and opening area at simulation node located in the opening. The conventional method in Equation (1) is thought to estimate this ventilation rate by approximating the flow introduction speed based on pressure differences. The net ventilation rates are thought to be more accurate than the conventional method; however, the wind speed averages in openings are needed for determination. The wind speed averages in openings were determined by simulated results for calculating the net ventilation rates for each case.
On the other hand, the gross ventilation rate
regards the reverse flow as positive ventilation:
Therefore, the gross ventilation rate can be defined individually for the windward and leeward sides. This treatment of reverse flow is the simplest and most often used method in coupled simulations of indoor and outdoor areas [
9]. Contrarily, spatial distributions of wind speeds in openings are necessary for determining the gross ventilation rate. Therefore, it might be difficult to determine the rate except for flow fields obtained by CFD. In the present estimation, spatial distributions of wind speeds were determined by CFD to derive the gross ventilation rates.
Although the present simulation is based on steady-state analysis, the effect of turbulent flow exchange on ventilation efficiency is another factor that must be considered. Therefore, the cumulative and average instantaneous ventilation rate
was estimated by assuming a normal distribution of fluctuating wind speeds through openings following Caciolo et al. [
7]. Instantaneous wind speeds
are approximated by 10
6 random numbers following a normal distribution with a standard deviation by eddy diffusivity and the mean velocity gradient (following eddy viscosity theory with eddy viscosity
, where
). The cumulative and average instantaneous ventilation rate
is determined by integrating the absolute values of estimated instantaneous wind speeds as follows:
Under the realistic condition of the urban boundary layer, exterior flow fields are turbulent in most cases; therefore, the instantaneous ventilation rate is the most realistic volume flow rate which can express air exchange between indoor and outdoor spaces. On the other hand, the ventilation rates require temporal and spatial distribution of wind speeds in openings, indicating assumptions are needed to determine the turbulent contribution on the ventilation as we followed from Cionco et al. [
7] when a RANS-type simulation was employed. According to the previous LES of flow and pressure field data for the same block arrangements conducted by Ikegaya et al. [
10], the probability density of the wall pressure coefficient on both windward and leeward walls can be regarded as a normal distribution because even the maximum skewness was only 0.17 near the opening position of A11. Therefore, the assumption of a normal distribution of fluctuating wind speed through openings is considered plausible for estimating the turbulent contribution to ventilation. We have to state that, however, the accuracy of standard deviation by present RANS simulation might be insufficient for quantifying instantaneous effect, as can be seen in the Reynolds stress profile in
Figure 3c. According to Xie and Castro [
30], the RANS simulation underestimates the TKE values in the canopy layer as compared with LES for an urban-like array, and the underestimation ratios of RANS to LES can be estimated as approximately 30% based on the vertical profiles of TKE [
30]. Therefore, the cumulative and averaged instantaneous ventilation rates with TKE values increased by 30% from the RANS simulation results are also calculated in order to quantify the expected ranges of instantaneous ventilation.
The ventilation rates of all 10 cases are compared in
Figure 9.
and
were individually determined at each windward and leeward opening because these values are unequal, depending on the flow distribution in the openings. The error bar on
indicates expected values when TKE increased 30%. In general, both windward and leeward openings in the higher row (A11C13, A11C11, and A12C12) produced high ventilation rates. In contrast, cases where the openings were in the bottom row on both windward and leeward walls (A32C32, A31C33, and A31C31) produced low ventilation rates. These results are considerably natural because inlets or outlets in upper positions can be exposed to flow with larger wind speed. Moreover, surprisingly, A31C33 and A31C31 showed reversed ventilation, indicating that the net flow direction was from the leeward opening to the windward opening (
and
are always positive, but are shown to be negative in comparison with their
and
magnitudes). These values are simply due to stronger positive or negative wind wall pressure near the higher opening position on the building’s windward or leeward wall.
In cases of one-way ventilation,
is larger than
, although
and
are almost comparable. Therefore, turbulent ventilation can enhance indoor ventilation even when one-way flow through openings is dominant, though increasing rates due to turbulence are not so significant as compared with the total ventilation rates. Even though the expected increase of TKE is considered, the qualitative tendency that
is slightly larger than
does not change because the mean flow mainly determined the ventilation rates, and the contribution of turbulent ventilation is significant for one-way ventilation. In contrast to the one-way types, the ratios of
to
are large in two-way types because inflow and outflow can occur simultaneously through the same openings. Moreover, the difference in the ratios of
and
for the two-way ventilation is especially large at the leeward opening. This may be due to upward shear flows, which can promote turbulence ventilation. The turbulence effect may have grate impact on two-way ventilations as can be seen the expected increasing rate of
due to 30% appreciation of TKE. Lastly, in the cases of one-way ventilation,
and
showed similar trends for each of the opening locations, whereas the
values were always 40% smaller than the
values. As shown in
Section 3.2, wind flows that are directed diagonally downward through an opening into the room result in decreased flow rate coefficients (Vickery and Karakatsanis [
31]). Therefore,
yields a smaller value than
. This result indicates that conventional methods are applicable when flows introduced through openings are one-way ventilation even for the building in the building array; however, estimating the effective opening area, including the flow rate coefficient, is important for the flow affected by the surrounding building arrays based on the wind direction through the openings. Conversely, such as the two-way ventilation cases, estimating the ventilation flow rates by the conventional method can be difficult because of reverse or turbulent flow near, and within, the openings.