Abstract

Previous studies showed that opening windows could help with kitchen ventilation in pollutant removal. However, no studies have systematically examined the impacts of window positions on kitchen hood performance, and there is insufficient information on indoor airflow characteristics and pollutant distribution from makeup air through open windows. Therefore, the objective of this study was to use a validated computational fluid dynamics approach with CO2 as an indoor air quality indicator (a surrogate for cooking emissions) to understand the impacts of exhaust flow rate and the window opening position on the flow characteristics, concentration distribution, and capture efficiency (CE) of the hood. We conducted four-point validation tests of the numerical models based on CO2 concentration and temperature measurements under steady-state conditions. The validated models were subsequently used in simulations to understand the effects of six different window opening positions and the two exhaust flow rates on exposure. We found that the CO2 concentration could be better reduced by having windows open at the higher location. Generally, the front windows were more effective with , followed by the back and the side windows, respectively. We also found that as the exhaust flow rate increased from 6.72 to 12.16 m3/min, CE reached >75% for all window positions, where the most significant increase was 1.58 times for the lower side window. To sum up, changing the relative position of the window and the exhaust hood could help disperse the incoming airflow from the window, improve the kitchen’s overall ventilation, and reduce pollutant concentration.

1. Introduction

Establishing a healthy indoor environment is essential because people spend most of their time indoors [1]. However, pollution from both indoor and outdoor sources has been linked to worsening indoor air quality (IAQ) and could pose various threats to human health [2]. Cooking activities were found to be a significant source of indoor pollution and could generate harmful substances such as particulate matter (PM), polycyclic aromatic hydrocarbons (PAHs), volatile organic compounds (VOCs), and other gaseous pollutants such as nitrogen dioxide (NO2) and carbon dioxide (CO2) [35]. For example, Ko et al. investigated the contributing factors of lung cancer in nonsmoking women in Taiwan and found that the risk was associated with an increased number of meals cooked per day [6]. Belanger et al. discovered that asthmatic children living in homes using gas stoves were more likely to have onsets of respiratory symptoms, and the association increased with rising indoor NO2 concentration [7]. In another study by Hu et al., a decline in indoor NO2 emissions was linked to a decreased incidence rate of pediatric asthma [8]. To reduce the concentration of cooking emissions and protect residents’ health, the use of a kitchen exhaust hood could be an effective intervention strategy [9, 10]. Kitchen hood performance has been characterized in experimental studies. Factors such as hood flow rates, burner positions, and hood geometry designs were found to have considerable impacts on the hood performance [1116]. However, these studies relied on measurements of cooking emissions at specific locations, and information on the indoor flow characteristics and how they influence the pollutant distribution in the kitchen was limited.

The computational fluid dynamics (CFD) approach has been widely used in IAQ studies to supplement findings from experimental studies through modeling and simulation. Compared to full-scale experiments, a validated CFD approach could be more cost-effective and provide additional information, including continuous spatial distributions of the flow field, temperature, and target air pollutants through visualization. This information could further assist in understanding the exposure pathway of air pollutants, which provides insight into improving hood performance and reducing human exposure. Several studies have used the CFD approach to simulate the flow characteristics in the kitchen environment under various hood and environmental settings. For example, Zhou and Kim simulated the spatial distribution of CO2 concentration and temperature in a residential open kitchen during cooking periods [17]. The results illustrated a better removal effect on CO2 concentration and heat when the angle between the air inlet and ceiling decreased and under the higher flow rate (16.67 m3/min). Zhou et al. investigated a new ventilation control system with the air curtain inlets surrounding the gas stove and the exhaust hood installed above in a residential kitchen [18]. The study demonstrated that the cooking emissions and the thermal environment could be well controlled and improved with the air curtain system. Chen et al. analyzed the distribution of PM-based cooking emissions and the thermal condition under different flow rates of the exhaust hood by using a validated kitchen environment simulation model. The optimal condition was achieved as the exhaust flow rate increased above 11 m3/min [19]. Le Hocine et al. assessed the impacts of the number of in-use burners and flow rate by modeling the exhaust hood performance in a residential kitchen [20]. They used CO2 as the tracer gas and found that the capture efficiency (CE) was 100% when the flow rate reached 8.5 m3/min. Yi et al. employed a concurrent supply and exhaust ventilation system in the kitchen and effectively improved the hood performance with the optimum flow rate of 1.6 m3/min [21]. In essence, most studies indicated that an increase in the hood flow rate or using an additional supply air system helped with air pollutants and heat removal. However, the former approach could also lead to excessive energy consumption and noise [22]. Therefore, setting the flow rate at a range that ensures effective pollutant removal while minimizing the undesirable effects is crucial.

The use of a kitchen hood could influence indoor ventilation in two ways: by providing local air exhaust in the proximity of the cooking area and general ventilation in the household with the introduction of ambient makeup air. For Chinese-style residential kitchens, ventilation using natural makeup air from window openings or cracks consisted of prevalent scenarios [23, 24]. However, very few studies have evaluated the influence of ventilation conditions on hood performance. Sun and Wallace investigated 132 Canadian homes and pointed out that exposure to cooking emissions could be well reduced with the window openings in the kitchen [25]. Wang et al. modeled the flow characteristics in the residential kitchen under nine scenarios with a combination of an open window and/or door under different exhaust flow rates [26]. The results showed that only opening the window would considerably improve IAQ around the cooking zone under the same flow rate, and the window location could be an influential factor. He et al. studied the effects of natural makeup air on hood performance by measuring time-elapsed particle concentrations during oil heating experiments in a simulated Chinese residential kitchen [27]. They set up five different openings (three window openings on the same side, a ceiling-mounted opening, and a floor-mounted opening) and found substantial leakage when opening the window closest to the stove. The authors concluded that the leakage could be alleviated by adjusting the window position or supplying makeup air through the mounted openings. In their study, the effects of the window’s spatial orientation and CE quantification remained unidentified.

Although opening windows could help with kitchen ventilation in pollutant removal, no studies have systematically examined the impacts of window positions on hood performance. Furthermore, there is insufficient information on indoor airflow characteristics and pollutant distribution from makeup air through open windows. Therefore, the objective of this study was to use a validated CFD approach with CO2 as an IAQ indicator (a surrogate for cooking emissions) to understand the impacts of exhaust flow rate and the window opening position on the flow characteristics, concentration distribution, and CE of the hood. A total of 12 scenarios (six window positions and two flow rates) were evaluated.

2. Materials and Methods

2.1. Study Zone

The experimental study was conducted in the laboratory where the dimensions of the study zone were 3.9 m in length, 3.1 m in width, and 3.6 m in height. Figure 1 shows the reconstructed schematic diagram of the study zone for the CFD simulation. The minor physical structures, such as pipes and cables, were simplified in the three-dimensional (3D) model for a more efficient computing process. A cooking station with an artificial pollutant generation system and a kitchen exhaust hood was set up to simulate hood operation during cooking activities. We used a heating chamber to provide a buoyant plume and a round opening on the top to approximate the cooking condition with a 28 cm diameter pan heated on an electric stove. Constant CO2 emission from the gas cylinder was directed into the heating chamber as a surrogate of gaseous pollutants through a duct equipped with a small mixing fan (). The chamber also contained particle emissions from the nebulizers at a constant rate; however, only CO2 measurements are presented in this study. An exhaust hood with two vertical side panels was installed 70 cm away from the opening of the heating chamber and 1.65 m from the ground level. The side panels were expected to ensure high removal efficiency for the cooking emissions, consistent with a similar design suggested by Zhao et al. The authors employed a CFD approach to evaluate the hood performance among 21 residential exhaust hoods with different geometric designs and found that the capture and containment efficiency could be enhanced by 20% with the installation of the side panels [28]. The dimensions of the hood were 100 cm in length, 50 cm in width, and 67 cm in height, with other detailed information described in the Supplemental Information (SI) (SI-1, Figure S1). Additionally, a 1.6 m tall human model was placed at a distance of 25 cm from the heating chamber to mimic the standing position of a person while cooking. The breathing zone was defined as a 30 cm radius from the nose of the human model, which was at the height of 1.5 m. Finally, there were two openings () in the wall and the ceiling, respectively. Under the operation of the exhaust hood, the fresh, ambient air was passively introduced into the study zone through these openings.

2.2. Experimental Procedure

Experiments were carried out to validate the predicted numerical results based on the steady-state conditions of CO2 concentration and temperature to ensure the reliability of the simulation in triplicate. In each experiment, we first started the constant CO2 emission at the flow rate of 0.00036 m3/min to introduce the pollutant into the heating chamber. The temperature of the coil in the heating chamber was set at 350°C to provide a buoyant plume to simulate that from cooking activities. We maintained the air exchange rate using the CO2 constant injection method at 2 per hour by manually adjusting the openings of the study zone [29]. Detailed measurement procedure of air exchange rate can be found in SI (SI-2). Based on the mass balance model, the concentration of the artificial pollutant followed an exponential increase and eventually reached the steady state, where the CO2 concentration and the temperature in the breathing zone of the human model were maintained at 3000 ppm and 32°C, respectively. We subsequently turned on the exhaust hood flow rate from low (6.72 m3/min) to high (12.16 m3/min) and waited for the CO2 concentration to reach the new steady states under the hood operation. We then measured the CO2 concentrations and the temperature concurrently at four vertically aligned locations ( m, 1.6 m, 2.4 m, and 3.2 m) on the right side of the human model (Figure 1), in conjunction with one measurement point outside of the study zone to monitor the ambient environmental condition. All instruments were calibrated, and their detailed information is listed in SI (SI-3, Table S1). The three repeated measurements from the four vertical locations were then used for the validation tests with the simulation results.

2.3. Numerical Method

The 3D computational model in the study was implemented in SolidWorks Flow Simulation [30]. The steady-state flow field in SolidWorks was solved based on the Favre-averaged Navier-Stokes equations, together with the modified turbulence model to close the equations. The Favre-averaged Navier-Stokes equations consist of the continuity, momentum, and energy equations, respectively (Equations (1)–(3)). They are defined as follows: where is the air density; and are the fluid velocity in -, -, and -directions (-, -, -directions); is the time; is the pressure; is the viscous shear stress tensor (Equation (5)); is the Reynolds stress tensor (Equation (6)); is the external force; is the diffusive heat flux; is the turbulence dissipation rate; is the internal heat generation; is the thermal enthalpy; is the viscosity; is the Kronecker delta function; represents the turbulent kinetic energy; and is the eddy viscosity. Additionally, the turbulence model proposed by Lam and Bremhorst was employed in the calculation [31]. Finally, to solve the complex algorithm, SolidWorks Flow Simulation used a pressure-based solver and SIMPLE algorithm to decouple the velocity and pressure to simulate the flow field [32].

2.3.1. Boundary Conditions

The boundary conditions were primarily based on the settings in the validation experiments and the ambient environmental conditions (Table 1). The wall materials were defined to account for the heat conduction in the laboratory environment. The initial condition inside the study zone for the CFD simulation was set to be 3000 ppm for the CO2 concentration and 32°C for the temperature, representing the steady-state condition when the exhaust hood was turned off. The study zone was divided by the cube shape and Cartesian-based grids for the grid generation, where the local grids were refined around the window, the hood, and the heating chamber regions. The detailed grid size is defined in Section 2.3.2.

2.3.2. Model Validation

We conducted the grid independence test before the CFD simulation results were valid for analysis. Four mesh cases were tested per number of total grids. They were named mesh cases 1, 2, 3, and 4, comprising , , , and grids, respectively. The simulation results of the optimal grid condition were subsequently compared with the averaged experimental results at the four measurement locations shown in Figure 1. The validation was considered acceptable if the calculated difference ratio (Equation (7)) between the simulation and experimental values was below 10% at each measurement location.

2.4. Simulation Cases

We modified the validated baseline CFD model to simulate a more representative kitchen layout. The cooking station was repositioned against the wall, and the ceiling opening was removed, as shown in Figure 2. Additionally, we increased the size of the window to in dimension. Six different window opening positions were selected, with two each on the front, side, and back orientation of the kitchen walls. The front side represents the wall where the human model was facing, whereas the back refers to the wall behind the human model. The side windows were installed in the wall on the right-hand side of the human model. Furthermore, considering the location of the transom openings commonly found in Taiwanese residences, we set the window positions at two different heights (1.2 and 2.4 m) at the same horizontal position of each wall. Information on the six window opening positions is summarized for the simulation cases in Table 2.

We presented the simulation results of the CO2 concentration and the flow field across the breathing plane defined as the vertical plane cut from the center of the measurement point in the breathing zone (Figure 3). We also determined the CO2-based CE using Equation (8) developed by Li and Delsante to characterize the hood performance for each simulation case under the steady-state conditions [33].

3. Results

3.1. Grid Independence Test and Model Validation

We compared the simulation results of the CO2 concentration and temperature by the height of the measurement locations for the four mesh cases (Figure 4). As the grid number increased, the simulation values showed a similar trend in mesh cases 3 and 4. In Figure 4, the values at the height of 2.5-3 m were slightly different between mesh cases 3 and 4; however, the values at the other sampling points were in good agreement. Overall, the grid test was considered to have reached independence. Using the experimental results as the baseline values, the difference ratios for mesh cases 3 and 4 were under 10% at all measurement locations and were considered acceptable for validation (Table 3). Mesh case 3 was chosen for the follow-up simulation cases (Section 3.2) because it had better computational efficiency due to the lower number of grids than mesh case 4. As a result, mesh case 3 could reasonably reflect the spatial distribution of CO2 concentration and temperature in the kitchen environment under hood operation.

3.2. Simulation Cases
3.2.1. Distribution of the Flow Field by Exhaust Flow Rate

We used mesh case 3 to simulate the cases for different window opening locations, where the maximum and minimum grid sizes were 60 mm and 3.75 mm, respectively (Figure 5). Figure 6 shows the distribution of the CO2 concentration at the steady state under the low (6.72 m3/min) and high (12.16 m3/min) flow rates for different window positions. When the exhaust hood flow rate was low, the hood demonstrated variable performance in removing the emitted CO2, which was leaked from the proximity of the exhaust hood in all six cases. The leaked CO2 accumulated in the upper space of the kitchen, following the vertical temperature gradient established by the buoyancy effect of the heat released from the simulated cooking activity (SI-4, Figure S2). When the exhaust flow rate increased to 12.16 m3/min, the hood showed improved performance in removing the pollutants, leading to a decreased level of escaped CO2. However, an accumulation of CO2 was observed in cases G, H, and L, suggesting that window positions could still be an influential factor in pollutant distribution under the high exhaust flow rate.

3.2.2. Distribution of the Flow Field by Window Height

For different window installation heights, CO2 concentration gradients with varying degrees were observed for cases with a lower window position (cases A, B, and C). The concentrations were higher in the space above the exhaust opening of the hood, compared to those in cases D, E, and F (Figure 6(a)). To assess the extent of the concentration gradient, we compared the maximum and minimum CO2 values from the plane of the measurement points used in the previous validation tests. The concentration difference ranged from 310 to 1080 ppm for the cases with a lower window position. In contrast, it decreased by approximately 50% for the cases with a higher window position and ranged from 160 to 510 ppm, suggesting better air mixing. Figure 7(a) shows the CO2 concentration and the airflow field under the low flow rate in the plane from the center of the window and the plane from the center of the heating chamber. As mentioned previously, there was a temperature gradient with higher CO2 concentration near the ceiling due to the buoyancy effects from simulated cooking activities. When the cooler ambient air entered from the upper window, it went downwards and carried part of the CO2 from the higher space to the bottom of the kitchen, which facilitated air mixing and the dilution of CO2 concentration near the ceiling. On the other hand, when the window was placed at lower positions, the cooler ambient air had the tendency to stay at the bottom of the space. The results showed that exhaust hood performance was better in cases with windows located at the height of 2.4 m than 1.2 m.

3.2.3. Distribution of the Flow Field by Window-Wall Locations

In addition to the window height, the window-wall locations (front, side, and back) were found to influence pollutant removal by the kitchen hood, and the effect was more pronounced under the low flow rate (Figure 7). When the window was in the front wall (cases A and D), the CO2 concentrations in the kitchen were lower than those in the back or side walls. We further used cases D, E, and F (the low flow rate and the high window position) as examples to discuss how the selection of the window-wall location affected the flow field inside the kitchen (Figure 8). For case D, as the fresh air entered from the upper front window, a counterclockwise vortex was generated in the kitchen. The vortex contributed to the pollutant removal in two aspects: (1) it enhanced the air mixing and (2) it carried the pollutant along the way to the hood at the end of the vortex and facilitated the pollutant exhaust. On the other hand, the flow field of case E showed that the fresh makeup air entered from the back window and flowed clockwise inside the kitchen, creating multiple vortexes as it reached the cooking station. The vortexes interfered with the pollutant capture at the hood, which subsequently led to the accumulation of CO2 concentration in the kitchen. As for case F, the makeup air moved towards and reached the front wall. As the air turned and circulated in the kitchen, it disturbed the thermal plume around the cooking station, leading to CO2 escape from the range of the hood and accumulation on the left side of the kitchen space.

3.2.4. CE of the Simulation Cases

We used CE as the hood performance metric to assess the impact of window position on CO2 removal under the two flow rates (Table 4). In general, CE was higher when the windows were in the front wall, followed by the back and side walls. Exceptions were found for the cases with high window positions under the high flow rate (cases K and L). Under the low flow rate, CE ranged from 33 to 83%, with the lowest value for the low window position in the side wall (case B) and the highest for the high window position in the front wall (case D). On the other hand, CE increased by flow rate and ranged from 77 to 97%, depending on the window positions. The most significant impact was observed for the case with a low window position in the side wall, where a 1.8 times increase in flow rate (from low to high) resulted in a 1.58 times increase in CE.

4. Discussion

We found that the CO2 concentration distribution and the hood performance were influenced by the window height, the relative position between the window and the hood, and the exhaust flow rate. We based the discussion on comparisons to findings from previous studies and additional simulation results.

4.1. Distribution of CO2 Concentration and CE by Window Height

We found that the exhaust hood generally performed better in removing CO2 for cases with windows at the height of 2.4 m than 1.2 m. However, for the back window cases I and L under the high flow rate, higher CE was observed for using the lower window. This could be explained by the CO2 leakage from the hood region due to the disturbance of buoyant plume by the high-velocity makeup air from the back window. For example, in case L, the contamination was pushed to the side of the hood, resulting in leakage (Figure 7(b)). When the rising thermal plume was disturbed by the makeup air, it would affect the capture area of the hood. Chen et al. examined the hood performance under the effects of the cross draft from lateral, oblique, and front directions [34]. They found that cross draft could lead to the spillage of cooking emissions, and the effect worsened as the airflow velocity increased. Similar findings were reported by Wang et al., who simulated the particle distribution under three different ventilation modes (only window open, only door open, and window and door open) with three different exhaust hood flow rates (5, 8.64, and 12.96 m3/min) during cooking processes [26]. They found that the incoming air from the window or the door (located at each side of the human model) would interfere with the thermal plume and reduce the hood performance around the cooking area. To confirm the mechanism mentioned above, we conducted an additional analysis to decrease the influence of the cross airflow by modifying the kitchen space to two times larger than the original space. This approach led to an elongated route for the makeup air to reach the cooking station and thus reduced the velocity of the incoming air. The findings from the additional analysis were consistent with those from the other scenarios. Further details are included in SI (SI-5, Figure S3 and S4). In sum, we found that the relative position between the window and the exhaust hood could affect the hood performance.

It is noteworthy that the variable vertical concentration gradient under different flow and window settings had a relatively mild impact on CE because it was determined based on CO2 concentration at the same location in the breathing zone. Additionally, the breathing zone was located in the middle region of the space and was less sensitive to the extreme values of the concentration gradient. As a result, CE could be considered a more robust index for hood performance and suitable for cross-study comparisons.

4.2. Comparison of CE to Other Studies

In this study, we assessed the impacts of six different window locations on CO2-based CE and found that the kitchen hood performed better under the high flow rate with the front, side, and upper back windows open. The performance decreased drastically when the flow rate decreased. The change in hood performance could be well explained by the presentation of the flow characteristics and the distribution of CO2 concentration in the study zone from the simulation results.

CE has been widely applied in previous studies to investigate the influential factors for kitchen hood performance, including hood design (coverage), the range of the flow rate, and burner type and position [13, 14, 22, 35, 36]. However, no experimental studies have investigated the effects of window openings on CE. Despite the differences in ventilation scenarios, CE values reported in this study were comparable to those from the laboratory and the household kitchen environment. Table 5 shows the comparison of CE for the back-burner positions from different studies under the flow rates comparable to the low flow scenarios in this study. In general, previous studies were conducted with windows closed. Their CE was overall >75%, except for the hoods from Singer et al. and was in a similar range to our findings with the front windows open [36]. However, we found that opening the back or the side windows for ventilation could impede the hood performance considerably. This suggested that CE from previous studies might have represented better-case scenarios regarding the protective effects of exposure. Additionally, the two-burner position in the present study is more similar to the back-burner positions in other studies by providing a well-covered range to capture the thermal plume and pollutant emissions. CE for the front-burner positions was found to be lower, with a deduction of 20-25%, as reported by Singer et al. [36].

4.3. Study Strengths and Limitations

In this study, we built the simulation models based on the boundary conditions from the validation tests in the laboratory space that shared common features of a realistic, closed-type kitchen environment. One strength of the study is that we conducted four-point validation tests under steady-state conditions consistent with the assumptions of the numerical models. By having concurrent, continuous measurements at four vertical points, the validation captured the changes in temperature and CO2 concentration gradients. Additionally, there was no interference of CO2 from human breathing in the experimental setting. Another strength lies in practicality. It was unrealistic to install windows in all six positions for the hood performance tests in the same kitchen. Alternatively, field experiments using different kitchen spaces with designated window positions would make the comparisons particularly challenging by introducing variability in additional parameters (e.g., kitchen layout, room volume, and environmental conditions) that could influence the simulation results. As a result, simulation using validated models could fill this gap and provide cost-effective, reliable results to characterize hood performance.

There are also limitations regarding the generalizability of our findings. First, the models should not be generalized to kitchens with distinctly different features, for example, kitchens with open space or single-sided natural ventilation in high-rise buildings. Based on our simulation, a kitchen of a smaller size with a high hood exhaust flow rate was found to have high-velocity makeup air which could disturb the thermal plume and lead to pollutant leakage into the kitchen space. As a result, simulations for open kitchens should be validated separately due to different ventilation scenarios and airflow characteristics. Second, we did not account for the effect of reentrainment of cooking emissions commonly seen in multifamily buildings, where the windows between flats were installed on the same side as the exhausted kitchen air [37]. Gao et al. simulated the upward pollutant transmission through the windows between flats of a four-floor building using CO2 as the tracer gas [38]. The authors found that around 7.5% of the exhausted polluted air from the second floor could reenter the upper flats (third floor) due to buoyancy effects under windless weather. In this case, opening the front window, as suggested by our findings, might not be the optimal approach to reduce human exposure to cooking emissions. Nevertheless, findings from this study provide valuable information on the choice of window opening positions to minimize residents’ exposure to cooking emissions for a typical closed-type kitchen in Taiwanese homes. Furthermore, the validated models could serve as a reference for future studies with more diverse kitchen layouts and environmental conditions.

5. Conclusion

We conducted four-point validation tests of the numerical models based on CO2 concentration and temperature measurements under steady-state conditions. The validated models were subsequently used in simulations to understand the effects of six different window opening positions and two exhaust flow rates on flow fields, CO2 concentration and temperature distributions, and CE for hood performance. We found that the CO2 concentration could be better reduced by having windows open at the higher location. Generally, the front windows were found to be more effective with , followed by the back and the side windows. We also found that as the exhaust flow rate increased from 6.72 to 12.16 m3/min, CE reached >75% for all window positions, where the most significant increase was 1.58 times for the lower side window. To sum up, changing the relative position of the window and the exhaust hood could help disperse the incoming airflow from the window, improve the kitchen’s overall ventilation, and reduce pollutant concentration. Findings from this study provide important implications that when the kitchen window is positioned and used correctly, the hood performance could be well improved to protect residents from exposure to cooking emissions, even under the low exhaust flow rate.

Abbreviations

3D:Three-dimensional
CE:Capture efficiency
CFD:Computational fluid dynamics
CO2:Carbon dioxide
HVAC:Heating, ventilation, and air conditioning
HVI:Home Ventilating Institute
IAQ:Indoor air quality
NO2:Nitrogen dioxide
PAHs:Polycyclic aromatic hydrocarbons
PM:Particulate matter
SIMPLE:Semi-implicit method for pressure-linked equations
Temp.:Temperature
VOCs:Volatile organic compounds.
Nomenclature
:Air density
and :Fluid velocity in -, -, and -directions
:Time
:Pressure
:Viscous shear stress tensor
:Reynolds stress tensor
:External force component in -, -, and -directions
:Diffusion heat flux component in -, -, and -directions
:Turbulence dissipation rate
:Internal heat generation
:Specific enthalpy
:Thermal enthalpy
:Viscosity
:Kronecker delta function
:Turbulent kinetic energy
:Eddy viscosity.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Disclosure

A preprint has previously been published (https://papers.ssrn.com/sol3/papers.cfm?abstract_id=4261535) [39].

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgments

This work was supported by Taiwan National Science and Technology Council (MOST 108-2314-B-002-101-MY3) and the Higher Education Sprout Project of National Taiwan University (NTU-111L7256).

Supplementary Materials

Supplemental Information SI-1: characteristics of the kitchen hood. SI-2: air exchange rate measurement. SI-3: specification of the instruments used in experimental measurement. SI-4: distribution of temperature by window position and hood flow rate. SI-5: additional analysis with two times the kitchen space in cases I and L. Table S1: specification of the instruments used in experimental measurement. Figure S1: front view of the kitchen exhaust hood. Figure S2: distribution of temperature in the plane in the breathing zone. Figure S3: distribution of the CO2 concentration in the plane of the breathing zone. Figure S4: distribution of the CO2 concentration with flow trajectories. (Supplementary Materials)