| Literature DB >> 34185796 |
Michael Saidani1, Harrison Kim1, Jinju Kim1.
Abstract
Providing sufficient testing capacities and accurate results in a time-efficient way are essential to prevent the spread and lower the curve of a health crisis, such as the COVID-19 pandemic. In line with recent research investigating how simulation-based models and tools could contribute to mitigating the impact of COVID-19, a discrete event simulation model is developed to design optimal saliva-based COVID-19 testing stations performing sensitive, non-invasive, and rapid-result RT-qPCR tests processing. This model aims to determine the adequate number of machines and operators required, as well as their allocation at different workstations, according to the resources available and the rate of samples to be tested per day. The model has been built and experienced using actual data and processes implemented on-campus at the University of Illinois at Urbana-Champaign, where an average of around 10,000 samples needed to be processed on a daily basis, representing at the end of August 2020 more than 2% of all the COVID-19 tests performed per day in the USA. It helped identify specific bottlenecks and associated areas of improvement in the process to save human resources and time. Practically, the overall approach, including the proposed modular discrete event simulation model, can easily be reused or modified to fit other contexts where local COVID-19 testing stations have to be implemented or optimized. It could notably support on-site managers and decision-makers in dimensioning testing stations by allocating the appropriate type and quantity of resources.Entities:
Mesh:
Substances:
Year: 2021 PMID: 34185796 PMCID: PMC8241042 DOI: 10.1371/journal.pone.0253869
Source DB: PubMed Journal: PLoS One ISSN: 1932-6203 Impact factor: 3.240
Fig 1COVID-19 testing workflow at the University of Illinois at Urbana-Champaign (UIUC).
Application of the STRESS checklist [38] to the present DES model.
| Category | Checklist item | Present simulation model |
|---|---|---|
| Objectives | Purpose of the model | Designing better COVID-19 testing stations |
| Model outputs | Number of vials being processed on a daily basis | |
| Experimentation aims | Testing different configuration (in terms of operators and machines number and allocation) | |
| Logic | Base model overview diagram | Gantt diagram of the testing process ( |
| Base model logic | Flowchart of the COVID-19 testing process ( | |
| Scenario logic | Based on the hotspots (bottlenecks) identified | |
| Algorithms | Not applicable (N/A) | |
| Components | Number and allocation of operators and machines | |
| Data | Data sources | The SHIELD team of the University of Illinois |
| Input parameters | Time distribution and resources allocation ( | |
| Pre-processing | N/A | |
| Assumptions | Provided with the initial data by the experts from the SHIELD team (see values in | |
| Experimentation | Initialisation | Initial configuration provided by the SHIELD team. See the initial transient regime for the first batch in |
| Run-length | Two consecutive processing days (10 to 12 working hours per day) | |
| Estimation approach | Multiple replications (and box plots) for each scenario | |
| Implementation | Software | AnyLogic PLE |
| Random sampling | Triangular distribution function in AnyLogic (Monte Carlo simulation) | |
| Model execution | AnyLogic simulation engine FIFO (first in, first out) | |
| System specification | Intel Core i7-8550U, 1.80Ghz, 8.0GB RAM (Windows 10 Enterprise environment 64-bit) | |
| Code access | Computer model | Supplementary digital file (DES_model.alp) |
Fig 2Overview of the modeling approach.
Fig 3Workflow model of the COVID-19 testing process, used for developing the DES model.
Fig 4Developed DES model, with pools of resources and parameters to optimize.
Time distribution and resources allocation in the DES process.
| Task | Time distribution (seconds) | Resources | In DES model |
|---|---|---|---|
| Open bag and extract vial | triangular(4,5,6) | Operator | OpPrep |
| Attach labels and scan vial | triangular(24,30,36) | Operator | OpPrep |
| Place vial in rack | triangular(8,10,12) | Operator | OpPrep |
| Transfer vial rack to tank | triangular(24,30,36) | Operator | OpTran |
| Heat to 95°C | triangular(1800,1830,1860) | Machine | EqHeat |
| Transfer vial rack | triangular(24,30,36) | Operator | OpTran |
| Open tube and pipett to PCR tube | triangular(24,30,36) | Operator | OpTran |
| Load test tube rack into Biomek | triangular(24,30,36) | Operator | OpColl |
| constant(30) | Machine | EqBio | |
| Load 96 well plate into Biomek | triangular(24,30,36) | Operator | OpColl |
| constant(30) | Machine | EqBio | |
| Transfer to 96 well plate | constant(10) | Machine | EqBio |
| Unload and store plate | triangular(24,30,36) | Operator | OpColl |
| constant(30) | Machine | EqBio | |
| Discard test tube | triangular(96,120,144) | Operator | OpColl |
| Load 4 96 well plate | triangular(24,30,36) | Operator | OpLoad |
| constant(30) | Machine | EqBio | |
| Load 384 well plate | triangular(12,15,18) | Operator | OpLoad |
| constant(30) | Machine | EqBio | |
| Transfer to 384 well plate | constant(1800) | Machine | EqBio |
| Transfer to Vortex | triangular(24,30,36) | Operator | OpLoad |
| Vortex | constant(180) | Machine | EqCent |
| Centrifuge | constant(180) | Machine | EqCent |
| Transfer plate to QuantiStudio | triangular(48,60,72) | Operator | OpTest |
| constant(60) | Machine | EqTest | |
| Test-RT-qPCR | constant(5400) | Machine | EqTest |
| Output results and prepare for next batch | triangular(240,300,360) | Operator | OpTes |
| constant(120) | Machine | EqTest |
Fig 5Evaluation of scenarios through simulation runs (in established and continuous regime).
Fig 6Detailed box and whisker chart for key configurations (Sim. #4 to Sim. #7).
Fig 7Testing time of vials batches in transient and continuous operation.
Optimal resources allocation, as a function of the number of samples to be tested.
| Samples | OpPrep | OpTran | OpColl | OpLoad | OpTest | EqHeat | EqBio | EqCent | EqTest |
|---|---|---|---|---|---|---|---|---|---|
| 2000 | 3 | 2 | 1 | 1 | 2 | 1 | 1 | 1 | 2 |
| 4000 | 6 | 4 | 1 | 1 | 2 | 1 | 1 | 1 | 2 |
| 6000 | 9 | 7 | 1 | 1 | 3 | 1 | 2 | 1 | 3 |
| 8000 | 11 | 10 | 1 | 2 | 3 | 2 | 2 | 1 | 3 |
| 10000 | 13 | 12 | 1 | 2 | 4 | 2 | 3 | 1 | 4 |
| 12000 | 15 | 12 | 1 | 2 | 5 | 2 | 3 | 1 | 5 |
| 14000 | 19 | 13 | 2 | 3 | 6 | 2 | 3 | 1 | 6 |
| 16000 | 21 | 14 | 2 | 3 | 7 | 2 | 3 | 1 | 7 |
| 18000 | 23 | 16 | 2 | 3 | 7 | 3 | 3 | 1 | 7 |
| 20000 | 26 | 19 | 3 | 3 | 8 | 3 | 3 | 1 | 8 |
Fig 8Variability of the number of persons being tested on a daily basis (source: [43]).