Physical and CFD model used in the analysis of particles dispersion.

Air pollution, representing one of the major environmental crises affecting us all, is responsible for half a million deaths each year in Europe. This research shows a numerical model based on Computational Fluid Dynamics (CFD) and a methodology for its validation allowing to know the evolution of particles in open urban environments. In this case, the model represents an area of Gijón (Asturias), specifically chosen as it serves to highlight a location at which the daily limit values of PM10 set by the regulations are most frequently exceeded. A numerical model has been developed at a scale of 1:10,000 to simulate the dispersion of pollutants, including particulate emissions. A physical model at the same scale was built using 3D printing. This model was tested in a wind tunnel and analysed in the four main wind directions. To carry out the tests, a particle generation system was designed and manufactured, and the necessary infrastructure was built to take measurements using an optical particle meter. Results show that the numerical model meets the expected objectives and is capable of predicting the behaviour of particle dispersion in the air. The numerical model produces results in the same order of magnitude as the physical model, although there is a tendency to underestimate the maximum values with respect to the measured ones.

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.

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.

Implementation of a specific boundary condition for a simplified symmetric single-path CFD lung model with OpenFOAM.

CFD modeling research about the lung airflow with a complete resolution and an adequate accuracy at all scales requires a great amount of computational resources due to the vast number of necessary grid elements. As a result, a common practice is to conduct simplifications that allows to manage it with ordinary computational power. In this study, the implementation of a special boundary condition in order to develop a simplified single conductive lung airway model, which exactly represents the effect of the removed airways, is presented. The boundary condition is programmed in the open-source software OpenFOAM®, and the developed source code is presented in the proper syntax. After this description, modeling accuracy is evaluated under different flow rate conditions typical of human breathing processes, including both inspiration and expiration movements. Afterward, a validation process is conducted using results of a Weibel's model (0-4 generations) simulation for a medium flow rate of 50 L/min. Finally, a comparison against the proposed boundary condition implemented in the commercial code ANSYS Fluent is made, which highlights the benefits of using the free code toolbox. The specific contribution of this paper will be to show that OpenFOAM® developed model can perform even better than other commercial codes due to a precise implementation and coupling of the default solver with the in-house functions by virtue of the open-source nature of the code.

Simulation of the human airways using virtual topology tools and meshing optimization.

A method is proposed to improve the quality of the three-dimensional airway geometric models using a commercial software, checking the number of elements, meshing time, and aspect ratio and skewness parameters. The use of real and virtual topologies combined with patch-conforming and patch-independent meshing algorithms results in four different models being the best solution the combination of virtual topology and patch-independent algorithm, due to an excellent aspect ratio and skewness of the elements, and minimum meshing time. The result is a reduction in the computational time required for both meshing and simulation due to a smaller number of cells. The use of virtual topologies combined with patch-independent meshing algorithms could be extended in bioengineering because the geometries handling is similar to this case. The method is applied to a healthy person using their computed tomography images. The resulting numerical models are able to simulate correctly a forced spirometry.