In silico prototype of a human lung with a single airway to predict particle deposition

A revision of the multiple-path particle dosimetry model focusing on tobacco product aerosol dynamics

To assess the health impact of inhaled aerosols, it is necessary to understand aerosol dynamics and the associated dosimetry in the human respiratory tract. Although several studies have measured or simulated the dosimetry of aerosol constituents, the respiratory tract focus areas have been limited. In particular, the aerosols generated from tobacco products are complex composites and simulating their dynamics in the respiratory tract is challenging. To assess the dosimetry of the aerosol constituents of tobacco products, we developed a revised version of the Multiple-Path Particle Dosimetry (MPPD) model, which employs (1) new geometry based on CT-scanned human respiratory tract data, (2) convective mixing in the oral cavity and deep lung, and (3) constituent partitioning between the tissue and air, and clearance. The sensitivity analysis was conducted using aerosols composed of four major constituents of electronic cigarette (EC) aerosols to investigate the parameters that have a significant impact on the results. In addition, the revised model was run with 4 and 10 constituents in ECs and conventional cigarettes (CCs), respectively. Sensitivity analysis revealed that the new modeling and the physicochemical properties of constituents had a considerable impact on the simulated aerosol concentration and dosimetry. The simulations could be carried out within 3 min even when 10 constituents of CC aerosols were analyzed simultaneously. The revised model based on MPPD is an efficient and easy-to-use tool for understanding the aerosol dynamics of CC and EC constituents and their effect on the human body.

A critical analysis of the CFD-DEM simulation of pharmaceutical aerosols deposition in upper intra-thoracic airways: Considerations on air flow

A well-corroborated numerical methodology ensuring reproducibility in the modeling of pharmaceutical aerosols deposition in the respiratory system via CFD-DEM simulations within the RANS framework is currently missing. Often, inadequately clarified assumptions and approximations and the lack of evidences on their quantitative impact on the simulated deposition phenomenology, make a direct comparison among the different theoretical studies and the limited number of experiments a very challenging task. Here, with the ultimate goal of providing a critical analysis of some crucial computational aspects of aerosols deposition, we address the issues of velocity fluctuations propagation in the upper intra-thoracic airways and of the persistence of secondary flows using the SimInhale reference benchmark. We complement the investigation by describing how methodologies used to drive the flow through a truncated lung model may affect numerical results and how small discrepancies are observed in velocity profiles when comparing simulations based on different meshing strategies.

Predicting particle deposition using a simplified 8-path in silico human lung prototype

Understanding particle deposition in the human lung is crucial for the assessment of environmental pollutants and the design of new drug delivery systems. Traditionally, research has been carried out by experimental analysis, but this generally requires expensive equipment and exposure of volunteers to radiation, resulting in limited data. To overcome these drawbacks, there is an emphasis on the development of numerical models capable of accurate predictive analysis. The most advanced of these computer simulations are based on three-dimensional computational fluid dynamics. Solving the flow equations in a complete, fully resolved lung airway model is currently not feasible due to the computational resources required. In the present work, a simplified lung model is presented and validated for accurate prediction of particle deposition. Simulations are performed for an 8-path approximation to a full lung airway model. A novel boundary condition method is used to ensure accurate results in truncated flow branches. Simulations are performed at a steady inhalation flow rate of 18 l min-1, corresponding to a low activity breathing rate, while the effects of particle size and density are investigated. Comparison of the simulation results with available experimental data shows that reasonably accurate results can be obtained at a small fraction of the cost of a full airway model. The simulations clearly evaluate the effect of both particle size and particle density. Most importantly, the results show an improvement over a previously documented single-path model, both in terms of accuracy and the ability to obtain regional deposition rates. The present model represents an improvement over previously used simplified models, including single-path models. The multi-path reduced airway approach described can be used by researchers for general and patient-specific analyses of particle deposition and for the design of effective drug delivery systems.

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.

Construction of a hybrid lung model by combining a real geometry of the upper airways and an idealized geometry of the lower airways

BACKGROUND AND OBJECTIVE

Health care costs represent a substantial an increasing percentage of global expenditures. One key component is treatment of respiratory diseases, which account for one in twelve deaths in Europe. Computational simulations of lung airflow have potential to provide considerable cost reduction and improved outcomes. Such simulations require accurate in silico modeling of the lung airway. The geometry of the lung is extremely complex and for this reason very simple morphologies have primarily been used to date. The objective of this work is to develop an effective methodology for the creation of hybrid pulmonary geometries combining patient-specific models obtained from CT images and idealized pulmonary models, for the purpose of carrying out experimental and numerical studies on aerosol/particle transport and deposition in inhaled drug delivery.

METHODS

For the construction of the hybrid numerical model, lung images obtained from computed tomography were exported to the DICOM format to be treated with a commercial software to build the patient-specific part of the model. At the distal terminus of each airway of this portion of the model, an idealization of a single airway path is connected, extending to the sixteenth generation. Because these two parts have different endings, it is necessary to create an intermediate solid to link them together. Physically realistic treatment of truncated airway boundaries in the model was accomplished by mapping of the flow velocity distribution from corresponding conducting airway segments.

RESULTS

The model was verified using two sets of simulations, steady inspiration/expiration and transient simulation of forced spirometry. The results showed that the hybrid model is capable of providing a realistic description of air flow dynamics in the lung while substantially reducing computational costs relative to models of the full airway tree.

CONCLUSIONS

The model development outlined here represents an important step toward computational simulation of lung dynamics for patient-specific applications. Further research work may consist of investigating specific diseases, such as chronic bronchitis and pulmonary emphysema, as well as the study of the deposition of pollutants or drugs in the airways.