Effect of resilient architecture in an ancient windmill in the Sistan region on natural ventilation enhancement

Scientific Reports volume 12, Article number: 18240 (2022) Cite this article

1414 Accesses

Metrics details

Over centuries different elements have been developed in architectures for ensuring adequate natural ventilation in residential units. This study assesses the different components of an ancient windmill in Sistan, Iran, on the structure's indoor air quality (IAQ) enhancement. Several climatic scenarios have been defined by the wind analysis of Sistan meteorological data and analyzed by CFD. The site measurements confirm the accuracy of the simulation results. In the windmill, two deflectors facing the prevailing wind are the significant elements which, in addition to directing wind toward the entrance, could form vortices near the east and west openings leading to suction ventilation. Alteration of the wind speed and angle from 10 to 15 m/s and 30° to 17° would increase the air change per hour (ACH) by 150% and 110%, respectively. Meanwhile, the ACHs were higher than the ASHRAE desired level (ACH > 0.35).

Nowadays, occupants' well-being in buildings places a strain on energy by increasing fossil fuel consumption. Zero net energy (ZNE) architecture would be a solution to reducing energy use. In ZNE, the natural ventilation1,2,3,4,5,6 in dwellings has long been common in arid7 and windy areas, such as some parts of Iran. Residents of these areas have utilized wind in and around buildings for centuries to improve indoor air quality without consuming energy. An examination of these works shows that they were undoubtedly aware of all aerodynamics principles and applied them in their designs.

The design parameters of openings and additional ventilating tools play an essential role in wind-induced ventilation performance. The effect of type and positions of transom windows8 (TWs) on ACH has been evaluated in high-rise buildings1. A combination of CFD analysis and the artificial neural network has indicated that the ACH would increase to 108.1% on average by proper design of TWs (the dimension and direction of windows in accordance with the velocity and direction of the wind).

Wind towers, also known as wind-catcher, Kolak et al.9, are typically vertical structures built on the roof of buildings with mainly rectangular openings facing the prevailing wind for cooling indoors10,11. Site measurements12,13,14,15 as well as numerical studies9,16,17,18,19,20,21,22 confirm the effectiveness and reliability of these structures, introducing them as new solutions to ventilate modern buildings naturally23,24,25,26.

Wing walls or deflectors27 are vertical28 or horizontal29 solid elements around and near windows or windcatchers30, enhancing natural single-side ventilation31,32 by creating extra pressure gradients among openings. It can induce natural suction ventilation by producing a low-pressure region on outlet windows33. This element is efficient when the wind flows by an oblique angle, and the performance diminishes if the wind direction tends to 90°32. Although the wing wall has been used in vernacular buildings (casement windows pivoting outward), modern architectures utilize its practical application in modern design; for example, the UMNO office corridors can be naturally ventilated by its wing walls34,35. The use of wing walls on the edge of a windcatcher revealed that this element would supply a maximum of 9.6 kW cooling capacity at 4 m/s wind speed, which could help to maintain the indoor temperature below the maximin adaptive temperature30.

The indoor air quality would diminish if the outside air is warm36, warm-humid, and polluted37. So, alternative methods are suggested to prevent direct outdoor–indoor air interaction under intensely hot/cold outdoor conditions. Concerning warm outdoor conditions, the combination of evaporation cooling and the solar chimney was suggested7,38,39,40,41 to decrease the indoor temperature by about 8 °C when the outdoor temperature was higher than 35 °C. In addition, as the applicability of evaporation cooling is not recommended in warm-humid regions, the venturi-shape windcatcher42 was designed which can enhance the ventilation rate by 8 times more than cross-ventilation43. This structure was designed based on the Bernoulli effect44 which explains how the negative pressure occurs when the speed increases on one side of the venturi-shape structure leading to suction-ventilation. Further, in air polluted conditions, the outdoors can affect negatively IAQ7,45,46 and the use of filters on entrances was recommended47. However, few studies have been performed on the couple effects of filters with windcatchers and solar chimneys as an effective way of inducing ventilation7,47.

Environmental pollution in Sistan (a region on the border of Iran and Afghanistan), where dust storms48,49,50 occasionally occur because of dry conditions and 120-day winds51, introduced this region as the most polluted city in 201652. In this region, ancient architectures utilized wind which, in addition to cooling, could control dust entrance to buildings15,19,21. Heidari has identified three elements regarding cross-ventilation in Sistan's buildings19, i.e., a windcatcher (called "Kolak"), openings (called "Daricheh"), and a wall groove28 (called "Surak"). Davtalab53 has identified a Sistani traditional element, an inlet opening called Kharkhona, which can affect human comfort by simultaneously utilizing wind and humidity in arid areas53. The results indicated that Kharkhona would transfer the indoors from "very hot" and "hot" to "warm" conditions by reducing Physical Equivalent Temperature (PET) and elevating humidity by 9 °C and 30%, respectively.

This study applied CFD analyses to assess the wind-driven ventilation mechanism employed in an ancient windmill in the Hozdar region, Sistan, Iran, at a longitude of \(61.27^\circ \) and latitude of \(30.55^\circ \). The buoyancy effect causes air change in the building due to the temperature difference between indoors and outdoors. However, in this study, the buoyancy ventilation was considered negligible since, in arid areas, the temperature gradient between indoors and outdoors is not significant6,43,54. Long-term meteorological data were analyzed to assess the wind speeds and directions considered for four scenarios of the CFD analyses, i.e., two different wind angles by two different wind speeds.

The windmill consists of three internal parts; the workplace (where occupants spent most of their time), the entrance, and the rear. The structure was orientated toward prevailing wind with two vertical deflectors constructed at the specified angles and lengths on the windward. Other ventilation elements included one opening as an inlet, six openings on two sidewalls (four out of six were located in two corridors built on two sidewalls of the workplace, and two remaining openings on sidewalls of the entrance), and two openings on the roof.

This paper aimed to assess how the ancients used elements in their architecture to enhance indoor air quality at a desirable level by natural ventilation. Indeed, the IAQ can be investigated by different criteria; in this research, the ACH was introduced as a well-known criterion based on international standards. The present study investigated how deflectors can cause the formation of vortices behind them and near the sidewalls, leading to natural suction ventilation in a wide range of incident angles.

Sistan is a region located in the eastern part of Iran and the western part of Afghanistan. Investigations about wind energy in Iran55,56,57,58,59,60,61 and the global wind atlas62 indicated that Sistan is a region with the highest wind energy potential in Iranian territory (Fig. 1).

The highest wind speed occurring in the Sistan region62.

The 56-year (from February 1963 to February 2020) wind data at Zabol airport meteorological station (Station ID = 40,829), located at 31° 5′ N longitude and 61° 32′ E latitude was used for wind statistical analyses. The wind speed and direction distribution have been analyzed and shown in the wind rose. The results indicated that the dominant wind directions were from northwest to southeast from \(315^\circ \) to \(360^\circ \). The angle of the resultant vector was \(343^\circ \) (or \(17^\circ \) counterclockwise from the North direction, Ω1), and the dominant wind direction was \(330^\circ \) (or \(30^\circ \) counterclockwise from the North direction, Ω2) (Fig. 2).

The wind rose of Sistan, as well as two dominant wind directions.

In addition, the analyses reveal that although the average wind speed increased to near 10 m/s (or 36 \(\mathrm{km}/\mathrm{h}\)) in July, the annual average of top wind speeds is approximately 15 m/s (or 54 \(\mathrm{km}/\mathrm{h}\)) (Fig. 3).

Average and Average of maximum wind speed (1963–2020).

Significant angles of the windmill are depicted in Fig. 4. Two deflectors were constructed on the structure's frontier wall to confront the prevailing wind. The deflectors had two main functions:

The most important function was to accumulate and direct the wind toward the propeller63,64,65,66,67,68,69 and wind entrance.

The second function was to reduce the turbulent flow (by forming vortices) on the east and west wall for enhancing indoor ventilation.

Dimension of elements in the windmill.

The west deflector was constructed by the angle of \({\alpha }_{1}=145^\circ \) and the length of \({L}_{1}=4.5 \mathrm{m}\), where the angle of the east deflector was \({\alpha }_{2}=119^\circ \) with the length of \({L}_{2}=4 \mathrm{m}\).

According to the deflectors and the angle of structure (\(17^\circ \) in the NE-SW direction), the dimensions of the building block in this windmill were 9.5 m (W)\(\times \) 22 m(length) \(\times \) 9 m(height), (Fig. 5). The CFD domain was allocated according to Autodesk guidance70. The dimensions of the domain were 5 W (97.5 m) \(\times \) 6d (156 m) \(\times \) 3 h (27 m). Regarding CFD analysis, four different boundary conditions were allocated to the domain: the steady-state velocity condition for the inlet (different for each scenario), slip/symmetry condition for lateral and top boundaries, zero pressure condition for the outlet, and non-slip condition for the ground.

Domain dimensions and boundary conditions.

In this windmill, nine openings have been considered for natural ventilation. Four vents have been embedded in two corridors (which were constructed for preventing indoor turbulent flow), i.e., three western corridor vents (annotated by WV1, WV2, and WV3) and one eastern corridor vent (annotated by EV1) (Fig. 6). Further, three vents were constructed in the middle room, one on the western side, WV4, one above the door as a transom windows1, EV2, and one on the roof acting as the chimney, RV1. The vent RV2 was an identical chimney on the roof of the rear room. A vent on the frontier wall (ENTRANCE) served as an air intake.

Elements on (a) eastern and (b) western walls of the windmill.

According to wind rose in Sistan (Fig. 2), the resultant wind vector and dominant wind vector are \(\Omega 1=17^\circ \) and \(\Omega 2=30^\circ \), respectively. Further, the maximum monthly wind speed (in July) and an average maximum wind speed are 10 m/s and 18 m/s (Fig. 3). Thus, this study considered four different scenarios based on a statistical analysis of wind in Sistan: two angles of \(17^\circ \) (the angle of resultant wind vector as well as the windmill) and 30°, and two different wind speeds, 10 \(\mathrm{m}/\mathrm{s}\) and 15 \(\mathrm{m}/\mathrm{s}\) (Fig. 7).

Scenarios for CFD of wind-driven ventilation.

The pseudo-steady-state incompressible Reynolds-average Navier–stokes (RANS) for the two-equation \(k-\varepsilon \) turbulent model30 has been chosen for considering the model accuracy and computational cost71. In other words, the ventilation model was carried out only by considering the wind. The governing equations have been listed in Eqs. (1)–(5)4,70.

Continuity:

Momentum conservation:

The Eddy viscosity (\({v}_{t}\)) is calculated in terms of \(k\) and \(\varepsilon \) by

Turbulent kinetic energy (TKE):

Turbulent dissipation rate (TKE):

where \(\overline{u}{ }_{i}\) and \(\overline{p }\) are the averaged velocity components and pressure, respectively. The five modeling constants \({C}_{1\varepsilon }\), \({C}_{2\varepsilon }\), \({C}_{\mu }\), \({\sigma }_{k}\), and \({\sigma }_{k\varepsilon }\) are 1.44, 1.92, 0.09, 1.0, and 1.3, respectively.

Wind direction and speed across three western and eastern corridor vents were measured by an anemometer in different sections near the walls. The anemometer used in this study was UNI-T UT363 with wind speed accuracy and resolution of (± 5%rdg + 0.5 m/s) and 0.1 m/s, respectively. The measurements were done in May 2020, where the average wind speed and direction were 10 m/s and \(20^\circ \), respectively (Fig. 8). The wind speed at specified locations was measured while the wind condition stayed identical around the windmill during measurement.

Field measurement of wind velocity on profiles perpendicular to vents.

The root means square error, RMSE, (Eq. 6) value was used to evaluate the validation of the CFD results.

where \({U}_{o}^{i}=({V}_{o}^{i}/{V}^{\infty })\) is the dimensionless measured average velocity in position i and \({U}_{s}^{i}=({V}_{s}^{i}/{V}^{\infty })\) represents the dimensionless simulated average velocity in the same position, and n denotes the number of measurements. The reference velocity in the measurement date was about \({V}^{\infty }=10\) m/s. According to the mentioned parameter, the average margin error of 17% was obtained (Table 1).

The unstructured tetrahedral mesh was generated by Autodesk CFD to perform the simulation. The computational domain was composed of ~ 5.6 million cells as the reference mesh (Fig. 9). Because of the complexity of geometry, the computational domain was refined in critical zones such as the vents and corridors.

Schematic view of the CFD mesh.

Grid independence analysis was conducted to ensure that fining the grid resolution would not significantly affect the results. Thus, coarse, medium, and fine meshes were generated for Sc17-10 in the CFD model. The average difference, \({R}_{a}\) (Eq. 7), between the velocity profiles for different mesh sizes was defined to study the mesh independence criteria.

where \({V}^{*}\) is the velocity value at finer mesh rather than \(V\) for coarser mesh. Table 2 presents the total number of cells and the percentage difference between the velocity profiles over five different lines being perpendicular to the vents EV1, WV1, WV2, WV3, and WV4, respectively. As seen in Table 1, the \({R}_{a}\) between the coarse and medium grid is close to 1.5% and between the medium and fine mesh is below 1%. Thus, the mesh with about 5.8 million cells was chosen as the main mesh for the present study.

To monitor the convergence, the automatic convergence assessment71 was selected. At first, the criteria were set to moderate, while the accuracy was increased by changing the criteria to tight. Also, the following solution convergence metrics were calculated during the convergence check:

Instantaneous Convergence Slope: This metric calculated the slopes of the convergence data from one iteration to the next. When the maximum instantaneous slope was below the set level, the solution was stopped. In the study, the set value was defined as \(1\times {10}^{-4}\).

Time Averaged Convergence Slope: This metric evaluated the slope of the convergence data over several iterations. The threshold value was set on \(1\times {10}^{-2}\).

Time Averaged Convergence Concavity: The derivation of maximum time-averaged convergence slope is a measure of whether the curve is flattening (slope is descending) or growing (slope is ascending). The solution will stop when the concavity is less than the set value. This value was set on \(1\times {10}^{-2}\) in tight criteria.

Field Variable Fluctuations: It calculates the standard deviation of the dependent variable where the solution will stop when the deviation is below the set level. This value was set on \(1\times {10}^{-5}\) in this study.

In this study, CFD analyses have revealed vortices' performance behind deflectors and assessed how they affect indoor wind-driven ventilation. The ventilation was based on the wind current in a boundary layer near the vents, which would lead to pressure drop and vacuum in them.

In this structure, deflectors play a crucial role in natural ventilation. They can provide air circulation due to the induced pressure gradient between the interior and exterior of the building. Thus, to understand the necessity of deflectors, the model's total pressure distribution is presented in Fig. 10. As depicted in Fig. 10, the circulation on the leeward of the deflectors provides a negative pressure gradient (the outdoor pressure is higher than indoor vents). The pressure gradient and air velocity through the vents in four scenarios have been shown in Fig. 11. This led to suction forces that ventilated the indoor air through corridor vents as well as WV4 and EV2.

Pressure distribution around and inside the windmill.

Pressure gradient in four scenarios for western and eastern vents.

Flow patterns around the windmill have been depicted in Fig. 12. The ventilation by western and eastern vents was directly related to the vortices formed behind the deflectors (Figs. 12, 13, 14). In general, when the wind direction was \(\omega =\beta =17^\circ \) (perpendicular to frontier wall with the angles of attack being \({\gamma }_{1}=55^\circ \) and \({\gamma }_{2}=29^\circ \)), i.e., Sc17-10 and Sc17-15, the length of vortex behind western deflector was larger than behind the eastern one and vice versa, increasing the flow angle to \(\omega =30^\circ \), i.e., Sc30-10 and Sc30-15 resulted in higher vortex length behind the eastern deflector.

Flow pattern around the building in different scenarios: (a) Sc17-10, (b) Sc17-15, (c) Sc30-10, and (d) Sc30-15.

\({V}_{y{^{\prime}}}\) on the profiles parallel and near the (a) western and (b) eastern wall, behind the deflectors.

\({\mathrm{V}}_{\mathrm{y{^{\prime}}}}\) on the perpendicular profiles to western and eastern vents in scenarios (a) Sc17-10, (b) Sc17-15, (c) Sc30-10 and (d) Sc30-15.

The parallel component of flow to the walls, \({V}_{y{^{\prime}}}\), would be an appropriate estimation for the shape of currents around the windmill. The positive values, \({V}_{y{^{\prime}}}>0\), showed reverse flow, while the negative values illustrated current in the same direction as the dominant wind in the space (ordinary current). Thus, the wind speed (\({V}_{y{^{\prime}}}\)) profiles were considered to assess the vortices' length and width.

The vortices' length has been assessed by two 25-m profiles from deflectors near the eastern and western walls (Fig. 13).

Furthermore, the vortices' widths were measured by six profiles, ten meters away from the vents (Fig. 14). These graphs also demonstrated how the wind angles and speeds affected the vortices behind deflectors.

At \(\omega =17^\circ \), the length of vortices was approximately 16.5 m and 8.5 m on the western and eastern wall, respectively, for both wind speeds, i.e., 10 and 15 m/s. Somehow, at \(\omega =30^\circ ,\) the length of vortices was 7 m and 9 m for the wind speed of 10 and 15 m/s, respectively, on the western wall. This value occurred at more than 25 m behind the eastern deflector (higher than the structure length) at both wind speeds.

In general, the widths of vortices were more expansive on the western wall when the wind angle was \(\omega =17^\circ \), while it was wider behind the eastern deflector at \(\omega =30^\circ \) (Fig. 14). In these profiles, on the vortices, the maximum \({V}_{y{^{\prime}}}\) occurred approximately 50 cm away from the walls and then diminished upon increasing the distances from vents. At \(\omega =17^\circ \), the flow currents were directed to a normal flow pattern at near 9 m on the western and 6.5 m on the eastern wall (Fig. 14). On the other hand, at \(\omega =30^\circ \), they tended to normal flow pattern at 5 m and more than 10 m for western and eastern walls, respectively.

In addition, the CFD results on sections A-A revealed that RV1 and RV2 contributed to natural ventilation properly (Fig. 15).

Flow pattern on the section A-A at different scenarios; (a) Sc17-10, (b) Sc17-15, (c) Sc30-10, and (d) Sc30-15.

The \({V}_{y\mathrm{^{\prime}}}\) near the western and eastern vents, and the pressure gradient in the vents has been depicted in Figs. 16 and 17, respectively. In general, the ventilation from vents was directly correlated with the wind current near the vents.

\({V}_{y{^{\prime}}}\) of flow near the (a) western and (b) eastern vents.

Pressure gradient in (a) western and (b) eastern vents.

The pressure gradient in vents was generally reduced upon increasing the distance from deflectors. The trend differed in WV3 and WV4 at \(\omega =30^\circ \) (Fig. 17a) as the absolute \({V}_{y\mathrm{^{\prime}}}\) was higher in WV4 than in WV3 (Fig. 16a).

When the wind angle was \(\omega =30^\circ ,\) no swirl formed on WV3 and WV4, and the flow currents were normal flow patterns on these vents (Fig. 13). This is due to the formation length of vortices behind deflectors. At this angle, the vortex formed before WV3 (Fig. 12), after which the flow reached zero at a specified point before this vent and then increased steadily on WV3 and WV4. Thus, the absolute wind speed in WV4 was higher than in WV3 (Fig. 16). Further, when the wind angle was \(\omega =17^\circ ,\) the vortices on eastern vents were small enough such that no vortices formed on EV2 (Fig. 13).

The flow rate (FR, \({\mathrm{m}}^{3}/\mathrm{h})\) in vents has been calculated and depicted in Fig. 18. Apart from the flow pattern near the walls (reverse or ordinary), the flow rate was outside in all vents (i.e., vacuum). The results indicated that the maximum flow rate in all scenarios occurred in WV4 followed by EV2. Furthermore, the western vents would contribute to ventilation more than the eastern vents do, by 343% at the angle of \(17^\circ \) and 169% at the angle of \(30^\circ \), respectively.

Contribution of different vents to ventilation of indoor.

The effect of windspeed on ventilation has been reported in Table 3. The results revealed that the flow rate has grown upon increasing wind speed from 10 to 15 m/s.

The effect of wind direction on ventilation has been outlined in Table 4. The results revealed that when the wind direction changed from \(17^\circ \) to \(30^\circ \), the ventilation increased dramatically in eastern vents (149% on average) and decreased in western vents (73% on average). The values dropped continually upon the increasing distance from the west deflector, except in WV4 because of different flow patterns on WV3 and WV4. The ventilation in roof vents did not change with the alteration of the wind direction.

The higher flow rate in WV4 and EV2 was related to the higher surface area of these two vents. Accordingly, the average speed has been calculated, with the results illustrated in Fig. 19.

Average speed (Flow rate per area), (\(\left({\mathrm{m}}^{3}/\mathrm{h}\right)/{\mathrm{cm}}^{2}\)).

By looking at each scenario, when \(\omega =17^\circ \), the highest and lowest average speeds occurred in WV1 and EV2, respectively. At \(\omega =30^\circ \), the maximum and minimum average speeds took place in EV2 and WV3, respectively. Further, the lowest average speed occurred in WV3, followed by WV4 in the scenario of Sc30-10, as expected.

The analyses revealed that although the ventilation was higher on the western wall at \(\omega =17^\circ \) by 126%, the eastern wall's contribution to ventilation was higher at \(\omega =30^\circ \) (Table 5).

Further, corridor vents (C.V.) were more effective in natural ventilation than other vents (O.V.) (Table 6). This effect was highlighted in western vents. The maximum ratio of mean average speed in three western corridor vents (WCV) to WV4 occurred at Sc30-10 because of the lowest average speed value in WV4. Furthermore, as expected, the minimum value of WV3/WV4 occurred in Sc30-15.

As stated earlier, near WV4, the flow pattern was normal in Sc30-15 and the absolute \({V}_{y{^{\prime}}}\) was higher than in Sc17-15 (on which the vortex has been formed), \( \left| {\underbrace {{V_{{yWV4}}^{{'30 - 15}} }}_{{ - 0.72}}} \right| > \left| {\underbrace {{V_{{yWV4}}^{{'17 - 15}} }}_{{0.41}}} \right| \); however, the pressure gradient and average speed had a reverse pattern, \( \underbrace {{\Delta P_{{WV4}}^{{30 - 15}} }}_{{3.77}} < \underbrace {{\Delta P_{{WV4}}^{{17 - 15}} }}_{{7.69}} \) and \( \underbrace {{FRPA_{{WV4}}^{{30 - 15}} }}_{{352}} < \underbrace {{FRPA_{{WV4}}^{{17 - 15}} }}_{{477}} \).

Hence, it could be concluded that the flow in vortices would affect ventilation more than ordinary currents.

Looking at Fig. 16, near the vents \({V}_{y{^{\prime}}}\) had the maximum rate on the west in Sc17-15 rather than in other scenarios, leading to the highest-pressure gradient and average speed. In addition, in Sc30-10 and Sc30-15, the pressure gradient and average speed changed with \({V}_{y{^{\prime}}}\) (Table 7).

This pattern was different in WV1 and WV2 in Sc17-10. In this scenario, although the average speed in WV2 was lower than in WV1 (\( \underbrace {{FRPA_{{WV2}}^{{17 - 10}} }}_{{355}} < \underbrace {{FRPA_{{WV1}}^{{17 - 10}} }}_{{393}} \)), \({V}_{y{^{\prime}}}\) near WV2 was higher than near WV1 (\( \underbrace {{V_{{{y'}} _{{WV2}}}^{{17 - 10}} }}_{{1.15~{\text{m}}/{\text{s}}}} > \underbrace {{V_{{{y'}} _{{WV1}}}^{{17 - 10}} }}_{{1.02~{\text{m}}/{\text{s}}}} \)). This could be due to microscale eddies formed near WV1 and WV2, while the reverse current in microscale eddy near WV1 was stronger than near WV2 (Figs. 20, 21). These microscale eddies formed just when the wind speed was 10 m/s.

Microscale eddies near WV1 and WV2 behind the western deflector.

Eddies near (a) WV1 and (b) WV2 in Sc17-10 as well as (c) WV1 and (d) WV2 in Sc30-10.

At \({V}_{\infty }=10 \)m/s, if the wind angle was \(\omega =17^\circ \), the micro eddy started at 35 cm away from the center of the WV1 (Fig. 21) and continued by 25 cm. In this eddy, the maximum absolute value of \({V}_{y{^{\prime}}}\) was \(\left|{{V}_{y{^{\prime}}}}_{WV1}^{Sc17-10}\right|=1.02 \mathrm{m}/\mathrm{s}\), which occurred at 15 cm away from the center of the vent. In addition, the eddy was started from 25 cm to near 10 cm from the center of WV2 (the length was 15 cm), and the maximum absolute value was \(\left|{{V}_{y{^{\prime}}}}_{WV2}^{Sc17-10}\right|=0.13 \mathrm{m}/\mathrm{s}\), occurring at a distance of 15 cm. When \({V}_{\infty }=10 \mathrm{m}/\mathrm{s}\) and \(\omega =30^\circ \), the eddy formed from 20 to 10 cm from WV1 (the length was 10 cm) where the absolute maximum value was \(\left|{{V}_{y{^{\prime}}}}_{WV1}^{Sc30-10}\right|=0.56 \mathrm{m}/\mathrm{s}\). Finally, in the same scenario and near WV2, the eddy formed 16 cm to 14 cm from the center of the vent (2 cm length) with the maximum value of \(\left|{{V}_{y{^{\prime}}}}_{WV2}^{Sc30-10}\right|=0.13 \mathrm{m}/\mathrm{s}\) 15 cm away from the vent center.

Thus, as the maximum reverse flow occurred at the eddy in Sc17-10 near WV1, \({V}_{y{^{\prime}}}\) has diminished more near this vent in this scenario than in others (Fig. 16).

In addition, according to Fig. 16, the minimum \({V}_{y{^{\prime}}}\) in WV1 occurred in Sc17-10, resulting from the described micro-scale eddy.

The results indicate that the mean value of indoor air velocity in the four considered scenarios was about 0.31 m/s. On the other hand, in naturally ventilated buildings, the occupants’ thermal comfort must be evaluated based on adaptive thermal comfort models. As such, in the present study, we use the adaptive thermal comfort model of ASHRAE standard 55-202072. Thus, by considering the summer outdoor mean temperature of the Sistan region (about 33 °C) and mean air velocity of 0.31 m/s, the adaptive comfort model reveals that without using any mechanical cooling systems for the indoor space of the windmill, the occupants’ thermal perception lies within the comfort range with 80% acceptability limit.

IAQ depends on various parameters such as air flow rate, air change rate, and local mean age of air (LMA) in a 1.5 m horizontal plane (level of standing human in activities). Since the air change per hour (ACH) parameter was extensively used in similar studies, this criterion has been utilized for evaluating the IAQ.

The air change per hour, ACH (/h), was calculated through dividing the integrated vent velocities by the room volume, V (Eq. 8)

where \({v}_{j}\) is velocity vector, \({n}_{j}\) is the ordinary vector to the vents surface, \({A}_{j}\) represents the area of jth cell, n is the total number of cells at the vents, and k shows the number of vents in the building.

The results (Fig. 22) showed that all scenarios were beyond the ASHRAE standard for residential buildings (ACH = 0.35)73,74. This figure indicates that air change increases when the wind blows faster. One of the interesting findings is that the air change diminished by increasing the wind angle.

Air change per hour (ACH) in the windmill in different scenarios.

The application of wind in the ventilation of vernacular buildings was explained in various literature. Indeed, the important issue is how we can enhance the ventilation rate in modern architecture with this knowledge. Projects such as Kidderminster College75, and Windcatcher Zero76 are the cases that utilize natural ventilation in modern architecture7,36. The characteristics of the studied structure, including its orientation, the use of deflectors, and the shape and placement of vents as well as corridors can be a suitable model for harnessing the wind in indoor air exchange. Further, each of the mentioned approaches (such as the effect of negative pressure43 behind deflectors) or a combination of these approaches can be proposed as an effective solution for cross or single-side ventilation strategies. For example, this system could be recommended as a suitable auxiliary solution for the proper ventilation process in large public spaces such as halls, museums77, warehouses, enclosed stadiums78, and factories, which on the one hand, overcrowding increases the indoor-induce pollutants79 in the space, and on the other hand, the cost of ventilation in these places is high.

This paper assessed wind-driven ventilation in an ancient structure, a windmill in Sistan, Iran. CFD analyses were applied to study the natural ventilation and validated by site measurement. The Sistan wind flows in the NW–SE direction by the resultant vector of \(17^\circ \), fluctuating from 0° to 60°. Accordingly, the windmill was built based on the resultant wind vector in the region. Two deflectors were constructed wide enough to collect all winds within the mentioned range to utilize the wind in the higher range. Two western and eastern corridors were built to enhance ventilation.

This study chose four different climatic scenarios according to the Sistan wind pattern, i.e., two significant wind speeds by two prevailing wind directions. The current approach illustrated deflectors' functions in the formation of vortices and the contribution of vents to indoor ventilation. An error of 17% was calculated for validation because of the transient condition in the natural wind pattern and the high uncertainty in site measurements. The results revealed that the functionality of vents on the eastern and western walls would change dramatically with changes in wind speed and direction. In general, the findings of this study are summarized as follows:

Regarding natural ventilation, the deflectors had two main functions:

The most important function was to accumulate and direct the wind toward the wind entrance.

The second function involved the diminution of turbulent flow (by forming vortices) on the eastern and western walls to enhance indoor ventilation.

According to the meteorological data, the wind resultant vector in Sistan is \(17^\circ \). In addition, regarding natural ventilation at \(\omega =17^\circ \), the west wall contributed more than the east wall by 126%. Thus, more vents were designed and constructed on the western wall to enhance ventilation.

By changing the wind direction from \(17^\circ \) to \(30^\circ \), the more considerable vortex length behind eastern deflectors led to higher ventilation by 149% on this wall.

Although the higher wind speed (15 m/s) enhanced the ventilation in roof vents, the wind direction did not affect this value in RV1 and RV2.

The microscale eddy near WV1 led to a reduction in \({V}_{y{^{\prime}}}\) on this vent, especially in Sc17-10. It suggests designing a vent with a rounded edge to reduce the flow loss coefficient near vents80.

When the wind speed increased from 10 m/s to 15 m/s, the ACHs dropped by 150%. Further, by changing the wind angle from \(17^\circ \) to \(30^\circ \), the ACH dropped by 10%. However, the ACHs were higher than the ASHRAE standard (ACH > 0.35).

Some aspects should be specified to utilize wind in residential regions, including the current dynamic around buildings, the angle and speed of wind hitting the building, low and high-pressure areas, and vortices around the buildings. In this way, the complete flow pattern around the buildings could be determined and provided for both engineers and architects.

It is necessary to consider some limitations in this paper, such as the deflectors' geometrical characteristics (i.e., length and angle) on vacuum rate. The human comfort analyses (i.e., thermal, humidity, and wind speed effects) should be performed in further investigations considering the effect of corridors on preventing turbulent indoor flow. The impact of moist and porous media on the entrance (Kharkhona53 in the native language of Sistan) should be analyzed on the ENTRANCE.

The underlying principles uncovered by this research and the design type of this structure can be adopted in many buildings and enclosed spaces with a high rate of indoor-induced pollution. However, it is worth noting that not all strategies in designing this structure apply to every location and it is recommended to use this under similar climatic conditions.

The author would like to thank the Global Wind Atlas for access to the wind map of Sistan. The datasets analyzed during the current study are available in the figshare repository, https://doi.org/10.6084/m9.figshare.19775980.

Liu, T. & Lee, W. L. Evaluating the influence of transom window designs on natural ventilation in high-rise residential buildings in Hong Kong. Sustain. Cities Soc. 62, 102406. https://doi.org/10.1016/j.scs.2020.102406 (2020).

Article Google Scholar

Gratia, E., Bruyère, I. & De Herde, A. How to use natural ventilation to cool narrow office buildings. Build. Environ. 39, 1157–1170. https://doi.org/10.1016/j.buildenv.2004.02.005 (2004).

Article Google Scholar

Rizk, A., El Morsi, M. & Elwan, M. A review on wind-driven cross-ventilation techniques inside single rooms. Int. J. Sci. Eng. Res. 6, 75 (2018).

Google Scholar

Cheung, J. O. P. & Liu, C.-H. CFD simulations of natural ventilation behaviour in high-rise buildings in regular and staggered arrangements at various spacings. Energy Build. 43, 1149–1158. https://doi.org/10.1016/j.enbuild.2010.11.024 (2011).

Article Google Scholar

Kosutova, K., van Hooff, T., Vanderwel, C., Blocken, B. & Hensen, J. Cross-ventilation in a generic isolated building equipped with louvers: Wind-tunnel experiments and CFD simulations. Build. Environ. 154, 263–280. https://doi.org/10.1016/j.buildenv.2019.03.019 (2019).

Article Google Scholar

Khan, N., Su, Y. & Riffat, S. B. A review on wind driven ventilation techniques. Energy Build. 40, 1586–1604. https://doi.org/10.1016/j.enbuild.2008.02.015 (2008).

Article Google Scholar

Ahmed, T., Kumar, P. & Mottet, L. Natural ventilation in warm climates: The challenges of thermal comfort, heatwave resilience and indoor air quality. Renew. Sustain. Energy Rev. 138, 110669. https://doi.org/10.1016/j.rser.2020.110669 (2021).

Article CAS Google Scholar

Aflaki, A., Hirbodi, K., Mahyuddin, N., Yaghoubi, M. & Esfandiari, M. Improving the air change rate in high-rise buildings through a transom ventilation panel: A case study. Build. Environ. 147, 35–49. https://doi.org/10.1016/j.buildenv.2018.10.011 (2019).

Article Google Scholar

Varela-Boydo, C. A. & Moya, S. L. Inlet extensions for wind towers to improve natural ventilation in buildings. Sustain. Cities Soc. 53, 101933. https://doi.org/10.1016/j.scs.2019.101933 (2020).

Article Google Scholar

Iyengar, K. Sustainable architectural Design: An overview. (2015).

Bahadori, M. N. Passive cooling systems in Iranian architecture. Sci. Am. 238, 144–155 (1978).

Article Google Scholar

Yaghoubi, M. A., Sabzevari, A. & Golneshan, A. A. Wind towers: Measurement and performance. Sol. Energy 47, 97–106. https://doi.org/10.1016/0038-092X(91)90040-4 (1991).

Article ADS Google Scholar

Sangdeh, P. K. & Nasrollahi, N. Windcatchers and their applications in contemporary architecture. Energy Built Environ. https://doi.org/10.1016/j.enbenv.2020.10.005 (2020).

Article Google Scholar

Elmualim, A. A. Verification of design calculations of a wind catcher/tower natural ventilation system with performance testing in a real building. Int. J. Vent. 4, 393–404. https://doi.org/10.1080/14733315.2005.11683717 (2006).

Article Google Scholar

Molanaei, S. & Soleimani, S. Insight into the Valuable Elements of Sistan Local Architecture in Relation to Climatic Factors of Sustainable Architecture. Vol. 13 (Bagh-E Nazar Garden of Vision), 2016).

Calautit, J. K., O’Connor, D. & Hughes, B. R. Determining the optimum spacing and arrangement for commercial wind towers for ventilation performance. Build. Environ. 82, 274–287. https://doi.org/10.1016/j.buildenv.2014.08.024 (2014).

Article Google Scholar

Montazeri, H. & Montazeri, F. CFD simulation of cross-ventilation in buildings using rooftop wind-catchers: Impact of outlet openings. Renew. Energy 118, 502–520. https://doi.org/10.1016/j.renene.2017.11.032 (2018).

Article Google Scholar

Liu, S., Mak, C. M. & Niu, J. Numerical evaluation of louver configuration and ventilation strategies for the windcatcher system. Build. Environ. 46, 1600–1616. https://doi.org/10.1016/j.buildenv.2011.01.025 (2011).

Article Google Scholar

Heidari, A., Sahebzadeh, S. & Dalvand, Z. Natural ventilation in vernacular architecture of Sistan, Iran; Classification and CFD study of compound rooms. Sustainability 9, 1048. https://doi.org/10.3390/su9061048 (2017).

Article Google Scholar

Elmualim, A. A. Dynamic modelling of a wind catcher/tower turret for natural ventilation. Build. Serv. Eng. Res. Technol. 27, 165–182. https://doi.org/10.1191/0143624406bse159oa (2006).

Article Google Scholar

Sahebzadeh, S., Heidari, A., Kamelnia, H. & Baghbani, A. Sustainability features of Iran's vernacular architecture: A comparative study between the architecture of hot–arid and hot–arid–windy regions. Sustainability 9, 749 (2017).

Article Google Scholar

Bahadori, M. N. An improved design of wind towers for natural ventilation and passive cooling. Sol. Energy 35, 119–129. https://doi.org/10.1016/0038-092X(85)90002-7 (1985).

Article ADS Google Scholar

Bahadori, M. N., Mazidi, M. & Dehghani, A. R. Experimental investigation of new designs of wind towers. Renew. Energy 33, 2273–2281. https://doi.org/10.1016/j.renene.2007.12.018 (2008).

Article Google Scholar

Alwetaishi, M. & Gadi, M. New and innovative wind catcher designs to improve indoor air quality in buildings. Energy Built Environ. 2, 337–344. https://doi.org/10.1016/j.enbenv.2020.06.009 (2021).

Article Google Scholar

Calautit, J. K., O’Connor, D. & Hughes, B. R. A natural ventilation wind tower with heat pipe heat recovery for cold climates. Renew. Energy 87, 1088–1104. https://doi.org/10.1016/j.renene.2015.08.026 (2016).

Article Google Scholar

Calautit, J. K. et al. Development of a natural ventilation windcatcher with passive heat recovery wheel for mild-cold climates: CFD and experimental analysis. Renew. Energy 160, 465–482. https://doi.org/10.1016/j.renene.2020.05.177 (2020).

Article Google Scholar

Che, W., Ding, J. & Li, L. Airflow deflectors of external windowsto induce ventilation: Towards COVID-19 prevention and control. Sustain Cities Soc 77, 103548. https://doi.org/10.1016/j.scs.2021.103548 (2022).

Article PubMed Google Scholar

Aflaki, A., Mahyuddin, N., Al-Cheikh Mahmoud, Z. & Baharum, M. R. A review on natural ventilation applications through building façade components and ventilation openings in tropical climates. Energy Build. 101, 153–162. https://doi.org/10.1016/j.enbuild.2015.04.033 (2015).

Article Google Scholar

Givoni, B. Man, Climate and Architecture (Applied Science Publishers, 1976).

Google Scholar

Nejat, P., Salim Ferwati, M., Calautit, J., Ghahramani, A. & Sheikhshahrokhdehkordi, M. Passive cooling and natural ventilation by the windcatcher (Badgir): An experimental and simulation study of indoor air quality, thermal comfort and passive cooling power. J. Build. Eng. 41, 102436. https://doi.org/10.1016/j.jobe.2021.102436 (2021).

Article Google Scholar

Mak, C. M., Cheng, C. & Niu, J. L. The application of computational fluid dynamics to the assessment of green features in buildings: Part 1: Wing walls. Archit. Sci. Rev. 48, 121–134. https://doi.org/10.3763/asre.2005.4816 (2005).

Article Google Scholar

Mozaffari Ghadikolaei, F., Ossen, D. R. & Mohamed, M. F. Effects of wing wall at the balcony on the natural ventilation performance in medium-rise residential buildings. J. Build. Eng. 31, 101316. https://doi.org/10.1016/j.jobe.2020.101316 (2020).

Article Google Scholar

Natural ventilation. https://www.briangwilliams.us/solar-energy/natural-ventilation.html (2020).

Wood, A. & Salib, R. CTBUH Technical Guide Natural Ventilation in High-Rise Office Buildings Principal Authors Contributors/Peer Review (Routledge, 2012).

Book Google Scholar

Gomaa, B. In PLEA.

Zafarmandi, S. & Mahdavinejad, M. Technology of modern windcatchers: A review. Int. J. Archit. Eng. Urban Plan. 31, 1–11. https://doi.org/10.22068/ijaup.31.3.549 (2021).

Article Google Scholar

Choi, Y. & Song, D. How to quantify natural ventilation rate of single-sided ventilation with trickle ventilator?. Build. Environ. 181, 107119. https://doi.org/10.1016/j.buildenv.2020.107119 (2020).

Article Google Scholar

Abdallah, A. S. H. et al. Integration of evaporative cooling technique with solar chimney to improve indoor thermal environment in the New Assiut City, Egypt. Int. J. Energy Environ. Eng. 4, 45. https://doi.org/10.1186/2251-6832-4-45 (2013).

Article CAS Google Scholar

Hassan Abdallah, A. S. Occupant comfort and indoor temperature reduction by using passive air conditioning system with solar chimney concept in hot arid climate. Proced. Eng. 205, 1100–1107. https://doi.org/10.1016/j.proeng.2017.10.178 (2017).

Article Google Scholar

Abdallah, A. S. H. et al. Parametric investigation of solar chimney with new cooling tower integrated in a single room for New Assiut city, Egypt climate. Int. J. Energy Environ. Eng. 5, 92. https://doi.org/10.1007/s40095-014-0092-6 (2014).

Article Google Scholar

Maerefat, M. & Haghighi, A. P. Natural cooling of stand-alone houses using solar chimney and evaporative cooling cavity. Renew. Energy 35, 2040–2052. https://doi.org/10.1016/j.renene.2010.02.005 (2010).

Article Google Scholar

Zhang, L. et al. Study on the influence and optimization of the venturi effect on the natural ventilation of buildings in the Xichang area. Energies 14, 5053 (2021).

Article CAS Google Scholar

Haw, L. C., Saadatian, O., Sulaiman, M. Y., Mat, S. & Sopian, K. Empirical study of a wind-induced natural ventilation tower under hot and humid climatic conditions. Energy Build. 52, 28–38. https://doi.org/10.1016/j.enbuild.2012.05.016 (2012).

Article Google Scholar

Lechner, N. Heating, Cooling, Lighting: Sustainable Design Methods for Architects (Wiley, 2014).

Google Scholar

Liu, C. et al. Influence of natural ventilation rate on indoor PM2.5 deposition. Build. Environ. 144, 357–364. https://doi.org/10.1016/j.buildenv.2018.08.039 (2018).

Article Google Scholar

Chen, J., Brager, G. S., Augenbroe, G. & Song, X. Impact of outdoor air quality on the natural ventilation usage of commercial buildings in the US. Appl. Energy 235, 673–684. https://doi.org/10.1016/j.apenergy.2018.11.020 (2019).

Article CAS Google Scholar

Kang, K. et al. Filtration efficiency and ventilation performance of window screen filters. Build. Environ. 178, 106878. https://doi.org/10.1016/j.buildenv.2020.106878 (2020).

Article Google Scholar

Kaskaoutis, D. G. et al. Meteorological aspects associated with dust storms in the Sistan region, southeastern Iran. Clim. Dyn. 45, 407–424. https://doi.org/10.1007/s00382-014-2208-3 (2015).

Article Google Scholar

Miri, A., Ahmadi, H., Ekhtesasi, M. R., Panjehkeh, N. & Ghanbari, A. Environmental and socio-economic impacts of dust storms in Sistan Region, Iran. Int. J. Environ. Stud. 66, 343–355. https://doi.org/10.1080/00207230902720170 (2009).

Article Google Scholar

Rashki, A., Kaskaoutis, D., Francois, P., Kosmopoulos, P. & Legrand, M. Dust-storm dynamics over Sistan region, Iran: Seasonality, transport characteristics and affected areas. Aeolian Res. 16, 25. https://doi.org/10.1016/j.aeolia.2014.10.003 (2015).

Article ADS Google Scholar

Alizadeh-Choobari, O., Zawar-Reza, P. & Sturman, A. The, "wind of 120days" and dust storm activity over the Sistan Basin. Atmos. Res. 143, 328–341. https://doi.org/10.1016/j.atmosres.2014.02.001 (2014).

Article CAS Google Scholar

Nelson, W. L. Zabol Has World's Most Polluted Air (The Free Library, 2016).

Google Scholar

Davtalab, J. & Heidari, A. The effect of kharkhona on outdoor thermal comfort in hot and dry climate: A case study of Sistan Region in Iran. Sustain. Cities Soc. 65, 102607. https://doi.org/10.1016/j.scs.2020.102607 (2021).

Article Google Scholar

Hughes, B. R. & Cheuk-Ming, M. A study of wind and buoyancy driven flows through commercial wind towers. Energy Build. 43, 1784–1791. https://doi.org/10.1016/j.enbuild.2011.03.022 (2011).

Article Google Scholar

Mostafaeipour, A., Jadidi, M., Mohammadi, K. & Sedaghat, A. An analysis of wind energy potential and economic evaluation in Zahedan, Iran. Renew. Sustain. Energy Rev. 30, 641–650. https://doi.org/10.1016/j.rser.2013.11.016 (2014).

Article Google Scholar

Fazelpour, F., Soltani, N., Soltani, S. & Rosen, M. A. Assessment of wind energy potential and economics in the north-western Iranian cities of Tabriz and Ardabil. Renew. Sustain. Energy Rev. 45, 87–99. https://doi.org/10.1016/j.rser.2015.01.045 (2015).

Article Google Scholar

Keyhani, A., Ghasemi-Varnamkhasti, M., Khanali, M. & Abbaszadeh, R. An assessment of wind energy potential as a power generation source in the capital of Iran, Tehran. Energy 35, 188–201. https://doi.org/10.1016/j.energy.2009.09.009 (2010).

Article Google Scholar

Kahkha Moghaddam, P. & Delbari, M. Evaluation of the feasibility of wind energy utilization in Sistan and Baluchestan Province. Phys. Geogr. Res. Q. 49, 441–455. https://doi.org/10.22059/jphgr.2017.218706.1006952 (2017).

Article Google Scholar

Fadai, D. The feasibility of manufacturing wind turbines in Iran. Renew. Sustain. Energy Rev. 11, 536–542. https://doi.org/10.1016/j.rser.2005.01.012 (2007).

Article Google Scholar

Fazelpour, F., Soltani, N. & Rosen, M. A. Wind resource assessment and wind power potential for the city of Ardabil, Iran. Int. J. Energy Environ. Eng. 6, 431–438 (2015).

Article Google Scholar

Fazelpour, F., Markarian, E. & Soltani, N. Wind energy potential and economic assessment of four locations in Sistan and Balouchestan province in Iran. Renew. Energy 109, 646–667. https://doi.org/10.1016/j.renene.2017.03.072 (2017).

Article Google Scholar

Atlas, G. W. https://globalwindatlas.info/en/area/Islamic%20Republic%20of%20Iran?print=true

Agha, A., Chaudhry, H. N. & Wang, F. Diffuser augmented wind turbine (DAWT) technologies: A review. Int. J. Renew. Energy Res. 8, 1369–1385 (2018).

Google Scholar

Altan, B. D., Atılgan, M. & Özdamar, A. An experimental study on improvement of a Savonius rotor performance with curtaining. Exp. Thermal Fluid Sci. 32, 1673–1678. https://doi.org/10.1016/j.expthermflusci.2008.06.006 (2008).

Article Google Scholar

Vaz, J. R. P. & Wood, D. H. Effect of the diffuser efficiency on wind turbine performance. Renew. Energy 126, 969–977. https://doi.org/10.1016/j.renene.2018.04.013 (2018).

Article Google Scholar

Khamlaj, T. A. & Rumpfkeil, M. P. Analysis and optimization of ducted wind turbines. Energy 162, 1234–1252. https://doi.org/10.1016/j.energy.2018.08.106 (2018).

Article Google Scholar

Bagheri-Sadeghi, N., Helenbrook, B. & Visser, K. Ducted wind turbine optimization and sensitivity to rotor position. Wind Energy Sci. 3, 221–229. https://doi.org/10.5194/wes-3-221-2018 (2018).

Article ADS Google Scholar

Şahin, B., İlhan, A. & Bilgili, M. Investigation of Diffuser Augmented Wind Turbine Technologies. Çukurova Üniversitesi Mühendislik-Mimarlik Fakültesi Dergisi (2017).

Kanya, B. & Visser, K. Experimental validation of a ducted wind turbine design strategy. Wind Energy Sci. Discuss. 20, 1–22. https://doi.org/10.5194/wes-2017-54 (2017).

Article Google Scholar

Autodesk. Natural Ventilation, https://knowledge.autodesk.com/support/cfd/learn-explore/caas/CloudHelp/cloudhelp/2014/ENU/SimCFD/files/GUID-F6011744-B69B-47D0-BB2B-50E1BA8DC32D-htm.html (2019).

Launder, B. E. & Spalding, D. B. The numerical computation of turbulent flows. Comput. Methods Appl. Mech. Eng. 3, 269–289. https://doi.org/10.1016/0045-7825(74)90029-2 (1974).

Article ADS MATH Google Scholar

American Society of Heating, R. a. A.-C. E. (2021).

Das, D. et al. Estimating residential air exchange rates in rural Bangladesh using a near field-far field model. Build. Environ. 206, 108325. https://doi.org/10.1016/j.buildenv.2021.108325 (2021).

Article Google Scholar

ASHRAE, A. A. S. (American Society of Heating, Refrigerating, and Air-Conditioning Engineers, Inc, Atlanta, GA, 2019).

Kidderminster College. https://www.monodraught.com/projects/kidderminster-college (2022).

Windcatcher Zero. https://www.monodraught.com/products/natural-ventilation/windcatcher-zero. (2022).

Litti, G., Audenaert, A. & Braet, J. Natural Ventilation as Passive Cooling Strategy Aimed at Summer Overheating Reduction in Heritage Buildings: The Case Study of Vleeshuis Museum in Antwerp (Belgium). (2014).

van Hooff, T. & Blocken, B. CFD evaluation of natural ventilation of indoor environments by the concentration decay method: CO2 gas dispersion from a semi-enclosed stadium. Build. Environ. 61, 1–17. https://doi.org/10.1016/j.buildenv.2012.11.021 (2013).

Article Google Scholar

Ali, M. Y., Hanafiah, M. M., Khan, M. F. & Latif, M. T. Quantitative source apportionment and human toxicity of indoor trace metals at university buildings. Build. Environ. 121, 238–246. https://doi.org/10.1016/j.buildenv.2017.05.032 (2017).

Article Google Scholar

Costa, N. P., Maia, R., Proença, M. F. & Pinho, F. T. Edge effects on the flow characteristics in a 90deg Tee junction. J. Fluids Eng. 128, 1204–1217. https://doi.org/10.1115/1.2354524 (2006).

Article Google Scholar

Download references

Tose Sazan Pahne Kiyansea (TSPK) Company, Incubator Center, University of Zabol, Zabol, Iran

Elham Mohammadi & Mohammadreza Jarkeh

Department of Mechanical Engineering, Faculty of Engineering, University of Birjand, Birjand, Iran

Seyed Alireza Zolfaghari & Vahid Arbabi

Department of Orthopedics, University Medical Center Utrecht, Utrecht, The Netherlands

Vahid Arbabi

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

M.J. and E.M. are responsible for developing methods, drawing charts, analyzing results, and writing the original draft. V.A. and S.A.Z. are responsible for data curation and methodology improvements.

Correspondence to Vahid Arbabi.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Mohammadi, E., Jarkeh, M., Zolfaghari, S.A. et al. Effect of resilient architecture in an ancient windmill in the Sistan region on natural ventilation enhancement. Sci Rep 12, 18240 (2022). https://doi.org/10.1038/s41598-022-23027-w

Download citation

Received: 09 May 2022

Accepted: 24 October 2022

Published: 29 October 2022

DOI: https://doi.org/10.1038/s41598-022-23027-w

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.