Abstract

A large amount of surgical smoke in electrosurgery seriously deteriorates the clean environment of the operating room and can potentially harm medical staff and patients. Exploring the distribution and removal of indoor particulate matter and selecting efficient ventilation patterns are effective ways to control harmful smoke. Therefore, in this study, we combined simulations and full-scale experiments to quantitatively explore the high-concentration spatial regions of particles and compared three ventilation patterns: vertical laminar airflow (VLAF), horizontal laminar airflow (HLAF), and hybrid ventilation, wherein unidirectional airflow (UDAF) was applied to the operating table along with peripheral mixing (UDAF + mixing). We found that simple laminar flow ventilation was significantly affected by the equipment layout and air change rate (air changes per hour; ACH), and the smoke particles were distributed in large amounts in the operating area and could not be removed completely. Conversely, hybrid ventilation can work effectively, and the optimal ACH is approximately 60, which can remove nearly 72% of smoke particles. The airflow distribution in the operating room is also an important factor affecting the distribution and removal of smoke particles. Therefore, medical staff should avoid prolonged exposure to areas with high particle concentrations and particle removal paths.

1. Introduction

Air pollution control in the operating room is a basic requirement to maintain the health of patients and medical staff [1, 2]. However, in many surgical procedures, high-frequency electric knives cut tissue by heating and burning, resulting in a large amount of surgical smoke [35]. These particles, which carry harmful substances, are often considered a major contamination source during surgery [6]. In several studies, electrosurgical particles have been found to be susceptible to airflow and quickly spread far from the wound [2, 7], spreading disease to healthcare workers through viral and bacterial pathogens in the smoke particles [8, 9]. Therefore, it is essential to study the ventilation patterns and the distribution and removal of smoke particles during the operation process to provide effective suggestions for maintaining a clean operating room environment.

Several standards and guidelines have addressed surgical smoke exposure. The Occupational Safety and Health Administration (OSHA) states that nearly half a million healthcare workers are exposed to electrosurgical smoke in operating rooms [10]. Although healthcare workers take precautions during surgery, regular masks can only protect against particles larger than 0.5 microns. A large number of smoke particles are fine particles with a diameter of about 0.3 microns or even smaller ultrafine particles [11]. Exposure and inhalation of surgical smoke can lead to acute adverse effects and severe health risks for patients and staff [3, 12, 13]. Many studies have shown the presence of dangerous substances in surgical smoke, including carcinogens, harmful chemicals, bacteria, and viruses, and long-term exposure can lead to eye irritation, headaches, nausea, weakness, and even the risk of cancer [14]. Currently, OSHA standards do not require ventilation and filtration systems, nor do they set specific particulate limits for controlling smoke from electrosurgery [10]. Increasing studies have focused on the properties and hazards of volatile organic compounds in smoke particles, but little attention has been paid to intraoperative particle distribution and concentration monitoring.

Ensuring high air quality in operating rooms is essential for reducing smoke particles and protecting the health of medical staff and patients. Among them, the use of ventilation patterns is the most effective way to control the risk of cross-infection [1519]. Laminar airflow (LAF) ventilation is currently widely used in operating rooms and other clean rooms to effectively create a sterile environment around patients and healthcare staff [2022]. Numerous studies have experimentally or numerically demonstrated that LAF is superior to mixed ventilation in reducing the risk of patient exposure [23]. Hansen et al. experimentally studied the effect of ventilation on the transport of fine and ultrafine particles in the operating room [15]. The results show that the LAF system is superior to the turbulent mixed airflow system in removing particles. Although LAF is used as a general design requirement, its operation remains flawed, and there is a risk of cross-infection between patients and healthcare staff [24]. Previous studies have shown that the direction of vertical laminar airflow (VLAF) is easily affected by surgical lamp blockage [25]. Additionally, the thermal plume and updraft caused by the temperature difference between the body surface and the surrounding environment greatly reduce the effect of the vertical airflow created by the VLAF on particle removal [23, 26]. As another form of laminar ventilation, horizontal laminar airflow (HLAF) greatly attenuates the influence of operating lights and medical staff on the concentration of particulate matter in the operating area. However, there are higher requirements in terms of the design parameters and operating room layout, and the stability is not guaranteed [27]. In recent years, owing to the complexity of the operating environment, more medical equipment has been used in the operating room, forming a hybrid between the ordinary operating room and the diagnostic room. Different forms of ventilation (unidirectional airflow (UDAF) on the operating table and a peripheral hybrid ventilation system) were used in each area. Romano et al. examined air pollution in operating rooms during surgery using this ventilation system. Their results show that operating rooms equipped with this system can better dilute and eliminate air pollution [3]. Different ventilation patterns have their own advantages and disadvantages in maintaining a clean environment in the operating room. Therefore, it is important to compare the spatial concentration distribution of smoke particles in an operating room under different ventilation patterns. Meanwhile, increasing research has focused on reducing the risk of wound infection during surgery [28, 29], and little attention has been paid to the effects of particulate matter on healthcare workers. Notably, cross-infection is bidirectional; thus, detailed comparisons and analyses of particulate concentrations near patients and healthcare workers are also of interest.

Air changes per hour (ACH) are an important factor affecting the ventilation effect of different ventilation systems [28, 30]. Romano et al. studied the relationship between ACH and particulate matter concentrations in the operating room, and an increase in ACH tended to reduce the ambient particulate matter concentrations [31]. A previous study showed that ACH in the operating room had a great impact on the indoor airflow distribution and the particulate matter distribution in the operating area and surrounding areas [32]. Another study showed that ACH had a greater impact on the distribution of bioaerosols in the surgical area than that in the temperature conditions [33]. Therefore, it is important to consider the influence of ACH on the distribution of intraoperative smoke particles.

The migration rule, spatial concentration distribution, and removal efficiency of smoke particles under the three ventilation patterns were studied in a standard class I operating room. In the experiment, electrosurgical equipment was used to cut through pork tissue to create a smoky environment. The operating area is an important part of the operating room. Therefore, the local pollutant exposure index and spatial concentration were used to quantitatively analyse the spatial concentration distribution and changes in particles near patients and medical staff in the operating area under different ACH values for the three ventilation patterns. The removal efficiency of particulate matter in the operating room was measured and compared. Our results can guide the prevention of surgical smoke and effectively reduce the harm caused by smoke particles to the surgical staff.

2. Experiment Settings and Methodology

2.1. Physical Model

As shown in Figure 1(d), a standard class I operating room environmental test chamber has dimensions of (). The operating room has three air supply forms: ceiling air supply, side-wall air supply, and slit air supply. In the middle of the ceiling, there were eight air outlets with a size of  m and a HEPA filter. Eight air supply outlets with a size of  m were present around the ceiling, and an entire air supply outlet with a size of  m was located on the side wall opposite the door. The left and right sides of the room were, respectively, equipped with four exhaust vents with a size of  m, and the lower edge of each exhaust outlet was 0.15 m away from the ground.

The basic facilities in the operating room adhere to the standard of a typical level I operating room, including a single-arm mechanical surgical tower, an operating table, two operating lights, an instrument table, electric surgical equipment, and some common equipment. Additionally, the operating room staff comprised one patient and four medical staff members, as shown in Figure 1. This arrangement is equivalent to that commonly used in abdominal or visceral surgeries. The height of the mannequin was 1.75 m, and the heat load of the mannequin was adapted to the level of individual activity.

2.2. Case Study and Instrumentation

The layout of the operating room equipped with the VLAF, HLAF, and UDAF + mixing ventilation systems is shown in Figures 1(a)1(c), and the specific design parameters of each ventilation system are listed in Table 1. The experimental equipment and parameters are presented in Table 2. Notably, the ACH in the UDAF + mixing system room was 60; however, the ACH in the surgical area was higher than that in the surrounding area.

2.3. Experimental Method

Studies using pig tissue models to assess tissue damage during electrosurgical procedures have found no differences between pig tissue models and human tissue samples [11, 34]. Therefore, in this study, pork tissue was used instead of human tissue for high-frequency electric knife-cutting experiments. To simulate most surgical situations, we used the cutting mode of monopole electrotomy, which is most commonly used in surgery, to study the distribution of smoke particles in the operating room. Figure 2 shows the layout of the equipment and the generation of smoke particles in the experiment.

Three particle monitoring points were set up in the operating room: S1, S2, and S3, as shown in Figure 3. S1 was placed 1.65 m from the floor to measure the change in the particulate matter concentration in the respiratory area of the doctor. Similarly, S2 was placed 0.8 m from the floor to measure the change in particulate matter concentration in the respiratory area of the patient. Point S3 is located near the outlet, 0.2 m above the ground. The particle concentration per second at each measurement point was determined using a handheld PGM-300 particle counter.

To reduce the influence of background concentration on the measurement results, the ventilation system of the operating room was opened for at least 2 h before each experiment until the concentration at each measurement point was stable. The prepared fresh pork was placed on the patient’s abdomen, which was slightly close to the attending physician. The patient lay flat on the operating table, 0.8 m above the ground. The particle counters were arranged at the smoke detection points, as shown in Figure 3. The simulation experiment was conducted using electric surgical equipment with an electric cutting power of 45 W. To measure the background concentration and eliminate environmental interference, the particle counter was turned on 5 min before the formal experiment began and then turned off 5 min after the experiment ended. The total duration of the formal surgical experiment was 15 min. The first 10 min was the surgical cutting time, and the last 5 min was the evacuation time of the surgical smoke. During the experiment, the electric knife was used to cut the tissue in a uniform linear manner. A TSI330 aerodynamic particle analyser was used to sample and record the source concentration and particle size distribution of the surgical smoke for subsequent analysis.

3. Computational Fluid Dynamic Methodology

3.1. Airflow Phase Simulation

In this study, the renormalization group (RNG) k-ε turbulence model in the analysis of systems (ANSYS) FLUENT 2020R2 was used to simulate the fluid flow. According to previous studies, the RNG k-ε model is suitable for indoor airflow simulation [3537]. The coupling of the pressure and velocity fields was realised using the Semi-Implicit Method for Pressure Linked Equations-Consistent (SIMPLEC) algorithm. The second-order upwind discrete scheme was used for the momentum, energy, turbulent kinetic energy (k), and dissipation rate of the turbulent kinetic energy (ε). The time discrete analysis of the transient term was conducted using the second-order implicit scheme, which ensures accuracy and improves calculation speed. An enhanced wall treatment was applied near the turbulent wall. In terms of the boundary conditions, the velocity inlet was adopted, and the temperature was 294.15 K. The outlet is the pressure outlet, whose pressure value is 10 Pa. The heating load of the room was primarily supplied by air supply at a certain temperature (21°C). The surface area of the medical staff was 1.49 m2. The patient was set at a constant temperature of 27.4°C, and the heat flux of the medical staff was 115 W/m2. The heat fluxes of the instrument table, surgical tower, and surgical lamp were set at 150, 200, and 200 W/m2. The remaining surfaces were set as adiabatic. The physical process of indoor radiation was not simulated in this study.

The RNG k-ε model is derived using a statistical method called RNG theory and is expressed as follows: where and , k is the turbulent kinetic energy per unit mass, and is the turbulent dissipation rate.

3.2. Particle Phase Simulation

The discrete phase model (DPM) was used to predict particle transport behavior in the air, which was essentially the simulation of continuous gas phase and dispersed particle two-phase flow in the air. This model used one-way coupling to calculate particle transport and was widely used to simulate particle dispersion in the air [38]. Firstly, the stable flow field is calculated, then introduced the discrete particle phase, tracked, and drew its motion trajectory. At the same time, the discrete random walk model was applied to simulate the dispersion of particles due to turbulence [29, 39]. Smoke particles were simulated as spherical particles with a density of μg/m3, a diameter of 0.3 microns, and a mass velocity of μg/s. The mass velocity was converted from the experimental measurements. The data measured by the TSI 3330 optical particle size analyser showed that more than 70% of the smoke particles had a particle size of approximately 0.3 microns, as shown in Figure 4. Notably, electrosurgical instruments can produce particles as small as 0.3 microns [3, 4, 40]. Therefore, the particle size in the simulation was set to 0.3 microns. In this study, the steady-state flow field of the three ventilation patterns was first calculated, and the transient calculation was then conducted by setting the particle phase. The unsteady state simulation calculation is 900 s in total, with particle release for the first 600 s. In order to ensure a stable state, we performed a comparative calculation of releasing 2000 particles, 2500 particles, and 3000 particles per second. Results show that the number of particle tracks of cannot reach the statistical result of convergence. The statistical effect of convergence can be achieved when 2,500 particles and 3,000 particles are released per second. Finally, we chose to release 2500 particles per second for subsequent research. It is closer to the experimental data and the particles can fill the room well. At the same time, the particle trace number of 2500 × 600 is enough to achieve the statistical result of convergence. The time step is set to 1 s to fully track the movement path of smoke particles. where is the particle velocity; is the air velocity; is the resistance per particle mass; and are the densities of the particles and air, respectively; is the additional force; and is the gravitational acceleration vector.

In calculating the motion of particles, additional forces depend on the properties of airflow and particles. The pressure gradient, Basset, and virtual mass forces were ignored because of the small ratio of air density to particle density [41, 42]. In this study, the Stokes-Cunningham drag law was applied to calculate the drag force. The additional forces included the Brownian, Saffman, and thermophoretic forces. We assumed that the particles have left the area when they reached the vent; therefore, we set the vent to the escape condition. Moreover, when particles encounter hard surfaces and walls, they generally do not possess enough rebound energy to resist adhesion; therefore, the surfaces and walls in the operating room are set as trap conditions. Owing to the constant temperature and humidity in the operating room, there was no significant difference in the temperature of the return air; therefore, the influence of the environmental temperature difference on airflow and particulate matter was ignored.

To compare the pollutant distribution in different regions, we introduce the local relative pollutant exposure index () [43]. This value was also used in similar studies and is defined as where is the average concentration of particulate matter at a certain point in the operating room, is the average concentration of particulate matter at the exhaust port, and is the concentration of particulate matter in the intake air. When is the indoor average concentration, is the spatial relative pollutant exposure index.

3.3. Numerical Model Grid Independence

The integrated computer-aided engineering and manufacturing (ICEM) computational fluid dynamics (CFD) ANSYS software was used to divide the physical model through structured grids of different sizes, and three types of grids, namely, 2569882, 3552736, and 5098912, were obtained. For the three different meshes, mesh sizes of 4, 2, and 2 cm were set for the surgical staff surfaces. The mesh sizes of the supply diffuser and exhaust were 5 cm. The mesh sizes at the surface of the equipment were 4, 2, and 2 cm. The global sizes were 12, 10, and 8 cm, respectively. The grid independence was verified by measuring the air velocities at points on two lines near the rear end of the operating table. The results are shown in Figure 5. The variation trends of the three grid speeds were consistent, and the difference was small, indicating that the independence of the grid was acceptable. Additionally, the sensitivity of the mesh size was tested by comparing the simulation results for different mesh numbers. The grid convergence index (GCI) [44] was calculated using Equation (4) to evaluate the quality of the selected grid. The results showed that all the GCI values were less than 6.2% for the selected mesh sizes, indicating that the mesh was adequate. Considering the complexity and accuracy of subsequent calculations, grid 3552736 was selected for subsequent calculations. where is the linear grid refinement factor, is the formal order of accuracy (second-order), and is the safety factor, as suggested by previous studies [44, 45], and is the difference of the variables between the described mesh and a refined mesh.

4. Results and Discussion

4.1. Numerical Model Validation

The airflow field in the simulation was verified to ensure the accuracy of numerical calculations. Figure 3 shows the positions of the six guides at the two heights. Figure 6 shows a comparison between the simulated velocity values at each point of the six lines and those at the experimental measurement points. The results showed that the simulated values were consistent with the experimental values. Owing to the error of the anemometer in the experiment and the simplification of the model in the numerical simulation, the difference between the experimental and simulated values was reasonable, and the simulation results were acceptable.

The simulation and experimental values of the smoke particle concentration were verified, and the results are shown in Figure 7. For the experimental data, we averaged the values of three repeated experiments at each sampling point at the determined time, subsequently averaged the values obtained at three sampling points, and standardised the concentration value at each sampling point with this value. The error value was then obtained by calculating the standard deviation. To process the simulation values, the smoke particle concentrations at the three sampling points at the determined time were averaged, and the simulated concentration values at each sampling point were standardised using this value. Although the experimental result is not completely consistent with the simulated value, considering the complexity of the experimental conditions and the error of the measuring instrument, the deviation does not affect the overall trend of smoke particle concentration at different sampling points at different times. Therefore, the result is acceptable, and the accuracy of the subsequent simulation can be guaranteed.

4.2. Spatial and Temporal Distribution of Airflow Field and Particles

At the designed air volume, the airflow fields of the three ventilation patterns and the particle positions at different times are shown in Figure 8. Figure 8(a) shows the velocity vector diagram, flow diagram, and the diffusion of smoke particles of and sections in the operating room under the VLAF system. The flow from the diffuser in the centre of the ceiling was accelerated by gravity and flowed vertically downward through the operating area. However, the high-speed airflow was forced to change direction as it was obstructed by the operating lamp, forming eddy currents under the operating lamp and resulting in no vertical airflow in front of the attending physician and assistant. As air flows over the staff and patients, it slows down near the body because of the body’s heat plume. Part of the airflow is deflected to both sides at a large angle after encountering the operating table, and a small part of the airflow forms a small eddy current under the operating table. When the diffused air hits the wall, it begins to move upward and forms a large eddy outside the surgical area with a lower velocity at its centre. As for the smoke particles, as the airflow is blocked by the operating lamp and its own thermal buoyancy, the smoke particles first move to the top of the patient’s head after generation, and most of the particles then quickly move down to the + and - corners under the influence of the vertical airflow. At this point, a small number of particles remained in the vortex region in the corner, and most of the particles changed direction and moved up the wall. When the particles reached the ceiling, they began to move diagonally downward and gradually diffused.

The airflow and particulate matter in the operating room equipped with the HLAF ventilation system are shown in Figure 8(b). The airflow begins from a wall diffuser on one side and blows horizontally across the surgical area. The airflow over the operating lamp is blocked, changing direction and forming a low-speed eddy zone behind the operating lamp. After the airflow reaches the wall through the surgical area, it is blocked from flowing into the ground. Finally, it rises near the surgical area, forming a large eddy current in the space. The airflow distribution on both sides of the operation area was chaotic. Except for some airflow to the exhaust port, most of the airflow was towards the roof and formed a low-speed vortex near the corner. Immediately, after the smoke is produced, it follows the horizontal airflow towards the leg of the patient. Subsequently, a large amount of deposition occurred along the direction of the vortex moving towards the wall. Simultaneously, some particles returned to the space with the airflow and gradually diffused.

The airflow and particulate matter in the operating room equipped with the UDAF + mixing ventilation system are shown in Figure 8(c). In this system, the airflow pattern in the operation area is similar to that of VLAF ventilation, but the vortex under the operating table is significant, which may be ascribed to the interference of the airflow outside the operation area. Moreover, external airflow affects the vertical airflow in the surgical area, forcing it to change its direction in advance. However, a relatively high-velocity horizontal airflow formed near the patient, which ensured the removal of smoke particles in the surgical area. Compared to that of traditional VLAF ventilation, the space outside the surgical area in the UDAF + mixing system produces a larger eddy area with a smaller but more stable airflow velocity. Owing to the obstruction of the operating lamp, the smoke particles first rose upward under the action of thermal buoyancy, and some of the particles remained in the eddy current near the operating lamp. With an increase in the number of particles, most particles moved towards the exhaust vents on both sides with the vertical airflow. However, owing to the obstruction from medical staff and equipment along the way, the direction of movement of the particles was deflected, and some particles diffused into space or were trapped in the eddy current area.

4.3. Distribution of Smoke Particles in the Operating Area

Additionally, a cube with dimensions of  m () was constructed in the entire area with the medical staff and patients, from the patient to the head of the medical staff, representing the working area. In the physical model, the length, width, and height coordinate ranges of this region are to 1.2 m, to 1 m, and to 1.8 m, respectively.

To further study the concentration distribution of particulate matter in the surgical area, we selected three areas on the upper side of the patient and six points in front of the medical staff (Figure 9). The six points were located on the same plane, near the mouth and nose of the medical workers, to monitor changes in the particle concentration in their respiratory areas. The four points P1, P2, P5, and P6 were located 0.2 m in front of the faces of the medical staff, and the distance between the two points on each side was 0.9 m. We used five different ACH values (40, 50, 60, 70, and 80) to study and evaluate the three different ventilation patterns.

Figure 10(a) shows that in the operating room equipped with the VLAF ventilation system, the area with the highest local relative pollutant exposure index near the patient was always located in the patient’s respiratory area, wherein the highest exposure index reached 15.73. This may be attributed to the obstruction by the operating lamp, which causes the accumulation of smoke particles in the patient’s breathing area. With an increase in the ACH value, the of each region changes significantly. When the ACH value was less than 60, the in each region increased with increasing ACH. The of each region decreased rapidly after the ACH value exceeded 60. The indoor relative pollutant exposure coefficient reached the lowest value of 0.68 when the ACH was 70. Beyond an ACH value of 70, did not decrease significantly. This may be because the vertical downward laminar airflow is blocked by obstacles midway, such as the surgical lamp and medical staff. This destroys the clearance effect of vertical airflow on particulate matter and generates many small vortices near the obstacles, thereby hindering the removal of accumulated particulate matter.

In the vicinity of the medical staff, as shown in Figure 10(b), the concentration at point P1 near the assistant’s face showed the most obvious range of change accompanied by a change in the ACH value and was much higher than that at the other points. When , the highest and lowest concentrations at P1 were and when the ACH values were 50 and 70, respectively. When the time reached 720 s (the end of surgery was after 120 s), the concentration at each point decreased to less than ; at this time, if the ACH is more than 60, the concentration at each point can be lower than . Notably, the particle concentration at point P2 near the attending doctor was the highest at this time, which may be because the airflow here was affected by the operating lamp, resulting in an eddy current, which led to the retention of particles that were difficult to eliminate.

Figure 11(a) shows that in the operating room equipped with a HLAF ventilation system, when the ACH value is less than 60, HLAF cannot remove smoke particles in the operating area, and the near the patient was up to 28.26. After more than 60 breaths, the horizontal airflow carries particles from the patient’s respiratory and wound areas to the patient’s legs. At this point, the maximum value in the leg zone was 12.27. When the ACH value increased to 70, the HLAF system had a washing effect on the entire surgical area; however, at this time, relatively high concentrations of particulate matter were observed in the leg and wound areas of the patient. There may be two reasons for the above: one is due to obstruction of laminar airflow by the patient and surgical lights that can generate the small changes and small eddy currents in the vicinity of a patient’s body resulting in particulate retention; the other is that the airflow hitting the wall forms a large vortex in the space, which affects the leg area of the patient.

In the vicinity of the medical staff, as shown in Figure 11(b), when the ACH value is less than 70, the particulate matter concentrations at P4, P5, and P6 are very high in the first 200 s due to the obstruction of horizontal airflow by the medical equipment. The smoke particles moved in the operating room with the air currents. At 200 s, particulate matter accumulated in large quantities at point P2 near the wound, resulting in a relatively high concentration. Driven by the horizontal airflow, particulate matter continuously moves towards the patient’s legs, resulting in an increasing concentration at point P4. Over time, more particles hit the wall with airflow and spread through the room, gathering near the assistant. Thus, at 600 s, the concentration at P1 is relatively high. At 720 s, most of the particles were discharged from the air with the airflow, but near P5, the obstruction of airflow by the equipment resulted in the accumulation of a large number of particles, which was difficult to remove. Therefore, the concentration of the particles at this point was relatively high.

As shown in Figure 12(a), compared to that with the other two systems, the operating room equipped with the UDAF + mixing system has little difference in in each area, and is less affected by the ACH, and the overall particulate matter level is more stable. This may be due to the reduced area of the vertical laminar air supply, more significant obstruction of the vertical airflow by the operating lamp, and the eddy current generated above the patient’s thermal plume, resulting in a uniform particle distribution. Compared to those in the other two conditions, the of each area above the patient in the UDAF + mixing system was below 5, the indoor relative pollutant exposure coefficient was always below 0.2, and the overall stability was at a lower level without a high peak. All the values were less than 1 when the ACH value was 60 and 400 s. This may be ascribed to the fact that the airflow provided by the surrounding air supply outlet not only dilutes the particle concentration in the surrounding area but also strengthens the airflow to the exhaust port and enhances the exclusion effect of particles in the surgical area and room. This is also observed in the values of each region after 720 s; at this time, the values of the region above the patient in the operating room equipped with the UDAF + mixing system basically tend to 0, lower than those of the other two conditions.

Figure 12(b) shows the respiratory surface of medical staff. The concentrations at P3 and P6 are always close to 0, and those at other points change significantly with changes in the ACH value; however, the concentrations at all the points are in a low range. When  s, the concentrations at point P2 near the attending doctor reached the highest and lowest values of and μg/m3 when the ACH values were 60 and 70, respectively. The remaining points showed a tendency to increase with the ACH. This may be because the airflow provided by the surrounding air supply port enhances the airflow to the exhaust port, which improves the removal effect of indoor particulate matter. Airflow carries more particles through medical staff, resulting in a large amount of particulate matter diffusion and retention. When the time reached 720 s, the particle concentration at all points in the operation area, except P2, decreased to μg/m3, which is lower than that of the other two ventilation patterns, indicating that the UDAF + mixing system is more effective in removing particulates in the operation area. At this time, the concentration at point P2 near the attending doctor is relatively high, which may be caused by the blocking effect of the surgical lamp on the vertical airflow, similar to that in the VLAF system.

In the VLAF system, the particulate matter concentration in the patient’s respiratory and wound areas should be noted during surgery, and local smoke exhaust devices should be strengthened in this area. The attending physician and assistant should pay attention to the effectiveness of protective measures to avoid harm caused by smoke particles. For the HLAF system, an adequate air volume should first be ensured in the operating room, followed by smoke removal measures below the patient’s wound. Medical staff near the patient’s leg should strengthen the protective measures. The peripheral areas on both sides of the exhaust vent have chaotic flow, and the residence time of particles is long; therefore, medical staff should avoid staying there for long periods. The UDAF + mixing system in the operating room had a low overall concentration distribution in the operating area during operation. The main focus should be on the protection from smoke particles at the wound site, and the medical staff on the path from the wound to the exhaust outlet should pay special attention to these particles. The HLAF system is effective in removing smoke particles in the surgical area when a sufficient air volume is maintained. The UDAF + mixing system is the most stable in the surgical area compared to the other two systems, which maintains the concentration in the area at a lower range and is more effective in removing smoke particles.

4.4. Removal Efficiency of Smoke Particles under Different Ventilation Patterns

Figure 13 shows the removal rate of the operating room under the three ventilation patterns and different ACH values. The removal rate here was defined as the ratio of the total emission of particulate matter in the operating room to the total production at the end of the operation [46]. Notably, large ACH values do not represent a better removal effect. The VLAF and UDAF + mixing systems had the highest removal rates of 58.76% and 71.39%, respectively, at an ACH value of 60. When the ACH value was 70, the removal rate of the HLAF system was the best, with a value of 61.54%. The UDAF + mixing system had a better removal efficiency than the other two systems.

5. Conclusions

In this study, experiments and CFD simulations were used to explore the spatial concentration distribution and variation of smoke particles under VLAF, HLAF, and UDAF + mixing ventilation systems in an operating room. The main conclusions are as follows: (1)For the VLAF system, the concentration of smoke particles in the operating area was high because of the blockage of the operating lights. Increasing the ACH did not significantly reduce the concentration of smoke particles. The HLAF system was significantly affected by the ACH, and once the ACH was below 70, the indoor air cleanliness deteriorated significantly(2)The UDAF + mixing system was a more effective ventilation pattern, and the optimal ACH value was approximately 60. At this time, the local pollutant exposure index near patients in the operating room was very small, the particle concentration near the medical staff was less than μg/m3, and the particle concentration in the operating area was stable at a low level(3)Increasing the ACH did not achieve a better removal effect. Compared to that in the other two cases, the airflow distribution of the UDAF + mixing system was more conducive to particle removal, and the highest removal rate was 71.39%. However, in this case, the medical staff on the path from the wound to the exhaust outlet should ensure protection from particles

Considering these results, this study shows that the UDAF + mixing system is an alternative to laminar ventilation systems. The system can stably control the concentration of smoke particles in the operating area within a low range while ensuring a high indoor particle removal efficiency. Moreover, the layout of the operating room equipment and airflow distribution has a significant influence on the distribution and removal of particles. A study on the optimal equipment layout scheme and airflow distribution optimisation for each ventilation pattern will be conducted in future studies.

Data Availability

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

Additional Points

Practical Implications. (i) Surgery smoke particle concentration is not to be ignored, and ventilation pattern, air change rate, and indoor layout seriously affect the indoor clean environment. (ii) The migration, distribution, and removal of smoke particle in the operating room under different ventilation patterns were quantitatively analysed, and the high concentration area was determined. The research results effectively support the prevention and control of smoke particles. (iii) Central unidirectional airflow and peripheral mixing are effective ventilation.

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.

Authors’ Contributions

Zhijian Liu conducted the conceptualization, methodology, funding acquisition, supervision, project administration, and writing of the review and editing. Zheng Zhang conducted the methodology, software, validation, formal analysis, data curation, and writing of the original draft. Jiabin Lv conducted the project administration, resources, and writing of the review and editing. Juntao Ma managed the software and investigation. Guangpeng Yao conducted the writing of the review and editing, software, and formal analysis. Junzhou He conducted the writing of the review and editing. Guoqing Cao managed the resources and supervision.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (nos. 42122058 and 41977368), the Natural Science Foundation of Hebei Province (no. E2021502046), the National Key R&D Program of China (2021YFF0604000), and the Fundamental Research Funds for the Central Universities (nos. 2021MS075 and 2020YJ007).