key: cord-0862131-wksg1wdp authors: Pan, Shi-yu; Ding, Ming; Huang, Jing; Cai, Yan; Huang, Ying-zi title: Airway resistance variation correlates with prognosis of critically ill COVID-19 patients: a computational fluid dynamics study date: 2021-06-25 journal: Comput Methods Programs Biomed DOI: 10.1016/j.cmpb.2021.106257 sha: 3c500e15a0ba5b2bb22f86cd2c494cf7b8c291e4 doc_id: 862131 cord_uid: wksg1wdp OBJECTIVE: To evaluate the quantitative changes of respiratory functions for critically ill COVID-19 patients with mechanical ventilation, computational fluid dynamics (CFD) analysis was performed based on patient-specific three-dimensional airway geometry. METHODS: 37 cases of critically ill patients with COVID-19 admitted to the ICU of Huangshi Traditional Chinese Medicine Hospital from February 1st to March 20th, 2020 were retrospectively analyzed. 5 patients whose clinical data met the specific criteria were finally cataloged into death group (2 patients) and survival group (3 patients). The patient-specific three-dimensional airways were reconstructed from the central airways down to the 4th-5th bifurcation of the tracheobronchial tree. The volume changes of bronchi were calculated during the disease progression according to the comparison of two CT scans. Additionally, the changes of air flow resistance were analyzed using numerical simulation of CFD. RESULTS: Pearson correlation analysis demonstrated that there was negative correlation between the change of volume (ΔV) and the change of resistance (ΔR) for all COVID-19 patients (r=-0.7025). For total airway volume, an average decrease of -11.41± 15.71% was observed in death group compared to an average increase of 1.86± 10.80% in survival group (p=0.0232). For air flow through airways in lower lobe, the resistance increases for death group by 10.97± 77.66% and decreases for survival group by -45.49± 42.04% (p=0.0246). CONCLUSION: The variation of flow resistance in the airway could be used as a non-invasive functional evaluation for the prognosis and outcome of critically ill patients with COVID-19. The ‘virtual’ pulmonary function test by integrating follow-up CT scans with patient-derived CFD analysis could be a potentially powerful way in improving the efficiency of treatment for critically ill patients with COVID-19. Track back to December 2019, an unexpected outbreak of a highly contagious new coronavirus pneumonia has rapidly swept around the globe [1] [2] [3] [4] [5] [6] . By November 2020, the outbreak had grown to infect more than 57 million people with more than 1.3 million deaths in more than 200 countries/territories 7 . It is estimated that approximately 13.8% of infections developed into severe disease and 6.1% were critically ill 8 . The mortality rate of critically ill patients is reported to be over 37% 9 . Unfortunately, there is no evidence of any effective treatment for COVID-19 to date. Clinical control in critically ill cases is mainly based on supportive care including invasive mechanical ventilation, which can significantly improve the severe respiratory distress of these patients 10, 11 . However, the quantitative and effective evaluation of functional changes in patients with mechanical ventilation is hardly be assessed since the regular pulmonary function testing is not applicable to these critically ill patients. Computational fluid dynamics (CFD) is a technique that allows the simulation of flow patterns through three-dimensional models using a computational grid 12 . Based on the three-dimensional reconstruction of patient's CT images, CFD can be applied in the field of respiratory modeling where the flows inside the virtual airways are quantitatively analyzed and the influences of different morphology of airway on the flow pattern are investigated [13] [14] [15] [16] . In previous studies, CFD analysis has demonstrated special advantages compared with traditional assessment in management of COPD patients 17, 18 , calculation of drug distribution in airways [19] [20] [21] and therapeutic effect evaluation of asthmatics 22, 23 . In this study, we hypothesize that the patient-specific CFD airway analysis shows airway resistance variations that correlate with the development of the disease, and subsequently provides a useful approach for evaluation of respiratory functions and further prognosis of critically ill COVID-19 patients. To this end, retrospective investigation is taken where the morphological and CFD analysis for patient-specific airways result in a novel non-invasive functional assessment. In this study, a total of 37 cases of critically ill patients with COVID-19 admitted to the ICU of Huangshi Traditional Chinese Medicine Hospital from February 1st to March 20th, 2020 were retrospectively analyzed. On admission, the nucleic acid test results of all these patients were positive, and their CT suggested pneumonia. According to the diagnostic criteria of Chinese Clinical Guidance for COVID-19 Pneumonia Diagnosis and Treatment (7th edition) issued by the National Health Council (China) 24 , all cases included in this study were critically ill patients, and all of them met at least one of the following criteria: 1) respiratory failure occurred, and required mechanical ventilation; 2) shock occurred; 3) combined with other organ failure, and admission to the ICU. This study was approved by the Ethics Committee of Huangshi Traditional Chinese Medicine Hospital (HSZYPJ-2020-018-01), and all patients provided written informed consent. This study collected all clinical data and CT data of critical patients with COVID-19 during their hospitalization. The specific criteria for inclusion are as follows: (1) all patients included in this study have their CT data reflecting the dynamic changes from the initial stage to the improvement of the condition or death; 2) the interval between CT scans was about 10 days, and at least 2 times of the CT scans were available for analysis; 3) the CT data were in DICOM format, with the layer thickness less than 1.5mm; The exclusion criteria of this study include: 1) with no CT examination; 2) with no CT data at the initial stage; 3) their CT layer thickness didn't meet the criterion of the analysis; 4) the interval between CT scans were too long that the CT data could not reflect the entire dynamic process. Based on the aforementioned criteria, 5 critical patients with COVID-19 were finally cataloged into two groups, i.e., death group (Patients #1 & #2) and survival group (Patients #3, #4 & #5), for the further morphological and CFD analysis on individual airways. The schematic representation for inclusion is shown in Figure 1 Table As a prerequisite for the assessment of changes in airway volumes and CFD-determined resistance, three-dimensional reconstruction of patient-specific airway geometries was carried out based on the CT images. All CT scans were taken with a multi-slice scanner with 64 receptors, and all images had a pixel size of approximately 0.7mm. The DICOM (Digital Imaging and Communications in Medicine) images obtained by CT scans were then imported to a commercially developed software Mimics (Materialise, Leuven, Belgium) to implement the semi-automatic segmentation and three-dimensional reconstruction. Limited by the resolution of the CT scans, the constructed airway models started from the central airways down to the 4th-5th bifurcation of the tracheobronchial tree, where the airway diameter is around 1-2mm. The same procedure was applied to the pre and post CT images for constructing three-dimensional model. Subsequently, the airways of pre and post models were made topology corresponding and the endings were pruned to the same length, which is necessary for a comparison of airway volume and resistance changes. The airway models were subdivided in accordance with Ikeda 25 , and the volume was measured before and after the treatment. Homogeneous, incompressible, and Newtonian flow was assumed in the models. Since the influence of viscous effects could be neglected in the regions with highest resistance, i.e., the 4th-7th bifurcation, the air in the airway was considered as a laminar flow. Air was modeled to have a density of 1.225kg/m3. The airway wall was assumed to be rigid without displacement. The governing equations were Navier-Stokes equations in three-dimension, which are time-averaged equations of motion for fluid flow. Therefore, the flow can be described by In these equations, t represents time, is the fluid density, ⃗ is the velocity vector, p is the fluid pressure, and the viscous stress tensor is given by ̿ . To solve the flow numerically, meshing procedure was performed by ANSYS software to generate a computational grid with average mesh size of 3×10 6 tetrahedral cells. Mesh independent analysis was implemented to ensure enough grids were 4 male 48 174 72 survived 13 27 5 male 48 175 80 survived 46 33 constructed to capture variations of CFD parameters. After that, transient CFD simulation was processed by Fluent (ANSYS, 19 .0) and the pressure-based coupled solver based on SIMPLE algorithms was used to the solve the equations. The respiratory cycle was set as a period of 5.1s, which represents the breath frequency of about 12 times per minute for a normal adult. Besides, the tidal volume is set as 500ml 26 based on physical facts. Building time defined function, the time-dependent velocity curve was specified as the inlet boundary condition as shown in Figure 2 . At the outlets, the static pressure was uniformly defined to be zero pascal using an iterative process to reflect the internal flow distribution. The time step size was set to 0.051s. The convergence criterion was satisfied when the residual of continuity was less than 10 -4 . CFD allows the calculation of flow properties, including pressure, velocity, and mass flow rate inside the entire flow domain, which enables the determination of airway resistance. The resistance of certain branches was related with the pressure drop over the branches ∆p and the mass flow rate through the branches , where denotes to resistances of separate bronchial section 27 . The statistical software was SPSS 20.0 (SPSS Inc., Chicago, IL, USA). Data were expressed as median±SEM unless stated otherwise. The student's t-test was used for comparison of changes of volumes and resistance between death and survival groups. p<0.05 was considered statistically significant. The Pearson correlation coefficient was calculated to investigate the correlation between the changes of volume and the changes of resistance. To investigate the changes in the air flow characteristics during the disease progression, the CFD analysis were performed on all patients based on 3D airway models reconstructed by their CT images at two time points. A movie of change of velocity magnitude at mid-plane for Patient #2 during one respiratory cycle was provided in Supplementary Material. The cross-section profiles of velocity magnitude and pressure for two cases (Patient #2 from the death group and Patient #5 from the survival group) at mid-plane at 0.85 s are presented in Figure 3 . In the main bronchi, a slight increase of 3.18% and 1.23% in average flow velocity was observed for the death and survival cases respectively. However, the drop of pressure showed a decrease of -4.07 % for the survival case. This indicates that the resistance in the main bronchi decreased for the survival case and that the improved respiratory functions can be expected. On the contrast, it is worth noting that the drop of pressure for the death case increases to 2 times during the disease progression, suggesting that the respiratory function was getting worsen for this patient. The change of volume (∆V) and the change of resistance (∆R) of each bronchial section were averaged in the survival group and the death group, respectively. Table 2 provides a summary of ∆V and ∆R in the main bronchi and several bronchial sections for all patients. The opposite trend between ∆V and ∆R was found in most bronchial sections as well as the main bronchi. A statistical analysis (Pearson correlation analysis) between the results yielded significant negative correlation between ∆V and ∆R for all these COVID-19 patients (Figure 4 ; r=-0.7025, p=0.002). Figure 5 presents the changes in bronchial volume for all patients. In the survival group and the death group, the mean of the volume changes and resistance changes were measured. Table 3 provides an overview of the quantitative comparison between the death group and the survival group. The volume changes in total bronchial (ΔV total ) showed a mean decrease of -11.41±15.71% for death group and an increase of 1.86±10.80% for survival group (p= 0.0232 To the best of our knowledge, this is the first study to investigate the air flow changes in the airway for critically ill patients with COVID-19 by using CFD analysis. The analysis of retrospective cases demonstrated that the patients' prognosis was significantly related to the changes of airway volume and resistance. Since there are no registered drugs to treat COVID-19 disease currently, management is based mainly on supportive care including mechanical ventilation especially in severe cases. However, the changes of respiratory function after ventilation can hardly be detected quantitatively for critical ill patients. To this end, the current study presented a novel non-invasive analysis by comparing morphological changes and performing CFD calculation on patient-derived 3D airway models. We believe the integration with computational techniques such as CFD analysis could play a valuable role in improving the efficiency of treatment for critically ill patients with COVID-19. Among all the patients included in this study, the age of the two cases in the death group was significantly higher than that of the survival group, which is consistent with previous studies reported that age is a warning factor with strong predictive value in COVID-19 disease 28, 29 . Figure 6 shows the chest CT images during the disease progression for Patient #2 in the death group. [30] [31] [32] . All these anatomical and pathological features during the disease progression contribute to the changes in bronchial volume and resistance in the death cases. To be specific, the volumes of most bronchial sections were decreasing, and consequently an increase of the airflow resistance was observed in the death group (Table 2 & 3) . In contrast to the death group, a general increase of airway volume and a reduction of flow resistance were found for the three cases in the survival group (Table 2) . This CFD result was demonstrated to be consistent with the prognosis of these patients. For example, the CT images of Patient #3 in survival group (Figure 7) indicated that the exudation from the lung was obviously absorbed after effective treatment. The disease progression of Patient #5 was much more complicate. This patient was transferred to the ICU with mechanical ventilation seven days after symptoms developed. After eight days, this patient received non-invasive ventilation. However, four days later, he received mechanical ventilation again for another nine days due to the worsened condition. Therefore, his lung showed obvious uneven lesions with fibrosis and exudation caused by the relapse of illness. Correspondingly, the CFD analysis for Patient #5 (Table 2) showed that the volumes of several bronchi (B6R, B9L and B12L) even decrease during the progression, which may produce bad effect on the prognosis and recovery of this patient. It is noteworthy that both the reduced volume of the lower lobe bronchi (B6R and B9L) and the significantly increased resistance in the airway could be observed for the patients who maintained mechanical ventilation for a long time and eventually died. In addition, the difference of flow resistance between the death group and the survival group showed significant at the lower lung lobe (Table 3 ). This can be attributed to the trapped airway due to the deposition of inflammatory secretion in response to gravity after a long period of time lying supine 33, 34 . To this end, prone ventilation may be suggested for critically ill COVID-19 patients who require long duration of mechanical ventilation, in order to improve their lung ventilation or even their prognosis. The CFD analysis for airway flow regarding to prone ventilation need further studies in the future work. First, one of the major limitations of the study is its small sample size. CFD calculations were applied on the limited retrospective cases, therefore the findings of this study may or may not be generalizable to other COVID-19 patients. Prospective studies should be performed to confirm the additive value of CFD parameters in respiratory function assessment for critically ill patients with mechanical ventilation. Second, CFD analysis based on individual CT images are limited by the resolution of CT scans, which only allows a reconstruction of 3D airway models up to the region of airway diameters around 1-2 mm. However, evidence has demonstrated that the main part of the airway resistance is caused by the region between the 4 th and 7 th bifurcation. In addition, CT scans is the most commonly used imaging modality for COVID-19 patients to date. Finally, the role of drug therapy was not included, which can be important in a study to investigate the improvement of CFD parameters caused by drug therapy. Combined CFD and drug therapy for COVID-19 patients or acute respiratory distress syndrome (ARDS) patients would be of huge interest. The correlation between ∆V and ∆R. CT scans for Patient #2 indicate typical changes caused by COVID-19 virus infection in the lung lesion, including typical ground glass opacity (GGO) in the upper and lower lung lobes at the early stage of the disease (A&C) and pulmonary fibrosis and its induced tractive pulmonary bullae at the middle and late stages (B&D). CT scans for Patient #3 indicate that the exudation from the lung was obviously absorbed after effective treatment. Clinical characteristics of coronavirus disease 2019 in China A novel coronavirus from patients with pneumonia in China Presumed asymptomatic carrier transmission of COVID-19 Importation and human-to-human transmission of a novel coronavirus in Vietnam The first two cases of 2019-nCoV in Italy: Where they come from Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. The lancet WHO Coronavirus Disease (COVID-19) Dashboard Information about the new coronavirus disease (COVID-19) Mortality rates of patients with COVID-19 in the intensive care unit: a systematic review of the emerging literature Di Napoli R. Features, evaluation and treatment coronavirus (COVID-19) The COVID-19 pandemic: a comprehensive review of taxonomy, genetics, epidemiology, diagnosis, treatment, and control A review of the implications of computational fluid dynamic studies on nasal airflow and physiology Validation of computational fluid dynamics methodology used for human upper airway flow simulations Anatomically based three-dimensional model of airways to simulate flow and particle transport using computational fluid dynamics Computational fluid dynamics simulations of particle deposition in large-scale, multigenerational lung models The acute effect of budesonide/formoterol in COPD: a multi-slice computed tomography and lung function study Computational modeling of the obstructive lung diseases asthma and COPD Large-scale CFD simulations of airflow and particle deposition in lung airway Computational fluid dynamics simulation of airflow and aerosol deposition in human lungs Magnetic drug targeting through a realistic model of human tracheobronchial airways using computational fluid and particle dynamics Computational fluid dynamics can detect changes in airway resistance in asthmatics after acute bronchodilation A case series on lung deposition analysis of inhaled medication using functional imaging based computational fluid dynamics in asthmatic patients: effect of upper airway morphology and comparison with in vivo data Chinese Clinical Guidance for COVID-19 Pneumonia Diagnosis and Treatment Pulmonary artery sling with tracheal stenosis Functional imaging using computer methods to compare the effect of salbutamol and ipratropium bromide in patient-specific airway models of COPD Age-dependent effects in the transmission and control of COVID-19 epidemics Features of anosmia in COVID-19. Médecine et Maladies infectieuses Pathogenesis of COVID-19 from a cell biology perspective Alveolar macrophage dysfunction and cytokine storm in the pathogenesis of two severe COVID-19 patients A pathological report of three COVID-19 cases by minimally invasive autopsies. Zhonghua bing li xue za zhi= Chinese journal of pathology Effects of the Lower Airway Secretions on Airway Opening Pressures and Suction Pressures in Critically Ill COVID-19 Patients: A Computational Simulation Anti-inflammatory therapy may ameliorate the clinical picture of COVID-19. Annals of the Rheumatic Diseases This retrospective study demonstrates the application of patient-derived CFD analysis in non-invasive assessment of the respiratory function for clinically ill COVID-19 patients with mechanical ventilation. It is found that the change of airway resistance not only correlated well with the volume changes of the bronchial, but also showed significantly increasing in the lower lobe for the death group. The 'virtual' pulmonary function test by integrating follow-up CT scans with patient-derived CFD analysis provides the quantitative evaluation of respiratory functional improvement, suggesting that CFD analysis may have a role in improving the efficiency of treatment and predicting the prognosis and outcome for clinically ill patients with COVID-19.