Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Computational simulation to assess patient safety of uncompensated COVID-19 two-patient ventilator sharing using the Pulse Physiology Engine

  • Jeffrey B. Webb ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Writing – original draft, Writing – review & editing

    jeff.webb@kitware.com

    Affiliation Kitware, Inc., Carrboro, North Carolina, United States of America

  • Aaron Bray,

    Roles Conceptualization, Data curation, Investigation, Methodology, Software, Visualization, Writing – original draft

    Affiliation Kitware, Inc., Carrboro, North Carolina, United States of America

  • Philip K. Asare,

    Roles Conceptualization, Data curation, Formal analysis, Methodology, Writing – original draft

    Current address: Institute for Studies in Transdisciplinary Engineering Education & Practice and Division of Engineering Science, University of Toronto, Toronto, Ontario, Canada

    Affiliation Electrical and Computer Engineering, Bucknell University, Lewisburg, Pennsylvania, United States of America

  • Rachel B. Clipp,

    Roles Conceptualization, Project administration, Supervision, Writing – original draft, Writing – review & editing

    Affiliation Kitware, Inc., Carrboro, North Carolina, United States of America

  • Yatin B. Mehta,

    Roles Conceptualization, Formal analysis, Methodology, Writing – original draft

    Affiliation Critical Care Medicine, Geisinger, Danville, Pennsylvania, United States of America

  • Sudheer Penupolu,

    Roles Conceptualization, Formal analysis, Methodology, Writing – original draft

    Affiliation Critical Care Medicine, Geisinger, Danville, Pennsylvania, United States of America

  • Aalpen A. Patel,

    Roles Conceptualization, Formal analysis, Methodology, Writing – original draft

    Affiliation Department of Radiology and Steele Institute for Health Innovation, Geisinger, Danville, Pennsylvania, United States of America

  • S. Mark Poler

    Roles Conceptualization, Formal analysis, Methodology, Project administration, Supervision, Writing – original draft, Writing – review & editing

    Affiliation Anesthesiology, Geisinger, Danville, Pennsylvania, United States of America

Abstract

Background

The COVID-19 pandemic is stretching medical resources internationally, sometimes creating ventilator shortages that complicate clinical and ethical situations. The possibility of needing to ventilate multiple patients with a single ventilator raises patient health and safety concerns in addition to clinical conditions needing treatment. Wherever ventilators are employed, additional tubing and splitting adaptors may be available. Adjustable flow-compensating resistance for differences in lung compliance on individual limbs may not be readily implementable. By exploring a number and range of possible contributing factors using computational simulation without risk of patient harm, this paper attempts to define useful bounds for ventilation parameters when compensatory resistance in limbs of a shared breathing circuit is not possible. This desperate approach to shared ventilation support would be a last resort when alternatives have been exhausted.

Methods

A whole-body computational physiology model (using lumped parameters) was used to simulate each patient being ventilated. The primary model of a single patient with a dedicated ventilator was augmented to model two patients sharing a single ventilator. In addition to lung mechanics or estimation of CO2 and pH expected for set ventilation parameters (considerations of lung physiology alone), full physiological simulation provides estimates of additional values for oxyhemoglobin saturation, arterial oxygen tension, and other patient parameters. A range of ventilator settings and patient characteristics were simulated for paired patients.

Findings

To be useful for clinicians, attention has been directed to clinically available parameters. These simulations show patient outcome during multi-patient ventilation is most closely correlated to lung compliance, oxygenation index, oxygen saturation index, and end-tidal carbon dioxide of individual patients. The simulated patient outcome metrics were satisfactory when the lung compliance difference between two patients was less than 12 mL/cmH2O, and the oxygen saturation index difference was less than 2 mmHg.

Interpretation

In resource-limited regions of the world, the COVID-19 pandemic will result in equipment shortages. While single-patient ventilation is preferable, if that option is unavailable and ventilator sharing using limbs without flow resistance compensation is the only available alternative, these simulations provide a conceptual framework and guidelines for clinical patient selection.

Introduction

Even in well-resourced healthcare systems the current COVID-19 (SARS-CoV-2 disease) pandemic has led to equipment shortages in Europe, the United States, and is a more severe challenge across the world. Inadequate availability of ventilators for patients in respiratory failure is especially challenging [13].

Ventilator sharing or splitting is one mitigation strategy that has been used in practice to rescue as many patients as possible when sufficient numbers of ventilators are not immediately available. Media reports from Las Vegas, Nevada, USA in October 2017 popularized this tactic employed during the surge of over 500 casualties after a mass shooting incident. Manual and improvised ventilation support was needed until stockpiled surge reserve ventilators could be delivered. Curiously, though news media reports are still readily discoverable, no citable healthcare literature directly addressing reported shared ventilator use after the Las Vegas event has been found using extensive searches of PubMed, Google Scholar, and OVID.

The simplest approach for multiplex ventilation is to connect branched tubing without any compensating resistance in any branch for limiting flow to the most compliant patient(s). This is the only way shared ventilation is likely to be implemented rapidly with available tubing and adapters, without ready availability of flow compensating devices.

Only limited testing and clinical experience are available to expose the potential pitfalls of simplistic implementation. Branched breathing circuits to ventilate more than one patient with a single ventilator or compressed gas sources were published by Sommer et al. in 1994 [4], after contemplation of needs for ventilation of mass casualties from gas warfare at the time of the 1990-1 Gulf War [5]. A simpler concept published in 2006 with rudimentary simulation was employed at Las Vegas [6, 7]. Tonetti et al. tested feasibility using test lungs of this simply assembled circuit without resistance compensation and tabulated procedures to guide use [8]. Chatburn et al. have studied ventilation of two lung simulators using pressure control and volume control modes of ventilation for six pairs of simulated patients without compensating flow resistance. Simulated carbon dioxide tensions were estimated and used to estimate pH for each modeled patient [9]. Hermann et al. have studied pressure and volume ventilation modes using a computational model of lung ventilation and by test lung experimental simulation [10].

Paladino et al. published a 12-hour proof-of-concept using a single ventilator for four healthy adult sheep, reporting some difficulties even with well-matched “patients” [11]. Letters in response to Paladino’s study raised concerns that the demonstration was not a sound foundation for translation to clinical use, especially for anticipated mass casualty scenarios producing a variety of lung pathophysiologies. They discussed anticipated issues including differing compliances, spontaneous breathing, dyssynchrony with the ventilator, expected shortages of pharmaceuticals and oxygen, monitors, and skilled personnel for vigilant management of patients on shared breathing circuits in a surge situation [5, 12].

Due to the perceived risks of having multiple patients per ventilator a consortium of North American professional societies issued a joint statement opposed to this practice, available from several internet sources [13]. The joint statement gives only slight acknowledgment to desperate triage situations whereby dedication of one ventilator exclusively to one patient may determine the demise of another if some means of ventilation cannot be provided. The enumerated considerations of the joint consensus statement are realistic and concerning. However, they abdicate the possibility that split ventilation could be safely and effectively implemented. The identified issues may be considered an agenda for further investigation.

Faced with mounting needs for ventilator support of COVID-19 patients and limited resources, the Federal Emergency Management Administration (FEMA) and healthcare systems prepared clinical guidelines for ventilator sharing using pressure support mode, pharmacological paralysis, and vaguely stated “similar mechanical support needs” [1416]. Potential risks are enumerated, including unintentional extubation of one patient, cross-contamination, and delayed observation of hypo- or hyper-ventilation. The US Department of Health and Human Services (HHS) issued a document containing a statement that the FDA would not object to split ventilation for the duration of the declared COVID-19 emergency when demand for invasive ventilation exceeds the supply of ventilators [17]. Protocols from a Federal Emergency Management Administration (FEMA) Washington DC COVID-19 Taskforce (Babcock et al. [14]) and the Columbia University College of Physicians & Surgeons collaboration with New York-Presbyterian Hospital (Beitler et al. [15]) are included within the attached HHS document.

Limited clinical demonstrations have been reported. Simple branched-flow circuits for a ventilator do not incorporate flow compensation to limit excessive distribution of ventilation to the most compliant lungs. Beitler et al. reported encouraging results for 3 pairs of carefully matched stable New York COVID-19 ARDS patients, each pair slowly and cautiously prepared for tandem ventilation under a detailed protocol. Each pair was ventilated for 2 days using different equipment or refinements to protocol. Limitations of an anesthesia machine as a ventilator for two patients and consequences of inadequate neuromuscular blockade are reported [18]. Another New York group at Mount Sinai empirically compensated for the maldistribution of ventilation by incorporating a needle valve to restrict breath flowing to the most compliant lungs. Initial simulation testing of shared ventilation using a single ventilator employed two high-fidelity lab manikins. Subsequently, two consented ICU patients were supported for one hour by a single ventilator with a split breathing circuit. At another time, two other consented ICU patients were supported for an hour by one ventilator. Some of the assertions of the previously cited joint statement are addressed by Levin et al. in an appendix [19]. Similar apparatus has been described by an on-line collaborative [20] and another pre-publication document from Yale [21]. Clarke et al. have reported using simple compression clamps on segments of tubing to increase resistance [22].

Han et al. have empirically demonstrated advantages of a somewhat more complicated bag-in-box design for parallel co-ventilation of test lungs driven by a single ICU ventilator, avoiding respired gas cross-contamination, and allowing somewhat independent control of VT, PEEP, and FiO2 for widely different lung compliances and single circuit faults [23].

The joint statement and limited empirical information suggest a range of inadequately addressed issues. These provide opportunities for future conventional research, modeling and simulation. However, the pace of conventional research is too slow when addressing a pandemic crisis. Characterizing “last resort” split-circuit ventilation raises complex considerations and questions needing urgent answers to provide “guard rails” parameters for reasonably safe use.

Computational simulation provides an opportunity to quickly develop guidance over a wide range of possible clinical scenarios, without incurring patient risk. While empirical reports have been published [19], we have found only limited reports of empirical human simulation (using masks on volunteers and hyperventilation to avoid spontaneous breathing dyssynchrony) [24], and no computational simulation literature to inform multiple patient shared ventilation, nor specific definition of “similar mechanical support needs” [17]. Thus, it is clear that a method is needed to study and navigate this multi-parameter space and multiple scenarios in this complex system. We undertook a computational simulation of dual patient ventilation with a single ventilator using a computational whole-body simulation model to produce informed clinical parameters addressing noted assumptions and assertions. Use of whole-body physiological simulation exposes secondary effects that would not be apparent when modeling mechanics of ventilation alone.

The Pulse Physiology Engine (Kitware, Inc., Clifton Park, NY) open source software was previously designed and validated to model a single patient-ventilator pairing [25]. We modified and employed this model to simulate two or more patients sharing a single ventilator resource. This computational model permits rapid evaluation of a wide range of possible scenarios to educate clinicians about the potential risks for mismatched patients ventilated through uncompensated branched breathing circuits attached to a common ventilator. As a computational model without risk to patients, clinically untenable possibilities can be explored to determine guidance or “guard rails” to inform clinical practice, hopefully averting risky clinical practices and potentially poor outcomes.

Methods

Pulse overview

Concurrent pressure controlled ventilation of two patients in a configuration without flow-compensating resistance is being modeled by computational simulation. The open source Pulse Physiology Engine (pulse.kitware.com) models are built on the pioneering work of Ty Smith [26], further developed by the Department of Defense as BioGears (W81XWH-13-2-0068), and finally forked, continually developed, and supported for diverging uses as Pulse. Pulse is comprised of lumped-parameter models, which use electrical circuit analogues (e.g., resistors and capacitors) to represent the behavior of physiological regions or systems of the human body [25]. Feedback mechanisms and interactions of systems are accomplished through circuit connections or scalable circuit elements [27]. These circuit models are then solved using transient circuit analysis. Substances can be circulated through the models using pharmacokinetic (PK) and pharmacodynamic (PD) models, which use the physiologic properties of the patient and the physicochemical properties of the drug in differential equations to represent drug diffusion and distribution in the body [28]. Disease and treatment models are designed with differential equations describing the effects of disease and treatment and then applied to the lumped-parameter models affecting the overall calculations. Thus, the Pulse Physiology Engine models patient physiology interacting with medical equipment. The Pulse engine has a rigorous validation standard driven by data from the literature, patient data, and our relationships with clinical experts in a number of medical fields. More details on the transient circuit analysis, substance transport, and validation process can be found in the supplemental data published on the Pulse website (https://pulse.kitware.com/_docs.html). For this study, the Pulse ventilator implementation has been improved. The Pulse lung models were also modified and extended to enable simultaneous simulation of more than one patient interacting with a single ventilator.

Patient model

Recent studies of laboratory-confirmed COVID-19 pneumonia patients that received invasive mechanical ventilation have shown the PaO2/FiO2 ratios were consistent with the Berlin criteria for moderate-to-severe Acute Respiratory Distress Syndrome (ARDS) [29, 30]. However, in the earliest reports some COVID-19 patients seemed to present an atypical form of ARDS with dissociation between well-preserved lung mechanics (i.e., compliance) and unusually severe pulmonary shunt fraction with consequent hypoxemia [31, 32]. Subsequently others have reported typical ARDS respiratory mechanics and gas diffusion parameters for COVID-19 ARDS patients [33, 34].

Based on these data, the existing Pulse ARDS pathophysiology model was used. Additionally, the patient airway resistance, compliance, alveolar surface area, and pulmonary shunt fraction can be modified individually or by specifying a lumped-parameter patient disease “severity.” The severity is defined between zero (no change from the healthy state), and one (life-threatening). Severity values of 0.3, 0.6, and 0.9 correspond to accepted criteria for mild, moderate, and severe ARDS, respectively, including PaO2/FiO2 and shunt fraction changes. To model COVID-19, the ARDS disease severity mappings that represent the combined effects of the reduction in alveolar surface area and increase in pulmonary shunt were used. The combined changes are defined here as the diffusion impairment factor (DIF). The DIF selectively hinders gas exchange and emulates hypoxemic respiratory failure.

Using Pulse’s whole-body modeling approach, the patient is pharmacologically paralyzed simulating use of a neuromuscular transmission blocking agent and endotracheally intubated, consistent with current clinical guidelines [17]. The patient’s computed diseased pulmonary steady state is achieved while also evaluating other effects on the Pulse whole-body model. Computational steps simulate every 2 milliseconds. If a steady state cannot be achieved, those model parameters are rejected.

Multi-patient ventilator model

The Pulse ventilator model was first developed and tested with single-patient simulations and has been validated through pure simulation and hybrid testing with physical systems [35]. For this study, the Pulse ventilator model was executed in Pressure Control Continuous Mandatory Ventilation (PC-CMV) mode, following existing multi-patient ventilation guidelines [17, 21]. A square driving pressure waveform was specified by user-defined Positive End Expiratory Pressure (PEEP), Peak Inspiratory Pressure (PIP), respiratory rate (RR), and Inspiratory: Expiratory (IE) ratio. The ventilator gas fractions, such as Fractional inspired oxygen (FiO2), were also specified. The inspiratory, expiratory, and endotracheal tubes are included as part of the lumped-parameter circuit in Pulse. These comprise a closed-loop circuit and a branch for each patient and the ventilator for pneumatic and substance transport analysis. All tubing was estimated to be 3 feet long with a 22 mm inner diameter. The associated tubing resistance is negligible compared to the patient’s tracheobronchial resistance. One-way (check) valves are present to prevent backflow.

Previously, the Pulse engine represented a single patient and the associated equipment. To simulate a single ventilator connected simultaneously to two or more Pulse patients, required the creation of a new multi-patient ventilation engine. The new Pulse multi-patient engine simulates multiple patient physiology engines in lockstep to compute the effects of unequal, and potentially dynamic, differences in patient breathing mechanics when connected in parallel branched breathing circuits subject to a single pressure-mode ventilator. Fig 1 shows an overview of the multi-patient simulation environment.

thumbnail
Fig 1. The approach for simulating multi-patient ventilation.

An entire system state is calculated every 2 ms. The Pulse dynamic circuit solver and transporter are leveraged to ensure sound physics-based results with conservation of energy and mass. Mechanistic interactions occur with all other Pulse physiological systems, most notably, the alveolar-capillary partial pressure gradient diffusion gas exchange with the cardiovascular system.

https://doi.org/10.1371/journal.pone.0242532.g001

Multi-patient ventilation simulation design

COVID-19 respiratory mechanics parameters were gleaned from the available literature. An initial in silico investigation (computational simulation) used a sparse sampling of the combined mechanical ventilator and COVID-19 patient parameter space to understand patient outcome patterns. These simulations combined various ventilator PEEP settings with patient compliance and disease severity values. Tidal Volume (TV) and arterial partial pressure of oxygen (PaO2) were evaluated. This simulation was implemented as a two-step process: 1) patients were simulated with various levels of disease states connected to an individually dedicated ventilator with increasing FiO2 values until a homeostatic point (pulse-oximetric oxygen saturation, SpO2 > 89%) was reached; then 2) simulated patients were paired for three separate simulations using each patient’s individual ventilator settings and average values for the two paired patients. In the initial model hyperoxia was considered an undesirable, exclusionary outcome, judged to be acceptable for subsequent simulation criteria.

A scoring methodology was developed to analyze the simulation results and to identify safe clinical “guard rails” for patient pairing and specifying ventilator settings. The patient pair outcome was scored as follows:

  • Positive (green)
    • TV (mL/kg)—both patients between 5.5 and 6.5
    • SpO2 (%)—both patients above 89
    • PaO2 (mmHg)—both patients below 150
  • Less Positive (yellow)
    • TV (mL/kg)—both patients between 4.5 and 7.5
    • SpO2 (%)—both patients above 89
    • PaO2 (mmHg)—both patients between 150 and 200 (hyperoxemia)
  • Negative (red)
    • TV (mL/kg)—either patient below 4.5 (underinflation) or above 7.5 (overinflation)
    • SpO2 (%)—either patient under 89 (severe hypoxia)
    • PaO2 (mmHg)—either patient above 200 (oxygen toxicity)

If all parameters are in the target range, the result is scored green. If a patient has a tidal volume or SpO2 parameter in the yellow or red criteria, the patient is scored in the lowest scoring category. Only in the initial analysis was hyperoxia scored yellow or red.

Results using the first model specification informed a refined discretization of the parameter space in the final simulation study. The initial experiments revealed that the PIP and FiO2 should be specified as values providing the best chance of positive outcomes for both patients, as opposed to using average values or those of a single patient. Therefore, the ideal drive pressure (PIP—PEEP) for each pair of patients was directly calculated based on a target tidal volume (6 mL/kg) and known average static respiratory compliance (Cstat). The FiO2 was determined through simulation by finding the lowest value between 0.21 and 1.0 that stabilized both patients to a SpO2 of at least 89%. In this final analysis, hyperoxia (oxygen toxicity risk) alone was not scored yellow or red, being considered a likely acceptable expression of good or improving gas exchange.

All possible combinations of patient parameters for the final model required evaluation of 12,642 unique patient and ventilator combinations defined by the following parameter space:

  • Mechanical ventilator (one per simulation):
    • Respiratory rate (bpm)—Fixed at 20
    • I:E ratio—Fixed at 1:2
    • PEEP (cmH2O)—range of 10 to 20 in increments of 5
    • PIP (cmH2O)—Derived from 6 mL/kg given PEEP
    • FiO2—Derived through simulation for both patients’ SpO2 > 90% (if possible)
  • Patient (two per simulation):
    • Total respiratory resistance (cmH2O-s/L)—Fixed at 5
    • Total respiratory static compliance (mL/cmH2O)—a range of 10 to 50 in increments of 1
    • DIF (severity 0 to 1)—0.3 to 0.9 in increment of 0.1

Two equal-sized 70 kg patients are modeled. This should not limit the applicability of these results because pressure mode ventilation does not enforce delivery of a set volume, which could be injurious. Only a single respiratory resistance per patient branch was modeled because clinical variations of airway resistance have not been found to be a prominent characteristic of COVID-19 ARDS [18, 29, 31, 33, 34].

Results

Patient simulation

This Pulse COVID-19 model is an extension of the existing ARDS methodology that has been validated with referenced empirical data and trends [36]. Scenarios with the Pulse standard, properly ventilated ARDS patient shows good agreement with expected outcomes [37, 38]: severities matching mild, moderate, and severe cases resulting in PaO2/FiO2 of 340, 130, and 50 mmHg, respectively, and shunt fractions of 8%, 27%, and 59%, respectively. Validated ARDS and mechanical ventilator models, with the added ability to specify the patient’s respiratory compliance, allows for reasonable COVID-19 pathophysiology simulations.

Each of the 287 unique COVID-19 patients created (based on DIF and compliance combinations) for the full study were simulated as single patients with their own ventilators to ensure they were capable of achieving 89% SpO2 (green outcomes) at maximum FiO2. Those patients that were unable to meet this criterion were excluded from further multi-patient ventilation analysis to prevent skewing resulting guidelines. S1 Table in S1 File shows that 71% of patients with 0.8 DIF and 100% of patients with 0.9 DIF were disqualified for shared ventilator simulation.

Multi-patient ventilation simulation

Results of an example paired patient simulation are shown in Fig 2. The TV is significantly different because of the Cstat mismatch (plot a). The PaO2 values and response to oxygenation differ largely from the DIF (plot b). Because CO2 transits tissue much more readily than oxygen, the DIF affects O2 more severely than CO2 (plot b).

thumbnail
Fig 2. Select outputs from an example multi-patient ventilation scenario from the initial investigation.

The outcome of this scenario was classified as negative (red) because both tidal volumes (as mL/kg) are outside the desired bounds for lung protective ventilation.

https://doi.org/10.1371/journal.pone.0242532.g002

Parameter differences between patient pairs were evaluated using the scoring methodology described previously. The absolute value of the difference of specific parameters from patient one and patient two was analyzed using a one-way analysis of variance (ANOVA) to determine relative correlations to outcomes. An F-test was used to statistically test the equality of means, with larger η2 values denoting higher correlation, as shown in Table 1. Five parameters were found to be most correlated to outcomes: Cstat, End-tidal CO2 (EtCO2), Alveolar-arterial (A-a) gradient, oxygen saturation index (OSI), and oxygen index (OI). Three of these parameters can be evaluated non-invasively for patients: Cstat, EtCO2, and OSI.

thumbnail
Table 1. The results of an outcome correlation statistical analysis.

https://doi.org/10.1371/journal.pone.0242532.t001

Parameter independence was assessed using the Pearson product-moment correlation coefficient (PPMCC) method. The results in Fig 3 show that two of the four correlated parameters are both non-invasive and independent, Cstat, and OSI. Therefore, these two parameters were plotted against simulated clinical outcomes in Fig 4.

thumbnail
Fig 3. The selected parameters for investigation are compared with each other to determine their dependence.

The PPMCC method is used to calculate a value between -1 (inversely correlated) and 1 (correlated). Those with low correction (close to 0) are more independent of each other and are therefore the best candidates for informed decision-making.

https://doi.org/10.1371/journal.pone.0242532.g003

thumbnail
Fig 4. Comparison of multi-patient ventilation simulated outcomes due to TV (plot a) and PaO2 (plot b) outcome bounds.

Each graphical dash is a full simulation. Included are univariate histogram plots for each axis using kernel density estimation to represent the distribution of all three outcomes described in two dimensions. The compliance (abscissa) has discrete values due to the chosen patient model parameter setting methodology and fluid mechanics. The OSI (ordinate) is dependent on all external settings, along with the complex interactions of internal mechanistic models. Note that while the OSI has units of mmHg (because it is a ratio of pressure divided by saturation), the interpretation is like the unitless OI value.

https://doi.org/10.1371/journal.pone.0242532.g004

While respiratory compliance is a direct indicator of outcome based on lung recruitment, Fig 5 shows that OSI is an important and effective measure of overall diffusion impairment.

thumbnail
Fig 5. Distribution of all simulated patient’s OSI grouped by the seven DIF settings used from mild to severe.

The OSI increases with DIF and is, therefore, useful for non-invasive clinical assessment of hypoxemic respiratory failure. The OSI increases with diffusion impairment because the SpO2, which plateaus at 100, is a proportionately larger fraction of PaO2 as diffusion impairment and shunt make PaO2 less than expected for a given FiO2.

https://doi.org/10.1371/journal.pone.0242532.g005

The combined outcomes as a function of the patient parameter mismatches are used to generate a complete decision matrix in Fig 6. The simulation results show that patients with similar respiratory compliances and comparable OSI are most likely to have satisfactory outcomes when paired to a single ventilator.

thumbnail
Fig 6. The simulations from Fig 4 were holistically taken into account to get a complete decision matrix.

Outcomes were assigned a normalized value of green = 1, yellow = 0.5, and red = 0 to encode a z-axis as colors or color gradient. The resulting three-dimensional scatter plot (plot a) was used to produce an interpolated surface using the first-order bivariate B-spline method (plot b).

https://doi.org/10.1371/journal.pone.0242532.g006

An interactive version of the model is available via a Jupyter notebook in the Pulse repository (https://gitlab.kitware.com/physiology/jupyter).

Discussion

The results demonstrate that lung compliance and OSI are useful clinical guidance parameters for multi-patient ventilation, as they are non-invasively measurable and highly correlated to outcomes. OSI has also been shown to be a significant indicator of clinical outcomes for ARDS [39]. Our simulations show that a difference up to approximately 12 mL/cmH2O in lung compliance and an OSI difference of less than 2 mmHg anticipate satisfactory patient outcomes. (While the OSI has nominal units of mmHg as an artifact of its calculation, it is really to be interpreted as a ratio of oxygen partial pressure to hemoglobin saturation, which is unitless. Conveniently, oxyhemoglobin saturation of 100% typically corresponds to about 100 mmHg oxygen partial pressure for normal lungs.) The range of simulated parameter variation is more liberal than the conservative protocol of Beitler et al. [18] With a simple circuit, there is less flexibility than may be achieved with the more complex implementation proposed by Han et al. [23] This simulation may inform clinical guidelines when pairing patients on a single ventilator. An obvious limitation of this conclusion is that this is a pure simulation study, and we do not have correlated patient data. Planned future work includes comparing our computational simulations to physical experiments using the simple lung simulators with discrete selectable settings for resistance, compliance, and effects of negative inspiratory pressure patient-initiated breath triggers [IngMar Medical, QuickLung® Adult Precision Test Lung with QuickTrigger, part number 15 00 100, www.ingmarmed.com]. Parameters of this simulation model were chosen to match nominal QuickLung® parameters to facilitate future bench testing. Clinical validation may never be possible, would likely be observational at best, and is very unlikely to undergo ethical randomized clinical trial. This is a value of simulation, that possibilities can be explored that are beyond the realm of clinical investigation, yet of clinical value.

Implementation of split breathing circuits without compensating resistances was the specific model addressed in the joint statement by professional societies recommending against this practice [13]. Various means have been employed empirically to limit flow to the most compliant lungs [1923]. Modeling flow restriction adds model parameters, greatly increasing the number of simulations needed. For those locales with the most rudimentary resources, flow restriction devices may not be possible. Future work will also investigate sharing ventilators by more than two patients, and the use of flow restrictors in circuit branches to facilitate separate control of tidal volume delivered to individual patients during multi-patient ventilation.

If there is an adequate supply of ventilators available, individual ventilators are preferable and optimal [5, 7, 12]. Numerous efforts to produce quantities of individual ventilators have been kindled by an anticipated surge of patients needing ventilator support. There has been a surge of creative projects, yet there are gaps to be bridged before mass production could be a reality [40]. The FDA has issued Emergency Use Authorizations (EUA) for some innovative ventilators not made by usual medical equipment vendors [41]. Whether a surge of supply will be needed or available remains to be seen.

However, even with numerous creative projects creating basic ventilators from readily available components, there may not be enough locally available in a surge emergency. The primary default is manual ventilation, usually with bag-valve-mask (BVM) devices. Prolonged manual ventilation is difficult to maintain. [42] Split breathing circuits can be quickly assembled from available breathing circuit components for a transition to mechanical ventilation. Ventilation equipment not usually employed for critical care may also be pressed into service in a crisis, e.g., CPAP and BiPAP as ventilators, anesthesia machines, and high-flow non-invasive ventilation devices.

Multi-patient ventilation may be the only available option for poorly resourced regions. Manufacturing of single patient ventilators, however primitive, is an unlikely prospect for countries with limited resources. A personal communication from a Sudanese engineer interested in producing ventilators in-country stated that the entire country has only 80 ventilators and very limited materiel or manufacturing resources. Whatever can be done with resources already present in the country is their only available option (Conversation with MS, 28 Mar 2020, SMP). Similar limitations also apply for very populous countries, such as India, though they may be able to rally more resources for high throughput manufacturing. For these critical scenarios, we intend that these simulations could inform clinical guidelines for patient pairing on a single ventilator.

Under extreme circumstances, and in resource-challenged regions limited by the availability of an adequate number of ventilators, multi-patient ventilation is a potentially viable option and can significantly increase the capacity to care for critically ill patients in surge scenarios. Our computational simulation study provides ranges of basic physiological parameters that could be used in patient selection for future studies evaluating multi-patient ventilation. As a trade-off under urgent conditions, this simulation suggests that less stringent matching criteria may be tolerable without imposing a very deliberate and conservative protocol [18]. However, implementation would be less tolerant of patient variations and faults than the more complicated arrangement of partially independent patient breathing circuits [23]. Our simulation could inform clinicians acutely needing to support two patients with a single ventilator. Further physical validation of these simulation results is needed.

Inevitably, ethical issues arise. Some will question whether it is ethical to undertake a desperate ventilator sharing tactic without a well defined exit strategy. This is not unlike many clinical interventions that secure time for informed triage and more deliberate consideration of longer term options. However, sharing a ventilator has implications for more than a single patient. Response to an epidemic challenge may foster innovation with broad implications, for instance extensive use of manual ventilation of polio patients in 1952 that kindled vigorous development of mechanical ventilators [43] and the discipline of critical care. We regard uncompensated ventilator sharing as a desperate “entrance strategy”. The underlying question is whether it is unethical to maintain ventilation with imperfect apparatus as an alternative to not providing ventilation at all, despite not having a long-term exit strategy. This conundrum is worthy of careful and nuanced discussion, yet beyond the scope of this technical study. As in Las Vegas, this approach might be mitigation while single-patient ventilators are in transit from a remote location. It could be an intermediate or longer term intervention. These are unexplored considerations forced upon us by advent of a global pandemic. The variety of possible surge scenarios defies adequate planning for every contingency. [42, 4446] Perfect can be the enemy of good. The best possible action, though not optimal, is not necessarily futile. [1, 2, 47]

The simulation software, data, and a UI for exploring the simulations is available open-source at the Kitware GitLab site (https://gitlab.kitware.com/physiology). For more details, see the supplement.

Conclusion

Evidence before this study

If numbers of patients requiring mechanical ventilation exceed the number of available ventilators in a surge, shared branched ventilator circuits have been proposed for sharing one ventilator by multiple patients. Under pressure of imminent pandemic surge, only small and rudimentary clinical studies were available. Testing over expected ranges of lung-chest wall compliance were not found. Increasing clinical experiences of mechanical ventilation parameters employed for COVID-19 patients are being reported.

Added value of this study

The number of possible combinations of ventilation and physiological parameters is very large. Time and resource constraints do not permit conventional research. Computational simulation provides rapid sensitivity evaluation of several factors over a wide range of hypothetical ventilation conditions. Envelopes of evaluated parameters may provide reasonably estimated safety boundaries for clinicians compelled in an emergency surge to employ a poorly characterized practice. A previously well-vetted computational model for ventilation of a single patient by a dedicated ventilator has been modified to model the sharing of a single ventilator by two or more patients. Only pairings of two equally sized 70 kg patients are modeled in this report. These simulations provide estimates of effects on ventilation and blood oxygenation by clinically measurable values using conceivable mismatched patient lung compliance and oxygenation (diffusion and shunt).

Implications of all the available evidence

These estimates are for pressure mode ventilation using a single ventilator shared by branched breathing apparatus for two patients. Individual patient flow restriction to compensate for compliance mismatch is not considered in this simulation. Reasonable though arbitrary bounds of acceptable parameters may guide clinicians when evaluating candidate pairings of patients with different physiological characteristics. Further laboratory testing and clinical experience will be needed to determine the validity or utility of these assessments. Different simulations will be needed for flow-compensated branches, more than two patients, and unmatched body habitus.

References

  1. 1. White DB, Lo B. A framework for rationing ventilators and critical care beds during the COVID-19 pandemic. Jama. 2020;. pmid:32219367
  2. 2. Truog RD, Mitchell C, Daley GQ. The toughest triage—allocating ventilators in a pandemic. New England Journal of Medicine. 2020;. pmid:32202721
  3. 3. Rubinson L, Hick JL, Curtis JR, Branson RD, Burns S, Christian MD, et al. Definitive care for the critically ill during a disaster: Medical resources for surge capacity: From a Task Force for Mass Critical Care summit meeting, January 26–27, 2007, Chicago, IL. Chest. 2008;133(5):32S–50S. pmid:18460505
  4. 4. Sommer DD, Fisher JA, Ramcharan V, Marshall S, Vidic DM. Improvised automatic lung ventilation for unanticipated emergencies. Critical care medicine. 1994;22(4):705–709.
  5. 5. Iscoe S, Fisher J. Still better multi-patient ventilation. Resuscitation. 2008;79(1):173–174.
  6. 6. Neyman G, Irvin CB. A single ventilator for multiple simulated patients to meet disaster surge. Academic emergency medicine. 2006;13(11):1246–1249.
  7. 7. Branson RD. A single ventilator for multiple simulated patients to meet disaster surge. Academic Emergency Medicine. 2006;13(12):1352–1353.
  8. 8. Tonetti T, Zanella A, Pizzilli G, Irvin Babcock C, Venturi S, Nava S, et al. One ventilator for two patients: feasibility and considerations of a last resort solution in case of equipment shortage. Thorax. 2020;75(6):517–519. pmid:32327566
  9. 9. Chatburn RL, Branson RD, Hatipoğlu U. Multiplex Ventilation: A Simulation-based Study of Ventilating Two Patients with One Ventilator. Respiratory Care. 2020;.
  10. 10. Herrmann J, da Cruz AF, Hawley ML, Branson RD, Kaczka DW. Shared Ventilation in the Era of COVID-19: A Theoretical Consideration of the Dangers and Potential Solutions. Respiratory Care. 2020;. pmid:32376612
  11. 11. Paladino L, Silverberg M, Charchaflieh JG, Eason JK, Wright BJ, Palamidessi N, et al. Increasing ventilator surge capacity in disasters: ventilation of four adult-human-sized sheep on a single ventilator with a modified circuit. Resuscitation. 2008;77(1):121–126. pmid:18164798
  12. 12. Branson RD, Rubinson L. One ventilator multiple patients—What the data really supports. Resuscitation. 2008;79(1):171–172.
  13. 13. SCCM, AARC, ASA, APSF, AACN, CHEST. Joint Statement on Multiple Patients Per Ventilator; Last accessed 2020-09-18. https://www.apsf.org/wp-content/uploads/news-updates/2020/Multiple-Patients-Single-Vent-Statement.pdf.
  14. 14. Babcock C, Franco R, Bunting L, Paladino L, Habashi N, Kaplan L, et al. Co-Ventilating Patients During a Critical Ventilator Shortage: A Method for Implementation; Last accessed 2020-09-18. https://sccm.org/getattachment/Disaster/FEMA-Report.pdf.
  15. 15. Beitler J, Kallet R, Kacmarek R, Branson R, Brodie D, Mittel A, et al. Ventilator Sharing Protocol: Dual-Patient Ventilation with a Single Mechanical Ventilator for Use during Critical Ventilator Shortages; Last accessed 2020-09-17. https://www.gnyha.org/wp-content/uploads/2020/03/Ventilator-Sharing-Protocol-Dual-Patient-Ventilation-with-a-Single-Mechanical-Ventilator-for-Use-during-Critical-Ventilator-Shortages.pdf.
  16. 16. Chen M, Shah A, Shah R, Zhou G, Kane E, Burnett G, et al. Mount Sinai Health System Emergency Ventilator Sharing Protocol and Instructions for Use; Last accessed 2020-09-18. https://docs.google.com/document/d/1eFxzKmt__0sbYAjFNiQxVRCuIv-a65mbs4QHzIu_1i0/edit.
  17. 17. Giroir B, Adams J. Optimizing Ventilator Use during the COVID-19 Pandemic; Last accessed 2020-09-18. https://www.hhs.gov/sites/default/files/optimizing-ventilator-use-during-covid19-pandemic.pdf.
  18. 18. Beitler J, Mittel A, Kallet R, Kacmarek R, Hess D, Branson R, et al. Ventilator Sharing during an Acute Shortage Caused by the COVID-19 Pandemic. Am J Respir Crit Care Med. 2020;202(4):600–604. pmid:32515988
  19. 19. Levin M, Chen MD, Shah A, Shah R, Zhou G, Kane E, et al. Differential ventilation using flow control valves as a potential bridge to full ventilatory support during the COVID-19 crisis. medRxiv. 2020;. pmid:32639236
  20. 20. Group DMIW. DIFFERENTIAL MULTIVENTILATION INTERNATIONAL WORKING GROUP OVERVIEW; Last accessed 2020-09-18. https://www.differentialmultivent.org/.
  21. 21. Raredon MSB, Fisher C, Heerdt P, Deshpande R, Nivison S, Fajardo E, et al. Pressure-Regulated Ventilator Splitting (PReVentS): A COVID-19 Response Paradigm from Yale University. medRxiv. 2020;.
  22. 22. Clarke A, Stephens A, Liao S, Byrne T, Gregory S. Coping with COVID-19: ventilator splitting with differential driving pressures using standard hospital equipment. Anaesthesia. 2020;. pmid:32271942
  23. 23. Han J, Mashari A, Singh D, Dianti J, Goligher E, Long M, et al. Personalized Ventilation to Multiple Patients Using a Single Ventilator: Description and Proof of Concept. Crit Care Explor. 2020;2(5):e0118. pmid:32671348
  24. 24. Smith R, Brown J. Simultaneous ventilation of two healthy subjects with a single ventilator. Resuscitation. 2009;80(9):1087.
  25. 25. Bray A, Webb JB, Enquobahrie A, Vicory J, Heneghan J, Hubal R, et al. Pulse physiology engine: an open-source software platform for computational modeling of human medical simulation. SN Comprehensive Clinical Medicine. 2019;1(5):362–377.
  26. 26. Fukui Y, Smith NT. Interactions among Ventilation, the Circulation and the Uptake and Distribution of Halothane—Use of a Hybrid Computer Multiple ModelI. The Basic Model. Anesthesiology. 1981;54(2):107–118.
  27. 27. Clipp R, Thames M, Webb J, Metoyer R, Swarm Z, Aaron B, et al. Integration of a baroreflex model into a whole body physiology engine. In: Summer Biomechanics, Bioengineering, and Biotransport Conference; 2016.
  28. 28. Clipp R, Bray A, Metoyer R, Thames M, Webb J. Pharmacokinetic and pharmacodynamic modeling in BioGears. In: Conf Proc IEEE Eng Med Biol Soc. vol. 2016; 2016. p. 1467–1470.
  29. 29. Bhatraju PK, Ghassemieh BJ, Nichols M, Kim R, Jerome KR, Nalla AK, et al. Covid-19 in critically ill patients in the Seattle region—case series. New England Journal of Medicine. 2020;. pmid:32227758
  30. 30. Yang X, Yu Y, Xu J, Shu H, Liu H, Wu Y, et al. Clinical course and outcomes of critically ill patients with SARS-CoV-2 pneumonia in Wuhan, China: a single-centered, retrospective, observational study. The Lancet Respiratory Medicine. 2020;. pmid:32105632
  31. 31. Gattinoni L, Chiumello D, Caironi P, Busana M, Romitti F, Brazzi L, et al. COVID-19 pneumonia: different respiratory treatments for different phenotypes? Intensive Care Medicine. 2020; p. 1.
  32. 32. Gattinoni L, Chiumello D, Rossi S. COVID-19 pneumonia: ARDS or not? Crit Care. 2020;24(1):154.
  33. 33. Ziehr D, Alladina J, Petri C, Maley J, Moskowitz A, Medoff B, et al. Respiratory Pathophysiology of Mechanically Ventilated Patients with COVID-19: A Cohort Study. Am J Respir Crit Care Med. 2020;201(12):1560–1564. pmid:32348678
  34. 34. Haudebourg A, Perier F, Tuffet S, de Prost N, Razazi K, Mekontso Dessap A, et al. Respiratory Mechanics of COVID-19- versus Non-COVID-19-associated Acute Respiratory Distress Syndrome. Am J Respir Crit Care Med. 2020;202(2):287–290. pmid:32479162
  35. 35. Rodriquez Jr D, Blakeman T, Gomaa D, Branson R. Advancing Mechanical Ventilation Management through Simulation; Last accessed accessed 2020-05-05. https://apps.dtic.mil/dtic/tr/fulltext/u2/1074448.pdf.
  36. 36. Webb JB, Bray A, Clipp RB. Parameterization of Respiratory Physiology and Pathophysiology for Real-Time Simulation. In: 2020 42nd Annual International Conference of the IEEE Engineering in Medicine & Biology Society (EMBC). Montreal, QC, Canada: IEEE; 2020. p. 2274–2278.
  37. 37. Force ADT, Ranieri V, Rubenfeld G, Thompson B, Ferguson N, Caldwell E, et al. Acute respiratory distress syndrome. Jama. 2012;307(23):2526–2533.
  38. 38. Radermacher P, Maggiore SM, Mercat A. Fifty years of research in ARDS. Gas exchange in acute respiratory distress syndrome. American journal of respiratory and critical care medicine. 2017;196(8):964–984.
  39. 39. DesPrez K, McNeil JB, Wang C, Bastarache JA, Shaver CM, Ware LB. Oxygenation saturation index predicts clinical outcomes in ARDS. Chest. 2017;152(6):1151–1158.
  40. 40. Pearce JM. A review of open source ventilators for COVID-19 and future pandemics. F1000Research. 2020;9.
  41. 41. FDA. Emergency Use Authorizations; 2020. https://www.fda.gov/medical-devices/emergencysituations-medical-devices/emergency-use-authorizations#covid19ventilators.
  42. 42. Branson R, Johannigman J, Daugherty E, Rubinson L. Surge capacity mechanical ventilation. Respir Care. 2008;53(1):78–88; discussion88.
  43. 43. Slutsky A. History of Mechanical Ventilation. From Vesalius to Ventilator-induced Lung Injury. Am J Respir Crit Care Med. 2015;191(10):1106–1115.
  44. 44. Hick J, Einav S, Hanfling D, Kissoon N, Dichter J, Devereaux A, et al. Surge capacity principles: care of the critically ill and injured during pandemics and disasters: CHEST consensus statement. Chest. 2014;146(4 Suppl):e1S–e16S. pmid:25144334
  45. 45. Christian M, Devereaux A, Dichter J, Rubinson L, Kissoon N, Task FfMCC, et al. Introduction and executive summary: care of the critically ill and injured during pandemics and disasters: CHEST consensus statement. Chest. 2014;146(4 Suppl):8S–34S. pmid:25144202
  46. 46. Einav S, Hick J, Hanfling D, Erstad B, Toner E, Branson R, et al. Surge capacity logistics: care of the critically ill and injured during pandemics and disasters: CHEST consensus statement. Chest. 2014;146(4 Suppl):e17S–43S. pmid:25144407
  47. 47. Biddison L, Berkowitz K, Courtney B, De Jong C, Devereaux A, Kissoon N, et al. Ethical considerations: care of the critically ill and injured during pandemics and disasters: CHEST consensus statement. Chest. 2014;146(4 Suppl):e145S–55S. pmid:25144262