No abstract

Regional deposition effects are important in the pulmonary delivery of drugs intended for the topical treatment of respiratory ailments. They also play a critical role in the systemic delivery of drugs with limited lung bioavailability. In recent years, significant improvements in the quality of pulmonary imaging have taken place, however the resolution of current imaging modalities remains inadequate for quantifying regional deposition. Computational Fluid-Particle Dynamics (CFPD) can fill this gap by providing detailed information about regional deposition in the extrathoracic and conducting airways. It is therefore not surprising that the last 15years have seen an exponential growth in the application of CFPD methods in this area. Survey of the recent literature however, reveals a wide variability in the range of modelling approaches used and in the assumptions made about important physical processes taking place during aerosol inhalation. The purpose of this work is to provide a concise critical review of the computational approaches used to date, and to present a benchmark case for validation of future studies in the upper airways. In the spirit of providing the wider community with a reference for quality assurance of CFPD studies, in vitro deposition measurements have been conducted in a human-based model of the upper airways, and several groups within MP1404 SimInhale have computed the same case using a variety of simulation and discretization approaches. Here, we report the results of this collaborative effort and provide a critical discussion of the performance of the various simulation methods. The benchmark case, in vitro deposition data and in silico results will be published online and made available to the wider community. Particle image velocimetry measurements of the flow, as well as additional numerical results from the community, will be appended to the online database as they become available in the future.

No abstract

Atherosclerosis is becoming the number one cause of death worldwide. In this study, three-dimensional computer model of plaque formation and development for human carotid artery is developed. The three-dimensional blood flow is described by the Navier-Stokes equation, together with the continuity equation. Mass transfer within the blood lumen and through the arterial wall is coupled with the blood flow and is modeled by a convection-diffusion equation. The low-density lipoproteins transports in lumen of the vessel and through the vessel tissue are coupled by Kedem-Katchalsky equations. The inflammatory process is modeled using three additional reaction-diffusion partial differential equations. Fluid-structure interaction is used to estimate effective wall stress analysis. Plaque growth functions for volume progression are correlated with shear stress and effective wall stress distribution. We choose two specific patients from MRI study with significant plaque progression. Plaque volume progression using three time points for baseline, 3- and 12-month follow up is fitted. Our results for plaque localization correspond to low shear stress zone and we fitted parameters from our model using nonlinear least-square method. Determination of plaque location and composition, and computer simulation of progression in time for a specific patient shows a potential benefit for the prediction of disease progression. The proof of validity of three-dimensional computer modeling in the evaluation of atherosclerotic plaque burden may shift the clinical information of MRI from morphological assessment toward a functional tool. Understanding and prediction of the evolution of atherosclerotic plaques either into vulnerable or stable plaques are major tasks for the medical community.

Vascular endothelial cells are continuously exposed to hemodynamic shear stress. Intensity and type of shear stress are highly relevant to vascular physiology and pathology. Here, we modeled shear stress distribution in a tissue culture well (R ¼ 17.5 mm, fill volume 2 ml) under orbital translation using computational fluid dynamics with the finite element method. Free surface distribution, wall shear stress, inclination angle, drag force, and oscillatory index on the bottom surface were modeled. Obtained results predict nonuniform shear stress distribution during cycle, with higher oscillatory shear index, higher drag force values, higher circular component, and larger inclination angle of the shear stress at the periphery of the well compared with the center of the well. The oscillatory index, inclination angle, and drag force are new quantitative parameters modeled in this system, which provide a better understanding of the hydrodynamic conditions experienced and reflect the pulsatile character of blood flow in vivo. Validation experiments revealed that endothelial cells at the well periphery aligned under flow and increased Kruppel-like Factor 4 (KLF-4), cyclooxygenase-2 (COX-2) expression and endothelial nitric oxide synthase (eNOS) phosphorylation. In contrast, endothelial cells at the center of the well did not show clear directional alignment, did not induce the expression of KLF-4 and COX-2 nor increased eNOS phosphorylation. In conclusion, this improved computational modeling predicts that the orbital shaker model generates different hydrodynamic conditions at the periphery versus the center of the well eliciting divergent endothelial cell responses. The possibility of generating different hydrodynamic conditions in the same well makes this model highly attractive to study responses of distinct regions of the same endothelial monolayer to different types of shear stresses thereby better reflecting in vivo conditions.