Three-Dimensional Convective Alveolar Flow Induced by Rhythmic Breathing Motion of the Pulmonary Acinus

2007 ◽  
Vol 129 (5) ◽  
pp. 658-665 ◽  
Author(s):  
Josué Sznitman ◽  
Fabian Heimsch ◽  
Thomas Heimsch ◽  
Daniel Rusch ◽  
Thomas Rösgen

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, Q̇a∕Q̇d. 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%⩽Ua∕Ud⩽9.1%, they may exhibit complex patterns linked to the three-dimensional nature of the flow and confirm findings from earlier three-dimensional simulations. Such patterns are largely determined by the interplay between recirculation in the cavity induced by ductal shear flow over the alveolar opening and radial flows induced by wall displacement. Furthermore, alveolar flow patterns under rhythmic wall motion contrast sharply with results obtained in a rigid alveolus, further confirming the importance of including inherent wall motion to understand realistic acinar flow phenomena. The present findings may give further insight into the role of convective alveolar flows in determining aerosol kinematics and deposition in the pulmonary acinus.

2009 ◽  
Vol 131 (3) ◽  
Author(s):  
Josué Sznitman ◽  
Thomas Heimsch ◽  
Johannes H. Wildhaber ◽  
Akira Tsuda ◽  
Thomas Rösgen

The inhalation of micron-sized aerosols into the lung’s acinar region may be recognized as a possible health risk or a therapeutic tool. In an effort to develop a deeper understanding of the mechanisms responsible for acinar deposition, we have numerically simulated the transport of nondiffusing fine inhaled particles (1 μm and 3 μm in diameter) in two acinar models of varying complexity: (i) a simple alveolated duct and (ii) a space-filling asymmetrical acinar branching tree following the description of lung structure by Fung (1988, “A Model of the Lung Structure and Its Validation,” J. Appl. Physiol., 64, pp. 2132–2141). Detailed particle trajectories and deposition efficiencies, as well as acinar flow structures, were investigated under different orientations of gravity, for tidal breathing motion in an average human adult. Trajectories and deposition efficiencies inside the alveolated duct are strongly related to gravity orientation. While the motion of larger particles (3 μm) is relatively insensitive to convective flows compared with the role of gravitational sedimentation, finer 1 μm aerosols may exhibit, in contrast, complex kinematics influenced by the coupling between (i) flow reversal due to oscillatory breathing, (ii) local alveolar flow structure, and (iii) streamline crossing due to gravity. These combined mechanisms may lead to twisting and undulating trajectories in the alveolus over multiple breathing cycles. The extension of our study to a space-filling acinar tree was well suited to investigate the influence of bulk kinematic interaction on aerosol transport between ductal and alveolar flows. We found the existence of intricate trajectories of fine 1 μm aerosols spanning over the entire acinar airway network, which cannot be captured by simple alveolar models. In contrast, heavier 3 μm aerosols yield trajectories characteristic of gravitational sedimentation, analogous to those observed in the simple alveolated duct. For both particle sizes, however, particle inhalation yields highly nonuniform deposition. While larger particles deposit within a single inhalation phase, finer 1 μm particles exhibit much longer residence times spanning multiple breathing cycles. With the ongoing development of more realistic models of the pulmonary acinus, we aim to capture some of the complex mechanisms leading to deposition of inhaled aerosols. Such models may lead to a better understanding toward the optimization of pulmonary drug delivery to target specific regions of the lung.


2014 ◽  
Vol 136 (8) ◽  
Author(s):  
N. K. C. Selvarasu ◽  
Danesh K. Tafti

The incidence of stent late restenosis is high (Zwart et al., 2010, “Coronary Stent Thrombosis in the Current Era: Challenges and Opportunities for Treatment,” Curr. Treat. Options Cardiovasc. Med., 12(1), pp. 46–57) despite the extensive use of stents, and is most prevalent at the proximal and distal ends of the stent. Elastic modulus change in stented coronary arteries subject to the motion of the myocardium is not studied extensively. It is our objective to understand and reveal the mechanism by which changes in elastic modulus and geometry contribute to the generation of nonphysiological wall shear stress (WSS). Such adverse hemodynamic conditions could have an effect on the onset of restenosis. Three-dimensional (3D), spatiotemporally resolved computational fluid dynamics (CFD) simulations of pulsatile flow with moving wall boundaries and fluid structure interaction (FSI) were carried out for a helical artery with physiologically relevant flow parameters. To study the effect of coronary artery (CA) geometry change on stent elastic modulus mismatch, models where the curvature, torsion and both curvature and torsion change were examined. The elastic modulus is increased by a factor of two, five, and ten in the stented section for all three modes of motion. The changes in elastic modulus and arterial geometry cause critical variations in the local pressure and velocity gradients and secondary flow patterns. The pressure gradient change is  47%, with respect to the unstented baseline when the elastic modulus is increased to 10. The corresponding WSS change is 15.4%. We demonstrate that these changes are attributed to the production of vorticity (vorticity flux) caused by the wall movement and elastic modulus discontinuity. The changes in curvature dominate torsion changes in terms of the effects to local hemodynamics. The elastic modulus discontinuities along with the dynamic change in geometry affected the secondary flow patterns and vorticity flux, which in turn affects the WSS.


Author(s):  
Oxana A. Tkachenko ◽  
Svetlana A. Tkachenko ◽  
Victoria Timchenko ◽  
John A. Reizes ◽  
Guan Heng Yeoh ◽  
...  

Choonpa Igaku ◽  
2014 ◽  
Vol 41 (2) ◽  
pp. 155-163
Author(s):  
Yoshihiro SEO ◽  
Tomoko ISHIZU ◽  
Akiko ATSUMI ◽  
Ryo KAWAMURA ◽  
Kazutaka AONUMA

2004 ◽  
Vol 126 (5) ◽  
pp. 861-870 ◽  
Author(s):  
A. Thakur ◽  
X. Liu ◽  
J. S. Marshall

An experimental and computational study is performed of the wake flow behind a single yawed cylinder and a pair of parallel yawed cylinders placed in tandem. The experiments are performed for a yawed cylinder and a pair of yawed cylinders towed in a tank. Laser-induced fluorescence is used for flow visualization and particle-image velocimetry is used for quantitative velocity and vorticity measurement. Computations are performed using a second-order accurate block-structured finite-volume method with periodic boundary conditions along the cylinder axis. Results are applied to assess the applicability of a quasi-two-dimensional approximation, which assumes that the flow field is the same for any slice of the flow over the cylinder cross section. For a single cylinder, it is found that the cylinder wake vortices approach a quasi-two-dimensional state away from the cylinder upstream end for all cases examined (in which the cylinder yaw angle covers the range 0⩽ϕ⩽60°). Within the upstream region, the vortex orientation is found to be influenced by the tank side-wall boundary condition relative to the cylinder. For the case of two parallel yawed cylinders, vortices shed from the upstream cylinder are found to remain nearly quasi-two-dimensional as they are advected back and reach within about a cylinder diameter from the face of the downstream cylinder. As the vortices advect closer to the cylinder, the vortex cores become highly deformed and wrap around the downstream cylinder face. Three-dimensional perturbations of the upstream vortices are amplified as the vortices impact upon the downstream cylinder, such that during the final stages of vortex impact the quasi-two-dimensional nature of the flow breaks down and the vorticity field for the impacting vortices acquire significant three-dimensional perturbations. Quasi-two-dimensional and fully three-dimensional computational results are compared to assess the accuracy of the quasi-two-dimensional approximation in prediction of drag and lift coefficients of the cylinders.


Crystals ◽  
2020 ◽  
Vol 10 (1) ◽  
pp. 19 ◽  
Author(s):  
Yaping Tao ◽  
Ligang Han ◽  
Andong Sun ◽  
Kexi Sun ◽  
Qian Zhang ◽  
...  

Methyl-3-aminothiophene-2-carboxylate (matc) is a key intermediate in organic synthesis, medicine, dyes, and pesticides. Single crystal X-ray diffraction analysis reveals that matc crystallizes in the monoclinic crystal system P21/c space group. Three matc molecules in the symmetric unit are crystallographically different and further linked through the N–H⋯O and N–H⋯N hydrogen bond interactions along with weak C–H⋯S and C–H⋯Cg interactions, which is verified by the three-dimensional Hirshfeld surface, two-dimensional fingerprint plot, and reduced density gradient (RDG) analysis. The interaction energies within crystal packing are visualized through dispersion, electrostatic, and total energies using three-dimensional energy-framework analyses. The dispersion energy dominates in crystal packing. To better understand the properties of matc, electrostatic potential (ESP) and frontier molecular orbitals (FMO) were also calculated and discussed. Experimental and calculation results suggested that amino and carboxyl groups can participate in various inter- and intra-interactions.


2008 ◽  
Vol 9 (2) ◽  
pp. 83-104 ◽  
Author(s):  
T. Bodnár ◽  
A. Sequeira

The process of platelet activation and blood coagulation is quite complex and not yet completely understood. Recently, a phenomenological meaningful model of blood coagulation and clot formation in flowing blood that extends existing models to integrate biochemical, physiological and rheological factors, has been developed. The aim of this paper is to present results from a computational study of a simplified version of this coupled fluid-biochemistry model. A generalized Newtonian model with shear-thinning viscosity has been adopted to describe the flow of blood. To simulate the biochemical changes and transport of various enzymes, proteins and platelets involved in the coagulation process, a set of coupled advection–diffusion–reaction equations is used. Three-dimensional numerical simulations are carried out for the whole model in a straight vessel with circular cross-section, using a finite volume semi-discretization in space, on structured grids, and a multistage scheme for time integration. Clot formation and growth are investigated in the vicinity of an injured region of the vessel wall. These are preliminary results aimed at showing the validation of the model and of the numerical code.


Sign in / Sign up

Export Citation Format

Share Document