Integrative Computational Models of Lung Structure-Function Interactions

Algebraic approximation of the distributed model for the pressure drop in the respiratory airways

The complexity of the human respiratory system causes that one of the main methods of analyzing the dynamic pulmonary phenomena and interpreting experimental results are simulations of its computational models. Among the most compound elements of these models, apart from the bronchial tree structure, is the phenomenon of flow limitation in flexible bronchi, which causes them to collapse with increasing flow, thus their properties, such as resistance, compliance and inertance, are highly nonlinear and time-varying. Commonly, this phenomenon is ignored, or a distributed model for the airway pressure drop is applied, simulated with a modified numerical solver of this differential equation (ODE). The disadvantages of this solution are the problems with taking into account the inherent singularity of the model and the long computation time due to iterative nature of the ODE procedure. The aim of the work was to derive an algebraic approximation of this distributed model, suitable for implementation in continuous dynamic models, to validate it by comparing the results of simulations with the respiratory system model including approximate and original (ODE solver) numerical procedures, as well as to evaluate possible acceleration of calculations. All simulations, including spontaneous breathing, mechanical ventilation with the optimal ventilatory waveform and forced expiration, proved that algebraic approximation yielded results negligibly differing from the ODE solution, and shortened the computation time by an order. The proposed approach is an attractive alternative in the case of computer implementations of pulmonary models, where simulations of flows and pressures in the complex respiratory system are of primary importance.