In this article, we develop a lung ventilation model. The parenchyma is described as an elastic homogenized media. It is irrigated by a space-filling dyadic resistive pipe network, which represents the tracheobronchial tree. In this model, the tree and the parenchyma are strongly coupled. The tree induces an extra viscous term in the system constitutive relation, which leads, in the finite element framework, to a full matrix. We consider an efficient algorithm that takes advantage of the tree structure to enable a fast matrix-vector product computation. This framework can be used to model both free and mechanically induced respiration, in health and disease. Patient-specific lung geometries acquired from computed tomography scans are considered. Realistic Dirichlet boundary conditions can be deduced from surface registration on computed tomography images. The model is compared to a more classical exit compartment approach. Results illustrate the coupling between the tree and the parenchyma, at global and regional levels, and how conditions for the purely 0D model can be inferred. Different types of boundary conditions are tested, including a nonlinear Robin model of the surrounding lung structures.

In spite of numerous clinical studies, there is no consensus on the benefit Heliox mixtures can bring to asthmatic patients in terms of work of breathing and ventilation distribution. In this article we use a 3D finite element mathematical model of the lung to study the impact of asthma on effort and ventilation distribution along with the effect of Heliox compared to air. Lung surface displacement fields extracted from computed tomography medical images are used to prescribe realistic boundary conditions to the model. Asthma is simulated by imposing bronchoconstrictions to some airways of the tracheo-bronchial tree based on statistical laws deduced from the literature. This study illuminates potential mechanisms for patient responsiveness to Heliox when affected by obstructive pulmonary diseases. Responsiveness appears to be function of the pathology severity, as well as its distal position in the tracheo-bronchial tree and geometrical position within the lung.