Three-dimensional convective alveolar flow induced by rhythmic breathing motion of the pulmonary acinus

Low Reynolds number flows (Re<1) in the human pulmonary acinus are often difficult to assess due to the submillimeter dimensions and accessibility of the region. In the present computational study, we simulated three-dimensional alveolar flows in an alveolated duct at each generation of the pulmonary acinar tree using recent morphometric data. Rhythmic lung expansion and contraction motion was modeled using moving wall boundary conditions to simulate realistic sedentary tidal breathing. The resulting alveolar flow patterns are largely time independent and governed by the ratio of the alveolar to ductal flow rates, Qa/Qd. This ratio depends uniquely on geometrical configuration such that alveolar flow patterns may be entirely determined by the location of the alveoli along the acinar tree. Although flows within alveoli travel very slowly relative to those in acinar ducts, 0.021%相似文献   

Flow field analysis in a compliant acinus replica model using particle image velocimetry (PIV)

Inhaled particles reaching the alveolar walls have the potential to cross the blood–gas barrier and enter the blood stream. Experimental evidence of pulmonary dosimetry, however, cannot be explained by current whole lung dosimetry models. Numerical and experimental studies shed some light on the mechanisms of particle transport, but realistic geometries have not been investigated. In this study, a three dimensional expanding model including two generations of respiratory bronchioles and five terminal alveolar sacs was created from a replica human lung cast. Flow visualization techniques were employed to quantify the fluid flow while utilizing streamlines to evaluate recirculation. Pathlines were plotted to track the fluid motion and estimate penetration depth of inhaled air. This study provides evidence that the two generations immediately proximal to the terminal alveolar sacs do not have recirculating eddies, even for intense breathing. Results of Peclet number calculations indicate that substantial convective motion is present in vivo for the case of deep breathing, which significantly increases particle penetration into the alveoli. However, particle diffusion remains the dominant mechanism of particle transport over convection, even for intense breathing because inhaled particles do not reach the alveolar wall in a single breath by convection alone. Examination of the velocity fields revealed significant uneven ventilation of the alveoli during a single breath, likely due to variations in size and location. This flow field data, obtained from replica model geometry with realistic breathing conditions, provides information to better understand fluid and particle behavior in the acinus region of the lung.     

Gravitational deposition in a rhythmically expanding and contracting alveolus.

In a previous simulation, our laboratory demonstrated that the flow induced by a rhythmically expanding and contracting alveolus is highly complex (Haber S, Butler JP, Brenner H, Emanuel I, and Tsuda A, J Fluid Mech 405: 243-268, 2000). Based on these earlier findings, we hypothesize that the trajectories and deposition of aerosols inside the alveoli differ substantially from those previously predicted. To test this hypothesis, trajectories of fine particles (0.5-2.5 microm in diameter) moving in the foregoing alveolar flow field and simultaneously subjected to the gravity field were simulated. The results show that alveolar wall motion is crucial in determining the enhancement of aerosol deposition inside the alveoli. In particular, 0.5- to 1-microm-diameter particles are sensitive to the detailed alveolar flow structure (e.g., recirculating flow), as they undergo gravity-induced convective mixing and deposition. Accordingly, deposition concentrations within each alveolus are nonuniform, with preferentially higher densities near the alveolar entrance ring, consistent with physiological observations. Deposition patterns along the acinar tree are also nonuniform, with higher deposition in the first half of the acinar generations. This is a result of the combined effects of enhanced alveolar deposition in the proximal region of the acinus due to alveoli expansion and contraction and reduction in the number of particles remaining in the gas phase down the acinar tree. We conclude that the cyclically expanding and contracting motion of alveoli plays an important role in determining gravitational deposition in the pulmonary acinus.     

Validation of CFD predictions of flow in a 3D alveolated bend with experimental data

Verifying numerical predictions with experimental data is an important aspect of any modeling studies. In the case of the lung, the absence of direct in vivo flow measurements makes such verification almost impossible. We performed computational fluid dynamics (CFD) simulations in a 3D scaled-up model of an alveolated bend with rigid walls that incorporated essential geometrical characteristics of human alveolar structures and compared numerical predictions with experimental flow measurements made in the same model by particle image velocimetry (PIV). Flow in both models was representative of acinar flow during normal breathing (0.82ml/s). The experimental model was built in silicone and silicone oil was used as the carrier fluid. Flow measurements were obtained by an ensemble averaging procedure. CFD simulation was performed with STAR-CCM+ (CD-Adapco) using a polyhedral unstructured mesh. Velocity profiles in the central duct were parabolic and no bulk convection existed between the central duct and the alveoli. Velocities inside the alveoli were approximately 2 orders of magnitude smaller than the mean velocity in the central duct. CFD data agreed well with those obtained by PIV. In the central duct, data agreed within 1%. The maximum simulated velocity along the centerline of the model was 0.5% larger than measured experimentally. In the alveolar cavities, data agreed within 15% on average. This suggests that CFD techniques can satisfactorily predict acinar-type flow. Such a validation ensure a great degree of confidence in the accuracy of predictions made in more complex models of the alveolar region of the lung using similar CFD techniques.     

Assessment of human left ventricle flow using statistical shape modelling and computational fluid dynamics

Blood flow patterns in the human left ventricle (LV) have shown relation to cardiac health. However, most studies in the literature are limited to a few patients and results are hard to generalize. This study aims to provide a new framework to generate more generalized insights into LV blood flow as a function of changes in anatomy and wall motion. In this framework, we studied the four-dimensional blood flow in LV via computational fluid dynamics (CFD) in conjunction with a statistical shape model (SSM), built from segmented LV shapes of 150 subjects. We validated results in an in-vitro dynamic phantom via time-resolved optical particle image velocimetry (PIV) measurements. This combination of CFD and the SSM may be useful for systematically assessing blood flow patterns in the LV as a function of varying anatomy and has the potential to provide valuable data for diagnosis of LV functionality.     

Validating CFD predictions of respiratory aerosol deposition: effects of upstream transition and turbulence

A number of computational fluid dynamics (CFD) studies have explored local deposition patterns of inhaled aerosols in the respiratory tract. These studies have highlighted the effects of multiple physiologic, geometric, and particle characteristics on deposition. However, very few studies have reported local or sub-branch quantitative comparisons to in vitro particle deposition data. The objective of this study is to numerically investigate the effects of transition and turbulence on highly localized particle deposition in a respiratory double bifurcation model in order to quantitatively validate CFD results. To perform the validations, local comparisons have been made to a specific in vitro case study of 10 microm particles depositing in a model of respiratory generations G3-G5. To achieve this objective, two geometric cases have been considered. The first case includes only the double bifurcation model. The second case includes a portion of the experimental particle delivery geometry, where transitional flow is expected. To evaluate the effectiveness of two-equation turbulence models in this system, the flow field solutions have been computed using laminar, standard k-omega, and low Reynolds number (LRN) k-omega approximations. Results indicate that even though the Reynolds number remained below the critical limit required for full turbulence, transition and turbulence have a significant impact on the flow field and local particle deposition patterns. For the experimental case considered, turbulence impacted the local deposition of 10 microm particles primarily by influencing the initial velocity and particle profiles. As such, both the laminar and LRN k-omega flow models provided good local quantitative matches to the in vitro deposition data, provided that the correct initial particle profile was specified. Implications of this study include the need for local quantitative validations of particle deposition results, the importance of correct inlet conditions, and the need to consider upstream effects in experimental and computational studies of the respiratory tract. Applications of these results to realistic respiratory geometries will require consideration on upstream flow conditions in the lung, transient flow, and intermittent turbulent structures.     

Simulation and visualization of flow pattern in microarrays for liquid phase oligonucleotide and peptide synthesis

Microfluidic microarrays have been developed for economical and rapid parallel synthesis of oligonucleotide and peptide libraries. For a synthesis system to be reproducible and uniform, it is crucial to have a uniform reagent delivery throughout the system. Computational fluid dynamics (CFD) is used to model and simulate the microfluidic microarrays to study geometrical effects on flow patterns. By proper design geometry, flow uniformity could be obtained in every microreactor in the microarrays.     

Radial transport along the human acinar tree

A numerical model of an expanding asymmetric alveolated duct was developed and used to investigate lateral transport between the central acinar channel and the surrounding alveoli along the acinar tree. Our results indicate that some degree of recirculation occurs in all but the terminal generations. We found that the rate of diffusional transport of axial momentum from the duct to the alveolus was by far the largest contributor to the resulting momentum in the alveolar flow but that the magnitude of the axial momentum is critical in determining the nature of the flow in the alveolus. Further, we found that alveolar flow rotation, and by implication chaotic mixing, is strongest in the entrance generations. We also found that the expanding alveolus provides a pathway by which particles with little intrinsic motion can enter the alveoli. Thus, our results offer a possible explanation for why submicron particles deposit preferentially in the acinar-entrance region.     

Streamline crossing: An essential mechanism for aerosol dispersion in the pulmonary acinus

The dispersion of inhaled microparticles in the pulmonary acinus of the lungs is often attributed to the complex interplay between convective mixing, due to irreversible flows, and intrinsic particle motion (i.e. gravity and diffusion). However, the role of each mechanism, the exact nature of such interplay between them and their relative importance still remain unclear. To gain insight into these dispersive mechanisms, we track liquid-suspended microparticles and extract their effective diffusivities inside an anatomically-inspired microfluidic acinar model. Such results are then compared to experiments and numerical simulations in a straight channel. While alveoli of the proximal acinar generations exhibit convective mixing characteristics that lead to irreversible particle trajectories, this local effect is overshadowed by a more dominant dispersion mechanism across the ductal branching network that arises from small but significant streamline crossing due to intrinsic diffusional motion in the presence of high velocity gradients. We anticipate that for true airborne particles, which exhibit much higher intrinsic motion, streamline crossing would be even more significant.     

Computational Fluid Dynamics Simulations at Micro-Scale Stenosis for Microfluidic Thrombosis Model Characterization

Platelet aggregation plays a central role in pathological thrombosis, preventing healthy physiological blood flow within the circulatory system. For decades, it was believed that platelet aggregation was primarily driven by soluble agonists such as thrombin, adenosine diphosphate and thromboxane A2. However, recent experimental findings have unveiled an intriguing but complementary biomechanical mechanism—the shear rate gradients generated from flow disturbance occurring at sites of blood vessel narrowing, otherwise known as stenosis, may rapidly trigger platelet recruitment and subsequent aggregation. In our Nature Materials 2019 paper [1], we employed microfluidic devices which incorporated micro-scale stenoses to elucidate the molecular insights underlying the prothrombotic effect of blood flow disturbance. Nevertheless, the rheological mechanisms associated with this stenotic microfluidic device are poorly characterized. To this end, we developed a computational fluid dynamics (CFD) simulation approach to systematically analyze the hemodynamic influence of bulk flow mechanics and flow medium. Grid sensitivity studies were performed to ensure accurate and reliable results. Interestingly, the peak shear rate was significantly reduced with the device thickness, suggesting that fabrication of microfluidic devices should retain thicknesses greater than 50 µm to avoid unexpected hemodynamic aberration, despite thicker devices raising the cost of materials and processing time of photolithography. Overall, as many groups in the field have designed microfluidic devices to recapitulate the effect of shear rate gradients and investigate platelet aggregation, our numerical simulation study serves as a guideline for rigorous design and fabrication of microfluidic thrombosis models.     

Numerical simulation of airflow and microparticle deposition in a synchrotron micro-CT-based pulmonary acinus model

The acinus consists of complex, branched alveolar ducts and numerous surrounding alveoli, and so in this study, we hypothesized that the particle deposition can be much influenced by the complex acinar geometry, and simulated the airflow and particle deposition (density = 1.0 g/cm³, diameter = 1 and 3 μm) numerically in a pulmonary acinar model based on synchrotron micro-CT of the mammalian lung. We assumed that the fluid–structure interaction was neglected and that alveolar flow was induced by the expansion and contraction of the acinar model with the volume changing sinusoidally with time as the moving boundary conditions. The alveolar flow was dominated by radial flows, and a weak recirculating flow was observed at the proximal side of alveoli during the entire respiratory cycle, despite the maximum Reynolds number at the inlet being 0.029. Under zero gravity, the particle deposition rate after single breathing was less than 0.01, although the particles were transported deeply into the acinus after inspiration. Under a gravitational field, the deposition rate and map were influenced strongly by gravity orientation. In the case of a particle diameter of 1 μm, the rate increased dramatically and mostly non-deposited particles remained in the model, indicating that the rate would increase further after repeated breathing. At a particle diameter of 3 μm, the rate was 1.0 and all particles were deposited during single breathing. Our results show that the particle deposition rate in realistic pulmonary acinar model is higher than in an idealized model.     

MRI and CFD studies of pulsatile flow in healthy and stenosed carotid bifurcation models

Pulsatile flow was studied in physiologically realistic models of a normal and a moderately stenosed (30% diameter reduction) human carotid bifurcation. Time-resolved velocity measurements were made using magnetic resonance imaging, from which wall shear stress (WSS) vectors were calculated. Velocity measurements in the inflow and outflow regions were also used as boundary conditions for a computational fluid dynamics (CFD) model. Experimental flow patterns and derived WSS vectors were compared qualitatively with the corresponding CFD predictions. In the stenosed phantom, flow in the bulb region of the "internal carotid artery" was concentrated along the outer wall, with a region of low and recirculating flow near the inner wall. In the normal phantom, the converse was found, with a low flow region near the outer wall of the bulb. Time-averaged WSS and oscillatory shear index were also markedly different for the two phantoms.     

Inertial particle dynamics in large artery flows – Implications for modeling arterial embolisms

The complexity of inertial particle dynamics through swirling chaotic flow structures characteristic of pulsatile large-artery hemodynamics renders significant challenges in predictive understanding of transport of such particles. This is specifically crucial for arterial embolisms, where knowledge of embolus transport to major vascular beds helps in disease diagnosis and surgical planning. Using a computational framework built upon image-based CFD and discrete particle dynamics modeling, a multi-parameter sampling-based study was conducted on embolic particle dynamics and transport. The results highlighted the strong influence of material properties, embolus size, release instance, and embolus source on embolus distribution to the cerebral, renal and mesenteric, and ilio-femoral vasculature beds. The study also isolated the importance of shear-gradient lift, and elastohydrodynamic contact, in affecting embolic particle transport. Near-wall particle re-suspension due to lift alters aortogenic embolic particle dynamics significantly as compared to cardiogenic. The observations collectively indicated the complex interplay of particle inertia, fluid–particle density ratio, and wall collisions, with chaotic flow structures, which render the overall motion of the particles to be non-trivially dispersive in nature.     

Effects of tumors on inhaled pharmacologic drugs

Lung carcinomas are now the most common form of cancer. Clinical data suggest that tumors are found preferentially in upper airways, perhaps specifically at carina within bifurcations. The disease can be treated by aerosolized pharmacologic drugs. To enhance their efficacies site-specific drugs must be deposited selectively. Since inhaled particles are transported by air, flow patterns will naturally affect their trajectories. Therefore, in Part I of a systematic investigation, we focused on tumor-induced effects on airstreams, in Part II (the following article [p. 245]), particle trajectories were determined. To facilitate the targeted delivery of inhaled drugs, we simulated bifurcations with tumors on carinas using a commercial computational fluid dynamics (CFD) software package (FIDAP) with a Cray T90 supercomputer and studied effects of tumor sizes and ventilatory parameters on localized flow patterns. Critical tumor sizes existed; e.g., tumors had dominant effects when r/R > or = 0.8 for bifurcation 3-4 and r/R > or = 0.6 for bifurcation 7-8 (r = tumor radius and R = airway radius). The findings suggest that computer modeling is a means to integrate alterations to airway structures caused by diseases into aerosol therapy protocols.     

Dependence of leukocyte capture on instantaneous pulsatile flow

Atherosclerosis, an artery disease, is currently the leading cause of death in the United States in both men and women. The first step in the development of atherosclerosis involves leukocyte adhesion to the arterial endothelium. It is broadly accepted that blood flow, more specifically wall shear stress (WSS), plays an important role in leukocyte capture and subsequent development of an atherosclerotic plaque. What is less known is how instantaneous WSS, which can vary by up to 5 Pa over one cardiac cycle, influences leukocyte capture. In this paper we use direct numerical simulations (DNS), performed using an in-house code, to illustrate that leukocyte capture is different whether as a function of instantaneous or time-averaged blood flow. Specifically, a stenotic plaque is modeled using a computational fluid dynamics (CFD) solver through fully three-dimensional Navier-Stokes equations and the immersed boundary method. Pulsatile triphasic inflow is used to simulate the cardiac cycle. The CFD is coupled with an agent-based leukocyte capture model to assess the impact of instantaneous hemodynamics on stenosis growth. The computed wall shear stress agrees well with the results obtained with a commercial software, as well as with theoretical results in the healthy region of the artery. The analysis emphasizes the importance of the instantaneous flow conditions in evaluating the leukocyte rate of capture. That is, the capture rate computed from mean flow field is generally underpredicted compared to the actual rate of capture. Thus, in order to obtain a reliable estimate, the flow unsteadiness during a cardiac cycle should be taken into account.     

Validation of CFD simulations of cerebral aneurysms with implication of geometric variations

BACKGROUND: Computational fluid dynamics (CFD) simulations using medical-image-based anatomical vascular geometry are now gaining clinical relevance. This study aimed at validating the CFD methodology for studying cerebral aneurysms by using particle image velocimetry (PIV) measurements, with a focus on the effects of small geometric variations in aneurysm models on the flow dynamics obtained with CFD. METHOD OF APPROACH: An experimental phantom was fabricated out of silicone elastomer to best mimic a spherical aneurysm model. PIV measurements were obtained from the phantom and compared with the CFD results from an ideal spherical aneurysm model (S1). These measurements were also compared with CFD results, based on the geometry reconstructed from three-dimensional images of the experimental phantom. We further performed CFD analysis on two geometric variations, S2 and S3, of the phantom to investigate the effects of small geometric variations on the aneurysmal flow field. Results. We found poor agreement between the CFD results from the ideal spherical aneurysm model and the PIV measurements from the phantom, including inconsistent secondary flow patterns. The CFD results based on the actual phantom geometry, however, matched well with the PIV measurements. CFD of models S2 and S3 produced qualitatively similar flow fields to that of the phantom but quantitatively significant changes in key hemodynamic parameters such as vorticity, positive circulation, and wall shear stress. CONCLUSION: CFD simulation results can closely match experimental measurements as long as both are performed on the same model geometry. Small geometric variations on the aneurysm model can significantly alter the flow-field and key hemodynamic parameters. Since medical images are subjected to geometric uncertainties, image-based patient-specific CFD results must be carefully scrutinized before providing clinical feedback.     

Computational fluid dynamics simulations of particle deposition in large-scale, multigenerational lung models

Computational fluid dynamics (CFD) has emerged as a useful tool for the prediction of airflow and particle transport within the human lung airway. Several published studies have demonstrated the use of Eulerian finite-volume CFD simulations coupled with Lagrangian particle tracking methods to determine local and regional particle deposition rates in small subsections of the bronchopulmonary tree. However, the simulation of particle transport and deposition in large-scale models encompassing more than a few generations is less common, due in part to the sheer size and complexity of the human lung airway. Highly resolved, fully coupled flowfield solution and particle tracking in the entire lung, for example, is currently an intractable problem and will remain so for the foreseeable future. This paper adopts a previously reported methodology for simulating large-scale regions of the lung airway (Walters, D. K., and Luke, W. H., 2010, "A Method for Three-Dimensional Navier-Stokes Simulations of Large-Scale Regions of the Human Lung Airway," ASME J. Fluids Eng., 132(5), p. 051101), which was shown to produce results similar to fully resolved geometries using approximate, reduced geometry models. The methodology is extended here to particle transport and deposition simulations. Lagrangian particle tracking simulations are performed in combination with Eulerian simulations of the airflow in an idealized representation of the human lung airway tree. Results using the reduced models are compared with those using the fully resolved models for an eight-generation region of the conducting zone. The agreement between fully resolved and reduced geometry simulations indicates that the new method can provide an accurate alternative for large-scale CFD simulations while potentially reducing the computational cost of these simulations by several orders of magnitude.     

Dependence of adhesive behavior of neutrophils on local fluid dynamics in a region with recirculating flow

We have recently described patterns of adhesion of different types of leukocytes downstream of a backward facing step. Here the predicted fluid dynamics in channels incorporating backward facing steps are described, and related to the measured velocities of flowing cells, patterns of attachment and characteristics of rolling adhesion for neutrophils perfused over P-selectin. Deeper (upstream depth 300 microm, downstream depth 600 microm, maximum wall shear stress approximately 0.1 Pa) and shallower (upstream depth 260 microm, downstream depth 450 microm, maximum wall shear stress approximately 0.3 Pa) channels were compared. Computational fluid dynamics (CFD) predicted the presence of vortices downstream of the steps, distances to reattachment of flow, local wall shear stresses and components of velocity parallel and perpendicular to the wall. Measurements of velocities of perfused neutrophils agreed well with predictions, and suggested that adhesion to P-selectin should be possible in the regions of recirculating flow, but not downstream in re-established flow in the high shear channel. When channels were coated with a P-selectin-Fc chimaera, neutrophils were captured from flow and immobilised. Capture showed local maxima around the reattachment points, but was absent elsewhere in the high shear chamber. In the low shear chamber there was depression of adhesion just beyond the reattachment point because of expansion of flow and depletion of neutrophils near the wall. Inside the recirculation zones, adhesion decreased approaching the step because of an increasing, vertically upward velocity component. When channels were coated with P-selectin, neutrophils rolled in all regions, but lifted off the surface as they rolled backwards into low shear regions near the step. Rolling velocity in the recirculation zone was independent of shear stress, possibly because of the effects of vertical lift. We conclude that while local wall shear stress influences adhesive behavior, delivery of cells to the wall and their behavior after capture also depend on components of flow perpendicular to the wall.     

Aerosol deposition characteristics in distal acinar airways under cyclic breathing conditions

Although the major mechanisms of aerosol deposition in the lung are known, detailed quantitative data in anatomically realistic models are still lacking, especially in the acinar airways. In this study, an algorithm was developed to build multigenerational three-dimensional models of alveolated airways with arbitrary bifurcation angles and spherical alveolar shape. Using computational fluid dynamics, the deposition of 1- and 3-μm aerosol particles was predicted in models of human alveolar sac and terminal acinar bifurcation under rhythmic wall motion for two breathing conditions (functional residual capacity = 3 liter, tidal volume = 0.5 and 0.9 liter, breathing period = 4 s). Particles entering the model during one inspiration period were tracked for multiple breathing cycles until all particles deposited or escaped from the model. Flow recirculation inside alveoli occurred only during transition between inspiration and expiration and accounted for no more than 1% of the whole cycle. Weak flow irreversibility and convective transport were observed in both models. The average deposition efficiency was similar for both breathing conditions and for both models. Under normal gravity, total deposition was ~33 and 75%, of which ~67 and 96% occurred during the first cycle, for 1- and 3-μm particles, respectively. Under zero gravity, total deposition was ~2-5% for both particle sizes. These results support previous findings that gravitational sedimentation is the dominant deposition mechanism for micrometer-sized aerosols in acinar airways. The results also showed that moving walls and multiple breathing cycles are needed for accurate estimation of aerosol deposition in acinar airways.     

Fabrication,Operation and Flow Visualization in Surface-acoustic-wave-driven Acoustic-counterflow Microfluidics

Surface acoustic waves (SAWs) can be used to drive liquids in portable microfluidic chips via the acoustic counterflow phenomenon. In this video we present the fabrication protocol for a multilayered SAW acoustic counterflow device. The device is fabricated starting from a lithium niobate (LN) substrate onto which two interdigital transducers (IDTs) and appropriate markers are patterned. A polydimethylsiloxane (PDMS) channel cast on an SU8 master mold is finally bonded on the patterned substrate. Following the fabrication procedure, we show the techniques that allow the characterization and operation of the acoustic counterflow device in order to pump fluids through the PDMS channel grid. We finally present the procedure to visualize liquid flow in the channels. The protocol is used to show on-chip fluid pumping under different flow regimes such as laminar flow and more complicated dynamics characterized by vortices and particle accumulation domains.