Indoor exposure to outdoor pollutants adversely affects health, varying with building dimensions and particularly ventilation that have critical role on their indoor dispersion. This study assesses the impact of outdoor air on indoor air quality in a child care center. Computational fluid dynamics was utilized to analyze the dispersion of particulate matter, with a specific focus on window screens featuring 6 distinct pore sizes ranging from 0.8 mm to 2 mm and 2 different thicknesses of 0.5 mm and 0.1 mm. Results indicate that the presence of a window screen offers significant advantages in controlling particle infiltration compared to scenarios without a screen, as larger particles tend to pass directly through the window within the breathing zone. The scenario without window screens minimizes pressure drop but lacks enhanced particle capture capabilities. However, for effective particle reduction, the window screen with a pore size of 0.8 mm (R0.8T2) and a thickness of 0.5 mm proves to be the most beneficial, achieving the particle filtering efficiency of approximately 54.16%, while the larger window screen with a pore size of 2 mm and a thickness of 1 mm exhibits the lowest efficiency at about 23.85%. Nonetheless, screens with very small sizes are associated with a high-pressure drop, impacting energy efficiency, and overall window performance. Larger pores with smaller thicknesses (0.5 mm) reduced particle count by approximately 45.97%. Therefore, the significance of window screen thickness beyond pore size for particle reduction efficiency is highlighted, emphasizing screens’ role in indoor air quality and health protection.
Introduction
The deterioration of the atmospheric environment in urban areas is a prominent environmental issue in numerous countries intensified as a consequence of accelerated industrialization, urbanization, and the proliferation of motor vehicles, all driven by human activities.1 This has led to human-caused air pollution becoming a significant global public health threat, responsible for approximately 4.2 million premature deaths worldwide per year in 2019.2 Among air pollutants, particulate matter (PM), specifically PM2.5, poses a significant health concern due to its small size and ability to penetrate human lungs and bronchial passages.3,4
In developed countries, people are dedicating an increasing amount of their time, typically ranging from 80% to 90%, to various indoor environments.5,6 This highlights the significance of indoor air quality on public health, particularly for vulnerable groups, such as children and the elderly.7-9 Given that most people prefer to spend their time indoors, it is important to understand the PM levels in the indoor air.
Preschools and childcare centers serve as the primary indoor environments where young children spend a significant portion of the day. The activities undertaken by children, such as playing with toys on carpets and bedcovers, have a considerable impact on the concentration of air pollutants.10,11 It is critical to increase the awareness that children are physiologically more prone to the negative health repercussions of air pollution compared to adults. This heightened susceptibility is due to their higher volume of air intake relative to their body mass and the ongoing development of their lungs.12,13
Indoor air quality (IAQ) is primarily influenced by building attributes, such as ventilation, building location, cleaning practices, and the operability of windows. Achieving improved indoor air quality is accomplished by ventilating rooms with external air.14-16 There are 3 main ventilation methods: natural, mechanical, and hybrid (mixed mode).17 Among these, natural ventilation is the most efficient and environmentally friendly due to its reliance on the elements and the presence of large openings.18 However, controlling the natural ventilation can be challenging and requires careful planning.
Air filtration methods have demonstrated both high particle removal efficiency and relatively low air resistance. In previous studies, researchers have noted that window screens, acting as supportive substrates, can effectively block outdoor PM2.5 from entering indoor spaces when windows are open, especially during natural ventilation.19-22 Notably, the high-efficiency particulate air filter, renowned for its ability to remove 99.97% of particles measuring ⩾0.3 µm in diameter, stands out as a highly effective technology for enhancing IAQ.23 Despite high particle removal efficiency achieved with these filters, a significant increase in energy consumption in heating, ventilation, and air-conditioning (HVAC) systems can be brought about by the large air resistance.24
In previous studies, the exploration of window screen properties involved a combination of experimental investigations and modeling analyses. It was reported that nanofiber filters were recommended as a viable option for use as window screens to prevent outdoor PM2.5 from entering indoor spaces through natural ventilation during the open window.24-26 An aspect to consider is that nanofiber filters may be less reusable due to the challenge of cleaning PM particles deeply embedded within the nanofiber structure.27Another type of air filtration system is the porous membrane filter, typically constructed by creating minute pores within a solid substrate or a thin film, meticulously designed to target larger PM particles.28 In addition to pore size, the filtration process benefits from various essential mechanisms, including inertial impaction, interception, Brownian motion, gravitational settling, and electrostatic forces.29 Previous research has demonstrated the impressive performance of isoporous through-hole membrane filters, with exceptional filtration efficiency maintained while achieving a low pressure drop.27 Furthermore, the ease of maintenance is provided by these filters, allowing for convenient cleaning through simple wipe with a damp cloth or a gentle spray with water. Additionally, they are designed for multiple uses, making them a cost-effective and eco-friendly choice.27
Computational fluid dynamics (CFD) is a well-established tool for assessing and understanding ventilation and pollutant dispersion patterns across various scales in the urban environment.30 CFD simulations have been extensively applied in various fields, including the food industry, biotechnology processes, agro-environmental applications, and the design of indoor and outdoor environments.31-35 The CFD models are widely used for analysis of airflow, temperature, and contaminant distributions in a diverse range of applications, providing detailed thermal and contaminant data through simulation, particularly emphasizing natural ventilation scenarios.36-38
In CFD models, particle dispersion can be approached via 2 methods: Eulerian and Lagrangian. Eulerian methods treat particles as a continuous medium, deriving particle concentration fields from conservation equations to understand dispersion broadly. In contrast, the Lagrangian discrete phase model meticulously tracks individual particle dynamics and is particularly effective for modeling dust dispersion.39,40 Therefore, in this study, the discrete phase model was selected to conduct the investigation of particle dispersion.
The primary objective of this study is to perform a numerical analysis exploring the impact of outdoor air on indoor air characteristics, examining the dispersion of particulate matter. Additionally, the study aims to assess the influence of window screens with different pore sizes and thicknesses on the air characteristics in a childcare center using CFD. The finding can assist childcare centers located near small scale industrial areas in formulating plans to implement suitable measures for controlling indoor air pollutant exposure, with a particular focus on the use of window screens.
Materials and Methods
Study area
The research site was located in Samut Sakhon province, approximately 48 km away from Bangkok, as illustrated in Figure 1a. This area is known for its substantial metal smelting industry, with the childcare center (Figure 1b) situated about 1 km from the small scale industrial area.
Data collection
The meteorological data, including various wind speeds used in this study, were collected from the King Mongkut’s University of Technology Thonburi meteorological station during the cool and hot dry seasons.41 During cool dry periods: November 17th, 2022. During this period, there was about 20-40 mm of rain with a relative humidity of 70%-75% and an average temperature of 29.0°C ± 5.8°C. In contrast, the hot dry season in February experienced 10-20 mm of rain, lower relative humidity (65%-70%), and a higher average temperature of 30.2°C ± 5.9°C.42 These measurements were recorded over an 8-hour (07.00 AM-03.00 PM) duration during each period.
Case description
The simulation was conducted using Autodesk inventor 2020 and ANSYS 2021 R2. The study primarily examined the horizontal sliding window type, a prevalent configuration in this study area where windows are commonly fully open. During this period, 2 scenarios were investigated: one with the window fully open without a window screen, and the other with the window fully open while using a window screen.
Geometry design
In CFD simulations, the generation of a highly detailed geometric model of the simulated area is essential before the simulation mesh is created. This study, specifically, focused on a selected room within the childcare center, as depicted in Figure 2a. The arrangement of window screen placement is illustrated in Figure 2b, presenting a side view of the window where the screen will be installed on the horizontal sliding window type. The window screen examined in this study closely mirrors an insect screen. The configuration of the air window screen is established through the porous jump method.43
In this model, the thin perforated plate is represented as a porous medium with a limited thickness (equation (1)). The decrease in pressure arises from a blend of Darcy’s law and an extra inertial loss factor.44
where, µ and ρ represent viscosity and density of air, C2 stands for the pressure-jump coefficient, υ denotes the velocity perpendicular to the porous plate, and ∆m represents the plate’s thickness.
Meshing strategies
The constructed room’s geometry was exported to ANSYS, where the model underwent numerical discretization through mesh generation. The computational domain had the dimensions of the classroom, serving as the primary focus of the study (3.75 m × 2.68 m × 2.96 m). Furthermore, the assessment of the mesh geometric quality in this study was achieved by examining parameters related to orthogonal quality and skewness (Figure 3) with mesh details presented in Table 1. The computation grid resolution, with varying element sizes set to 0.03 m for the filter, 0.065 m for the environment room, and 0.06 m for the window zone, resulted in mesh quality showing acceptable skewness (0.25-0.50) and very good orthogonality (0.70-0.95). In this study, the grid is first constructed in a horizontal plane before being extended vertically. This method provides a high level of control over the shapes and sizes of the individual cells.45 The structured grid was generated in ANSYS Mesh using the sweep method, along with edge sizing applied to parallel corresponding edges. These parameters were rated on a scale ranging from 0 to 1, where values higher values closer to 1 indicate better orthogonal quality and lower values closer to 0 suggest improved skewness.46
Table 1.
The details of the mesh structure.
Boundary condition
The accuracy of CFD simulations is dependent on the correct configuration of boundary conditions and numerical simulation parameters.47 Figure 4 depicts the boundary conditions, consisting of a mass flow inlet (Figure 4a) where the primary wind enters perpendicular to the inflow area, with its direction and velocity determined by the inflow area’s boundary conditions. The outlet and wall conditions are illustrated in Figure 4b. For the outflow area, boundary conditions must align with those of the inflow area to ensure conservation of flux within the computational domain. The walls are characterized by the fluid velocity profile near the contact area, determined by y+, a dimensionless value representing the distance between the fluid and the wall and assessing turbulent effects according to the law of walls.48,49
In Table 2, the numerical parameters utilized in this study were drawn from field data and existing literature,50-52 with the assumption that the wind velocity and particle velocity were equal. The wind velocity was determined based on the average wind speed over an 8-hour period. In this study, 3 different groups were formed, consisting of round perforated sheets: R0.8T2, R1T2, and R2T3, with each assigned a thickness of 0.5 mm and 1 mm, as shown in Table 3. The open area was reported as a percentage that reflects the extent to which the sheet is occupied by holes. Furthermore, the range of particle diameters employed in this study was set as 1 to 100 µm.
Table 2.
Summary of numerical parameters used in this study.50-52
Table 3.
Window screen properties.
Numerical setup
The numerical setup assumed the continuous phase characterized as air with its flow simulated based on the realizable k − ε model, while the dust particles were regarded as the discrete phase, and the trajectory of each particle was computed utilizing the discrete phase model (DPM).
Governing equation
The movement of air is governed by the Navier-Stokes equation, which are given in equation (2) to equation (5)53,54:
where, p, ρg, ug are static pressure, gas density, gas volume fraction, gas velocity, and gravity acceleration, while is the viscous stress tensor, respectively. µ, µt, k, I, and F are the molecular viscosity, turbulent viscosity, turbulent kinetic energy unit tensor, and volumetric momentum transfer rate between particles and gas phases. Furthermore, in this study, the standard k − ε model is widely employed as the predominant turbulence model (equations (6)-(9)).
where ε, Gk are turbulent dissipation rate and turbulent production with the following values for the constant: Cµ, Clε, C2ε, σk, σε are set to 0.09, 1.44, 19.92, 1, and 1.3.
Discrete Phase Model (DPM)
The Euler-Lagrange model approach within the DPM is utilized for tracking particle movements. The continuous phase is coupled with the k-ε turbulence model to facilitate calculations alongside the discrete phase.55 Based on the DPM each particle was traced during the whole dispersion process.
In accordance with Newton’s second law of motion, the movement of particles within the airflow field can be described as in equation (10):
where, mp represents the particle mass; up signifies particle velocity; Fg is the gravitational force, while FD represents drag force as in equation (11):
where, is particle velocity and dp signifies particle diameter. Additionally, CD represents the drag coefficient that can be considered as in equation (12):
where, a1, a2 and a3 represent constants relevant to smooth spherical particles within multiple Reynolds number (Re) ranges while µ stands for dynamic viscous as in equation (13).51,52
Mass imbalance
Due to the limited availability of experimental data, the mass imbalance method was employed in this study. In a previous study, the mass imbalance method was considered an indicator of result accuracy.56-59 The mass imbalance in this study serves as a measure of the simulation error, determined by comparing the mass flow rates entering and leaving the computational domain. Subsequently, the results were examined and compared based on the mass imbalance factor at the conclusion of the modeling practice. A recommended practice is to target solution imbalances of less than 1%, as indicated by equation (14).60
Window screen efficiency
The window screen efficiency is defined in equation (15)61
where, η represents the window screen efficiency and Nin and Nout represent the number of particles at the inlet and outlet.
Result and Discussion
Comparison of pressure distribution between without and with a window screen
In this study, a pressure outlet boundary condition was specified at the air outlet. Within this configuration, static pressure is set as a constant at this boundary, and all other flow properties are derived from the internal flow. The careful management of pressure is considered important for regulating the transfer of airborne contaminants between the indoor and outdoor spaces of the building.62 The pressure distribution during the full open window condition without the presence of a window screen is depicted in Figure 5. Notably, high pressure is observed at the inlet, where it is recorded at 7.184 Pa with a pressure drop to 0.176 Pa (Figure 5a). However, during the hot dry season (Figure 5b), the pressure value significantly, reaching 20.287 Pa. A more substantial pressure drop of 0.481 Pa is also noted during this season. The fluctuations in indoor pressure, after air passes through natural ventilation in this study, are linked to the principles of the Bernoulli equation which suggests that a pressure disparity emerges when the airflow channel undergoes narrowing or experiences alterations in its geometry.63 Furthermore, Window Screen E was selected to show an overview of pressure distribution (Figure 5c) during cool dry season. Its inclusion offers a visual representation of the associated pressure values, allowing for a comparison of different pressure values with and without the window screen. In the realm of natural ventilation, these pressure variations play a pivotal role in shaping the broader dynamics of airflow within the indoor environment.64 Significantly, this investigation reveals a dynamic in which the indoor air pressure is reduced, resulting in the creation of a negative pressure differential compared to the outdoor air pressure.65 This imbalance, as unveiled in this study, leads to the infiltration of outdoor air into the building through the window. It should be noted that the high pressure drop across a window screen is directly related to 2 key factors: velocity profile and filter thickness.25 With an increased pressure drop, greater removal efficiency is achieved,66 reinforcing the significance of the filter’s role in maintaining indoor air quality.
Comparison of velocity distribution between without and with a window screen
The velocity distribution across the domain is presented in Figure 6. A peak velocity of around 3.219 m/s was documented (Figure 6a), signifying the airflow speed in the absence of a window screen. In contrast, a notably higher maximum velocity of around 5.395 m/s was observed during the hot dry season (Figure 6b). This significant difference between the 2 seasons suggests a pronounced variation in airflow characteristics, with swifter air movement featured during the hot dry season, surpassing that of the cool dry season. The introduction of Window Screen E during cool dry season resulted in a higher maximum velocity of 6.092 m/s (Figure 6c). Notably, both cases experienced a considerable reduction in velocity, potentially attributed to friction between air molecules and the complete opening of a horizontally sliding window. This particular window type, when fully opened, can impede the airflow, leading to a reduction in velocity. As demonstrated by Heiselberg et al.,67 a direct correlation is found between the angle at which a window is opened and the ensuing impact on airflow and velocity levels. The more open the window, the greater the airflow and velocity it accommodates.
In Table 4, the summary of pressure and velocity results in an open window condition without a window screen is presented. These findings are in agreement with the results observed in the previous study, where a decrease in pressure and an increase in velocity were observed.68
Table 4.
Summary results of inlet and outlet in open window condition without window screen.
Mass imbalance without a window screen
The summary of the mass imbalance of a window without a screen is presented in Table 5. The numerical model results consistently indicate the presence of mass imbalance, a critical aspect to consider. It is noteworthy mentioning that previous research has emphasized the importance of achieving the lowest possible mass imbalance59 as it often leads to the most favorable outcomes in numerical simulations. Throughout this study, it was determined that the percentage of mass imbalance is 0.05%, which confirms reliability of the numerical model applied in this study.
Table 5.
The summary details of mass imbalance in an open window condition without window screen.
The result of inlet 10.76 kg/s represents the mass flow rate of the air. The cross-sectional area equals the surface area of the classroom, and in this study, the surface area plays a crucial role in determining whether a given mass flow rate is appropriate. A larger surface area can disperse and accommodate a larger mass flow rate more effectively.
Comparison of particle dispersion between without and with a window screen
Figure 7 illustrates the dispersion of outdoor aerosol particles indoors, influenced by natural ventilation. The particle range considered in this study ranges from 1 µm to 100 µm, assumed to be transported indoors via a widely opened horizontally sliding window. The results show that without a window screen (Figure 7a), the particles pass directly through the window within the height of the breathing zone with some of larger particles found to be deposited quickly on the floor. Previous research has reported a symmetrical distribution of velocity profiles and pollutant concentrations, particularly for particles equal to or smaller than 0.8 µm, which tend to exhibit higher concentrations under lower wind speeds, primarily due to the accumulation process.69,70 Figure 7b illustrates the dispersion of particles with the addition of Window Screen F. Despite the largest pore size among the window screens considered in the study, an elevated dispersion profile emerges, indicating fewer particles passing through when using the window screens, suggesting on the reduced distribution of the outdoor air particles in the indoor environment.
Window screen efficiency
The effectiveness of the 6 different types of window screens on limiting the entry of particles into a childcare center was assessed and presented in Table 6. The pressure jump coefficient (C2) is a parameter that accounts for the depth of the porous medium to estimate the pressure drop across the porous jump zone, allowing for the determination of the pressure drop across the study area. The investigation revealed a significant range in window efficiency, influencing the potential particle entry into the childcare facility, with efficiencies ranging from 23.85% to 54.16%, respectively. Furthermore, the mass imbalances of the window screen simulations were recorded at less than 0.1% across the 6 window screen types. The study revealed that the smaller the pore size in the window screen the higher the efficiency in preventing outdoor atmospheric particles from entering indoors, with window screen A, featuring a 0.8 mm pore size and 0.5 mm thickness, emerging as the most efficient at 54.16% with the highest pressure drop. It has been reported that there is a relationship between pressure drop and particle removal.71 When there is a higher pressure drop in a filtration system, it becomes more difficult for air to move through the filter. Consequently, this resistance prevents airborne particles from entering easily, acting as a protective measure to prevent the infiltration of contaminants. Conversely, the Window Screen F with the largest pore size (2 mm) and thickness (1 mm) displayed the lowest efficiency at 23.85%. The study also showed that screen thickness of 0.5 mm has higher efficiencies than screens with larger thickness of 1 mm. This might be due to the higher pressure across the window screen, indicating improved airflow dynamics and enhanced performance. This study aligns with previous findings that indicate higher capture efficiency in screens with lower thickness.72 Additionally, it is important to note that the high filtration efficiency, despite the larger pore size compared to the particle diameters, can be attributed to various filtration mechanisms such as interception, diffusion, electrostatic attraction, inertial impaction, and Brownian motion.73
Table 6.
Summary details of mass imbalance, pressure, velocity, and window screen efficiency in filtering outdoor particles from entering indoors.
Conclusion
Indoor air quality holds significant importance, particularly as people spend more time indoors, and it relies heavily on building attributes, such as ventilation. Natural ventilation can inadvertently introduce contaminants from outdoors to indoors, potentially causing adverse health effects on individuals. In this study, a numerical analysis was conducted to explore the impact of outdoor air on indoor air quality within a childcare center located near a small-scale industrial area. The findings revealed that the presence of a window screen provides notable advantages in controlling particle infiltration when compared to scenarios without a screen. Specifically, larger particles tend to pass directly through the window within the height of the breathing zone. On the other hand, even with a screen featuring larger pores, there is a discernible shift in dispersion patterns, indicating a reduction in particle penetration through the window. In a comparative study of window screens with diameters ranging from 0.8 mm to 2 mm and varying thickness, this research examined pressure and velocity distributions, revealing distinctions in maximum velocities and pressures. The introduction of window screens resulted in a notable reduction in velocity, with the smallest screen size (R0.8T2) exhibiting the highest velocities. Pressure variations were associated with screen size effects. Nevertheless, introducing a smaller window screen in the childcare center led to a significant increase in pressure drop compared to configurations without any window screen and with a larger pore size window screen. Prioritizing particle reduction, it was found that the window screen, R0.8T2 with 0.5 mm thickness had the highest efficiency (54.16%) in particle removal, while the larger pore size (R2T3) with a larger thickness of 1 mm showed the lowest efficiency (23.85%). However, beyond pore size, the study emphasizes the critical role played by thickness in determining window screen efficiency. This study recommends considering window screens for particle reduction in the childcare center, for the protection of human health.
Author Contribution Statement
Cathleen Ariella Simatupang: Conceptualization, Methodology, Formal analysis, Investigation, Writing - Original Draft. Vladimir Strezov: Conceptualization, Writing – Review & Editing, Supervision. Suwanna Kitpati Boontanon: Conceptualization, Writing – Review & Editing, Funding acquisition, Visualization, Supervision. Prapat Pongkiatkul: Conceptualization, Methodology, Investigation, Formal analysis, Writing – Review & Editing. Narin Boontanon: Writing – Review & Editing, Visualization, Supervision. Ranjna Jindal: Writing – Review & Editing, Visualization, Supervision.
Ethics Approval and Consent to Participate
Not applicable.
Consent for Publication
All authors approved the manuscript for publication.
Availability of Data and Material
The authors confirm that the data supporting the findings of this study are available within the article. Data will be made available on request.
© The Author(s) 2024
This article is distributed under the terms of the Creative Commons Attribution-NonCommercial 4.0 License ( https://creativecommons.org/licenses/by-nc/4.0/) which permits non-commercial use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages ( https://us.sagepub.com/en-us/nam/open-access-at-sage).