Abstract
We have developed a computational model of gas mixing and ventilation in the human lung represented as a bifurcating network. We have simulated multiple-breath washout (MBW), a clinical test for measuring ventilation heterogeneity (VH) in patients with obstructive lung conditions. By applying airway constrictions inter-regionally, we have predicted the response of MBW indices to obstructions and found that they detect a narrow range of severe constrictions that reduce airway radius to 10%-30% of healthy values. These results help to explain the success of the MBW test to distinguish obstructive lung conditions from healthy controls. Further, we have used a perturbative approach to account for intra-regional airway heterogeneity that avoids modelling each airway individually. We have found, for random airway heterogeneity, that the variance in MBW indices is greater when indices are already elevated due to constrictions. By quantifying this effect, we have shown that variability in lung structure and mechanical properties alone can lead to clinically significant variability in MBW indices (specifically the Lung Clearance Index-LCI, and the gradient of phase-III slopes-Scond), but only in cases simulating obstructive lung conditions. This method is a computationally efficient way to probe the lung's sensitivity to structural changes, and to quantify uncertainty in predictions due to random variations in lung mechanical and structural properties.