Comparing Design Schemes and Infection Risk Assessment of Negative Pressure Isolation Cabin

: At present, various public health emergencies have forced a deeper study of measures to prevent infectious diseases. To prevent the spread of infectious diseases on large cruise ships, the use of negative pressure isolation cabins is an eﬀective method. However, existing cruise ships rarely use negative pressure isolation cabins or use them with shortcomings such as unreasonable layout of air inlets, which do not eﬀectively reduce the risk of infection for medical staﬀ while isolating patients. To solve this problem, ﬁrst, the method of active air interference and the principle of proximity of the fresh air outlet were examined in this paper, and four groups of optimization schemes were designed. Second, by analyzing the diﬀusion of droplets in three breathing modes (coughing, sneezing, and talking while wearing or not wearing a mask), the direction of droplets and the eﬃciency of pollutant discharge under the condition of coughing were compared. Finally, in this paper, the infection risk of contact transmission and aerosol transmission was optimized by using the linear quantitative evaluation method and MSDR method, respectively. The results showed that the auxiliary air intake layout of optimal scheme 2 can eﬀectively reduce the risk of infection for medical personnel in negative pressure isolation cabins. This study provides a useful reference for the design and optimization of negative pressure isolation cabins in future cruise ships.


Introduction
In recent years, the development of the transportation industry has made the global movement of people and goods more convenient and frequent, and ocean cruise tourism has become one of the fastest growing travel modes. Sea travel, including cruise tourism, has increased the spread of various pathogens, particularly the spread of infectious diseases [1,2]. Studies have shown that respiratory infections on cruise ships are the most common cause of travel and crew visits [3]. Particles larger than 5 µm expelled by patients while speaking, coughing, or sneezing are generally considered to be droplet transmission, and particles smaller than 5 µm are classified as airborne transmission. The airborne transmission of viruses refers to the spread of infectious pathogens caused by droplets (aerosols) suspended in the air, allowing them to maintain infectivity over long distances and for a long time [4]. Recent outbreaks of SARS, MERS, and COVID-19 were caused by the spread of droplet aerosols, and prior to the outbreak of COVID-19, such infectious diseases frequently occurred on cruise ships, causing serious public health security incidents [5,6].
Due to the large capacity of passengers, the small internal space, the poor effect of ventilation, intensive personnel activities, and the majority of passengers being elderly, once infectious diseases occur on cruise ships, they can easily cause the spread of cascading epidemics [7,8]. For example, in the summer of 2009, the outbreak of influenza A(H1N1) on the large luxury cruise ship "Pacific Dawn" and the COVID-19 outbreak on the Diamond Princess in 2020 were both caused by inadequate protection and the failure to take timely quarantine measures [9,10]. In the past three years, the Center for Disease Control and Prevention in Morgantown, West Virginia, has used a self-designed cough aerosol simulator to analyze the effect of patients wearing face masks on reducing the excretion of COVID-19 virus aerosols during coughing and has designed the corresponding experimental devices for verification [11]. In addition to personal precautions such as wearing masks, negative pressure isolation wards are the most effective facilities to control the source of infection and cut off the transmission routes. Diseases transmi ed on cruise ships are mainly respiratory infections caused by droplets or aerosols, and every cruise ship should be equipped with medical isolation facilities [12].
The negative pressure isolation cabin has good sealing. When there is an epidemic of infectious diseases on the cruise ship, confirmed patients can be quickly isolated to prevent the further spread of the epidemic and ensure the safety of the cruise ship personnel. On 14 September 2014, the newly constructed negative pressure isolation room was officially opened at the home port of China's Tianjin Cruise Line. This negative pressure isolation room has functions such as rapid screening, X-ray detection, quarantine detention, and sample processing, which can minimize the spread of infectious diseases [13]. With the rapid development of computer hardware and software, more and more researchers are now using computational fluid dynamics (CFD) methods to analyze biological pollutants [14][15][16][17]. After the outbreak of SARS in 2003, many nonvirological and biological experts and scholars adopted numerical simulation methods to simulate and analyze the organization of viral transmission airflow and evaluate regional airflow [18][19][20][21][22][23][24].
The existing negative pressure isolation cabin was designed by adding a negative pressure ventilation system to an ordinary cabin. However, the complexity of the flow field in the negative pressure isolation cabin was not considered in the arrangement of the air inlets, resulting in the risk of infection still existing during the process of treating patients. To solve the above issues, the following contributions are made in this paper by using numerical simulations. (1) Based on the active interference airflow method, combined with the fresh air outlet proximity principle, four optimization schemes are proposed to reduce the infection risk of medical staff. (2) The droplet diffusion in the negative pressure isolation cabin of a large cruise ship was compared and analyzed under three conditions: coughing, talking, and sneezing, and the effectiveness of wearing masks was theoretically analyzed. (3) The risk of infection was assessed from both contact transmission and aerosol transmission. Finally, the feasibility of the optimization scheme is verified.
This paper comprises six sections. Section 2 describes the numerical simulation and boundary conditions. Section 3 compares and analyzes the effectiveness of wearing a mask under three working conditions. In Section 4, the optimized model is compared and analyzed. In Section 5, the risk of infection is assessed. Section 6 gives the conclusion of this paper.

Building of Optimization Models
Due to the limited space available on cruise ships, negative pressure isolation cabins are designed to take up less space than negative pressure wards. Therefore, the negative pressure isolation cabin of the cruise ship is designed and analyzed by reference to the compartment of negative pressure ambulances. This design not only saves space but also achieves negative pressure isolation. The isolation cabin is basically the same structure as the ambulance cabin, including basic medical equipment, stretchers, seats, etc. Not only are two individuals, a caregiver and a patient accommodated, but the patient must be kept in isolation until the ship is docked without developing any other symptoms. This is important for older people with underlying medical conditions, as most cruise passengers are elderly. A large amount of medical equipment is arranged on the right side of the patient, including oxygen cylinders, ECG monitors, first aid kits, infusion hangers, and more. Moreover, after a reasonable layout, it has li le effect on the fluid and can be reasonably simplified.
Depending on the difficulty of the model and the characteristics of surface modeling, SolidWorks was selected to complete solid modeling, and the model was modified in combination with the modeling software Design Modeler of Ansys Workbench. The mesh division method is the Fluent Meshing mixed mesh division method, the surface mesh is a polyhedron, including boundary layer division, and the body mesh is a hexahedral mesh structure. Due to the complex internal structure of the negative pressure isolation cabin, a reasonable simplification was made during modeling. Considering that the structure of the mannequin has a certain influence on the local fluid, a curved solid model with a height of 170 mm was selected without gender being ignored. Regardless of the effect of the holes and grilles on the fluid, the air inlet adopts an approximate rectangular port, which is divided into groups of two; the size of a single port is 140 mm × 60 mm, the interval within the group is 80 mm, and the interval between the groups is 400 mm. To generate the wall backflow, the influence of the grating on the fluid is considered in the exhaust outlet. The single mouth size is 55 mm × 5 mm matrix grid, and the patient's mouth is elliptical: the long axis is 40 mm and the short axis is 12 mm. The fluid simulation analysis is completed by Fluent 2020R2. The established fluid domain model is shown in Figure 1. Based on the characteristics of air flow and active interference methods, this paper proposes increasing the auxiliary inlet of the fresh air system to protect the individual, so that the breathing area of the individual is always in a clean, fresh air environment. Air intake above or behind the head of the individual and the two sets of seats form four combinations. The four schemes are air intake above the head of individual, air intake behind the head of individual, air intake above the head of caregiver + air intake behind the head of individual, and air intake behind the head of caregiver + air intake above the head of individual, as shown in Figure 2.

Basic Control Equations
(1) Conservation of mass equations Any flow problem satisfies the conservation law: the increase in mass per unit time in a fluid cell is the same as the net mass flowing into the cell in the same time interval [25]. Following this law, the general form of the mass conservation equation shown in Equation (1) can be obtained: where ρ is the density of the fluid; t is time; u is the velocity vector of the fluid; Sm is the source term.
(2) Momentum conservation equation The flow of any flow system complies with the law of momentum conservation. According to Newton's second law, the external cooperative force on a fluid element is the rate of change of its momentum with time [25].
From this, the expression of the momentum conservation equation that can be obtained in the three directions of x, y, and z is: where p represents the pressure on the fluid microbody; τij indicates surface stress; u, v, and w represent the velocity in the three directions of x, y, and z, respectively; Fx, Fy, and Fz refer to the physical strength on the microbody.

(3) Equation of energy conservation
When heat exchange is included in a fluid system, it must comply with the equation of energy conservation. This law is the first law of thermodynamics: the work done by the fluid microelement by physical or surface force and its own net heat flux is equal to the increase in the internal energy of the microelement. The equation is expressed as follows: where T represents temperature; ST represents the internal heat source of the fluid and the part of the mechanical energy of the fluid that is converted into heat energy due to viscous action; cp indicates the specific heat capacity; k represents the heat transfer coefficient of the fluid.

Boundary Conditions
In this study, the DPM model is mainly used to simulate the movement of droplet aerosol particles generated by respiratory methods such as coughing by patients, and the simulation is a transient flow field. However, a steady-state result should be obtained before analysis of the transient flow field. The parameter of the steady-state calculation is shown in Table 1. When using Fluent software for the simulation calculation, the same boundary conditions need to be added for each working condition as follows: the air inlet of the isolation cabin of a large cruise ship is the inlet boundary condition. The number of air changes in the isolation cabin of the cruise ship is set to 20 times /h, the volume of the fluid domain measured by Fluent is 7.208255 m 3 , the size of the four air intake port of the air conditioner is 140 mm × 60 mm, and the minimum inlet wind speed can be obtained, and the inlet wind speed v = 1.5 m/s is finally determined. In the case of low temperature influence, the temperature is set at 21 °C combined with the air conditioning gear and human comfort. Hydraulic diameter D = 84 mm. Turbulence intensity I = 5.20%. At the same time, the temperature of the medical staff is taken into account, the temperature parameter for all human bodies in the isolation cabin is set to 28 °C, the patient's mouth temperature at 37 °C, and other default parameters can be selected. Set the outlet to the pressure outlet. The total pressure parameter is set to −30 Pa, and the direction is the vertical outlet surface.
The DPM model is used to simulate the movement of aerosol particles produced by respiratory pa erns such as coughing. Common respiratory behaviors such as sneezing, coughing, and talking are the main ways of virus transmission and the main source of aerosol particles from droplets. According to the data of droplet particles generated by different breathing modes, the fixed values within each parameter interval are selected and the Rosin-Rammler [26] distribution is used for simulation. In the numerical calculation, two conditions are considered: without a mask and wearing a mask, and the particle temperature is lowered to 30 °C when wearing a mask. The particle velocity is automatically corrected by the resistance of DEFINE_DPM_DRAG macro in the unmasked condition. Other important parameters related to droplets are shown in Table 2 [27][28][29][30][31]. In order to ensure the effectiveness of the simulation of wearing masks, a doublelayer non-woven mask was selected for calculation [28]. The fiber diameter of the mask material is 227 µm, the layer spacing is 1 mm, and the fiber spacing is about 350 µm. In the Fluent DPM model, to refine the initial ejection of droplet aerosol particles, the UDF macro DEFINE_DPM_INJECTION_INIT is used to specify the change law of incident source position with time, including the change law of the initial state of the particle (such as particle size, speed, flow, temperature, etc.). When simulating wearing a mask, resistance is loaded with functional formulas through the macro DEFINE_DPM_DRAG. The air flow ejected by the patient is entered directly using the if condition function. In the Fluent simulation calculation, to be er reflect the relationship between the discrete phase and the watershed model boundary, the four boundary conditions of patients, caregiver, individuals, and ground are set to Trap, the exit is set to Escape, and other boundaries such as walls are set to Reflect.

Grid Independence Verification
The data obtained in this paper comes from Fluent simulation software, so the validity of the data needs to be further studied. Tsz-Wun Tsang et al. [32] mentioned three common ways to validate a model: (1) using experimental or measurement results obtained in a real simulated environment; (2) making use of experimental or measurement results in the literature; (3) adopting the benchmark model. In this study, the method of grid independence verification [33] is used to illustrate the accuracy of the CFD model to ensure that the mesh size is reasonable and does not affect the calculation accuracy.
In order to ensure that there are no gaps, holes, or other problems in the grid, the free plane, multiple planes, overlapping plane, and skew rate (>0.85) are all 0, indicating that the model has a good structure. Then, according to the size of the special plane and the use occasion, the opposite grid is refined, the specific size is shown in Table 3.

Parameter
Value The volume mesh size was selected as 100 mm, 80 mm, 60 mm, 50 mm, 40 mm, 30 mm, 25 mm, 20 mm, and 15 mm, respectively. In order to be er explain the influence of the volume mesh size on the data results, the speed and pressure of the reference plane V plane at 15 mm above the patient were used for data analysis, as shown in Figure 3. After reasonable se ing, simple boundary conditions and methods were used to obtain the velocity and pressure values of the reference plane V plane, and the results of grid independence verification were shown in Figure 4.  According to the result data of grid independence verification, combined with the selection principle of small data fluctuation and large mesh size, the volume mesh size is finally selected as 40 mm.

Cough
As shown in Figure 5, at 10 s, the droplets have moved to the rear of the isolation cabin. They move down from the top of the back of the isolation cabin under the influence of airflow and are mainly distributed around the patient near the medical equipment side. At this time, there are a small number of droplets in front of the caregiver and a very small number of droplets around the individual. At 120 s, the number of droplets suspended in the entire fluid domain is significantly reduced. Although more than half of the droplets have been exhausted or se led on the ground through the exhaust port, some droplets are evenly distributed throughout the fluid domain, indicating that the individuals and caregiver are still in a dangerous environment at this time. (2) Wearing a mask As shown in Figure 6, at 10 s, a large number of droplets at the exit have been discharged. Some droplets are still in the drainage basin above the head of the patient and tend to discharge from the outlet. Nearly half of the droplets spread out on the right side of the caregiver and are evenly dispersed around the corner of the isolation cabin. A small fraction of droplets have reached the chest cavity of the caregiver, and the fastest-moving particles have reached the top of the caregiver's head. At 120 s, the droplets have not been significantly reduced, and there is a certain amount of particles around the caregiver and individuals. This indicates that droplets cannot be discharged quickly in a short time and the medical staff in the cabin are still at risk of infection.

Distribution of Droplets of Different Sizes
The diameter range of droplets produced by coughing is 1~100 µm, and the main particle size value is 90 µm.
(1) No masks As shown in Figure 7, the droplets at 1 s are still moving due to the force of the airflow generated by coughing, but due to the small airflow velocity near the patient's mouth and nose, and the large-size droplets of 63.2~90.1 µm are greatly affected by gravity, so no large displacement is produced. The droplets at 10 s are diffused in the watershed, and the large droplets with a large particle size of 16.4~27.1 µm in the remaining particles are mainly distributed in the front of the isolation cabin, and most of them are located at the exit and on the ground. (2) Wearing a mask As shown in Figure 8, the droplets at 1 s have an obvious stratification phenomenon, the larger droplets of 63.2~90.1 µm slide from the patient's face to the outlet direction, and the small droplets of small particle size of 0.271~27.2 µm float in the watershed. The large droplets at 10 s are basically in a stable state, and the main distribution in the watershed is the small droplets, and their diameter is less than 20 µm.

Discrete Droplets Particle Residence Time
(1) No masks As shown in Figure 9, at 10 s, droplets are distributed throughout the watershed, some of which are discharged under the action of the flow field in the cabin, and some particles finally se le on the ground of the isolation cabin by gravity. Under the action of the mainstream, the droplets in the watershed move mainly in the area above the patient, near the wall of the isolation cabin, and the number of droplets near the medical staff is relatively small. At 120 s, most of the droplets have been discharged or se led on the ground, but a small number of particles are still floating in the basin inside the isolation cabin. It appears to be in a state of more before and less after, and continues to float near medical staff, so that medical staff is at risk of infection. (2) Wearing a mask As shown in Figure 10, at 10 s, most of the droplets have been discharged from the outlet, and there is no obvious se lement on the stretcher of the patient's head. The remaining droplets eventually accumulate on the right side of the caregiver and move upward to the facial area. The risk of infection in caregiver is higher than that in other workers. At 120 s, only a few particles float throughout the fluid domain, droplets are present near caregiver and individuals, and most of the particles have been discharged.

Distribution of Droplets of Different Sizes
Sneezing will produce many droplets, and compared to coughing and talking, droplets are much larger. Research has concluded that their diameter range is 1~340 µm, with most of the droplets being about 210 µm in diameter.
(1) No masks As shown in Figure 11, the distribution of droplets at 1 s is clear. Droplets with a large particle size of 232 ~332 µm begin to fall near the patient's head, indicating that the mainstream does not have a significant effect on the droplets of large size. The droplets with a small particle size of 0.271~99.7 µm are located in the top area of the isolation cabin. At 10 s, the watershed is full small of droplets, all of a diameter of less than 80 µm, and they float in the watershed under the action of air flow. (2) Wearing a mask As shown in Figure 12, the diameter of the droplets at 1 s is less than 7.45 µm, and are distributed in a lump shape with the movement of the air flow to the outlet. The droplets of 10 s do not appear to be distributed differently due to the different sizes of the droplets, and all of the droplets are mixed together and spread to the isolation cabin.

Discrete Droplets Particle Residence Time
(1) No masks As shown in Figure 13, at 10 s, the two-second speaking process produces a very large number of droplets, and speaking generates airflow at a relatively low speed. Droplets do not reach the top of the isolation cabin and are mainly distributed around the caregiver under the influence of the mainstream. At this time, the caregiver is at high risk of infection. At 120 s, there are still many suspended droplets in the air, and fewer particles settle on the ground, indicating that most droplets can have been discharged. But the suspended particles are evenly distributed in the watershed, which poses a greater threat to medical staff. (2) Wearing a mask As shown in Figure 14, the droplets at the exit decreased significantly at 10 s. Some of the droplets are near the floor of the isolation cabin. A small amount of droplets are dispersed around the caregiver and have reached the caregiver's face as far as possible. At this time, the risk of infection of the caregiver is increased and the individuals are not threatened. At 120 s, the droplets remaining in the air were not completely discharged. Although the amount is very small and far from the facial area of the medical staff, the effect of the air flow makes it possible to reach the breathing area of the medical staff at any time. The risk of infection cannot be ignored.

Distribution of Droplets of Different Sizes
Speaking, unlike coughing and sneezing, produces very small air currents and droplets. The diameter range of the selected droplets is 1~50 µm. Most droplets are about 30 µm in diameter.

(1) No masks
As shown in Figure 15, at 1 s, most droplets with a diameter below 20 µm move rapidly and dispersed. Droplets with larger particle sizes of 33.9~48.3 µm are near the mouth and nose of the patient. Droplets at 10 s are distributed near the outlet and caregiver, and droplets are mixed relatively evenly. (2) Wearing a mask As shown in Figure 16, the droplets move in the outlet direction at 1 s, and most of the droplets have a particle size value of less than 20 µm. There is a tendency to spread upward. At 10 s, some droplets are discharged from the outlet, resulting in a significant decrease in the number of droplets. At this time, a small number of large particle size droplets of 10.9~13.6 µm remain in the facial area of the patients. The particle size of the droplets in the basin is less than 13.6 µm.

Comparison Analysis of the Optimized Model
In this paper, the four sets of optimization models adopt the active interference airflow method and are designed based on the principle of fresh air inlet. The optimization design only uses the most common cough in the new coronavirus for simulation. To be er obtain the features of each scheme under the condition of wearing a mask or not, the two working conditions of wearing a mask and no mask are still used to compare and analyze the location of droplets and the efficiency of sewage discharge, including being captured by medical staff and other walls, floating in the watershed, and discharged. According to the characteristics of the droplet particle size range produced under each working condition, the relationship between the droplets of different sizes and the five boundary conditions of patient, caregiver, individuals, ground, and outlet in the time dimension is determined.

Comparative Analysis of No Masks
As shown in Figure 17, a large number of droplets in each scheme are captured by each boundary. The large particle size of 50~70 µm has almost been discharged from the outlet in 10 s. Except for the second scheme, the patient boundary did not capture droplets in a short period of time, and except for the particles discharged from the outlet, the particle size captured at the remaining boundaries does not exceed 20 µm. Scheme 2 does not capture droplets after 60 s; the droplets of the four groups of schemes se led to the ground in the early stage and the outlet boundary is continuously discharged droplets, although the droplets se led in scheme 3 are relatively small, but not as good as scheme 2, and the number of droplets se led in the later stage of scheme 2 gradually decreased. Caregiver and individuals in the four groups catch droplets throughout the process, and although in small numbers, they are already a threat to medical staff.

Comparative Analysis of Mask Use
As shown in Figure 18, affected by the mask, patients in all protocols instantly capture droplets and the particle size exceeds 50 µm. Schemes 1 and 4 capture relatively few droplets on the ground, and the most droplets adsorb on the surface of the caregiver in scheme 3. When the droplets are captured by the individual boundary, scheme 3 captures relatively more after 30 s, and the particle size range is 0~30 µm, while scheme 1 and scheme 4 do not capture droplets within 120 s. The export boundary can reflect the number and time point of droplets discharged by the isolation cabin within 120 s, scheme 1 and scheme 4 basically discharge a large number of droplets within 10 s, and particles larger than 20 µm have been discharged within 5 s, while scheme 2 and scheme 3 continue to discharge particles in the range of 0~20 µm for the entire time period.

Comparative Analysis of Five Optimization Solutions
The purpose of the four optimization schemes is to actively interfere with the intensity and direction of airflow in the isolation cabin to affect the patient's breathing pa ern and the movement of particles. In the general quarantine cabin model, seemingly stable and regular watersheds are not optimal.
As shown in Figure 19, the simulation result data from the four sets of optimization schemes and the general model simulation data are compared, and from the number of droplets captured on the body surface of medical personnel, it can be judged that schemes 1 and 4 perform well in mask mode, and both schemes reduce the probability of medical staff carring the virus on the body surface. In the absence of masks, both schemes 2 and 3 can reduce the probability that medical staff carry the virus on the surface of their bodies. Patients can capture droplets under all working conditions, and droplets that se led on the floor of the negative pressure isolation cabin are also dangerous. Schemes 1, 3 and 4 for the mask-wearing condition are all reasonable, and schemes 2, 3 and 4 for the maskfree scenario are reasonable. From the perspective of the number of droplets discharged, all schemes in the mask-wearing mode are effective, and the maskless mode schemes 2 and 3 are effective.

Comparative Analysis of Sewage Efficiency
Discharge efficiency refers to the ability to discharge pollutants in a specific area through the input and outlet equipment, and the calculation formula of the discharge efficiency of pollutants at a certain point in time is as follows [34]: where C0 is the initial concentration of pollutants in the region, kg/m 3 . The Ct is the concentration of pollutants at the time of t, kg/m 3 .
Using the above method, the concentration of the droplets dispersed in the total volume of t = 0.50 s and t = 120 s is calculated, and the efficiency of sewage discharge under mask-free and mask-wearing conditions is shown in Tables 4 and 5.  The sewage efficiency needs to be higher than 93% and 86%, respectively, in the case of wearing a mask and without a mask to be called superior. In Figure 20a, schemes 1, 2 and 4 can all outperform the sewage efficiency of the general model. Only scheme 3 is lower than the sewage discharge efficiency of the general model. At this time, the concentration of droplets in the negative pressure isolation cabin is 10.0306 × 10 −12 kg/m 3 , which is higher than the simulation results of the other four groups. The simulation results of the four sets of optimization schemes in Figure 20b are be er than the simulation results of the general model, among which the simulation results of scheme 1 are basically the same as those of the general model. The sewage efficiency of scheme 2 is the best under this working condition, which can reach 92%.

Infection Risk Assessment
In addition to judging the flow of droplets from the patient's cough and the discharge efficiency of the entire watershed system, another important index is the assessment of infection risk. Although the negative pressure isolation cabin is be er at protecting medical staff and the external environment than a general isolation cabin, its internal structure is complicated and the airflow direction is more chaotic. The medical staff in the cabin remain in a prolonged and highly hazardous environment. Currently, there are four known transmission routes of droplets, among which contact transmission and aerosol transmission are the most common in the process of facilitating the isolation of patients in the isolation cabin [35]. Therefore, this study divided the risk area into two parts: human surface and respiratory area, and analyzed the infection risk from two aspects: contact transmission and aerosol transmission.

Risk Assessment of Contact-Borne Infections
We analyzed the infection risk to the medical staff accompanying the patient in the negative pressure isolation cabin, including one caregiver and two individuals. The collected droplet particle data included the number of droplets contained on the body surface of the three medical staff, so that the dose relationship of the simulation results could be accurately obtained. After analysis, a linear quantitative evaluation model is adopted to assess the infection risk of possible contact transmission pa erns on the body surface of medical personnel [36]: where Rj is the total infection risk of the medical staff j. Ri,j are single risks, which are multiplied by the number of viruses collected that may cause infection Ni,j, and the risk factor Ci. N0 is the baseline number of infection risk and the value is 500. In order to accurately calculate the risk of infection on the human surface, the droplet virus particles trapped on the human surface are divided into two regions, which are sorted according to the variable i of Equation (5) [36]: When I = 1, the value is the number of droplets captured by the face of the medical staff. Since the droplets are very close to the mouth and nose of medical staff, the value is 1.
When I = 2, the value is the number of droplets captured on the surface of the medical staff. At this time, medical staff will be infected with the virus through contact with clothing or skin, and there is a risk of sensitization, the value is 0.1.
When I = 3, the value is the number of virus particles on the floor of the carriage. Although the simulation process ignores the factor of human movement, the droplets settled on the ground will still blow up again due to the air flow and the bumps of the car's driving process, with a value of 0.01.
According to Equation (5), the risk of human surface infection is calculated in five groups of models, including the general model, as shown in Figure 21. In Figure 21a, a patient without a mask can produce a very high number of virus-carrying droplets, and the risk of infection does not exceed 7% in all protocols. In the first scheme, the risk of infection by caregiver and individuals carrying the virus through the human surface is the largest, both greater than that of the general model. The total risk of infection can reach 6.59%, and the risk of infection of caregiver is the highest. The overall infection risk of scheme 2 is lower than that of the general model and the total infection risk value was 3.75%. Scheme 3 has the lowest total infection risk, only 3.99%. The overall infection risk in scheme 4 is smaller than that of the general model, but the total individual infection risk is higher than that of the general model. In Figure 21b, the patient wearing a mask can effectively reduce the number of droplets produced, greatly reduce the intensity of airflow generated by the patient's breathing mode, and also greatly reduce the chance of infection of medical staff. Both scheme 1 and scheme 4 reduce the chance of infection to 0% for all medical staff under these conditions.

Risk Assessment of Aerosol-Borne Infections
The dose-response model was first used to assess the risk of tuberculosis infection. The EPA's four-step model of population health risk assessment defines the model as describing the likelihood and severity of adverse health effects at a certain dose and exposure to a chemical. In 2008, Sze To et al. proved the potential application of the viability function of airborne viral pathogens in exposure assessment and infection risk analysis, and obtained the exposure estimation method [37], which is expressed as: where c is the initial pathogen concentration; p is the respiration volume, m 3 /h; v (xi, ti) is the concentration of virus volume that varies over time in xi, the volume fraction of the pathogen m particle size pathogen; and f(t) is the pathogen's survivability parameter. Finally, based on the exponential dose-response model, a risk model containing estimated pathogen exposure is proposed. The sorted formula of the random non-threshold dose-response model is expressed as follows [37]: cpq v x t f t dt where PI (xi, ti) is the probability of infection or the risk of infection at the xi position within ti time; ti is the exposure time, h; m is the total number of pathogens; rm is the infectious parameter of the m particle size pathogen; βm is the percentage of exposure of the m particle size pathogen at the target site; and q is the frequency of coughing. From Equation (7) above, it can be seen that the model is used to evaluate the virus collected over time at xi position, while the aerosol droplets mentioned in this paper are dispersed in the drainage basin of the negative pressure isolation cabin. Although many droplets are produced when the patient does not wear a mask, the concentration of aerosol particles in the entire basin is not very high, so analyzing the risk value of infection at one point cannot fully reflect the risk value of infection of medical staff. After comprehensive consideration, in order to obtain more accurate and effective risk assessment result data, it is necessary to improve the theoretical basis of the dose-response model.
The MSDR model (Medical Staff Dose-response infection Risk Assessment Model) is improved on the basis of dose-response model, which can further analyze the influence of air flow organization factors on the infection risk of medical staff. It is a quantitative evaluation mechanism for the infection risk of medical staff [38,39].
According to the characteristics of MSDR model, Equation (6) is modified as follows: where E (z, tl) is the exposure evaluation model of a single droplet particle in the respiratory region within the time step l; d0k is the diameter of the k droplets, m; ∆ti is the i time step; H is the vertical height of the respiratory area, 0.2 m [40]; and S is the horizontal area of the medical personnel breathing in the isolation cabin, S0.60 = 4.58 m 2 , S1.05 = 4.87 m 2 , S1.45 = 4.66 m 2 . The MSDR model is obtained by modifying on the basis of Equation (7), and the expression is: In the formula, when d0k ≤ 3 µm, the infection parameter of the j pathogen rj = 0.387, otherwise rj = 0.0031. The percentage of the j particle size pathogen exposed at the target site βj = 2.587 .
. Pathogen survival parameters f(t) = 0.2 (0.975) t 60 . The initial concentration c was 1 × 10 5 TCID50/mL [41]. The respiratory volume p was 8.6 L/min. (1) No mask infection risk assessment As shown in Figure 22, the result of the infection risk when medical staff is standing is shown in Figure 22a. The infection risk of the five groups increased rapidly before 12 s, and the value was in the range of 0.25963~0.07923. The risk value of infection in scheme 1 and scheme 4 is 2~3 higher than that of the general model. The infection risk value of scheme 2 and the general model is basically equal, and the infection risk value of scheme 3 is the lowest. The infection risk of medical staff while si ing is shown in Figure 22b (2) Infection risk assessment by wearing a mask This is shown in Figure 23. The infection risk result when the caregiver is standing is shown Figure 23a. The infection risk value of the general model is 0.00590. Except for the third scheme, the infection risk value of other optimized schemes is lower than that of the general model. The result of the infection risk of the medical staff si ing is Figure 23b. The overall infection risk value is much smaller than that of Figure 23a, but the change trend is basically the same. The results of the infection risk at work for medical staff is shown in Figure 23c. The infection risk value in the general model is 1.11982 × 10 −4 . Although the effect of schemes 2 and 3 is worse than that of the general model, their infection risk values are of the same order of magnitude as that of the general model, and the change is not great, and the best effect is scheme 4. The infection risk value is only 1.54398 × 10 −9 . Since it has been mentioned in this paper that it is possible for patients to wear masks or not, from the perspective of infection risk assessment, the optimal scheme is finally selected as scheme 2. With normal si ing posture, the risk of infection is reduced by 71.77%.

Conclusions
As an important tool to prevent aerosol-borne diseases and to protect medical staff, the complex flow field structure in the negative pressure isolation cabin cannot be ignored. Optimizing the internal structure of the isolation cabin is of great significance for the prevention and treatment of infectious diseases. According to the research content of this paper, the main conclusions are as follows: (1) Based on the principle of active interference, auxiliary air inlets are added above or behind the faces of caregiver and individuals, and four sets of optimization schemes are proposed. The diffusion of droplets generated by the three respiratory modes of coughing, sneezing, and talking by patients in the general model is simulated, and the residence time and particle size distribution of droplets in the basin are compared and analyzed with or without the patient wearing a mask. The results show that the medical staff in the negative pressure isolation cabin have different degrees of exposure risk. (2) The common cough breathing mode in COVID-19 is analyzed and the superiority of optimal scheme 2 is demonstrated by comparing its droplet particle destination and pollutant discharge efficiency. Compared to the general model, the sewage efficiency of the scheme is increased by 3% when the patient wears a mask and 6% when the patient does not wear a mask. The sewage efficiency, as a reference data, is 96% and 92% without and with masks, respectively. (3) The risk of infection is evaluated from the aspects of contact transmission and aerosol transmission. The linear quantitative evaluation model and MSDR model are used to evaluate these transmission risks, respectively, and the results again show that optimization scheme 2 is the best. Compared with the general model, the risk of infection without mask is reduced by 34.85% and the risk of infection with mask is reduced by 71.77%. These results indicate the importance of structural optimization of the negative pressure isolation cabin of large cruise ships.
The results of this paper have limitations to a certain extent because the droplets produced by human coughing and other breathing methods contain sodium, potassium, and other elements, and the failure to refine the droplet parameters will cause deviations in the calculated results. Due to various factors, the relevant experiments have not been carried out, resulting in the simulation results not being compared and analyzed with the experimental data, and the parameter correction link is missing. Under the influence of many factors, it is still necessary to explore new methods and new ideas to study and improve the special isolation tool of the negative pressure isolation cabin in cruise ships.