Modelling levels of nitrous oxide exposure for healthcare professionals during EMONO usage
Annals of Occupational and Environmental Medicine volume 28, Article number: 30 (2016)
Computational fluid dynamics (CFD) has been used to compute nitrous oxide (N2O) levels within a room during the administration of an equimolar mix of N2O/oxygen (EMONO) in the clinical setting. This study modelled realistic scenarios of EMONO usage in hospital or primary care, in order to estimate the potential N2O exposure of healthcare professionals (HCP) with routine EMONO use and to provide guidance for EMONO users.
Sixteen scenarios were defined by carrying out a survey of practitioners. CFD simulations were performed for each scenario and N2O concentrations over time were calculated. N2O exposures (time-weighted average of concentration over 8 h [TWA-8 h]) were calculated at the HCPs’ mouth to be compared with a predefined occupational exposure limit (OEL).
Administration duration and ventilation type were the main factors influencing N2O levels; ventilation type also influenced wash-out time between EMONO administrations. N2O concentration showed a plume distribution towards the ceiling and was highly heterogeneous, highlighting the importance of measurement location. Although estimated TWA-8 h varied widely, 13 of the 16 scenarios had an N2O TWA-8 h of <100 parts per million.
Data demonstrate that EMONO usage in well ventilated rooms – as recommended – helps to ensure that N2O exposure does not exceed the OEL and does not signal any major risks for HCPs when recommendations are followed. Although these data are numerical simulations and should be considered as such, they can provide guidance for EMONO users.
The N2O/O2 50 %/50 % admixture, also known as Equimolar Mixture of Oxygen and Nitrous Oxide (EMONO, Kalinox™ Air Liquide Santé International, France, Equanox™, Alnox™, Nitronox™), is indicated for the treatment of short-term pain conditions of mild to moderate intensity when rapid-onset powerful short-term relief of pain is required [1, 2]. The risks associated with chronic occupational N2O exposure for healthcare professionals (HCP) are still controversial [3–5]. Hence, for several decades, some cases of repeated and prolonged exposure have been reported to be linked to medical conditions including reproductive, neurologic, hepatic, and renal disorders [5, 6].
The National Institute for Occupational Safety and Health (NIOSH) recommended exposure limit (REL) is 25 parts per million (ppm) as a time-weighted average (TWA) in operating room during the period of general anesthesia and the American Conference of Governmental Industrial Hygienists (ACGIH) threshold limit value (TLV) for N2O is 50 ppm as an 8-h TWA (TWA-8 h) . For HCPs administrating EMONO, specific occupational exposure limits (OELs) have been defined in some countries, for example in the US, the N2O OEL TWA-8 h is 50 ppm (90 mg/m3)  and in the UK and some European countries it is 100 ppm (180 mg/m3) [8–10].
These data are not easily measurable or predictable and experimental airflow studies are difficult and expensive due to their complexity ; therefore, Computational Fluid Dynamics (CFD) software is a useful tool that can be used to generate 3-dimensional (3D) numerical simulations to predict gas flow and gas concentrations in specified enclosed environments. It has been used in a wide range of applications such as investigating ventilation effectiveness in clinical settings such as operating rooms [11–18] and birthing rooms . In other healthcare applications, CFD has been used successfully to assess the benefits of a laminar flow screen to reduce bacterial contamination in hospitals , and to investigate airborne contaminants in medical [18, 20–24] and dental clinics . Other medical applications of CFD encompass analysis of blood flow through arteries [26, 27], assessing rupture risk in cerebral aneurysms , and analysis of exhaled respiratory gases [29, 30].
Ideally, adequate ventilation systems, room design, equipment maintenance, and clinical practices should ensure that staff are never subjected to excessive levels of N2O OEL . In view of the difficulties and cost of measuring N2O distribution in real-life situations, the feasibility of using 3D numerical simulation with CFD software to predict the distribution of N2O concentrations with EMONO use within a room has been assessed. The comparison of computed data with experimental measures of N2O concentration within a laboratory showed good agreement, both in terms of maximal concentrations achieved during EMONO administration and in terms of the estimated TWA-8 h. The changes in N2O concentration over time during EMONO administration and wash-out phases were also correctly described in the simulation, supporting the relevance of using 3D numerical simulation to compute the N2O TWA-8 h for HCPs.
Hence, the objective of our work was to use a numerical modelling approach to compute the N2O concentration distribution in a room during various realistic scenarios of EMONO administration, to predict the associated TWA-8 h ppm N2O exposure for a range of HCPs using EMONO and assess the pertinence of the current EMONO usage recommendations.
Scenarios and parameters used: In order to perform CFD simulation, detailed modelling of the room features such as dimensions, patient’s location, ventilation, door and window, the initial conditions within the room, and various input parameters such as EMONO flow rate and administration time are required. As an infinite number of scenarios are possible, for practical reasons, we defined 16 scenarios most representative of real life EMONO administration by surveying various HCPs using EMONO for different indications, such as painful procedures in children and adults emergency rooms, or for dermatology, dentistry, and obstetrics procedures.
For each scenario, a full set of environmental and EMONO administration parameters was specified. Environmental parameters included room size, the type of ventilation present (natural via windows and doors or a controlled mechanical ventilation system or extraction system), the number of subjects located within the room in addition to the HCP (an adult patient or a child patient and his/her parent), the temperature of the room (and outside temperature if windows and doors were open), the number of windows and doors and the location of the air inlet/outlet. The EMONO administration parameters included EMONO flow rate and duration, continuous or intermittent administration mode, the use of an on-demand valve or not, and the number of successive repetitions of EMONO administrations (Table 1).
Numerical simulations of these 16 scenarios based on real conditions of EMONO in hospital or in primary care setting were performed using CFD software. A virtual representation of a room was modelled in 3D (Fig. 1) in accordance with each scenario, using seven geometries with variations in room size, presence of window, single or double doors, mechanical ventilation systems, child or adult patient (and presence of parent): two different geometries for dentistry, one for dermatology, two for emergency treatment rooms (child and adult), one for obstetrics and one for radiology.
The patient was represented as having a supine position on a bed at the center of the room, and the HCP was represented as standing close to the patient’s bed. The parent, if featured in the scenario, was represented as standing close to the child patient on the opposite side from the HCP.
Human bodies (HCP + patient + parent if child patient) were represented as being realistic human body shapes (Fig. 1). Indeed, the shape of the body has been shown to affect gas distribution and concentrations around a human body with a rounded realistic body shape compared to oversimplified geometrical shape to represent the body [22, 24]. Any motion of practitioners or parent was not included in the model.
A leak was assumed to exist at the mask level – defined as a cylindrical surface 10 cm in diameter and 1 mm in width at the edge of the mask – and was the source of N2O presence in the room. The leak was assumed to be 10 % of the EMONO flow rate (i.e. 10 % of 6 L/min for a child and 10 % of 12 L/min for an adult patient). The species model assumed 0.79 N2 + 0.21 O2 mole fractions in the room initially and the leak at the mask level was defined to be 0.5 N2O + 0.5 O2 mole fractions. For scenarios using continuous EMONO administration (i.e., without on-demand valve), the leak was assumed to be continuous, whereas with an on-demand valve, the leak was assumed to occur only during the inspiration phase.
Turbulence modelling was included, as well as appropriate models to take into account thermal effects, e.g. temperature of the walls and room, and natural convection due to heat flux from the bodies. Turbulence has been observed to be an important phenomenon for accuracy of CFD modelling and for the optimum location of sampling points [14, 22, 24, 31–33]. The temperature of the room and EMONO were included in the model. The temperature of primary care offices (dentistry and dermatology) was assumed to be 20 °C with corridor 18 °C and outside 15 °C, and a hospital room (emergency, radiology, and obstetrics) was assumed to be 23 °C, with corridor 21 °C. The temperature of EMONO was assumed to be equal to the ambient room temperature. To model the heat released by people present in the room (and model natural convection), the following values were used: supine child patient = 55 W (heat flux Q = 74 W/m2); supine adult patient = 80 W (Q = 65.5 W/m2); relaxed parent = 110 W (Q = 63.7 W/m2); and working medical staff = 143 W (Q = 82.6 W/m2). Since the air initially present in the room and the EMONO administered to the patient can be considered as incompressible gases, the volume of gas entering the room due to EMONO inflow has to be equal to the volume of gas leaving the room to ensure mass conservation. Therefore, in the model an outflow surface at the bottom of the door (a door leak) was included to permit some air flow even when closed.
Fluent CFD software (ANSYS, PA, USA) was used to compute the N2O concentration in each part of the room (millions of points), at every second of the administration process, for each scenario. Because the location of the measurement point is known to influence the TWA-8 h value [22, 33, 34], a virtual sensor (probe) was positioned by the HCP’s mouth, to be representative of the inhaled exposure. The N2O exposure level for the HCP (TWA-8 h ppm) was computed using the simulated values of N2O concentration at this probe throughout the simulated scenario.
A wash-out phase was included in each scenario, following the EMONO administration phase. The wash-out time is defined as the time needed to evacuate at least 95 % of the N2O mass present in the room at the end of the EMONO administration phase. When there is no mechanical ventilation system in place, the (single or double) door or window is supposed to be closed during administration and open during the wash-out phase. Hence, when this natural ventilation starts, the N2O concentration within the room starts decreasing. For scenarios considering an active ventilation system in place, this system is supposed to run during both administration and wash-out phases.
Data on N2O concentration (ppm), temperature (°C), and gas velocity (m/s) were gathered for each scenario throughout the EMONO administration and wash-out phases. Quantitative data were visualized in 2D and 3D; an example of the computed N2O concentration values represented in two 2D slices intersecting the mask for one scenario is shown in Fig. 1.
Computed temperature and gas motion and velocity within the room were modelled to demonstrate the influence of natural convection (Fig. 2).
As the N2O concentration is numerically computed throughout the room and over time, the N2O concentration at the HCP’s mouth was calculated at each second of the administration and wash-out phases. The N2O exposure in TWA-8 h ppm at the HCP’s mouth probe was calculated for each scenario to be compared with country-dependent OELs. In addition, the evolution of N2O TWA-8 h ppm over time can be computed and gives useful insight on exposure levels.
N2O did not accumulate at floor level, instead the N2O concentration tends to show a plume distribution, orientated toward the ceiling (Figs. 1 and 2). However, the global N2O concentration distribution is very different for the various scenarios tested (affected by duration, use of mechanical ventilation, on-demand valve, etc.).
The change in N2O concentration in the room calculated for each scenario showed a clear influence of the specific ventilation system used. For example, in the dental treatment room, the N2O concentration level was higher without an active ventilation system, particularly with longer administration durations (Fig. 3). With active ventilation, a plateau effect is induced over time, at which point, the N2O global concentration in the room remains constant to the end of the current EMONO administration. Similarly, administration without an on-demand valve resulted in considerably higher N2O concentrations compared with administration using an on-demand valve, particularly with longer administration durations (Fig. 4).
It has to be noted that computed wash out times varied from 20 to 40 min, highlighting the role of a delay between two EMONO administrations. The decrease in N2O concentration during the wash-out phase was quicker when a double-door was open (compared with a single door or window), but is highly dependent on the temperature gradients used in the model. Therefore, the duration of the wash-out phase would be different if other temperatures were selected for the model.
The estimated TWA-8 h values varied widely (from 2 to 371 ppm), illustrating the influence of the conditions of EMONO administration. Thirteen of the 16 scenarios had TWA-8 h <100 ppm, including eight scenarios with TWA-8 h <25 ppm. The changes in TWA-8 h values over time obtained by numerical simulations were computed (Fig. 5). The influence of the administration duration is clearly shown in panels A and B, whereas the importance of the ventilation system is highlighted in panel C. Shorter procedures are presented in panel D.
Our computed numerical results of N2O levels in common clinical situations tend to demonstrate that administration of EMONO in well ventilated rooms, and/or with on-demand valve, and/or for short duration procedures, as recommended, would ensure that the N2O exposure remains below the OEL. These findings do not signal any major risks for the HCP when these recommendations are followed.
The simulated data show very high N2O concentration near the mask and considerable differences among the various scenarios tested, as expected. The data also show substantial heterogeneity in N2O concentration throughout the room in many scenarios, highlighting the importance of the location of the probe used to determine exposure and showing that using one point of measurement may not be representative of HCP N2O exposure.
Importantly, N2O does not accumulate at floor level, as might be expected due to the higher density of N2O compared with that of air, instead the N2O concentration tends to show a plume distribution, oriented toward the ceiling. Plume formation is due to the uneven temperature distribution throughout the room inducing movement of gases via natural convection. The temperature gradients observed in the room are mainly due to the heat generated by the human bodies (HCP, patient, and parent). Numerical simulation with CFD has been used to predict heat release from a human body using modelling of airflow and thermal radiation [14, 35] and CFD has also shown that using realistic human shapes is important in accurate modelling of thermal radiation and impact on air flow and turbulence .
Our data also demonstrate the importance of the recommended wash-out phase between successive N2O administrations to minimize the overall cumulative N2O exposure over time. Wash-out times were similar for children and adults, with and without an on-demand valve, and with and without active ventilation.
In this study, it was calculated that the N2O TWA-8 h was >100 ppm (i.e. higher than the most relaxed OEL) at the HCP’s mouth in only three scenarios – all of which were related to using a long EMONO administration (at least 59 min) with no active ventilation system. This highlights the importance of active ventilation particularly for longer durations of EMONO use.
Limitations of this study are those inherent to any model in that the data are hypothetical and only pertain to the 16 scenarios tested (e.g. size of the window, outside temperature). Beyond the influences on N2O concentration distribution described herein, there are several others that are beyond the scope of this study. For example, the motion of the practitioners is not included in the simulations (i.e., the HCP stays in the same room and same position for eight hours). This is unlikely in real life situations, but using this hypothesis leads to an overestimation of the theoretical N2O exposure of the HCP. However, in terms of CFD, computational aspects of the simulations have been checked for accuracy using standard practices in the field.
In conclusion, the current EMONO usage recommendations state that the product should be used in airy, well ventilated rooms. Our simulated data demonstrate that following these recommendations, together with use of an on-demand valve and/or for short procedures, would ensure that the N2O exposure for HCPs does not exceed the TWA-8 h OEL defined for each country. It should be noted that these results are derived purely from numerical simulations and should be considered as such; however, earlier research has shown the validity of using CFD to simulate gas levels in treatment rooms and so they can be considered a good estimation of N2O levels likely with EMONO used in the scenarios studied.
It is hoped that these data will provide useful information and guidance for EMONO users to ensure exposure recommendations are followed.
3D, 3-dimensional; ACGIH, American conference of governmental industrial hygienists; CFD, computational fluid dynamics; EMONO, equimolar mix of N2O/oxygen; HCP, healthcare professionals; N2O, nitrous oxide; ODV, on-demand valve; OEL, occupational exposure limit; ppm, parts per million; TWA-8 h, time-weighted average of concentration over 8 h.
Air Liquide Healthcare. Equanox: Summary of Product Characteristics. 2011. http://www.uk.airliquide.com/file/otherelement/pj/fd/14/21/29/equanox%20spc8535403540716584597.pdf. Accessed 1 July 2016.
Air Liquide Healthcare. Kalinox 50 %/50% medicinal gas, compressed: Summary of Product Characteristics. 2014. https://docetp.mpa.se/LMF/Kalinox%20medicinal%20gas%20compressed%20ENG%20. Accessed 1 July 2016.
Howard WR. Nitrous oxide in the dental environment: assessing the risk, reducing the exposure. J Am Dent Assoc. 1997;128(3):356–60.
Krajewski W, Kucharska M, Wesolowski W, Stetkiewicz J, Wronska-Nofer T. Occupational exposure to nitrous oxide - The role of scavenging and ventilation systems in reducing the exposure level in operating rooms. Int J Hyg Environ Health. 2007;210:133–8.
Clark MS, Brunick AL. Potential biohazards for health personnel associated with chronic exposure to N2O. Handbook of nitrous oxide and oxygen sedation. Thirdth ed. St Louis, MI: Mosby, Elsevier; 2008.
Kugel G, Norris LH, Zive MA. Nitrous oxide and occupational exposure: it’s time to stop laughing. Anesth Prog. 1989;36(6):252–7.
NIOSH. Controlling Exposures to Nitrous Oxide During Anesthetic Administration: DHHS (NIOSH) Publication Number 94–100. US Department of Health and Human Services, Public Health Service, Centers for Disease Control and Prevention, National Institute for Occupational Safety and Health. 1994. http://www.cdc.gov/niosh/docs/94-100/pdfs/94-100sum.pdf. Accessed 1 July 2016.
HMSO. Control of Substances Hazardous to Health (COSHH) Regulations: Approved Codes of Practice. London: HMSO; 1999. Accessed January 13 2015.
Health and Safety Commission. Table 1: List of approved workplace exposure limits (as consolidated with amendments October 2007). 2007. http://epsassets.manchester.ac.uk/medialand/psi/formsandguidance/WorkplaceExposureLimits.pdf. Accessed 13 Jan 2015.
Henderson KA, Matthews IP, Adisesh A, Hutchings AD. Occupational exposure of midwives to nitrous oxide on delivery suites. Occup Environ Med. 2003;60(12):958–61.
Sadrizadeh S, Tammelin A, Nielsen PV, Holmberg S. Does a mobile laminar airflow screen reduce bacterial contamination in the operating room? A numerical study using computational fluid dynamics technique. Patient Saf Surg. 2014;8:27.
Chow TT, Yang XY. Ventilation performance in the operating theatre against airborne infection: numerical study on an ultra-clean system. J Hosp Infect. 2005;59(2):138–47. doi:10.1016/j.jhin.2004.09.006.
Cheong KWD, Phua SY. Development of ventilation design strategy for effective removal of pollutant in the isolation room of a hospital. Build Environ. 2006;41:1161–70. doi:10.1016/j.buildenv.2005.05.007.
Qian H, Li Y, Nielsen PV, Hyldgaard CE. Dispersion of exhalation pollutants in a two-bed hospital ward with a downward ventilation system. Build Environ. 2008;43:344–54. doi:10.1016/j.buildenv.2006.03.025.
Yam R, Yuen PL, Yung R, Choy T. Rethinking hospital general ward ventilation design using computational fluid dynamics. J Hosp Infect. 2011;77:31–6. doi:10.1016/j.jhin.2010.08.010.
Yuen PL, Yam R, Yung R, Choy KL. Fast-track ventilation strategy to cater for pandemic patient isolation surges. J Hosp Infect. 2012;81:246–50. doi:10.1016/j.jhin.2012.04.013.
Beggs CB, Kerr KG, Noakes CJ, Hathway EA, Sleigh PA. The ventilation of multiple-bed hospital wards: review and analysis. Am J Infect Control. 2008;36(4):250–9. doi:10.1016/j.ajic.2007.07.012.
Sadrizadeh S, Holmberg S. Surgical clothing systems in laminar airflow operating room: a numerical assessment. J Infect Public Health. 2014;7:508–16. doi:10.1016/j.jiph.2014.07.011.
Nedin P, Coppel A. Priorities for birthing rooms set out. Health Estate. 2014;68(7):61–6.
Chau OKY, Liu C-H, Leung MKH. CFD Analysis of the Performance of a Local Exhaust Ventilation System in a Hospital Ward. Indoor Built Environ. 2006;15:257–71. doi:10.1177/1420326X06066123.
Inthavong K, Ge QJ, Li XD, Tu JY. Detailed predictions of particle aspiration affected by respiratory inhalation and airflow. Atmos Environ. 2012;62:107–17. doi:10.1016/j.atmosenv.2012.07.071.
Li J, Yavuz I, Celik IB, Guffey SE. A numerical study of worker exposure to a gaseous contaminant: variations on body shape and scalar transport model. J Occup Environ Hyg. 2005;2:323–34. doi:10.1080/15459620590961684.
Richmond-Bryant J, Eisner AD, Brixey LA, Wiener RW. Transport of airborne particles within a room. Indoor Air. 2006;16:48–55. doi:10.1111/j.1600-0668.2005.00398.x.
Russo JS, Ezzat Khalifa H. CFD assessment of intake fraction in the indoor environment. Build Environ. 2010;45:1968–75. doi:10.1016/j.buildenv.2010.01.017.
Chen C, Zhao B, Cui W, Dong L, An N, Ouyang X. The effectiveness of an air cleaner in controlling droplet/aerosol particle dispersion emitted from a patient’s mouth in the indoor environment of dental clinics. J R Soc Interface. 2010;7(48):1105–18. doi:10.1098/rsif.2009.0516.
Lee J, Smith NP. The multi-scale modelling of coronary blood flow. Ann Biomed Eng. 2012;40(11):2399–413. doi:10.1007/s10439-012-0583-7.
Shim EB, Kamm RD, Heldt T, Mark RG. Numerical analysis of blood flow through a stenosed artery using a coupled, multiscale simulation method. Comput Cardiol. 2000;27:219–22.
Cebral JR, Castro MA, Burgess JE, Pergolizzi RS, Sheridan MJ, Putman CM. Characterization of cerebral aneurysms for assessing risk of rupture by using patient-specific computational hemodynamics models. Am J Neuroradiol. 2005;26(10):2550–9.
Gao N, Niu J. Transient CFD simulation of the respiration process and inter-person exposure assessment. Build Environ. 2006;41:1214–22. doi:10.1016/j.buildenv.2005.05.014.
Zhang T, Yin S, Wang S. Quantify impacted scope of human expired air under different head postures and varying exhalation rates. Build Environ. 2011;46:1928–36. doi:10.1016/j.buildenv.2011.03.014.
Deevy M, Sinai Y, Everitt P, Voigt L, Gobeau N. Modelling the effect of an occupant on displacement ventilation with computational fluid dynamics. Energy Build. 2008;40:255–64. doi:10.1016/j.enbuild.2007.02.021.
Dai Y-T, Juang Y-J, Wu Y, Breysse PN, Hsu D-J. In vivo measurements of inhalability of ultralarge aerosol particles in calm air by humans. J Aerosol Sci. 2006;37:967–73. doi:10.1016/j.jaerosci.2005.10.002.
Ojima J. Gaseous contaminant distribution in the breathing zone. Ind Health. 2012;50:236–8.
Rim D, Novoselac A. Occupational exposure to hazardous airborne pollutants: effects of air mixing and source location. J Occup Environ Hyg. 2010;7:683–92. doi:10.1080/15459624.2010.526894.
Murakami S, Kato S, Zeng J. Combined simulation of airflow, radiation and moisture transport for heat release from a human body. Build Environ. 2000;35:489–500. doi:10.1016/S0360-1323(99)00033-5.
Medical writing assistance with the preparation of this manuscript was provided by Mary Hines of Springer Healthcare Communications. Funding for this assistance was provided by Air Liquide Santé International.
This work was supported by Air Liquide Santé International.
MP conceptualized/designed the study, collected the data, conducted analysis and drafted and reviewed the draft manuscript. CB conceptualized/designed the study, drafted and reviewed the draft manuscript. GC conceptualized/designed the study, drafted and reviewed the draft manuscript. All authors read and approved the final manuscript.
Marine Pichelin PhD, Catherine Billoet MD, Georges Caillibotte PhD are employees of Air Liquide Santé International, Medical R&D entity. They hold stocks in Air Liquide.
About this article
Cite this article
Pichelin, M., Billoet, C. & Caillibotte, G. Modelling levels of nitrous oxide exposure for healthcare professionals during EMONO usage. Ann of Occup and Environ Med 28, 30 (2016). https://doi.org/10.1186/s40557-016-0116-1
- 3D numerical simulation
- Equimolar mix of nitrous oxide (N2O)/oxygen (O2)
- Occupational exposure limit