Dynamics and Numerical Simulation of Contaminant Diffusion for a Non-Flushing Ecological Toilet

The poor indoor air quality (IAQ) of severely polluted toilets is associated with increased risk of severe disease. This study aimed to evaluate the overall IAQ according to the contaminant removal efficiency, volume average concentration, and breathing zone control level. The characteristics of contaminant transmission in a non-flushing ecological toilet (NFET) were analyzed using the computational fluid dynamics (CFD) methodology, and the proposed model was further validated based on experimental measurements. Both an orthogonal experimental design and CFD were used to analyze factors such as exhaust fan position (EFP), air change rate per hour (ACH), natural vent location (NVL), and grid height (G-h). The EFP and ACH were demonstrated to be the dominant factors affecting the IAQ, whereas NVL and G-h were found to play key roles. Single-factor analysis based on the significance levels of the ACH, EFP, and NVL was conducted using the CFD methodology to define three exhaust behaviors—namely, “ineffective”, “enhanced”, and “excessive”. These results provide key insights that may be used to improve the IAQ of NFETs.


Introduction
Restrooms, which are among the most frequently used rooms in daily life, emit large amounts of undesirable odors [1]. According to statistics from the World Health Organization, approximately 2.3 billion people worldwide lack sanitary restroom facilities, of which approximately 946 million people live without conventional restroom facilities.
Toilets contain numerous contaminants, and feces and urine produce odors as they ferment. Sato et al. [2] analyzed the malodorous volatile substances in human waste (feces and urine), and reported that approximately 90% of the malodor-causing substances are fatty acids, of which approximately 8% are N-containing compounds-including ammonia (6.5%), pyridine, pyrrole, indole, skatole, and trimethylamine. In addition, S-containing compounds (2%) such as hydrogen sulfide (1.6%) and methyl mercaptan were detected.
Various deodorization methods have been developed to deal with this problem: chemical deodorization [3], biological deodorization [4], activated sludge aeration [5], activated sludge scrubbing [6], physicochemical deodorization [7,8], and sensory deodorization, which employs aromatic and masking chemicals [9]. However, each of these solutions presents both advantages and disadvantages. Although some public restrooms have adopted air fresheners to cover odors, air fresheners have been identified as the primary example, Tung et al. [31] studied a new negative-pressure wall-type bathroom exhaust ventilation system, and performed an experiment-based analysis of the effects of the air change rate per hour (ACH) and the distance between the pollution source and exhaust vent on the distribution of the gaseous pollutant concentration in the restroom.
The main volatile substances in human waste are water soluble, and other than ammonia, these are denser than air when in a gaseous state. Tung et al. [32] introduced a local air quality index and exposure index for medical staff to evaluate the ventilation effectiveness. They performed experiments to determine the influences of various negative pressure differences on the distribution of contaminants in toilets, using sulfur hexafluoride (SF 6 ) tracer gas instead of pollutants to study the ventilation efficiency and air exchange rate of the test room. The results showed that the ventilation system performed best when the room was maintained at −15 Pa. In a separate study, Tung et al. [33] replaced restroom contaminants with SF 6 and measured the concentrations of contaminant gases in the restroom under three different ventilation modes to compare the ventilation efficiencies of these modes. The results showed that indoor contaminant concentration was lowest when the ventilation system that combined air supplies from the ceiling of the house and side wall exhaust was adopted.
Non-flushing ecological toilets (NFETs) are environmentally friendly toilets that do not require a sewer but directly collect feces and urine to a collection tank and conduct harmless treatment [34]. The benefits of non-sewered onsite sanitation systems are increasingly being recognized [35]. These are in great demand at construction sites, in remote areas, at Fangcang Hospital, in drought-affected areas, and in other areas where sewer systems are not available [36]. It can achieve the main aim of sanitation: protecting public health and the environment [37].
In addition, because NFETs reduce the water and wastewater flows within buildings, the United States Green Building Council has used NFETs to achieve leadership in energy and environmental design certifications [38]. Although several studies have been conducted on the diffusion of contaminants in private restrooms, hospital wards, and flushing toilets, limited studies have been conducted on the odor diffusion dynamics of NFETs.
Conventional ventilation systems in public restroom feature small exhaust fans, typically mounted in the ceiling or wall of the restroom, which extract unpleasant odors through the exhaust ducts. However, a greater ACH is required for systems with longer contaminant diffusion paths; in such cases, the concentration of contaminants in the immediate vicinity of restroom users is higher, and the wind speed in the restroom becomes too high, with a noticeable blowing sensation affecting the indoor comfort levels.
In this study, we used a combination of numerical simulations and experiments to study the odor concentration reduction in NFETs, to improve the air quality. The water-free ecological toilets considered in this research represent a departure from traditional toilets in several aspects. For example, although the pollutants in traditional toilets diffuse from the toilet location, those in the NFETs considered in this study diffuse from the surface of the fermentation tank. The manure produces pollutants, such as ammonia and hydrogen sulfide, during the fermentation process. Several previous studies have considered the influence of the ACH on the pollutant removal efficiency; however, the effects of air inlet and outlet positions on the pollutant diffusion dynamics are often overlooked. Therefore, we performed orthogonal experiments to study the effects of different factors on the IAQ of NFETs. In particular, we focused on the influences of the ACH, exhaust fan installation position, and air inlet position and height. Figure 1 illustrates the configuration of the test bench adopted in this experimental study, which measured 0.96 m (length) × 0.96 m (width) × 2.2 m (height), and the pollution source measuring 0.5 m (length) × 0.2 m (width), where the length, width, and height correspond to the x-, y-, and z-directions, respectively. Representative odor-release experiments were performed on the pollution source. The assumptions and simplifications involved in this study are as follows. (1) The effects of differences between indoor and outdoor temperatures on the contaminant diffusion process were not considered. (2) Pollutants diffuse uniformly from the surface of the pollution source at a constant rate. Monitoring points for collecting the experimental data were placed above the center of the cubicle (see Figure 1a).

Tracer Gas Release and Measurement
We selected SF 6 as the tracer gas because it has properties similar to those of the main volatile substances in human waste, being colorless, soluble in water, and denser than air. Under normal conditions, SF 6 is non-toxic and non-flammable, rendering it perfectly safe. In addition, it is commonly used as a tracer gas in short-term experiments to determine the ventilation efficiencies of buildings and indoor enclosures.

Data Acquisition Equipment
An XLA-BX-SF 6 portable SF 6 gas detector was selected to monitor the behavior of the tracer gas during the experiments (Figure 1d). The detector was equipped with new original imported sensors and a 32-bit low-power processor, which included a circuit based on four-layer wiring to enhance its performance against weak signals and interference. In addition, it featured a built-in high-performance suction pump, stable flow rate, and switchable control based on an independent high-precision imported chip with a unique signal processing algorithm, which resulted in accurate and stable measurements. The XLA-BX-SF 6 detector (Shenzhen Bridgecom Electronic Technology Co., Ltd., Shenzhen, China) had an error of ±0.01 ppm, with data recorded every 30 s for a duration of 45 min at each measurement point. A PNT400-H2S/O3/NH3-G portable ammonia and hydrogen sulfide gas detector was used to detect the concentration of pollutants in public toilets with a resolution of 0.001 ppm (Figure 1e).

Validation of the CFD Simulation Method
CFD modeling enables efficient analysis of different configurations under various conditions. However, CFD modeling cannot replace experimental methods, as experiments are necessary to support the CFD modeling as benchmark tests [39,40]. Experiments were conducted in the school laboratory of Tongji University, using the measurement point locations indicated in Figure 1a. The laboratory was ventilated prior to the experiment to ensure that the results would not be distorted by any residual SF 6 gas.
The fluid dynamics of the system were simulated using ANSYS Fluent, to determine details such as the concentrations of the various contaminants. For the simulations, the component transport model was adopted for the contaminant inlet. The standard, RNG, realizable, and low-Re k-ε models were used to simulate the gas turbulence [28,[41][42][43][44]. These models consider the natural buoyancy-driven convection during the diffusion process and use the finite volume method and pressure-velocity coupling format of the SIMPLEC algorithm to obtain the discrete solution of the control equation. The CFD model was built according to the conditions applied in the experiments (Figure 1). An unstructured computational grid of polygon elements with a mesh quantity of 1.09 million, which was refined at the relevant zone, was used; it was chosen after independent verification (see Figure 2a). As shown in Figure 2b, the results of the CFD simulation were compared with the experimental data collected under identical conditions at each monitoring point. The results obtained using the standard k-ε model were rational, revealing maximum and minimum disparities of 10% and 4%, respectively. Overall, Figure 2b shows a close agreement between the experimental and simulation data, indicating that the proposed model is suitable for predicting the distribution of indoor contaminants.
The standard k-ε model agrees well with the experimental data and has been widely adopted [45][46][47]. Therefore, considering the model performance, the standard k-ε model was selected for use in this study.

Boundary Conditions
The following boundary conditions were adopted in our simulations [11,48]: the exhaust outlet was set in accordance with the velocity outlet boundary conditions (velocity outlet); the grid and window were set as pressure inlets (pressure inlet); the walls and partitions comprising the restroom were defined as wall; the pollution source was set as mass flow outlet; and the ammonia release concentration and hydrogen sulfide release concentration were set to 2.5 × 10 −7 kg/s and 6 × 10 −8 kg/s, respectively [11,48].
A polygon element grid, which is an unstructured computational grid, was applied. Grid independent verification was conducted for different mesh numbers, namely, 1,840,000, 2,330,000, 3,130,000, 3,930,000, and 7,080,000. There was no significant change in the pollutant concentration for mesh numbers above 2,330,000. To enhance calculation efficiency and accuracy, a mesh number of 2,330,000 was therefore chosen (Figure 2c). The iterative calculations continued until a prescribed relative convergence of 10 −6 was satisfied by all the field variables.  Airflow movements in NFETs are so complex that it is difficult to study a number of influential factors simultaneously. Thus, a preliminary investigation is necessary.

Preliminary Study
As demonstrated in a previous study [31], the exhaust fan position (EFP) is an important factor with regard to contaminant removal. To obtain a deep insight into the influence of the EFP, two sets of cases were designed with different EFPs (see Figure S2: Case 1-F (0,0.62,0.59), Case 2-F (0.07,0.62,0.85)), unit: m. The pollutant concentration variation with height was studied at five points A, B, C, D, and E (see Figure 3c). National standards stipulate that the ACH should be greater than 5 h −1 (GB/T 17217−1998 for urban public toilets, CJJ 14−2016) [49,50], and 10−15 h −1 (National Civil Building Engineering Design Technical Measures HVAC Power) [51]. Therefore, in our preliminary research, air supply rates of 5 h −1 , 10 h −1 , 15 h −1 , and 20 h −1 were selected. The concentration standards for ammonia and hydrogen sulfide in public toilets are shown in Table 1 [49]. The ACH clearly affects the contaminant concentrations, which gradually decrease with increasing ACH. The contaminant concentrations along the vertical lines are shown in Figure 4, and Figures S3-S6. Figure 4a,b show the changes in ammonia concentration as functions of the ACH and height, respectively, whereas Figure 4c,d depict the corresponding plots for varying hydrogen sulfide concentration. As shown in Figure 4a, at an ACH of 5 h −1 , the ammonia concentration decreases slowly as the height increases, before stabilizing. At an ACH of 10 h −1 , it is concentrated between 0 m and 0.5 m. Within this zone, the ammonia concentration decreases rapidly with increasing height. In contrast, between 0.5 m and 2.2 m, the ammonia concentration first increases, stabilizes, and then decreases.
As the ACH is increased to 15 h −1 and 20 h −1 , the ammonia concentration approaches 0 between the heights of 0.5 m and 2.2 m. The variation of the hydrogen sulfide concentration follows similar trends, as depicted in Figure 4c.
As shown in Figure 4b, at an ACH of 5 h −1 , the ammonia concentration decreases rapidly as the height increases from 0 m to 0.1 m. Above this height, the ammonia concentration tends to 0 mg × m −3 . Moreover, when the ACH increases to 10 h −1 , 15 h −1 , and 20 h −1 , the contaminant concentration tends to 0; as shown in Figure 5d, the maximum hydrogen sulfide concentrations for these ACH values at the measurement locations are less than 8 × 10 −5 mg × m −3 . The ammonia concentrations along the vertical lines normal to the corners of the rectangle ABCD in Cases 1 and 2 are presented in Figures S3 and S4, respectively, and the corresponding hydrogen sulfide concentration measurements are provided in Figures S5 and S6, respectively. Along these vertical lines, the contaminant concentration varies with the ACH in a similar manner to that observed along the vertical line at point E. At an ACH of 5 h −1 , the contaminant concentration is relatively high in both cases. In response to increasing the ACH to 10 h −1 , the ammonia concentration increases only in Case 1. For ACH values exceeding 10 h −1 , the contaminant concentrations are relatively low in both cases. As shown in Figures 5 and 6, for the equivalent ACH values, the contaminant concentrations in Case 2 are much lower than those in Case 1, indicating that the CRE is higher in Case 2. Figures 7 and 8 also indicate that Case 2 results in more effective contaminant removal, showing that vortices are generated in the restroom when the ACH is 10 h −1 .    When the exhaust fan was installed under the side of the toilet, the pollutant concentrations were the lowest. Figure 8 exhibits the contours representing the mean age of air at different cross-sections. In general, the age of air is defined as the elapsed time since an air element entered the room. The distribution of the mean air age can reveal the airflow direction on the basis that air always moves from locations with lower air ages to locations with higher air ages. Thus, the air age can be used to evaluate an entire room. In addition, as the air freshness reflects the ability of the room to remove pollutants, it is widely used as an evaluation index for measuring IAQ.
At an ACH of 5 h −1 , the air age in Case 1 is higher, meaning that the contaminants exist for a longer time in the restroom. In Case 2, the air age increases with increasing ACH, indicating that the ACH is larger than necessary, inducing eddy currents that affect the contaminant diffusion. The research results in this section show that the pollutant removal effect is good when the exhaust fan is installed in the lower part, which paves the way for further research.

Orthogonal Experimental Design
For maximum efficiency, we designed our study to provide comprehensive testing while minimizing the number of experiments required. Typically, experiments are designed according to either orthogonal or uniform methodologies, with the former based on the equilibrium distribution theory.
In recent years, orthogonal experimental design (OED) has gained popularity in many research areas as an approach for optimizing the system performance [42], and is a reliable and representative method. Specifically, researchers have used the OED to explore the relationships between the experimental factors and related indicators, such as the CRE and volume average concentration (VAC). In this study, we adopted the OED to analyze the ACH, EFP, natural vent location (NVL), and grid height (G-h). As shown in Figure 7, when the number of air changes increases from 3 h −1 to 4 h −1 , the ammonia concentration drops sharply. Based on the results presented in Section 3.1, the ACH values of 3.5 h −1 , 4 h −1 , 5 h −1 , and 6 h −1 were selected, and exhaust fans were positioned below the fermentation tanks F1, F2, F3, and F4. Furthermore, the values of 50, 100, 200, and 300 mm were chosen for G-h, and natural vents were positioned at L 1 , L 2 , L 3 , and L 4 (see Figures S1b and S3a). An L 16 (4 5 ) orthogonal experiment was performed to determine the effects of various factors on the IAQ by enabling multifactor testing with the minimum number of experiments. The orthogonal test schemes are listed in Table S1.
Simulations were designed to consider the CRE and VAC with various factors and levels. Analyses of range and variance were performed on the collected data to identify the orders of magnitude of these factors. The full details of the calculation process employed are provided elsewhere [52].

Evaluation Index
To conduct a comparative analysis of the different ventilation strategies, the normalized values of contaminant concentration and ventilation effectiveness were used in the analysis.
The CRE evaluation index is expressed as where C e , C in , and C represent the contaminant concentrations at the exhaust outlet and air supply outlet and the average concentration in the immediate environment, respectively.

Orthogonal Significance Analysis of Relevant Factors
As shown in Figure 9, the orthogonal experiments indicate that the ACH, EFP, and NVL are statistically significant factors affecting the CRE and VAC, whereas G-h is less influential. The F-value calculation results are as follows: The results of the orthogonal experiments clearly demonstrate that ACH is the dominant element determining the VAC: when the ACH increases from 3.5 h −1 to 6 h −1 , the VAC values for ammonia and hydrogen sulfide decrease by 41% ( Figure S7) and 37% ( Figure S8), respectively. The second most influential factor is the EFP, with the VAC of ammonia decreasing by 14.7% from F4 to F1 ( Figure S7), whereas the hydrogen sulfide VAC decreases by 8% ( Figure S8). However, the significances of the ACH and EFP are reversed with respect to the CRE, with the EFP emerging as the dominant element. In response to changing the exhaust fan from F4 to F1, the ammonia removal efficiency increases by 1.3% (Figure 10), whereas the hydrogen sulfide removal efficiency decreases by 7.2% ( Figure S9). As the next most influential factor, increasing the ACH from 3.5 h −1 to 6 h −1 results in the ammonia and hydrogen removal efficiencies decreasing by 1.2% ( Figure 10) and 6.4% ( Figure S9), respectively. Therefore, contaminant type should be considered when positioning the exhaust fan.   Figure 9 shows the significance of a few elements, including the NVL and G-h. The NVL is statistically significant, with a p-value < 0.05 (see Figure 9a,b): changing the NVL from L 1 to L 2 results in the VAC and CRE of ammonia decreasing by 12.7% and 0.6%, respectively; and changing the NVL from L 3 to L 2 causes the VAC and CRE of ammonia to decrease by 11.6% and 1.9%, respectively. The G-h is statistically significant (p-value < 0.05) in Figure 9a, but not in Figure 9b. In addition, whereas the NVL is statistically significant (p-value < 0.1) in Figure 9c,d, the G-h is not. As demonstrated in Figure 9, for ammonia, the factors can be ordered from most to least influential as EFP > ACH > NVL > G-h.
Regarding hydrogen sulfide, the order of significance is ACH > EFP > NVL > G-h. According to previous research, the sequence of these factors regarding contaminants is EFP > ACH > NVL > G-h.
The factors examined in this study have different effects on ammonia and hydrogen sulfide because these two pollutants have different concentrations and obey different diffusion characteristics. Ammonia is less dense than air and has a higher diffusion rate, meaning that it diffuses to the upper part of the toilet more easily. In contrast, hydrogen sulfide gas is denser than air and has a slower diffusion rate, meaning that it tends to be deposited in the lower part of the toilet. The farther the exhaust fan is from the pollution source, the easier it is to extract and eject ammonia; however, the result differs in the case of hydrogen sulfide, with the ACH having a greater impact on its removal.
In short, the ACH and EFP are the most influential factors affecting the distributions of ammonia and hydrogen sulfide, respectively. Therefore, we focused on the effects of these two factors on the pollutant diffusion.

Single-Factor Analysis
After the ACH and EFP, the orthogonal experiments identified the NVL as the third most important parameter affecting the CRE and VAC. As each of these factors affects the CRE and VAC in different ways, identifying the interactions among the ACH, EFP, and NVL is key to understanding the effective use of exhaust fans for contaminant removal. Therefore, we conducted single-factor analyses for the ACH and EFP, with the G-h fixed at 300 mm and the natural vent height fixed at 2050 mm. The observation plane was located in the middle of the restroom (see Figure S1b in the Supplementary Materials). Figure 11 shows the changes in contaminant concentration and removal efficiency as functions of the ACH for ammonia and hydrogen sulfide. At an ACH of 0 h −1 (i.e., the exhaust fan is turned off), the contaminant concentration is high. As the ACH increases, the VAC of hydrogen sulfide decreases (Figure 11c).

Effect of the ACH
In contrast, in the case of ammonia, the VAC does not change between the ACH values of 0 h −1 and 3 h −1 , although when the ACH exceeds 3 h −1 , the VAC decreases rapidly. Furthermore, the ammonia concentration distribution results shown in Figure 12 indicate that the ammonia originates from the surface of the fermentation tank and diffuses through the toilet to the restroom. When the ACH is low, the ammonia is distributed uniformly in the space, as shown in Figure 12b. However, as the ACH increases, the contaminant concentrations in the space below decrease, and fewer contaminants spread to the top. As shown in Figure 13, the concentrations of hydrogen sulfide and ammonia gas exhibit similar trends in response to the ACH changes. At the same ACH, the concentration of hydrogen sulfide is lower than that of ammonia because hydrogen sulfide is denser than air, which restricts its upward diffusion.    The velocity contours in Figure 14 help illustrate the concentration distribution dynamics. When the ACH is 0 h −1 , the ammonia gas diffuses upward freely, forming a vortex in the upper space, and it exists in the space for a longer time. Increasing the ACH to 2 h −1 results in a reduction in the upward emission rate of ammonia; however, the ammonia mixes fully with air. Therefore, the contaminant concentrations in the upper space are extremely high. When the ACH increases to 4 h −1 , the air enters the toilet from the upper discharge port and is discharged together with the contaminants through the exhaust fan. Consequently, the airflow direction changes to inhibit the upward diffusion of ammonia. Therefore, the ammonia concentration in the upper part of the toilet is extremely low.
Compared with an exhaust without a fan, at an ACH of 4 h −1 , the VAC of ammonia is reduced by 71.7% on average, with this reduction increasing to 85% and 88% on average at ACHs of 8 h −1 and 10 h −1 , respectively. As in the case of ammonia, hydrogen sulfide diffuses upward freely and is distributed uniformly in the space at ACH = 0 h −1 . At an ACH of 2 h −1 , the hydrogen sulfide concentration in the upper space is extremely low, although this remains high in the lower space. However, increasing the ACH to 4 h −1 also results in a reduction in the hydrogen sulfide concentration in the lower space. Hydrogen sulfide and ammonia exhibit similar diffusion behaviors. When utilizing an exhaust without a fan, the VAC is reduced by 61.1%, 82%, 89.4%, and 89.8% on average for the ACH values of 2 h −1 , 4 h −1 , 8 h −1 , and 10 h −1 , respectively. The concentrations of indoor pollutants in large spaces were explored by Duci et al. [53] using CFD modeling, with the results showing that reasonable mechanical ventilation can reduce the concentrations of indoor pollutants. Specifically, at the ventilation speeds of 1 m/s and 2 m/s, the pollutant concentrations in the local respiratory area of restroom users were reduced by 30% and 50%, respectively.
As shown in Figure S10, when the volume of air in the exhaust is relatively small, the pollutants diffuse into the upper space and are discharged from the grille and the upper air outlet. When the ACH exceeds a certain value, the grille and upper air outlet become air inlets, allowing fresh outdoor air to enter the room and improve the IAQ. As the ACH increases, an increasing amount of air enters the room, forming a vortex above the inlet and a significant sense of air blowing is felt by occupants. Therefore, the ACH should be maintained within a certain range to avoid causing discomfort to restroom users. Figures 15 and 16 show that the B-Z control level, that is, the average contaminant concentration in the B-Z, and VAC exhibit similar behaviors. For example, Figure 15 demonstrates that at an ACH of 4 h −1 , a small amount of contaminant accumulates at the corners of the B-Z, whereas the average concentration approaches 0 mg × m −3 . Moreover, Figure 16 indicates that an ACH exceeding 3 h −1 is sufficient to achieve a first-class standard regarding average contaminant concentration in the B-Z.  As the ACH increases, the CRE first increases, then decreases, and finally stabilizes, as shown in Figure 16. Hydrogen sulfide and ammonia removal efficiencies reach their peaks when the ACH is 3 h −1 and 4 h −1 , respectively. The contaminant concentrations inside the toilet and the exhaust fan outlet gradually decrease with increasing ACH. These findings indicate that the ACH and CRE follow a nonlinear relationship in which the gain ratio is reduced when the ACH exceeds a certain value. Accordingly, the effect of the ACH can be divided into three stages. The first of these is the 'ineffective exhaust' stage, which is defined as the stage in which the contaminant cannot be removed, owing to an insufficient air supply velocity. At an ACH of 3 h −1 , the concentration of ammonia gas is equal to that when there is no exhaust fan, and the ammonia gas is not fully discharged. Second, the 'enhanced exhaust' phase is defined as the stage in which the air supply velocity exceeds the critical velocity required to maintain a stable discharge of contaminants. Additional amounts of contaminants may be discharged by increasing the ACH, whereas the concentration of contaminants in the B-Z can be reduced to below the standard concentration required for an acceptable IAQ in restrooms.
Finally, 'excessive ventilation' is defined as the stage in which the average contaminant concentration in the B-Z inside the restroom does not decrease as the ACH increases because the concentration is close to 0 mg × m −3 and no longer changes. When the ACH is greater than 4 h −1 , contaminant concentrations in the B-Z cease to change; therefore, increasing the volume of air in the exhaust has little effect on the contaminant discharge. Figure 17a,b show that the ammonia gas concentration and CRE decrease and increase (nonlinearly), respectively, as the ACH increases. To some extent, the CRE is influenced by the EFP. When the ACH is between 3.1 h −1 and 3.3 h −1 and the exhaust fan is located at F1 or F2, the VAC is reduced, and CRE is maximized. A slight increase in the ACH to between 3.4 h −1 and 3.5 h −1 and relocating the exhaust fan to F3 also maximizes the CRE. At an ACH of 3.6 h −1 , the concentrations of all contaminants are approximately equal. Furthermore, Figure 17c,d show that the contaminant concentrations in the B-Z change similarly with the ACH. When the ACH is less than 3.5 h −1 , the substances in the exhaust do not satisfy the first-class standards. At an ACH of 3.5 h −1 and with the exhaust fan positioned at F3 rather than F4, the first-class standards may be achieved. By increasing the ACH slightly to 3.6 h −1 , the first-class standards are satisfied for all the configurations examined in this study.

Effect of the Exhaust Fan Position
As shown in Figures S11 and S12, at an ACH of 3.3 h −1 , there are differences in the distribution of pollutant concentrations in the upper part of the toilet. The closer the exhaust fan is to the pollution source, the shorter the time spent by the pollutant in the toilet and the lower the pollutant concentration in the toilet. When the exhaust fan is located at F1 or F2, fewer vortices are generated close to the pollution source, and the pollutants are discharged more easily. As shown in Figures S13 and S14, when the ACH is 3.6 h −1 , the upper space of the toilet has very low pollutant concentrations. There is almost no difference between the restroom configurations. The pollutants exhibit minimal upward diffusion, with the majority of pollutants discharged through the exhaust fan. However, at lower pollutant concentrations, there is a certain difference: when the exhaust fan is positioned at F 3 or F 4 , the high-concentration area is relatively large.
In summary, increasing the ACH causes the vortex in the upper space to grow, whereas part of the vortex is generated around the fermentation tank, which is not conducive to pollutant discharge. Figure S15a,b show that the average hydrogen sulfide concentration decreases with increasing ACH. Specifically, when the ACH increases from 2 h −1 to 2.1 h −1 , the concentration of hydrogen sulfide decreases rapidly, whereas the ventilation efficiency increases rapidly. When the ACH is equal to at least 2.1 h −1 and the exhaust fan is in the F1 position, the pollutant concentration in the B-Z is at its lowest and the ventilation efficiency is at its highest. Conversely, relocating the exhaust fan to F4 results in the lowest observed ventilation efficiency. Figure S15c,d show that when the ACH is at least 2.2 h −1 , the hydrogen sulfide concentration in the B-Z for each configuration meets the requirement for the first-class standard. Chung et al. [54] investigated the ventilation efficiencies of different ventilation configurations involving two inlet and two outlet diffusers at different locations. The results indicated that floor-based exhaust ventilation systems are effective in removing pollutant-related odors, with only 1.2% and 1.5% of pollutant particles remaining in the male and female restrooms, respectively. Positioning the exhaust fan closer to the pollution source facilitates contaminant removal. However, as ammonia gas is less dense than air, it diffuses upward easily. In contrast, as hydrogen sulfide is denser than air, it accumulates below. The different spatial distributions of the ammonia and hydrogen sulfide arise from their different concentrations. The fermentation tank obstructs air circulation and changes the airflow direction, which causes vortices in some areas; consequently, the ventilation efficiency does not change linearly with increase in the ACH. In addition, the vortices mean that the optimal exhaust position depends on the ACH value. According to our rigorous analysis, an ACH of 3.5 h −1 provides the best solution overall. At this value, the average concentrations of ammonia and hydrogen sulfide satisfy the first-class standards, with the highest ventilation efficiency when the exhaust fan is placed at F 4 .

Effects of the Natural Vent Location
The positions of natural vents (L 1 , L 2 , L 3 , and L 4 ) also influence the pollutant diffusion. To analyze this influence, the position of the exhaust fan was fixed at F 4 , and the VAC and CRE were calculated for ammonia and hydrogen sulfide, with the results shown in Figures S16 and S17. First, Figure S16 shows that when the ACH is between 3 h −1 and 3.6 h −1 , the VAC is at its minimum (and the CRE is at the maximum) when the natural discharge port is located at L 1 . However, when the ACH exceeds 4 h −1 , placing the natural discharge port at L 3 yields the highest CRE. There is little variation in VAC between the various configurations. As shown in Figure S17, when the ACH is between 2 h −1 and 2.5 h −1 , the VAC and CRE of each configuration are similar; however, when the ACH exceeds 3 h −1 , positioning the natural discharge port at L 3 results in the highest CRE and lowest VAC. When the exhaust fan discharge is low, the natural air intake in the upper part of the toilet has a low air volume and the airflow disturbance is less pronounced, reducing any negative effects on restroom user comfort.

Conclusions
CFD simulations were designed and performed to determine the CRE to investigate the effects and weights of various factors that influence the CRE and B-Z control level in orthogonal experiments on NFETs. Single-factor analysis was conducted based on the results obtained from the orthogonal tests. The main conclusions of this study are as follows: 1.
Common toilet ventilation factors (e.g., the ACH) and factors specific to NFETs (e.g., the EFP, NVL, and G-h) were ranked based on their statistical significance as follows: EFP > ACH > NVL > G-h, with the EFP achieving statistical significance (p-value < 0.01) in the case of ammonia. The CRE is high when the exhaust fan is installed in the lower part of the toilet (near the fermentation tank).

2.
Contaminant concentration distributions were simulated for different exhaust positions and air changes, with the results showing that both the ACH and exhaust fan location must be considered in toilet exhaust design. Toilet ventilation efficiency may be optimized by installing the exhaust at the optimal location, thus maximizing the air quality improvement inside the toilet. 3.
The VAC decreases with increasing ACH, with the rate of decrease decelerating gradually. In contrast, the CRE first increases and then decreases before finally stabilizing in response to increasing ACH. The CRE varies between 1 and 3 but does not exceed 3, which may be related to the small size of the toilet investigated, and the relatively low pollutant concentrations contained within it. 4.
Single-factor analysis revealed three stages of exhaust behavior, namely, "ineffective", "enhanced", and "excessive". Beyond guaranteeing sufficient air quality, an appropriate number of air changes should be selected to reduce the energy consumption and indoor air speed, thus reducing the blowing sensation and improving user comfort.
The findings presented in this study provide a theoretical basis for the design of toilets with effective ventilation systems.