Next Article in Journal
Shaking Table Test and Numerical Simulation Study of the Reinforcement Strengthening of a Dam
Previous Article in Journal
Experimental and Numerical Study of the Behavior of Concrete-Filled High-Strength Steel Tube Columns with Large D/t Ratio under Axial Compression
Previous Article in Special Issue
How the Energy Price Variability in Italy Affects the Cost of Building Heating: A Trnsys-Guided Comparison between Air-Source Heat Pumps and Gas Boilers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Determination of Window Dimensions Using Mathematical and Simulation Methods to Improve Window Performance: A Case Study on the Jatinegara Barat Flats, Jakarta

by
Sally Septania Napitupulu
1,*,
Gagoek Hardiman
2 and
Rumiati Rosaline Tobing
1
1
Architecture Postgraduate Program, Faculty of Engineering, Universitas Katolik Parahyangan, Bandung 40141, Indonesia
2
Architecture Postgraduate Program, Faculty of Engineering, Universitas Diponegoro, Semarang 50275, Indonesia
*
Author to whom correspondence should be addressed.
Buildings 2022, 12(11), 1954; https://doi.org/10.3390/buildings12111954
Submission received: 10 October 2022 / Revised: 2 November 2022 / Accepted: 4 November 2022 / Published: 11 November 2022

Abstract

:
The construction of flats is one of the steps toward fulfilling the need for housing in Jakarta and other big cities in Indonesia. This study investigated the thermal problems that focus on air movement in residential units of flats due to window dimensions that cannot accommodate the air velocity that surrounds the buildings because of monotonous window dimensions and the high elevation of the residential units. The position of the interest of this research is on the search for window dimensions that allow comfortable air movement in the residential unit. Based on these problems and interest, the pertinent issue is the design of window dimensions within the facades of the buildings according to the actual air velocity and the elevation of the residential units. The selected object of study was the Jatinegara Barat flats, a block of high-rise flats located in the middle of East Jakarta, Indonesia, which is adjacent to the river. We used a quantitative deductive methodology for the problem analysis via mathematical and simulation methods that use the Ansys R1 2020 software. The final results showed that the relationships between window dimensions, air velocity, and the elevation of residential units can be mapped horizontally and vertically by paying attention to the locations of the window openings with respect to the wind direction and building layout. The horizontal and vertical maps showed repetition of the window dimensions every four floors, with the four floors in the center of the building (read vertically) having window dimensions one-quarter smaller than the four floors above and four floors below.

1. Introduction

Vertical housing in Indonesia consists of four types: public flats, particular needs flats, state flats, and commercial flats. Public flats are buildings that function as housing needs for low-income people; particular needs flats are buildings that function to meet particular needs (i.e., flats for disaster victims, fisherman’s flats, and flats for elderly residents); state flats are state-owned flats for civil servants or officials; and commercial flats aim to make a profit (in Indonesia, commercial flats are also called apartments) [1,2]. Public and commercial flats are the most common types of flats built in Indonesia. Based on the Central Statistics Agency, by 2021, the number of flats in Jakarta will reach 51 towers with 28,766 units (there are no definite records regarding the number of apartments, but according to local reports issued by the local daily news, until 2018, the supply of apartments was 228,000 units). The increase in the construction of the two types of flats is due to the need for residential houses being greater than the vacant land for residential functions.
The difference between public and commercial flats lies in the person in charge of development and the limitations in planning. Commercial flats are private parties with various residential units, namely studio, one-bedroom, two-bedroom, three-bedroom, and penthouse units while public flats built by the government have one residential unit type. Limitations regarding planning of public flats makes the planning of public flats complex.
Public flats are subsidized housing for low-income communities and urbanized communities. Based on the Regulation of the Minister of Public Works No. 05/PRT/M/2007, the planning for residential units of these flats focuses on human activities, consisting of a living room, two bedrooms, bathroom, kitchen, and laundry room with a room size of 30 to 36 m2 [3,4] in a prototype form that is used at all locations for flats in Indonesia. Apart from human activities, the construction volume determines the use, with sizes of 30 and 36 m2 for each flat. Based on the Law of the Republic of Indonesia No. 20 of 2011, residential units must reach at least 20% of the construction volume of the flats to reach the maximum number determined by the building height coefficient and the number of people who will live in the flats (in contrast to apartments, the government determines the occupants of each flats, such as the eviction community in Kampung Pulo, who live in the Jatinegara Barat flats).
The arrangement for the size of flats is that flats only consist of one type of dwelling. If the analysis of the flats shows that a size of 30 m2 is the right size in terms of the fulfilment of residential units, the basic coefficient of the building, and the coefficient of the building height, then the flats will only use the 30 m2 type in all residential units. Each tower of flats will have the same unit size. The purpose of using this prototype is to speed up development and reduce development costs. However, this prototype can create problems for the residents of the flats, especially in terms of fulfilling human thermal requirements.
The problems begin with the basic design, which is not focused on the environment and is unresponsive to the climate. Factors related to the environment and climate are assessed in softscape planning and ventilation calculation based on the Indonesian National Standard 03–6572-2001 concerning the Procedures for Planning for Ventilation and Air Conditioning Systems in Buildings [5,6]. This prototype has characteristics such as precast walls on the facade, monotonous window dimensions, and a limited size. The use of this prototype results in the performance of the building being poor, especially the performance of the windows on the facades of the buildings. The window dimensions, which have a monotonous size from the lower floors to the upper floors, cannot accommodate pressure differences, which affects the air velocity [7,8].
Window performance is a complex process that is affected by several factors and is often associated with the passive performance of low-rise and high-rise buildings. The performance of windows as inlets is one of the main factors affecting cross ventilation, the quality of air movement, and energy use [9,10,11]. For flats users, window performance is essential because it is related to the expenses that they must spend every month. Natural ventilation helps residents to reduce the cost of using cooling technologies, such as air conditioning and fans. In addition, natural ventilation supports air exchange, which helps to improve occupants’ health.
Window performance also relates to Indonesia’s emission reduction goals. Iwan Prijanto (chairman of the Green Building Council Indonesia), in the Sustainable Finance for Net Zero Green Building seminar, stated that the housing sector emits approximately 80% of carbon emissions [12]. High carbon emissions from residential houses occurs due to their daily use. In addition, high emissions are released by residential homes because there is a variety of energy being used in one place.
The window dimensions in high-rise flats are determined using the wind force ventilation formula listed in the Indonesian National Standard 03–6572-2001 concerning the Procedures for Planning for Ventilation and Air Conditioning Systems in Buildings, which calculates the airflow rates using area variables, inlet-free areas, air velocity, and the effectiveness of openings [13]. However, this calculation for airflow rates cannot establish the airflow rates within the comfort zone and instead provides new problems that threaten the safety of residents in the form of loose window connections due to a high air velocity. The main reason that the wind force ventilation formula cannot determine the window dimensions is that there is no variable height of the buildings (or elevation of the room), which affects the air velocity.
The wind is observed through pressure in both pure and building physics. Gravity affects the pressure and at different heights, gravity is different. The taller the height, the lower the pressure and the higher the air velocity. Another thing that affects the relationships between wind and buildings is temperature. Temperature affects the volume of air molecules and cause different speeds of wind under different temperature conditions [14,15,16,17,18]. Especially, in high-rise flats, the elevation of the residential units represents the height.
This study aimed to investigate the relationships between these four factors and their influence on the planning design of high-rise flats, especially in terms of the dimensions of the windows on building facades. Theoretically, windows on buildings allow wind and solar radiation to enter the residential units. Therefore, a poor window performance results in poor thermals in the rooms. Bernoulli’s principle of wind and height states that windows on the facades of high-rise flats do not function optimally because the window dimensions are the same from the lower floors to the upper floors even though the pressure from gravity affects the air velocity at different magnitudes [19].
Several studies on the performance of windows on the facades of flats in Indonesia have shown that four factors prevent windows from working optimally: the window location, wind direction, planning location, and elevation of residential units. Window location and wind direction are interrelated factors; for windows to have maximum wind suction capabilities, the location and placement of the windows must be in a positive field (i.e., facing the wind direction) [20,21]. The third and fourth factors indicate that the location and elevation of the residential units affect the responses of the flats’ design through air velocity changes, which affect the building facades [22]. Globally, studies on the relationship between air velocity and building height are still lacking. The relationship between air velocity and building height is a component of thermal comfort. Studies on thermal compliance, wind tunnels, and air movement have focused on heat transfer, represented by temperature differences [23], surface cooling [24], airflow [25,26], window location and position [27,28,29], and ventilation strategies [30], but have ignored the relationship between air velocity and building height.
Other research on windows as ventilation has shown other factors affecting ventilation performance: climate and form [31,32]. The climate factor that most affects the performance of windows in buildings is wind. Temperature and pressure are the main factors that affect the speed of the wind that hits buildings while altitude is the main factor affecting temperature and pressure. The factor that affects window performance the most is the angle of the opening, the type of the window, and the building form. The angle and type of window affect the size of the inlet-free area, which is how wind enters the room. Meanwhile, the building form and the placement of the buildings within one area affect the air velocity around the buildings.
Based on the existing phenomena, problems, and studies, the correct hypothesis is that the higher the building, the smaller the window dimensions, depending on the locations of the windows within the layout of the building, the direction of the wind, the actual wind speed, and the elevation of the residential units. This study aimed to determine the relationship between the air velocity and building height by comparing theories, mathematical calculations, and simulation results. This study specifically aimed to determine the best dimensions of windows in flats, in particular the Jatinegara Barat flats, Jakarta. The results provide preliminary ideas for further research on the same topic but in different locations and on different building forms. Long-term research aims to obtain new foundations that could become a theory for academics and practitioners.
Research on window dimensions has used mathematical and simulation methods. The mathematical methods use the power law, Torricelli, and MacFarlane formulae while the simulations use computational fluid dynamics (CFD). CFD simulation programs can describe the velocity of air movement and fluid flow around the surfaces of buildings, which are created by the actual boundaries in plans [33,34] and [35] (pp. 1–12). The variables used in both methods are the air velocity, size of existing residential units, dimensions of existing windows, temperature, humidity, and experimental windows. The novelty of this research lies in the horizontal and vertical mapping of window dimensions based on the relationships between air velocity, wind direction, the location of residential units, and the elevation of residential units.

2. Materials and Methods

The window in the Jatinegara Barat flats is a tipping window with a one-way opening (outward). The maximum window opening angle is 10 degrees, with the airflow entering through the window following the window opening angle. Based on Boutet’s theory, the tipping window allows the wind to move to the upper area of the room when it enters the room and split when it collides with the ceiling area (Figure 1).
The material object used in this study was the Jatinegara Barat flats, which has 16 floors with a unit size of 30 m2 (Figure 2a), tipping windows, and an existing inlet-free area of 0.941750378 m2 (Figure 2b). The Jatinegara Barat flats building has a zigzag shape with a floor finish level (FFL) height of 2.65 m for each residential unit (Figure 2c). The zigzag shape of the building allows for various kinds of wind motion (parallel, oblique, and turbulent) around the building (Figure 2d).
The variables used in this study of window dimensions were the actual air velocity, measured residential unit elevation, existing residential units, existing window dimensions, alternative window dimensions, dominant wind direction over 3 years (Figure 3), existing temperature of the residential units, and existing humidity of the residential units. The mathematical calculations used all variables, divided into two parts (as discussed in the method description). In contrast, the simulation experiments used the variables of the actual air velocity, dominant wind direction, existing residential units, dimensions of existing windows, and dimensions of alternative windows.

2.1. Basic Theory and Research Data

The basic factor that affects the determination of window dimensions is pressure. Pressure affects air velocity by differences in the gravity (in this study, represented by the elevation of the residential units), temperature, humidity, and roughness length [39] (Table 1). Theoretically, the comfortable air velocity for humans is 0.4 to 3.1 m/s (Beaufort theory) [40,41]; however, in tropical areas, the comfortable air velocity for humans is 0.25 to 1.5 m/s [42]. Air velocity measurements need to be carried out at a minimum height of 30 feet or ±9.4 m [43] above the ground; so, in the context of this study, these measurements started from the fourth floor of the flats (Figure 4). The control variables consisted of the comfortable limits for air movement, temperature, and humidity (Table 2).
The initial data for this study came from field measurements based on data from the Meteorology, Climatology, and Geophysics Agency. The determination of the elevation measurement points for the field measurements was based on horizontal and vertical simulations of the area using a maximum average air velocity of 5.6 m/s and the conditions around the Jatinegara Barat flats while the timing of the field measurements used the culmination time of 2021, which was 20 March 2021 and 23 September 2021. The field measurements used a hot wire anemometer to measure the air movement in the residential unit (Figure 5a) and a vane anemometer (multifunction) to measure the air velocity outside the building and the indoor temperature (Figure 5b).
The area outside the window, the window sill, the two bedrooms, and the living room were used as the measuring points for the field measurement (Figure 6a). The fourth floor is the first floor that has no obstructions (trees, buildings, etc.), the tenth floor is one of the floors in the middle of the building (read vertically from building section), and the fifteenth floor is on the top floor (Figure 6b), with the metric units facing the wind tunnel area and the urban area (Figure 6c).

2.2. Research Methods

This study of window dimensions based on the air velocity and elevation of residential units used two methods: mathematical and simulation methods. The mathematics method was divided into two stages while the simulation consisted of only one stage and this is a method that bridges the stages of mathematics 1 and mathematics 2.
Mathematics 1 consisted of a search for air velocity references and air movement references while mathematics 2 consisted of a search for the window dimension change range, determination of the vertical window dimension (section), and validation of the window dimension mapping generated by stage 1 and the simulation stage. The simulation process consisted of a search for air movement in the residential unit based on the alternative window. The simulation used a residential unit plan (2D drawing), with the final result of the simulation being the determination of the horizontal dimensions of the window (plan) (Figure 7).

2.2.1. Mathematical Calculation

The mathematical calculations consisted of two stages: before the simulations and after the simulations. The first stage of the mathematical calculations used the wind force ventilation from the Indonesian National Standard 03–6572-2001, power law formula, Torricelli formula, Bernoulli formula, and Macfarlane formula. The wind force ventilation used the free-inlet area, air velocity, and opening effectiveness (Equation (1)). The wind force ventilation formula was used to show the inability of the existing window to provide comfortable air movement:
Q = C v × A . V
where Q is the airflow rate, Cv is the aperture effectiveness (0.5–0.6 for perpendicular wind and 0.25–0.35 for skewed wind), A is the free-inlet area, and V is the air velocity.
The power law formula uses air movement variables at low elevations, air movement at reference elevations, low elevations, comparison elevations, and wind shear exponents [49]. The first step of the power law formula (Equation (2a)) aims to determine the 100% speed and the second step (Equation (2b)) aims to show the ideal reference air movement at each measured height, thereby becoming a validation tool for existing conditions:
V z = V g ( Z Z g ) α
where Vz is the mean air velocity at the height Z, Vg is the mean air velocity at the height Zg, Z is the height, Zg is the gradient height, and α is the wind shear exponent:
V r e f = V g ( Z r e f Z g ) α
where Vref is the reference speed, Vg is the mean of an area’s 100% air velocity, Zref is the reference height, Zg is the gradient height, and α is the wind shear exponent.
The following formulae are the Torricelli and Bernoulli formulae, which were used to determine the pressure and velocity, respectively. The Torricelli formula (Equation (3)) uses the residential unit elevation variable and aims to find the pressure at the measured height. In contrast, the Bernoulli formula uses the variables of gravity acceleration with a magnitude of 9.8067 m/s2 [50,51], air movement at a measured elevation, height at a measured elevation, and enthalpy (Equation (4)) to calculate the difference between the field measurement and theoretical calculations:
P = 76   c m H g h 100
where P is the air pressure and h is the height (in this study, it is the elevation of the measured residential unit):
v 2 2 + g h + ω
where g is acceleration of gravity, ϑ is the air velocity, h is the height, and ω is the enthalpy, with a value of 77.43 kJ/Kg (at 30 °C).
The Macfarlane formula (Equation (5)) calculates the reference air movement. The Macfarlane formula uses the variables of the temperature of the existing residential units and the humidity of the existing residential units. In the final results, the reference air movement is a control variable that validates the simulation results and data interpolation:
C v = 0.15   { D B T 27.2 + ( [ R H 60 ] / 10 ) 0.56 }
where Cv is the air movement reference, DBT is the dry-bulb temperature, and RH is the relative humidity.
The second stage of the mathematical calculations was data interpolation, which determined the air motion in the residential units on unmeasured floors. The data interpolation formula (Equation (6)) uses the variables of the von Karman constant, in the range of 0.4 to ± 0.02, the elevation of residential units, and the roughness length:
U ( z ) = U * K × l n   ( Z + Z 0 Z 0 )
where Uz is the friction velocity, K is the von Karman constant, z is the elevation of the occupancy unit, and Z0 is the roughness length.

2.2.2. CFD Ansys Simulation

In determining the window dimensions using the air velocity and elevation of residential units, simulation analysis has limited scope because the study object has already been built, and it is not possible to use simulations with total changes in its interior. The simulation focused on the window size and air velocity to enhance the finish details. The wall and door layouts were considered according to the existing conditions and did not change in the simulation. The limitation of the simulation is that the final results show the difference in the airflow that occurs when the size of the window dimensions are reduced while the interior conditions are the same as the existing conditions.
Ansys CFD is a simulation program that shows air movement and ventilation performance [35] (pp. 1–2) and [52,53]. This study used Ansys R1 2020 to measure the air motion formed by the window dimensions. The simulations used set equations for fluid dynamics, namely the continuity equation (Equation (7a,b), the x-y-z momentum equation (Equation (8a–c), and the energy equation (Equation (9)), which were solved automatically by the system in the simulation program. The simulations used polyhedral meshing and focused on local meshing to accommodate the flow around the model. The mesh around the building was denser compared to that in areas that were far from the building to increase the computational efficiency. In the residential unit model of the Jatinegara Barat flats, the model produced a mesh with 104,997 nodes and 476,716 elements, as shown in Figure 8:
t V ρ d V + A ρ V · d A = 0
ρ t + ρ · V = 0
( ρ u ) t + · ( ρ u V ) = p x + τ x x x + τ y x y + τ z x z + ρ f x
( ρ v ) t + · ( ρ v V ) = p y + τ x y x + τ y y y + τ z y z + ρ f y
( ρ w ) t + · ( ρ w V ) = p z + τ x z x + τ y z y + τ z z z + ρ f z
t [ ρ ( e + V 2 2 ) ] + · [ ρ ( e + V 2 2 ) V ] = ρ q ˙ ( ρ p ) x ( v p ) y ( w p ) z + ρ f · V
The computational settings of Ansys R1 2020 consisted of the following turbulence, pressure, and boundary conditions, and hybrid initialization settings:
  • Turbulence: k-omega SST (Share Stress Transport)
This study used k-omega modeling to model narrow slits, details, and high-pressure gradients. K-omega SST is a refinement of the k-epsilon model. The k-omega model is combined with k-epsilon to achieve good capability in the area around the wall, a low Reynold number value, flow with a low adverse pressure gradient, and insensitivity in the free stream area. The combination of k-omega and k-epsilon is then constructed b considering the turbulent shear stress transport equation, thereby increasing the predictive ability of flow separation in the adverse pressure gradient. This model is called k-omega SST [54]. This modification of SST served to accommodate free streams that were far away from the object;
  • Pressure-based
This study on the window dimensions in the Jatinegara Barat flats building used a couple models, so the approach to setting the proper pressure settings was pressure based. The pressure-based approach was the result of solving the continuity and momentum equations and could simultaneously solve all object coordinates;
  • Boundary conditions
Boundary conditions consist of two types: inlets and walls. For inlets, the input data were the air velocity outside the building, with the wind direction settings based on the wind rose analysis. The wind direction settings used vector settings. For the walls, we used a no-slip condition to represent the friction between the flow and the wall;
  • Hybrid initialization
The hybrid initialization method technically interpolates the values of boundary conditions by solving the Laplace equation to produce a velocity field that follows the geometry and a pressure field that smoothly corresponds.
The simulation data consisted of the existing residential units, existing windows, two alternative windows, elevation of the existing residential units (Table 3), actual air velocity, and wind direction (Table 4). On each measured floor, there were 9 simulated units in tower A and 8 simulated units in tower B, providing a total of 27 simulated units in tower A and 24 simulated units in tower B. All simulated residential units used for the experiments included three dimensions of the windows and the dominant wind direction.

3. Results

This study on the window dimensions in the facades of the Jatinegara Barat flats consisted of three parts: the air velocity patterns outside the building, the horizontal mapping of the window dimensions, and the vertical mapping of the window dimensions. The air velocity patterns outside the building were calculated from the results of the initial area simulations, initial simulations with the existing conditions, field measurements, power law calculations, Torricelli calculations, and Bernoulli calculations. The horizontal window dimension mapping was calculated from the simulations of the measured floors and the Macfarlane calculations. The vertical window dimension mapping was calculated from the data interpolation, which aimed to determine the air movement that was created in the residential units on the unmeasured floors.

3.1. Air Velocity Patterns

The initial simulations were area simulations using the average air velocity data from the Meteorology, Climatology, and Geophysics Agency, with a focus on the direction of the wind coming through the wind tunnel, namely west and east (Table 5). The regional simulation results showed that when the wind came from the east, the streamlined area on the 10th floor is longer than that on the 4th and 15th floors. The streamlined area on the 15th floor widened and produced the following two options: either the air velocity on the 15th floor was lower than that on the 10th floor or the air velocity on the 15th floor was faster than that on the 4th and 10th floors. According to Bernoulli’s pressure theory, both possibilities can occur in field measurements. This theory states that when streamlined areas are close together, the pressure decreases and the air velocity increases.
The Torricelli and Bernoulli calculations showed that the pressure and air velocity at the elevation of the measured residential units changed following the theory, namely, the higher the elevation of the measuring unit, the lower the pressure and the higher the air velocity. This result was different from the field measurements, which showed that the 10th floor had the highest air velocity (Table 6) and proved that the streamline areas were close to each other in the wind tunnel on the 10th floor (Figure 9).

3.2. Existing Window and Calculation of the Indonesian National Standard

The calculations using the Indonesian National Standard show that the existing window can allow comfortable air movement when the wind speed reaches 7.5 m/s in the perpendicular wind condition (Table 7). The field conditions show that the residents close the windows at a wind speed of 7.5 m/s; so, the existing windows are not functioning optimally. Another thing that makes the calculation using the wind force ventilation formula in the Indonesian National Standard less precise is that the wind does not continuously have the same direction of motion; so, it is difficult to determine how often the wind is perpendicular wind or skewed wind.

3.3. Window Dimensions: Horizontal Mapping

The comparison of the existing conditions to the power law calculations was the first step in the simulation experiments, which aimed to model the worst conditions. The worst conditions were modeled by the initial data from the experiment simulations of the residential units and alternative window dimensions in the building. The power law calculations showed that one of the problems is that the window performance of the building was not optimal. The existing air velocity is above the reference air velocity. The calculations of the 100% air velocity are shown in Table 8 for the following conditions:
  • 4.9–8.6 m/s on the 4th floor of tower A;
  • 8.6–9.3 m/s on the 10th floor of tower A;
  • 8.3–9.4 m/s on the 15th floor of tower A;
  • 5–8.2 m/s on the 4th floor of tower B;
  • 7.5–10 m/s on the 10th floor of tower B;
  • 7.6–10.2 m/s on the 15th floor of tower B.
The 100% air velocity produced a reference air velocity that was lower than the field measurements, namely 2.1–3.8 m/s on the 4th floor of tower A, 4.6–5 m/s on the 10th floor of tower A, 4.9–5.6 m/s on the 15th floor of tower A, 2.1–3.6 m/s on the 4th floor of tower B, 4–5.4 m/s on the 10th floor of tower B, and 4.5–6.1 m/s on the 15th floor of tower B as shown in Table 9.
Examination of the existing conditions using the power law calculations shows that the existing conditions were worse than the power law calculations. Therefore, the data used in the simulations consisted of the existing conditions that occurred in the field. According to the existing temperature and humidity conditions in the Macfarlane calculations, the appropriate dimensions for windows on the facades of high-rise buildings require data on the reference air movement. The Macfarlane calculation results show that the units located on the 4th floor of tower A required air movement of 0.41–0.45 m/s, those on the 10th floor of tower A required air movement of 0.46–0.54 m/s, those on the 15th floor of tower A required air movement of 0.29–0.46 m/s, those on the 4th floor of tower B required air movement of 0.29–0.54 m/s, those on the 10th floor of tower B required air movement of 0.36–0.39 m/s, and those on the 15th floor of tower B required air movement of 0.41–0.57 m/s, as shown in Table 10.
The simulations showed two results: the maximum and minimum air movement formed through the dimensions of the window in the facades of the residential units. The control variable used was the academic standard and the results of the Macfarlane calculations with air movement coupled with the control variable showed the minimum air movement. The main reason for only using the minimum air movement is that maximal air movement never occurs under actual conditions. In air movement conditions of more than 0.5–0.55 m/s, occupants tend to close the windows (actual conditions are conditions that occur in the field with non-verbal data because the data come from the experience of researchers after repeatedly observing the objects under study for one year). In the simulations, there were four simulated units on each floor in each tower. Each unit was simulated with the highest air velocity based on the wind direction as shown in Figure 10.
The number of residential units was chosen according to the order of the field measurements to ensure that the data and simulations produced followed the field measurement process. The simulation results showed the final results, as shown in Table 11.
Based on the simulation analysis of the measured units on the 4th floor, 10th floor, and 15th floor, the properties of the dimensions of the inlet-free area were (Figure 11):
  • The window dimensions on the 4th floor were the same as those on the 15th floor, except for unit 2 in tower A, which faced away from the wind direction;
  • The window dimensions of unit 1 in tower A were the same as those of unit 4 in tower B;
  • The window dimensions of unit 3 and unit 4 in tower A were the same as those of unit 2 in tower B;
  • The windows dimensions of unit 1 in tower B were the same as those of unit 3 in tower B.
The simulation results indicate that the window dimensions could be determined horizontally based on the unit position within the building layout, wind direction, air velocity, and the elevation of the residential unit. Data interpolation was carried out in the next step to validate the horizontal mapping and the repetition of the window dimensions on the 4th and 15th floors.

3.4. Window Dimensions: Vertical Mapping

The data interpolation aimed to validate the results of the horizontal mapping and determine the air movement that occurred on all floors of the residential units. The data interpolation consisted of two stages. The first stage aimed to determine the range of the repetition of the vertical window dimensions. In contrast, the second stage aimed to validate the repetition range of the window dimensions generated by the first stage of data interpolation. The first stage of data interpolation used data on the height of each residential floor (in the case of the Jatinegara Barat flats, the initial residential floor is the 3rd floor) and the air movement generated by the simulations on the measured floors (4th floor, 10th floor, and 15th floor).
The first stage showed that the window dimensions were repeated every four floors, i.e., floors 3 to 6 had larger window dimensions than floors 7 to 10 and had the same window dimensions as floors 11 to 16. Floors 15 to 16 had the same window dimensions as floors 7 to 10, as shown in Table 12; however, due to the limitation of the construction costs of the flats, the window dimensions on floors 15 and 16 were the same as those on floors 11 to 14, with the condition of air movement within the comfort zone.
The results from the first-phase data interpolation showed the occurrence of air movement that was not within the comfort zone due to the repetition of the window dimensions every four floors, namely, unit 1 in tower A with a southerly wind direction, unit 2 in tower A with a southeasterly wind direction, unit 3 in tower A with a southeasterly wind direction, unit 1 in tower B with a southwesterly wind direction, and unit 3 in tower B with a southwesterly or southeasterly wind direction. Therefore, to validate the horizontal mapping and repetition generated by the first-phase data interpolation, a second interpolation was carried out using the measured floor height data and the simulated air movement generated by the first-phase data interpolation on the measured floors and selected window dimensions, as shown in Table 13.
The second-phase data interpolation shows that the air movement that occurred on the 3rd floor to the 16th floor was within the comfort zone according to academic standards and the Macfarlane calculations, except for unit 3 in both tower A and tower B with a southeasterly wind direction, which showed air movement of 0 m/s. The second-phase data interpolation results show that overall, the air movement in the units was in the range of 0.2 to 0.6 m/s, with air movement in tower A ranging from 0.2 to 0.52 m/s and air movement in tower B ranging from 0.2 to 0.62 m/s, as shown in Figure 12.

4. Discussion

The final results from each analysis process, i.e., both the mathematical and simulation methods, showed consistent results and supported each other. The final results showed that the window dimensions were repeated every four floors, with the window dimensions of the top floor being repeated for six floors. The consistency of the results from each method was as follows:
  • The results of the comparison of the first-phase data interpolation (with simulated air movement) in the southwesterly direction with the control variable showed that the air movement was predominantly within the comfortable limit, with the lowest air movement being 0.2 m/s and the highest being 0.6 m/s. In the simulated sections of units 1 and 3 in tower B on floors 7 to 10, the air movement exceeded the standard; however, the second-phase data interpolation showed that the speed reduced to 0.39–0.63 m/s;
  • The results of the comparison of the first-phase data interpolation in the southerly direction with the control variable showed that the air movement was predominantly within the comfortable limit, with the lowest air movement being 0.2 m/s and the highest being 0.61 m/s. The second-phase data interpolation showed that the air movement stayed within the same range, namely 0.2 to 0.61 m/s;
  • The comparison results of the first-phase data interpolation in the southeasterly direction with the control variable showed that the air movement was predominantly within the comfortable limit, with the lowest air movement being 0 m/s and the highest being 0.6 m/s. The second-phase data interpolation showed that the air movement was reduced and stayed within the comfort zone, i.e., 0.2 to 0.505 m/s;
  • The comparison results of the first-phase data interpolation in the easterly direction with the control variable showed that the air movement was predominantly within the comfortable limit, with the lowest air movement being 0.2 m/s and the highest being 0.6 m/s; however, the second-phase data interpolation showed that the air movement was lower, i.e., 0.2 to 0.534 m/s.
The comparison data for the southwest, south, southeast, and east directions showed that the air movement produced by the experiments (simulations and mathematical calculations) was better than the existing conditions (average air movement for March and September) as shown in Figure 13.
Based on the results of the experiments using mathematical methods and CFD Ansys, the design for window dimensions on the facades of high-rise buildings (in this study, the Jatinegara Barat flats) must pay attention to several things, namely:
  • Building placement;
  • Wind direction (records are needed for at least the last three years to establish the dominant wind direction);
  • The actual air velocity at measured altitudes (not only at an altitude of ±9.14 m);
  • The elevation of the residential units.
  • For the vacant land, the layout of buildings around the planned location provides the existing conditions. The layout of buildings is a factor that aims to create a comfortable wind flow around the buildings so that the air movement inside the buildings is within the comfort zone. For land that already has buildings, the layout of the existing buildings provides the existing conditions. Buildings affect the direction of wind that hits the buildings. The building placement also needs to consider the orientation of the buildings with respect to the direction of wind and sun, which affects the thermal conditions of the residential units [31,55,56,57].
The building placement, direction of the wind, actual wind speed, and measured height are the four main factors because they are directly related to the pressure that affects the wind around buildings. In addition, these four factors affect the formation of negative and positive pressure areas around the inlets and outlets. This pressure affects the wind density and suction areas. The final results of this study show that the window dimensions on the 10th floor were one-quarter smaller than those on the 4th and 15th floors. The window dimensions were repeated every four floors. In addition, the window dimensions of the units that had inlets facing away from the direction of the wind were different; namely, the window dimensions on the 4th floor were one-quarter larger than those on the 10th floor and the window dimensions on the 10th floor were one-quarter smaller than those on the 15th floor as the window dimensions were repeated every 4 floors (10th floor < 4th floor = 15th floor), as shown in Figure 14. Based on the final results, the initial hypothesis stating that window dimensions are reduced depending on the position of the window with respect to the direction of the wind, the layout of the building, the actual wind speed, and the elevation of the residential unit was proven correct.
Determining the window dimensions on flats’ facade (which is the function of the object of study used) in the planning of new buildings can increase the quality of the room through the layout of the room. This increase occurs because the dimensions of the window affect the location of the walls and doors in the inner room so that each part of the building (windows, doors, and walls) supports each other to form of a thermal comfort zone in the interior through cross ventilation (this statement is the basis for hypotheses in further research). In contrast to the new building, the determination of the window dimensions on the facade of the old building can improve the quality of the room by adjusting the window dimensions to the existing layout. The layout of the walls and doors cannot be changed.

5. Conclusions

The problem with the Jatinegara Barat high-rise flats is that the performance of the windows on the building’s facade is not optimal. The leading cause of this sub-optimal window performance is the window dimensions, which cannot accommodate the air velocity around the building. The main concern with high-rise buildings is the gradation of air velocity around the buildings due to different pressures at specific heights.
In the case of the Jatinegara Barat flats, the location of the building forms a wind tunnel between the two towers, meaning that pressure differences occur in the wind tunnel area. In the field measurements, the air velocity in the wind tunnel area was the highest at the 10th floor. The height of the 10th floor was used as the specific height for the center area of the building when viewed vertically. The highest speed at the 10th floor showed that the streamline was closer and provided a pressure difference, with the air velocity at the 10th floor being higher than that at the 15th floor.
The mathematical methods and simulations showed that the window dimensions could be read horizontally and vertically. The horizontal mapping showed a relationship between the locations of the residential units with respect to the wind direction and the building layout. The horizontal mapping of the window dimensions also showed the following results:
  • Units located at the end of the wind tunnel had the same window dimensions (unit 1 in tower A and unit 4 in tower B);
  • Units located in the wind tunnel area had the same window dimensions (unit 3 in tower A, unit 4 in tower A, and unit 2 in tower B);
  • Units facing the urban area and the direction of the wind had the same window dimensions (unit 1 and unit 3 in tower B);
  • Units facing the urban area and facing away from the direction of the wind had the same window dimensions but were different from those of the units in the wind tunnel and the units facing the urban area and the direction of the wind (unit 2 in tower A).
The vertical mapping showed that the horizontal mapping is the basis for vertical mapping. The vertical mapping showed that the window dimensions were repeated every four floors, i.e., the window dimensions on the 4th floor were the same as the window dimensions on the 3rd floor to 6th floor and the 11th floor to 16th floor, and the window dimensions on the 10th floor were the same as the window dimensions on the 7th floor to the 9th floor. The vertical mapping showed a relationship between the locations of the residential units, the elevation of the residential units, and the actual wind speed that occurred around the facade of the residential units.
The final results showed that the windows on the 10th floor were one-quarter smaller than those on the 4th and 15th floors. Knowledge of this tendency to reduce window dimensions could be used for further research that is more specific to the pressure that occurs in the inlet free-area of the window dimensions used. Further research should be carried out to validate the horizontal and vertical window dimension mapping of the Jatinegara Barat flats by mapping other flats with different building shapes and locations. Long-term research aims to obtain a new formula that can be used as a calculation tool and a planning concept for flats that accounts for the surrounding environmental conditions. In addition, long-term research can also be carried out with different objects of study such as:
  • Office buildings, malls, and other public buildings (objects of study have been built): analysis of window dimensions based on building shape, building thickness, building height, and building user activities;
  • Urban (study object has been built): analysis of the existing condition of an area based on the layout of the building part of each existing building function, changes in the thermal condition of the area (urban heat island issue) through changes in window dimensions, and the window layout of each existing building function;
  • Vernacular: window dimension analysis based on the shape and position of the window (height and location). Several traditional Indonesian houses have sloping walls, so the window angle is also crucial in this case.
This research will provide benefits to several related parties, such as:
  • Researchers, as a starting point for knowledge about the pattern of window dimensions in high-rise buildings that is focused on wind, especially high-rise flats;
  • Academics, as a basic idea that can be used for the development of science, especially in the fields of building comfort science, building physics, building technology, and environmental engineering;
  • Practitioners, as a reference and basic conception in the design of a high-rise residential building;
  • Stakeholders, as a basis of reference, which can be used in developing standards used in high-rise buildings.

Author Contributions

Conceptualization, S.S.N., G.H. and R.R.T.; methodology, S.S.N. and G.H.; software, S.S.N.; validation, G.H. and R.R.T.; formal analysis, S.S.N.; investigation, S.S.N. and G.H.; resources, S.S.N. and R.R.T.; data curation, S.S.N.; writing—original draft preparation, S.S.N., G.H. and R.R.T.; writing—review and editing, S.S.N.; visualization, S.S.N.; supervision, G.H. and R.R.T.; project administration, S.S.N.; funding acquisition, S.S.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Not applicable.

Acknowledgments

I would like to thank the Dinas Pekerjaan Umum dan Penataan Ruang for its support for the fulfillment of data as built drawings for study objects and field survey permits, Gagoek Hardiman and Rumiati Rosaline Tobing for they guidance during research and journal writing, as well as residents of the flats for their positive response during the field survey.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Peraturan Pemerintah Republik Indonesia No. 13 Tahun 2021; Penyelenggaraan Rumah Susun; Presiden Republik Indonesia: Jakarta, Indonesia, 2021; p. 2.
  2. Undang-Undang Republik Indonesia No. 20 Tahun 2011; Rumah Susun; Presiden Republik Indonesia: Jakarta, Indonesia, 2011; p. 2.
  3. Nur Fajri Alfata, M.; Rani, W.; Fanny, K.; Ridlo, H.M.; Emir, H.R.M.; Ifta, I.U.; Arif, N.; Pandu, W.Y.; Resha, F. Rumah Susun di Indonesia. In Pedoman Desain Pasif Rumah Susun: Untuk Iklim Tropis Panas Lembap di Indonesia, 1st ed.; Prasasto, S., Agung, M.N., Eds.; Balai Sains Bangunan, Direktorat Bina Teknik Permukiman dan Perumahan Kementerian Pekerjaan Umum dan Perumahan Rakyat: Jakarta, Indonesia, 2020; Chapter 2; pp. 13–14. [Google Scholar]
  4. Dadang, R.; Pangihutan, M.; Tjuk, K.; Maliki, M.M.; Dedy, P.; Dwitri, W.; Nazir, A.; Bayu, W.; Niken, P. Rumah Sederhana Sehat. In Jejak Langkah Hunian Layak Indonesia, 1st ed.; Imam, S.E., Ed.; PT. Mediatama Saptakarya: Jakarta, Indonesia, 2019; Chapter 3; pp. 89–91. [Google Scholar]
  5. SNI 03-7013-2004; Tata Perencanaan Fasilitas Lingkungan Rumah Susun Sederhana. BSN: Jakarta, Indonesia, 2004.
  6. Peraturan Menteri Pekerjaan Umum No. 05/PRT/M/2007; Pedoman Teknis Pembangunan Rumah Susun Sederhana Bertingkat Tinggi; Direktorat Penataan Bangunan dan Lingkungan Direktorat Jenderal Cipta Karya: Jakarta, Indonesia, 2007; pp. 8–11.
  7. Terry, S.B. Meso-climatic Air Movement. In Controlling Air Movement: A Manual for Architects and Builders, 1st ed.; Nadine, M., Galen, H.F., Eds.; McGraw-Hill Book Company: New York, NY, USA, 1987; Chapter 5; pp. 43–52. [Google Scholar]
  8. Pasquill, P. Aerodynamics Characteristic of Atmospheric Boundary Layers. J. Fluid Dyn. 1972, 51, 619–623. [Google Scholar]
  9. Amy, K.; Shuoqi, W.; Ji-Eun, K.; Dorothy, R. Indoor/Outdoor Environmental Parameters and Window-Opening Behavior: A Structural Equation Modeling Analysis. Buildings 2019, 9, 94. [Google Scholar] [CrossRef] [Green Version]
  10. Rijal, H.B.; Tuohy, P.; Humphreys, M.A.; Nicol, J.F.; Samuel, A.; Clarke, J. Using Results from Field Surveys to Predict the Effect Open Windows on Thermal Comfort and Energy Use in Buildings. Energy Build. 2007, 39, 823–836. [Google Scholar] [CrossRef] [Green Version]
  11. Choi, J.H.; Loftness, V.; Aziz, A. Post-occupancy Evaluation of 20 Office Buildings as Basis for Future IEQ Standard and Guidelines. Energy Build. 2012, 46, 167–175. [Google Scholar] [CrossRef]
  12. Iwan, P. Alasan Netzero pada Piramida Terbawah: Mengapa Kita Semua harus Berjuang untuk Netzero? In Proceedings of the Seminar Permukiman Sehat Netzero, Universitas Katolik Parahyangan, Bandung, Indonesia, 27 August 2022. [Google Scholar]
  13. SNI 03-6572-2001; BSN: Jakarta, Indonesia, 2001; pp. 5–6.
  14. Julian, A. Relations Among Wind, Temperature, Pressure, and Density with Particular Reference to Monthly Average. Am. Meteorol. Soc. 1967, 95, 531–539. [Google Scholar]
  15. Daniel, J.J. Vertical Profiles of Pressure and Temperature. In Introduction to Atmospheric Chemistry, 1st ed.; Princeton University Press: Princeton, NJ, USA, 1999; Chapter 2; pp. 14–15. [Google Scholar]
  16. David, H.; Robert, R.; Jearl, W. Measurement, Motion Along a Straight Line, Vectors, Motion in Two and Three Dimension, Force and Motion I, Force and Motion II, Center of Mass and Linear Momentum, Gravitation, Fluid. In Fundamental of Physic Extended, 11th ed.; Ayantika, C., Alden, F., Eds.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2018; Volume 1, pp. 1–148, 214–256, 354–412. [Google Scholar]
  17. Klaus, W. Misinterpretation of Bernoulli’s Law. Researchgate. 2015. Available online: https://www.researchgate.net/publication/303974495 (accessed on 25 August 2022).
  18. Airfoils and Airflow. Available online: https://www.av8n.com/how/htm/airfoils.html (accessed on 27 August 2022).
  19. John, D.A. Incompressible Flow over Airfoils. In Fundamentals of Aerodynamics, 6th ed.; McGraw-Hill Education: Hightstown, NY, USA, 2017; Chapter 4; pp. 321–421. [Google Scholar]
  20. Hedy, C.I. Kinerja Ventilasi Pada Hunian Rumah Susun Dupak, Bangunrejo. J. Dimens. Inter. 2008, 6, 9–23. [Google Scholar]
  21. Rahmawati, R.; Arif, K.F.A.; Fathia, K.A. Penghawaan Alami Terkait Sistem Ventilasi Terhadap Kenyamanan Termal Rumah Susun Industri Dalam. Reka Karsa 2016, 4, 1–12. [Google Scholar]
  22. Sally, S.N. Pengaruh Orientasi Bangunan dan Kecepatan Angin Terhadap Bentuk dan Dimensi Filter Pada Fasad Bangunan Rumah Susun: Studi Kasus Rumah Susun Marunda. E-J. Grad. Unpar 2014, 1, 75–89. [Google Scholar]
  23. Zixuan, C.; Ahmed, W.A.H.; Imriyas, K.; Haddad, A. Optimizing Window Design on Residential Building Facades by Considering Heat Transfer and Natural Lighting in Non-tropical Regions of Australia. Buildings 2020, 10, 1–27. [Google Scholar]
  24. Qin, S.Y.; Cui, X.; Yang, C.; Jin, L.W. Thermal Comfort Analysis of Radiant Cooling Panels with Dedicated Fresh-Air System. Indoor Built Environ. 2020, 30, 1–13. [Google Scholar] [CrossRef]
  25. Per, H.; Kjeld, S.; Peter, V.N. Characteristic of Airflow from Open Windows. Build. Environ. 2001, 26, 859–869. [Google Scholar]
  26. Daniel, S.; Tobias, G.E.; Maximilian, H.; Ursula, E. Experimental Validation of Different Airflow Correlations for Natural Single Sided Ventilation. Energy Procedia 2015, 78, 2838–2843. [Google Scholar]
  27. Gao, C.F.; Lee, W.L. Evaluating the Influence of Window Type on the Natural Ventilation Performance of Residential Buildings in Hong Kong. Int. J. Vent. 2011, 10, 227–238. [Google Scholar] [CrossRef]
  28. Jinkyun, C.; Changwoo, Y.; Yundeok, K. Effective Opening Area and Installation Location of Windows for Single Sided Natural Ventilation in High-rise Residences. J. Asian Archit. Build. Eng. 2012, 11, 391–398. [Google Scholar]
  29. Bangalee, M.Z.I.; Miau, J.J.; San-Yih, L.; Mohammad, F. Effects of Lateral Window Position and Wind Direction on Wind-Driven Natural Cross Ventilation of a Building: A Computational Approach. J. Comput. Eng. 2014, 1, 310358. [Google Scholar] [CrossRef] [Green Version]
  30. Babak, R.; Martin, J.T.; Regina, B.; Andy, V.D.D. Natural Summer Ventilation Strategies for Energy-saving in High-rise Buildings: A Case Study in the Netherlands. Int. J. Vent. 2020, 19, 25–48. [Google Scholar]
  31. Zhaoming, D.; Weihong, G.; Weicong, L.; Xuyi, G. A Study on the Optimization of Wind Environment of Existing Villa Buildings in Lingnan Area: A Case Study of Jiangmen’s “Yunshan Poetic” Moon Island Houses. Buildings 2022, 12, 1304. [Google Scholar] [CrossRef]
  32. Mohamed, A.; Yassine, K.; Abdelaziz, M. Optimization of Passive Design Features for a Naturally Ventilated Residential Building According to the Bioclimatic Architecture Concept and Considering the Northern Morocco Climate. Build. Simul. 2020, 1, 677–689. [Google Scholar] [CrossRef]
  33. Sawako, K.; Roland, B.; Karen, W.; Suresh, N. Computational Fluid Dynamics for Architectural Design. Archit. Des. 2013, 83, 118–123. [Google Scholar]
  34. Francis, H.H.; Jacob, E.F. Computer Experiments in Fluid Dynamics. Sci. Am. 1965, 212, 104–111. [Google Scholar]
  35. Caesar, W. Teori Aliran Fluida dan Komputasi Finite Volume Method. In Teori dan Best Practices Computational Fluid Dynamics (CFD), 1st ed.; AE Publisher: Yogyakarta, Indonesia, 2021; pp. 1–12. [Google Scholar]
  36. Steven, V.S. Climatic Desing Archetypes. In Introduction to Architectural Science the Basis of Sustainable Design, 1st ed.; Elsevier: Oxford, UK, 2004; Part 1; pp. 66–68. [Google Scholar]
  37. DeKay, M.; Brown, G.Z. Favorite Design Strategies. In Sun, Wind, & Light: Architectural Design Strategies, 3rd ed.; Bennett, S., Ed.; Wiley & Sons, Inc.: Hoboken, NJ, USA, 2018; Part VI; pp. 238–262. [Google Scholar]
  38. Terry, S.B. Micro-climatic Air Movement. In Controlling Air Movement: A Manual for Architects and Builders, 1st ed.; Nadine, M., Galen, H.F., Eds.; McGraw-Hill Book Company: New York, NY, USA, 1987; Chapter 6; pp. 55–84. [Google Scholar]
  39. Richard, M.A. Architectural Aerodynamics; Applied Science Ltd.: London, UK, 1997. [Google Scholar]
  40. Fransesco, C.; Lovatelli, A. Site Selection. In Aquaculture Operations in Floating HDPE Cages. A Field Handbook, 1st ed.; FAO Fisheries and Aquaculture Technical Paper; FAO: Rome, Italy, 2015; Chapter 1; p. 12. [Google Scholar]
  41. Mangunwijaya, Y.B. Pengaruh Angin. In Pengantar Fisika Bangunan, 3rd ed.; Djambatan: Yogyakarta, Indonesia, 1988; Part 3; pp. 62–63. [Google Scholar]
  42. Lippsmeier, G. Persyaratan-persyaratan Kenyamanan. In Bangunan Tropis, 2nd ed.; Kiran, M., Paulhans, P., Johan, B., Paula, C., Rita, F., Kazuo., O.H.S., Eds.; Erlangga: Jakarta, Indonesia, 1980; Chapter 3; pp. 36–38. [Google Scholar]
  43. Brown, G.Z. Angin. In Matahari, Angin, dan Cahaya, 1st ed.; Aris, K.O., Ed.; Intermatra: Daerah Khusus Ibukota Jakarta, Indonesia, 1987; Chapter A; p. 26. [Google Scholar]
  44. Itard, L.C.M. Thermal Comfort: Air Temperature and Velocity. Available online: https://learning.edx.org/course/course-v1:DelftX+ECObuild3x+2T2020/block-v1:DelftX+ECObuild3x+2T2020+type@sequential+block@b9a3a1858aa54479a319288437e206e8/block-v1:DelftX+ECObuild3x+2T2020+type@vertical+block@2a3745f8dd3447e8b39ae880484b8e60 (accessed on 21 April 2021).
  45. BCA Green Mark. Green Mark for Residential Buildings; Building Construction Authority: Singapore, 2017; p. 8.
  46. Prasasto, S. Istilah-istilah dan Pengertian dalam Ventilasi. In Fisika Bangunan, 1st ed.; Sigit, S., Ed.; Andi Offset: Yogyakarta, Indonesia, 2009; Chapter 1; p. 9. [Google Scholar]
  47. Heinz, F.; Antonius, D.; Darmawan, A.M.S. Nilai Lembap Udara. In Ilmu Fisika Bangunan: Pengantar Pemahaman Cahaya, Kalor, Kelembapan, Iklim, Gempa Bumi, Bunyi, dan Kebakaran, 1st ed.; Kanisius: Yogyakarta, Indonesia, 2008; Part 3; pp. 52–56. [Google Scholar]
  48. Otto, H.K.; Ingersoll, T.G.; Alan, M.; Steven, V.S. Composite or Monsoon Climate. In Manual of Tropical Housing and Building: Climatic Design, 1st ed.; Universities Press: Himayatnagar, India, 1973; Chapter 1. [Google Scholar]
  49. Richard, M.A. Architectural Aerodynamics: Handbook of Architectural; Applied Science Ltd.: London, UK, 1995. [Google Scholar]
  50. Barry, N.T.; Ambler, T. The International System of Units. NIST Spec. Publ. 2008, 330, 52. [Google Scholar]
  51. Bureau International des Poids et Mesures. The International System of Units; BIPM: Paris, France, 2006; pp. 142–143. ISBN 92-822-2213-6. [Google Scholar]
  52. Ladeinde, F.; Michelle, D.N. CFD Application in the HVAC & R Industry. ASHRAE J. 1997, 39, 44–48. [Google Scholar]
  53. ANSYS Inc. ANSYS Fluent Theory Guide; ANSYS: Canonsburg, PA, USA, 2008. [Google Scholar]
  54. Menter, F.R. Improved Two-equation K-ω Turbulence Model for Aerodynamic Flows; NASA TM-103975; NASA AMES Research Center: Moffett Field, CA, USA, 1992.
  55. Rania, E.; Hamdy, H. Impact of Window Parameters on the Building Envelope on the Thermal Comfort, Energy Consumption and Cost and Environment. Int. J. Vent. 2019, 19, 233–259. [Google Scholar] [CrossRef]
  56. Weixin, Z.; Panu, M.; Sami, L.; Simo, K.; Juha, J.; Risto, K. Numerical and Experimental Study on the Indoor Climate in a Classroom with Mixing and Displacement Air. Buildings 2022, 22, 1314. [Google Scholar] [CrossRef]
  57. Simos, Y. Physics and Architecture: Issues of Knowledge Transfer and Translation to Design. Sol. Wind. Technol. 1989, 6, 301–308. [Google Scholar]
Figure 1. Airflow on the tipping window.
Figure 1. Airflow on the tipping window.
Buildings 12 01954 g001
Figure 2. General data for the Jatinegara Barat flats: (a) residential unit plan (30 m2/unit); (b) existing windows and inlet-free areas; (c) residential unit sections; (d) illustration of the air velocity around the Jatinegara Barat flats (inspired by the basic theories of Szokolay, DeKay, and Boutet [36,37,38]).
Figure 2. General data for the Jatinegara Barat flats: (a) residential unit plan (30 m2/unit); (b) existing windows and inlet-free areas; (c) residential unit sections; (d) illustration of the air velocity around the Jatinegara Barat flats (inspired by the basic theories of Szokolay, DeKay, and Boutet [36,37,38]).
Buildings 12 01954 g002aBuildings 12 01954 g002b
Figure 3. Wind rose data for the last 3 years.
Figure 3. Wind rose data for the last 3 years.
Buildings 12 01954 g003
Figure 4. Experiment floors and starting point of measurements.
Figure 4. Experiment floors and starting point of measurements.
Buildings 12 01954 g004
Figure 5. Measurement tools: (a) Hot wire anemometer; (b) vane anemometer.
Figure 5. Measurement tools: (a) Hot wire anemometer; (b) vane anemometer.
Buildings 12 01954 g005
Figure 6. General data of the field measurement: (a) Field measurement points on residential units; (b) field-measured floor position; (c) field-measured residential unit.
Figure 6. General data of the field measurement: (a) Field measurement points on residential units; (b) field-measured floor position; (c) field-measured residential unit.
Buildings 12 01954 g006aBuildings 12 01954 g006b
Figure 7. Research workflow.
Figure 7. Research workflow.
Buildings 12 01954 g007
Figure 8. Local meshing at tower B as a depiction of the local meshing in the simulation.
Figure 8. Local meshing at tower B as a depiction of the local meshing in the simulation.
Buildings 12 01954 g008
Figure 9. Comparison of the field calculations and measurements: (a) Air velocity from the field measurement; (b) Torricelli equation; and (c) Bernoulli equation.
Figure 9. Comparison of the field calculations and measurements: (a) Air velocity from the field measurement; (b) Torricelli equation; and (c) Bernoulli equation.
Buildings 12 01954 g009
Figure 10. Residential unit numbers in the simulation.
Figure 10. Residential unit numbers in the simulation.
Buildings 12 01954 g010
Figure 11. The horizontal mapping of the window dimensions.
Figure 11. The horizontal mapping of the window dimensions.
Buildings 12 01954 g011
Figure 12. Second-phase data interpolation results: (a) unit 1—tower A; (b) unit 2—tower A; (c) unit 3—tower A; (d) unit 4—tower A; (e) unit 1—tower B; (f) unit 2—tower B; (g) unit 3—tower B; (h) unit 4—tower B.
Figure 12. Second-phase data interpolation results: (a) unit 1—tower A; (b) unit 2—tower A; (c) unit 3—tower A; (d) unit 4—tower A; (e) unit 1—tower B; (f) unit 2—tower B; (g) unit 3—tower B; (h) unit 4—tower B.
Buildings 12 01954 g012
Figure 13. Comparison of the field measurement and experimental results: (a) 4th floor—tower A; (b) 10th floor—tower A; (c) 15th floor—tower A; (d) 4th floor—façade B; (e) 10th floor—tower B; (f) 15th floor—tower B.
Figure 13. Comparison of the field measurement and experimental results: (a) 4th floor—tower A; (b) 10th floor—tower A; (c) 15th floor—tower A; (d) 4th floor—façade B; (e) 10th floor—tower B; (f) 15th floor—tower B.
Buildings 12 01954 g013aBuildings 12 01954 g013b
Figure 14. The vertical mapping of the window dimension: (a) Unit 1 in tower A; (b) Unit 2, 3 in tower A and unit 1 in tower B; (c) unit 3 in tower B, unit 2 in tower B, and unit 4 in tower A; (d) unit 4 in tower B.
Figure 14. The vertical mapping of the window dimension: (a) Unit 1 in tower A; (b) Unit 2, 3 in tower A and unit 1 in tower B; (c) unit 3 in tower B, unit 2 in tower B, and unit 4 in tower A; (d) unit 4 in tower B.
Buildings 12 01954 g014
Table 1. Values for the roughness length [39].
Table 1. Values for the roughness length [39].
VariableZg (m)Zo (m) α
Open sea, ice, tundra, desert2500.0010.11
Open country with low bushes2000.030.15
Suburban areas, small towns, well-wooded areas4000.30.25
Numerous tall buildings, city centers, dense industrial developments50030.36
Table 2. Research control variables.
Table 2. Research control variables.
VariableComfort Zone
Air movement [44,45]0.2–0.6 m/s
Temperature [42]21–27 °C
Humidity [46,47,48]17–80%
Table 3. Simulation data.
Table 3. Simulation data.
VariableData
Residential unitBuildings 12 01954 i001Buildings 12 01954 i002
Field measurementSimulated unit
KeyplanBuildings 12 01954 i003Buildings 12 01954 i004Buildings 12 01954 i005
Buildings 12 01954 i006
Inlet-free area dimensionsWindow 1
(Existing)
Window 2Window 3
0.9417503780.58292142280.3355409635
Elevation4th floor10th floor15th floor
14.65 m33.55 m49.3 m
Table 4. Field measurement data used in the simulation.
Table 4. Field measurement data used in the simulation.
ElevationAir Velocity (m/s)
SouthwestSouthSoutheastEast
4th floor4.93.946.577.37
10th floor10.96.727.2612.18
15th floor7.595.927.257.94
Table 5. Regional simulation results.
Table 5. Regional simulation results.
ElevationAir Velocity (m/s)
EastWest
4th floorBuildings 12 01954 i007Buildings 12 01954 i008
10th floorBuildings 12 01954 i009Buildings 12 01954 i010
15th floorBuildings 12 01954 i011Buildings 12 01954 i012
Table 6. Results of the Torricelli and Bernoulli equations.
Table 6. Results of the Torricelli and Bernoulli equations.
FloorHeightTorricelli (cmHg)Bernoulli
414.6575.8535228.5173
1033.5575.6645440.038
1549.375.507578.1401
Table 7. Existing window: Indonesian National Standard Calculations.
Table 7. Existing window: Indonesian National Standard Calculations.
A (m2)FloorWind DirectionV (m/s)CvQStandard
0.9417503784thSouthwest4.90.62.768Min:
4.250 cfs
Max:
13.034 cfs
0.351.615
South3.940.62.226
0.351.298
Southeast6.570.63.712
0.352.165
East7.370.64.164
0.352.429
10thSouthwest10.90.66.159
0.353.592
South6.720.63.797
0.352.214
Southeast7.260.64.102
0.352.392
East12.180.66.882
0.354.014
15thSouthwest7.590.64.288
0.352.501
South5.920.63.345
0.351.951
Southeast7.250.64.096
0.352.389
East7.940.64.486
0.352.617
Table 8. The 100% air velocity values at an altitude of 400 m.
Table 8. The 100% air velocity values at an altitude of 400 m.
TowerFloorUnitVz (m/s)Vg (m/s)
A4th12.1854.994
22.1774.9
32.6636.087
43.7618.597
10th14.6988.636
24.1357.683
34.9019.107
44.999.272
15th14.9438.342
23.7636.350
34.8538.1905
45.5689.397
B4th12.1915.008
22.8676.553
32.134.868
43.5558.126
10th14.0587.540
25.80110.779
33.4456.401
45.38310.002
15th14.5517.6808
25.1438.68
34.4457.501
46.05510.219
Table 9. Reference air velocity at the measured elevations.
Table 9. Reference air velocity at the measured elevations.
TowerFloorUnitVref (m/s)Actual Air Velocity (m/s)
A4th12.1843.81
22.1433.86
32.6624.90
43.76086.57
10th14.6478.34
24.1346.76
34.900910.90
44.9897.26
15th14.9437.21
23.7627.08
34.8537.59
45.5677.25
B4th12.19084.31
22.8664.90
32.1293.94
43.5547.37
10th14.0576.72
25.800712.18
33.4445.12
45.3827.39
15th14.55095.92
25.14297.56
34.445.92
46.0547.94
Table 10. Reference air movement at the measured elevation.
Table 10. Reference air movement at the measured elevation.
TowerFloorUnitCv (m/s)
A4th10.41
20.44
30.44
40.45
10th10.46
20.54
30.53
40.49
15th10.4
20.46
30.29
40.46
B4th10.29
20.36
30.39
40.54
10th10.38
20.39
30.37
40.36
15th10.46
20.43
30.41
40.57
Table 11. CFD simulation results.
Table 11. CFD simulation results.
TowerFloorUnitSimulationResultStandard
A 4th1Buildings 12 01954 i013Buildings 12 01954 i014Window 1
Buildings 12 01954 i015
Theory:
0.2–0.6 m/s

Macfarlane:
0.41–0.45 m/s
SouthwestSouth
Buildings 12 01954 i016Buildings 12 01954 i017
SoutheastEast
2Buildings 12 01954 i018Buildings 12 01954 i019Window 1
Buildings 12 01954 i020
Theory:
0.2–0.6 m/s

Macfarlane:
0.41–0.45 m/s
SouthwestSouth
Buildings 12 01954 i021Buildings 12 01954 i022
SoutheastEast
3Buildings 12 01954 i023Buildings 12 01954 i024Window 2
Buildings 12 01954 i025
Theory:
0.2–0.6 m/s

Macfarlane:
0.41–0.45 m/s
SouthwestSouth
Buildings 12 01954 i026Buildings 12 01954 i027
SoutheastEast
4Buildings 12 01954 i028Buildings 12 01954 i029Window 2
Buildings 12 01954 i030
Theory:
0.2–0.6 m/s

Macfarlane:
0.41–0.45 m/s
SouthwestSouth
Buildings 12 01954 i031Buildings 12 01954 i032
SoutheastEast
10th1Buildings 12 01954 i033Buildings 12 01954 i034Window 2
Buildings 12 01954 i035
Theory:
0.2–0.6 m/s

Macfarlane:
0.46–0.54 m/s
SouthwestSouth
Buildings 12 01954 i036Buildings 12 01954 i037
SoutheastEast
2Buildings 12 01954 i038Buildings 12 01954 i039Window 2
Buildings 12 01954 i040
Theory:
0.2–0.6 m/s

Macfarlane:
0.46–0.54 m/s
SouthwestSouth
Buildings 12 01954 i041Buildings 12 01954 i042
SoutheastEast
3Buildings 12 01954 i043Buildings 12 01954 i044Window 3
Buildings 12 01954 i045
Theory:
0.2–0.6 m/s

Macfarlane:
0.46–0.54 m/s
SouthwestSouth
Buildings 12 01954 i046Buildings 12 01954 i047
SoutheastEast
4Buildings 12 01954 i048Buildings 12 01954 i049Window 3
Buildings 12 01954 i050
Theory:
0.2–0.6 m/s

Macfarlane:
0.46–0.54 m/s
SouthwestSouth
Buildings 12 01954 i051Buildings 12 01954 i052
SoutheastEast
15th1Buildings 12 01954 i053Buildings 12 01954 i054Window 1
Buildings 12 01954 i055
Theory:
0.2–0.6 m/s

Macfarlane:
0.29–0.46 m/s
SouthwestSouth
Buildings 12 01954 i056Buildings 12 01954 i057
SoutheastEast
2Buildings 12 01954 i058Buildings 12 01954 i059Window 3
Buildings 12 01954 i060
Theory:
0.2–0.6 m/s

Macfarlane:
0.29–0.46 m/s
SouthwestSouth
Buildings 12 01954 i061Buildings 12 01954 i062
SoutheastEast
3Buildings 12 01954 i063Buildings 12 01954 i064Window 2
Buildings 12 01954 i065
Theory:
0.2–0.6 m/s

Macfarlane:
0.29–0.46 m/s
SouthwestSouth
Buildings 12 01954 i066Buildings 12 01954 i067
SoutheastEast
4Buildings 12 01954 i068Buildings 12 01954 i069Window 2
Buildings 12 01954 i070
Theory:
0.2–0.6 m/s

Macfarlane:
0.29–0.46 m/s
SouthwestSouth
Buildings 12 01954 i071Buildings 12 01954 i072
SoutheastEast
B4th1Buildings 12 01954 i073Buildings 12 01954 i074Window 2
Buildings 12 01954 i075
Theory:
0.2–0.6 m/s

Macfarlane:
0.29–0.54 m/s
SouthwestSouth
Buildings 12 01954 i076Buildings 12 01954 i077
SoutheastEast
2Buildings 12 01954 i078Buildings 12 01954 i079Window 2
Buildings 12 01954 i080
Theory:
0.2–0.6 m/s

Macfarlane:
0.29–0.54 m/s
SouthwestSouth
Buildings 12 01954 i081Buildings 12 01954 i082
SoutheastEast
3Buildings 12 01954 i083Buildings 12 01954 i084Window 2
Buildings 12 01954 i085
Theory:
0.2–0.6 m/s

Macfarlane:
0.29–0.54 m/s
SouthwestSouth
Buildings 12 01954 i086Buildings 12 01954 i087
SoutheastEast
4Buildings 12 01954 i088Buildings 12 01954 i089Window 1
Buildings 12 01954 i090
Theory:
0.2–0.6 m/s

Macfarlane:
0.29–0.54 m/s
SouthwestSouth
Buildings 12 01954 i091Buildings 12 01954 i092
SoutheastEast
10th1Buildings 12 01954 i093Buildings 12 01954 i094Window 3
Buildings 12 01954 i095
Theory:
0.2–0.6 m/s

Macfarlane:
0.36–0.38 m/s
SouthwestSouth
Buildings 12 01954 i096Buildings 12 01954 i097
SoutheastEast
2Buildings 12 01954 i098Buildings 12 01954 i099Window 3
Buildings 12 01954 i100
Theory:
0.2–0.6 m/s

Macfarlane:
0.36–0.38 m/s
SouthwestSouth
Buildings 12 01954 i101Buildings 12 01954 i102
SoutheastEast
3Buildings 12 01954 i103Buildings 12 01954 i104Window 3
Buildings 12 01954 i105
Theory:
0.2–0.6 m/s

Macfarlane:
0.36–0.38 m/s
SouthwestSouth
Buildings 12 01954 i106Buildings 12 01954 i107
SoutheastEast
4Buildings 12 01954 i108Buildings 12 01954 i109Window 2
Buildings 12 01954 i110
Theory:
0.2–0.6 m/s

Macfarlane:
0.36–0.38 m/s
SouthwestSouth
Buildings 12 01954 i111Buildings 12 01954 i112
SoutheastEast
B15th1Buildings 12 01954 i113Buildings 12 01954 i114Window 2
Buildings 12 01954 i115
Theory:
0.2–0.6 m/s

Macfarlane:
0.41–0.57 m/s
SouthwestSouth
Buildings 12 01954 i116Buildings 12 01954 i117
SoutheastEast
2Buildings 12 01954 i118Buildings 12 01954 i119Window 2
Buildings 12 01954 i120
Theory:
0.2–0.6 m/s

Macfarlane:
0.41–0.57 m/s
SouthwestSouth
Buildings 12 01954 i121Buildings 12 01954 i122
SoutheastEast
3Buildings 12 01954 i123Buildings 12 01954 i124Window 2
Buildings 12 01954 i125
Theory:
0.2–0.6 m/s

Macfarlane:
0.41–0.57 m/s
SouthwestSouth
Buildings 12 01954 i126Buildings 12 01954 i127
SoutheastEast
4Buildings 12 01954 i128Buildings 12 01954 i129Window 1
Buildings 12 01954 i130
Theory:
0.2–0.6 m/s

Macfarlane:
0.41–0.57 m/s
SouthwestSouth
Buildings 12 01954 i131Buildings 12 01954 i132
SoutheastEast
Table 12. First-phase data interpolation results.
Table 12. First-phase data interpolation results.
TowerUnit 1Window 1Window 2Window 3Air Movement (m/s)
ASouthwestFloors 3–6 and 11–16Floors 7–10--0.46–0.6
SouthFloors 3–6 and 11–16Floors 7–10--0–0.3
SoutheastFloors 3–6 and 11–16Floors 7–10--0.2–0.3
EastFloors 3–6 and 11–16Floors 7–10--0.33–0.41
Unit 2Window 1Window 2Window 3Air Movement (m/s)
SouthwestFloors 3–6Floors 7–10Floors 11–160.4–0.5
SouthFloors 3–6 Floors 7–10Floors11–160.2–0.4
SoutheastFloors 3–6 Floors 7–10Floors 11–160–0.2
EastFloors 3–6 Floors 7–10Floors 11–160.2–0.4
Unit 3Window 1Window 2Window 3Air Movement (m/s)
Southwest--Floors 3–6 and 11–16Floors 7–100.2–0.5
South-- Floors 3–6 and 11–16Floors 7–100.29–0.37
Southeast-- Floors 3–6 and 11–16Floors 7–100
East-- Floors 3–6 and 11–16Floors 7–100.23–0.45
Unit 4Window 1Window 2Window 3Air Movement (m/s)
Southwest--Floors 3–6 and 11–16Floors 7–100.3–0.6
South-- Floors 3–6 and 11–16Floors 7–100.2–0.3
Southeast-- Floors 3–6 and 11–16Floors 7–100.2–0.6
East-- Floors 3–6 and 11–16Floors 7–100.2–0.54
TowerUnit 1Window 1Window 2Window 3Air Movement (m/s)
BSouthwest--Floors 3–6 and 11–16Floors 7–100.4–0.8
South-- Floors 3–6 and 11–16Floors 7–100.5–0.6
Southeast-- Floors 3–6 and 11–16Floors 7–100.2–0.3
East-- Floors 3–6 and 11–16Floors 7–100.25–0.48
Unit 2Window 1Window 2Window 3Air Movement (m/s)
Southwest--Floors 3–6 and 11–16Floors 7–100.3–0.5
South-- Floors 3–6 and 11–16Floors 7–100.3–0.5
Southeast-- Floors 3–6 and 11–16Floors 7–100.3–0.4
East-- Floors 3–6 and 11–16Floors 7–100.3–0.6
Unit 3Window 1Window 2Window 3Air Movement (m/s)
Southwest--Floors 3–6 and 11–16Floors 7–100.4–0.9
South-- Floors 3–6 and 11–16Floors 7–100.5–0.6
Southeast-- Floors 3–6 and 11–16Floors 7–100
East-- Floors 3–6 and 11–16Floors 7–100.28–0.44
Unit 4Window 1Window 2Window 3Air Movement (m/s)
SouthwestFloors 3–6 and 11–16Floors 7–10--0.4–0.5
SouthFloors 3–6 and 11–16Floors 7–10-- 0.3–0.5
SoutheastFloors 3–6 and 11–16Floors 7–10-- 0.2–0.33
EastFloors 3–6 and 11–16Floors 7–10-- 0.4–0.54
Table 13. Second-phase data interpolation: basic data.
Table 13. Second-phase data interpolation: basic data.
TowerFloorUnitWindowAir Movement (m/s)
SouthwestSouthSoutheastEast
A4th110.470.20.10.36
210.40.200.4
320.30.2900.38
420.30.20.260.29
10th120.600.20.41
220.50.30.20.3
330.50.3300.23
430.60.20.240.23
15th110.450.30.30.47
230.50.400.2
320.40.3700.45
420.50.30.60.52
B4th120.480.50.20.42
220.40.40.30.5
320.40.4900.44
410.40.30.190.42
10th130.80.60.20.29
230.30.30.30.3
330.90.6100.3
420.50.50.330.49
15th120.480.50.30.48
220.50.50.40.5
320.50.5800.44
410.50.50.310.53
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Napitupulu, S.S.; Hardiman, G.; Tobing, R.R. Determination of Window Dimensions Using Mathematical and Simulation Methods to Improve Window Performance: A Case Study on the Jatinegara Barat Flats, Jakarta. Buildings 2022, 12, 1954. https://doi.org/10.3390/buildings12111954

AMA Style

Napitupulu SS, Hardiman G, Tobing RR. Determination of Window Dimensions Using Mathematical and Simulation Methods to Improve Window Performance: A Case Study on the Jatinegara Barat Flats, Jakarta. Buildings. 2022; 12(11):1954. https://doi.org/10.3390/buildings12111954

Chicago/Turabian Style

Napitupulu, Sally Septania, Gagoek Hardiman, and Rumiati Rosaline Tobing. 2022. "Determination of Window Dimensions Using Mathematical and Simulation Methods to Improve Window Performance: A Case Study on the Jatinegara Barat Flats, Jakarta" Buildings 12, no. 11: 1954. https://doi.org/10.3390/buildings12111954

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop