Effects of extracellular fiber architecture on cell membrane shear stress in a 3D fibrous matrix

Interstitial fluid flow has been shown to affect the organization and behavior of cells in 3D environments in vivo and in vitro, yet the forces driving such responses are not clear. Due to the complex architecture of the extracellular matrix (ECM) and the difficulty of measuring fluid flow near cells embedded in it, the levels of shear stress experienced by cells in this environment are typically estimated using bulk-averaged matrix parameters such as hydraulic permeability. While this is useful for estimating average stresses, it cannot yield insight into how local matrix fiber architecture-which is cell-controlled in the immediate pericellular environment-affects the local stresses imposed on the cell surface. To address this, we used computational fluid dynamics to study flow through an idealized mesh constructed of a cubic lattice of fibers simulating a typical in vitro collagen gel. We found that, in such high porosity matrices, the fibers strongly affect the flow fields near the cell, with peak shear stresses up to five times higher than those predicted by the Brinkman equation. We also found that minor remodeling of the fibers near the cell surface had major effects on the shear stress profile on the cell. These findings demonstrate the importance of fiber architecture to the fluid forces on a cell embedded in a 3D matrix, and also show how small modifications in the local ECM can lead to large changes in the mechanical environment of the cell.     

Computational predictions of the embolus-trapping performance of an IVC filter in patient-specific and idealized IVC geometries

Embolus transport simulations are performed to investigate the dependence of inferior vena cava (IVC) filter embolus-trapping performance on IVC anatomy. Simulations are performed using a resolved two-way coupled computational fluid dynamics/six-degree-of-freedom approach. Three IVC geometries are studied: a straight-tube IVC, a patient-averaged IVC, and a patient-specific IVC reconstructed from medical imaging data. Additionally, two sizes of spherical emboli (3 and 5 mm in diameter) and two IVC orientations (supine and upright) are considered. The embolus-trapping efficiency of the IVC filter is quantified for each combination of IVC geometry, embolus size, and IVC orientation by performing 2560 individual simulations. The predicted embolus-trapping efficiencies of the IVC filter range from 10 to 100%, and IVC anatomy is found to have a significant influence on the efficiency results (\(P < 0.0001\)). In the upright IVC orientation, greater secondary flow in the patient-specific IVC geometry decreases the filter embolus-trapping efficiency by 22–30 percentage points compared with the efficiencies predicted in the idealized straight-tube or patient-averaged IVCs. In a supine orientation, the embolus-trapping efficiency of the filter in the idealized IVCs decreases by 21–90 percentage points compared with the upright orientation. In contrast, the embolus-trapping efficiency is insensitive to IVC orientation in the patient-specific IVC. In summary, simulations predict that anatomical features of the IVC that are often neglected in the idealized models used for benchtop testing, such as iliac vein compression and anteroposterior curvature, generate secondary flow and mixing in the IVC and influence the embolus-trapping efficiency of IVC filters. Accordingly, inter-subject variability studies and additional embolus transport investigations that consider patient-specific IVC anatomy are recommended for future work.     

Depth-dependent analysis of the role of collagen fibrils, fixed charges and fluid in the pericellular matrix of articular cartilage on chondrocyte mechanics

The pericellular matrix of articular cartilage has been shown to regulate the mechanical environment of chondrocytes. However, little is known about the mechanical role of collagen fibrils in the pericellular matrix, and how fibrils might help modulate strains acting on chondrocytes when cartilage is loaded. The primary objective was to clarify the effect of pericellular collagen fibrils on cell volume changes and strains during cartilage loading. Secondary objectives were to investigate the effects of pericellular fixed charges and fluid on cell responses. A microstructural model of articular cartilage, in which chondrocytes and pericellular matrices were represented with depth-dependent structural and morphological properties, was created. The extracellular matrix and pericellular matrices were modeled as fibril-reinforced, biphasic materials with swelling capabilities, while chondrocytes were assumed to be isotropic and biphasic with swelling properties. Collagen fibrils in the extracellular matrix were represented with an arcade-like architecture, whereas pericellular fibrils were assumed to run tangential to the cell surface. In the early stages of a stress-relaxation test, pericellular fibrils were found to sensitively affect cell volume changes, even producing a reversal from increasing to decreasing cell volume with increasing fibril stiffness in the superficial zone. Consequently, steady-state volume of the superficial zone cell decreased with increasing pericellular fibril stiffness. Volume changes in the middle and deep zone chondrocytes were smaller and opposite to those observed in the superficial zone chondrocyte. An increase in the pericellular fixed charge density reduced cell volumes substantially in every zone. The sensitivity of cell volume changes to pericellular fibril stiffness suggests that pericellular fibrils play an important, and as of yet largely neglected, role in regulating the mechanical environment of chondrocytes, possibly affecting matrix synthesis during cartilage development and degeneration, and affecting biosynthetic responses associated with articular cartilage loading.     

Examination of continuum and micro-structural properties of human vertebral cancellous bone using combined cellular solid models

Two- and three-dimensional structural models of the vertebral body have been used to estimate the mechanical importance of parameters that are difficult to quantify experimentally such as lattice disorder, trabecular thickness, trabecular spacing, connectivity, and fabric. Many of the models that investigate structure–function relationships of the vertebral body focus only on the trabecular architecture and neglect solid–fluid interactions. We developed a cellular solid model composed of two idealized unit cell geometries to investigate the continuum and micro-structural properties of human vertebral cancellous bone in a mathematically tractable model. Using existing histomorphological data we developed structure–function relationships for the mechanical properties of the solid phase, estimated the micro-structural strains, and predicted the fluid flow characteristics. We found that the micro-structural strains may be 1.7 to 2.2 times higher than the continuum level strains between the ages of 40 and 80. In addition, the predicted permeability agrees well with the experimental data.     

Influence of stent configuration on cerebral aneurysm fluid dynamics

Embolic coiling is the most popular endovascular treatment available for cerebral aneurysms. Nevertheless, the embolic coiling of wide-neck aneurysms is challenging and, in many cases, ineffective. Use of highly porous stents to support coiling of wide-neck aneurysms has become a common procedure in recent years. Several studies have also demonstrated that high porosity stents alone can significantly alter aneurysmal hemodynamics, but differences among different stent configurations have not been fully characterized. As a result, it is usually unclear which stent configuration is optimal for treatment. In this paper, we present a flow study that elucidates the influence of stent configuration on cerebral aneurysm fluid dynamics in an idealized wide-neck basilar tip aneurysm model. Aneurysmal fluid dynamics for three different stent configurations (half-Y, Y and, cross-bar) were first quantified using particle image velocimetry and then compared. Computational fluid dynamics (CFD) simulations were also conducted for selected stent configurations to facilitate validation and provide more detailed characterizations of the fluid dynamics promoted by different stent configurations. In vitro results showed that the Y stent configuration reduced cross-neck flow most significantly, while the cross-bar configuration reduced velocity magnitudes within the aneurysmal sac most significantly. The half-Y configuration led to increased velocity magnitudes within the aneurysmal sac at high parent-vessel flow rates. Experimental results were in strong agreement with CFD simulations. Simulated results indicated that differences in fluid dynamic performance among the different stent configurations can be attributed primarily to protruding struts within the bifurcation region.     

FRAP analysis of membrane-associated proteins: lateral diffusion and membrane-cytoplasmic exchange

Obtaining quantitative kinetic parameters from fluorescence recovery after photobleaching (FRAP) experiments generally requires a theoretical analysis of protein mobility and appropriate solutions for FRAP recovery derived for a given geometry. Here we provide a treatment of FRAP recovery for a molecule undergoing a combined process of reversible membrane association and lateral diffusion on the plasma membrane for two commonly used bleach geometries: stripes, and boxes. Such analysis is complicated by the fact that diffusion of a molecule during photobleaching can lead to broadening of the bleach area, resulting in significant deviations of the actual bleach shape from the desired bleach geometry, which creates difficulty in accurately measuring kinetic parameters. Here we overcome the problem of deviations between actual and idealized bleach geometries by parameterizing, more accurately, the initial postbleach state. This allows for reconstruction of an accurate and analytically tractable approximation of the actual fluorescence distribution. Through simulated FRAP experiments, we demonstrate that this method can be used to accurately measure a broad range of combinations of diffusion constants and exchange rates. Use of this method to analyze the plextrin homology domain of PLC-δ1 in Caenorhabditis elegans results in quantitative agreement with prior analysis of this domain in other cells using other methods. Because of the flexibility, relative ease of implementation, and its use of standard, easily obtainable bleach geometries, this method should be broadly applicable to investigation of protein dynamics at the plasma membrane.     

The influence of cell mechanics, cell-cell interactions, and proliferation on epithelial packing

BACKGROUND: Epithelial junctional networks assume packing geometries characterized by different cell shapes, neighbor number distributions and areas. The development of specific packing geometries is tightly controlled; in the Drosophila wing epithelium, cells convert from an irregular to a hexagonal array shortly before hair formation. Packing geometry is determined by developmental mechanisms that likely control the biophysical properties of cells and their interactions. RESULTS: To understand how physical cellular properties and proliferation determine cell-packing geometries, we use a vertex model for the epithelial junctional network in which cell packing geometries correspond to stable and stationary network configurations. The model takes into account cell elasticity and junctional forces arising from cortical contractility and adhesion. By numerically simulating proliferation, we generate different network morphologies that depend on physical parameters. These networks differ in polygon class distribution, cell area variation, and the rate of T1 and T2 transitions during growth. Comparing theoretical results to observed cell morphologies reveals regions of parameter space where calculated network morphologies match observed ones. We independently estimate parameter values by quantifying network deformations caused by laser ablating individual cell boundaries. CONCLUSIONS: The vertex model accounts qualitatively and quantitatively for the observed packing geometry in the wing disc and its response to perturbation by laser ablation. Epithelial packing geometry is a consequence of both physical cellular properties and the disordering influence of proliferation. The occurrence of T2 transitions during network growth suggests that elimination of cells from the proliferating disc epithelium may be the result of junctional force balances.     

Metabolic pathways of the fossil dinosaur bones. Part V. Morphological differentiation of osteocyte lacunae and bone canaliculi and their significance in the system of extracellular communication

The study was carried out on dinosaur bones nearly 80 million years old. Samples for examination were prepared with specially elaborated methods. The light and transmission electron microscopic images permitted two kinds of bone lacunae and two types of paralacunar canalicular endings to be distinguished. The lacunae of the first kind were characterized by their elongated shape, their length exceeding their width several times, their dimensions being 31.2/9.4 microns. The lacunae of the other kind were not so long, their mean measurements amounting to 21.32/9.7 microns. Among the paralacunar canalicular endings those of small diameter were more numerous. The canaliculi of wider, funnel-shaped endings amounted to two or three, they were usually localized in the polar part of the lacuna, and were defined as the axial canaliculi. These were canaliculi of a large diameter. The canalicular wall was constructed of collagen fibres. The same fibres were found in the lacunar wall. Also a relationship between the structure of the lacunar wall and the localization of an osteocyte in the lacuna was analysed in the light and electron microscopes. In regard to the structure of the bone lacuna and the localization of an osteocyte in it, zones A and B were distinguished. Zone A had a characteristic loose and disorderly system of collagen fibres building the lacunar wall. The fibres in this area were by nature open to view. Besides, this region of the lacunar wall revealed specific terraced hollows. Zone B was distinguished by a compact system of parallelly arranged collagen fibres, which formed characteristic ridges in the lacunar wall. The localization of the osteocyte in the lacuna was irregular, the pericellular space around it being of variable width. This space was shown to contain mucopolysaccharides. The images obtained from dinosaur bone were compared with those already known for modern bone. These comparisons permitted it to be ascertained that zone A corresponds to a spot in the lacuna in which the osteocyte exhibits a decreased activity. Zone B is the area of the actual direction of the osteocyte's activity aiming at the shaping of the wall of its lacuna. It can be supposed that the widened endings of the paralacunar canaliculi perform more important functions in conveyance, this being evident from comparisons of analogous areas in modern bone.     

A prediction of cell differentiation and proliferation within a collagen–glycosaminoglycan scaffold subjected to mechanical strain and perfusive fluid flow

Mesenchymal stem cell (MSC) differentiation can be influenced by biophysical stimuli imparted by the host scaffold. Yet, causal relationships linking scaffold strain magnitudes and inlet fluid velocities to specific cell responses are thus far underdeveloped. This investigation attempted to simulate cell responses in a collagen–glycosaminoglycan (CG) scaffold within a bioreactor. CG scaffold deformation was simulated using μ-computed tomography (CT) and an in-house finite element solver (FEEBE/linear). Similarly, the internal fluid velocities were simulated using the afore-mentioned μCT dataset with a computational fluid dynamics solver (ANSYS/CFX). From the ensuing cell-level mechanics, albeit octahedral shear strain or fluid velocity, the proliferation and differentiation of the representative cells were predicted from deterministic functions. Cell proliferation patterns concurred with previous experiments. MSC differentiation was dependent on the level of CG scaffold strain and the inlet fluid velocity. Furthermore, MSC differentiation patterns indicated that specific combinations of scaffold strains and inlet fluid flows cause phenotype assemblies dominated by single cell types. Further to typical laboratory procedures, this predictive methodology demonstrated loading-specific differentiation lineages and proliferation patterns. It is hoped these results will enhance in-vitro tissue engineering procedures by providing a platform from which the scaffold loading applications can be tailored to suit the desired tissue.     

A new view of mechanotransduction and strain amplification in cells with microvilli and cell processes

In this paper we shall describe new mechanical models for the deformation of the actin filament bundles in kidney microvilli and osteocytic cell processes to see whether these cellular extensions, like the stereocilia on hair cells in the inner ear, can function as mechanotransducers when subject to physiological flow. In the case of kidney microvilli we show that the hydrodynamic drag forces at the microvilli tip are <0.01 pN, but there is a 38-fold force amplification on the actin filaments at the base of the microvilli due to the resisting moment in its terminal web. This leads to forces that are more than sufficient to deform the terminal web complex of the microvillus where ezrin has been shown to couple the actin cytoskeleton to the Na(+)/H(+) exchanger. In the case of bone cell processes we show that the actin filament bundles have an effective Young's modulus that is 200 times > the measured modulus for the actin gel in the cell body. It is, therefore, unlikely that bone cell processes respond in vivo to fluid shear stress, as proposed in [59]. However, we show that the fluid drag forces on the pericellular matrix which tethers the cell processes to the canalicular wall can produce a 20-100 fold amplification of bone tissue strains in the actin filament bundle of the cell process.     

Diagnostic imaging of small amounts of pleural fluid: pleural effusion vs. physiologic pleural fluid

The aim of this article is to present an overview of our 10 years clinical research work and early clinical experience with small pleural effusions. Small amounts of pleural fluid are severely difficult to identify with imaging methods (chest x-rays and ultrasound). Nevertheless, it may be an important finding, sometimes leading to a definitive diagnosis of pleural carcinomatosis, infection or other pathologic condition. Chest x-rays were used for many years for the diagnosis of small pleural effusions. Lateral decubitus chest radiographs represented a gold standard for imaging of small amounts of plural fluid for more than 80 years. In the last two decades, ultrasonography of pleural space became a leading real-time method for demonstrating small pleural effusions. Furthermore, the advent of sonographic technology actually enables detection of physiologic pleural fluid in some otherwise healthy individuals. In conclusion, new definitions of the key terms in the field of diagnostic imaging of small amounts of pleural fluid seem to be justified. We suggest that the term pleural fluid should determine physiologic pleural space condition while the term pleural effusion should only be used in the cases of pleural involvement or pleural illness.     

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.     

Effect of inlet velocity profiles on patient-specific computational fluid dynamics simulations of the carotid bifurcation

Patient-specific computational fluid dynamics (CFD) is a powerful tool for researching the role of blood flow in disease processes. Modern clinical imaging technology such as MRI and CT can provide high resolution information about vessel geometry, but in many situations, patient-specific inlet velocity information is not available. In these situations, a simplified velocity profile must be selected. We studied how idealized inlet velocity profiles (blunt, parabolic, and Womersley flow) affect patient-specific CFD results when compared to simulations employing a "reference standard" of the patient's own measured velocity profile in the carotid bifurcation. To place the magnitude of these effects in context, we also investigated the effect of geometry and the use of subject-specific flow waveform on the CFD results. We quantified these differences by examining the pointwise percent error of the mean wall shear stress (WSS) and the oscillatory shear index (OSI) and by computing the intra-class correlation coefficient (ICC) between axial profiles of the mean WSS and OSI in the internal carotid artery bulb. The parabolic inlet velocity profile produced the most similar mean WSS and OSI to simulations employing the real patient-specific inlet velocity profile. However, anatomic variation in vessel geometry and the use of a nonpatient-specific flow waveform both affected the WSS and OSI results more than did the choice of inlet velocity profile. Although careful selection of boundary conditions is essential for all CFD analysis, accurate patient-specific geometry reconstruction and measurement of vessel flow rate waveform are more important than the choice of velocity profile. A parabolic velocity profile provided results most similar to the patient-specific velocity profile.     

Image-based axon model highlights heterogeneity in initiation of damage

Head injury simulations predict the occurrence of traumatic brain injury by placing a threshold on the calculated strains for axon tracts within the brain. However, a current roadblock to accurate injury prediction is the selection of an appropriate axon damage threshold. While several computational studies have used models of the axon cytoskeleton to investigate damage initiation, these models all employ an idealized, homogeneous axonal geometry. This homogeneous geometry with regularly spaced microtubules, evenly distributed throughout the model, overestimates axon strength because, in reality, the axon cytoskeleton is heterogeneous. In the heterogeneous cytoskeleton, the weakest cross section determines the initiation of failure, but these weak spots are not present in a homogeneous model. Addressing one source of heterogeneity in the axon cytoskeleton, we present a new semiautomated image analysis pipeline for using serial-section transmission electron micrographs to reconstruct the microtubule geometry of an axon. The image analysis procedure locates microtubules within the images, traces them throughout the image stack, and reconstructs the microtubule structure as a finite element mesh. We demonstrate the image analysis approach using a C. elegans touch receptor neuron due to the availability of high-quality serial-section transmission electron micrograph data sets. The results of the analysis highlight the heterogeneity of the microtubule structure in the spatial variation of both microtubule number and length. Simulations comparing this image-based geometry with homogeneous geometries show that structural heterogeneity in the image-based model creates significant spatial variation in deformation. The homogeneous geometries, on the other hand, deform more uniformly. Since no single homogeneous model can replicate the mechanical behavior of the image-based model, our results argue that heterogeneity in axon microtubule geometry should be considered in determining accurate axon failure thresholds.     

Modeling the large-scale geometry of human coronary arteries

Two principles suffice to model the large-scale geometry of normal human coronary arterial networks. The first principle states that artery diameters are set to minimize the power required to distribute blood through the network. The second principle states that arterial tree geometries are set to globally minimize the lumen volume. Given only the coordinates of an arterial tree's source and "leaves", the model predicts the nature of the network connecting the source to the leaves. Measurements were made of the actual geometries of arterial trees from postmortem healthy human coronary arteriograms. The tree geometries predicted by the model look qualitatively similar to the actual tree geometries and have volumes that are within a few percent of those of the actual tree geometries. Human coronary arteries are therefore within a few percent of perfect global volume optimality. A possible mechanism for this near-perfect global volume optimality is suggested. Also, the model performs best under the assumption that the flow is not entirely steady and laminar.     

A novel in vitro loading system to produce supraphysiologic oscillatory fluid shear stress

A multi-well fluid loading (MFL) system was developed to deliver oscillatory subphysiologic to supraphysiologic fluid shear stresses to cell monolayers in vitro using standard multi-well culture plates. Computational fluid dynamics modeling with fluid-structure interactions was used to quantify the squeeze film fluid flow between an axially displaced piston and the well plate surface. Adjusting the cone angle of the piston base modulated the fluid pressure, velocity, and shear stress magnitudes. Modeling results showed that there was near uniform fluid shear stress across the well with a linear drop in pressure across the radius of the well. Using the MFL system, RAW 264.7 osteoclastic cells were exposed to oscillatory fluid shear stresses of 0, 0.5, 1.5, 4, 6, and 17 Pa. Cells were loaded 1 h per day at 1 Hz for two days. Compared to sub-physiologic and physiologic levels, supraphysiologic oscillatory fluid shear induced upregulation of osteoclastic activity as measured by tartrate-resistant acid phosphatase activity and formation of mineral resorption pits. Cell number remained constant across all treatment groups.     

Interstitial fluid flow within bone canaliculi and electro-chemo-mechanical features of the canalicular milieu

Canalicular fluid flow is acknowledged to play a major role in bone functioning, allowing bone cells’ metabolism and activity and providing an efficient way for cell-to-cell communication. Bone canaliculi are small canals running through the bone solid matrix, hosting osteocyte’s dendrites, and saturated by an interstitial fluid rich in ions. Because of the small size of these canals (few hundred nanometers in diameter), fluid flow is coupled with electrochemical phenomena. In our previous works, we developed a multi-scale model accounting for coupled hydraulic and chemical transport in the canalicular network. Unfortunately, most of the physical and geometrical information required by the model is hardly accessible by nowadays experimental techniques. The goal of this study was to numerically assess the influence of the physical and material parameters involved in the canalicular fluid flow. The focus was set on the electro-chemo-mechanical features of the canalicular milieu, hopefully covering any in vivo scenario. Two main results were obtained. First, the most relevant parameters affecting the canalicular fluid flow were identified and their effects quantified. Second, these findings were given a larger scope to cover also scenarios not considered in this study. Therefore, this study gives insight into the potential interactions between electrochemistry and mechanics in bone and provides the rational for further theoretical and experimental investigations.     

Three-dimensional finite element modeling of pericellular matrix and cell mechanics in the nucleus pulposus of the intervertebral disk based on in situ morphology

Nucleus pulposus (NP) cells of the intervertebral disk (IVD) have unique morphological characteristics and biologic responses to mechanical stimuli that may regulate maintenance and health of the IVD. NP cells reside as single cell, paired or multiple cells in a contiguous pericellular matrix (PCM), whose structure and properties may significantly influence cell and extracellular matrix mechanics. In this study, a computational model was developed to predict the stress–strain, fluid pressure and flow fields for cells and their surrounding PCM in the NP using three-dimensional (3D) finite element models based on the in situ morphology of cell–PCM regions of the mature rat NP, measured using confocal microscopy. Three-dimensional geometries of the extracellular matrix and representative cell–matrix units were used to construct 3D finite element models of the structures as isotropic and biphasic materials. In response to compressive strain of the extracellular matrix, NP cells and PCM regions were predicted to experience volumetric strains that were 1.9–3.7 and 1.4–2.1 times greater than the extracellular matrix, respectively. Volumetric and deviatoric strain concentrations were generally found at the cell/PCM interface, while von Mises stress concentrations were associated with the PCM/extracellular matrix interface. Cell–matrix units containing greater cell numbers were associated with higher peak cell strains and lower rates of fluid pressurization upon loading. These studies provide new model predictions for micromechanics of NP cells that can contribute to an understanding of mechanotransduction in the IVD and its changes with aging and degeneration.     

Stereoscopic particle image velocimetry analysis of healthy and emphysemic alveolar sac models

Emphysema is a progressive lung disease that involves permanent destruction of the alveolar walls. Fluid mechanics in the pulmonary region and how they are altered with the presence of emphysema are not well understood. Much of our understanding of the flow fields occurring in the healthy pulmonary region is based on idealized geometries, and little attention has been paid to emphysemic geometries. The goal of this research was to utilize actual replica lung geometries to gain a better understanding of the mechanisms that govern fluid motion and particle transport in the most distal regions of the lung and to compare the differences that exist between healthy and emphysematous lungs. Excised human healthy and emphysemic lungs were cast, scanned, graphically reconstructed, and used to fabricate clear, hollow, compliant models. Three dimensional flow fields were obtained experimentally using stereoscopic particle image velocimetry techniques for healthy and emphysematic breathing conditions. Measured alveolar velocities ranged over two orders of magnitude from the duct entrance to the wall in both models. Recirculating flow was not found in either the healthy or the emphysematic model, while the average flow rate was three times larger in emphysema as compared to healthy. Diffusion dominated particle flow, which is characteristic in the pulmonary region of the healthy lung, was not seen for emphysema, except for very small particle sizes. Flow speeds dissipated quickly in the healthy lung (60% reduction in 0.25 mm) but not in the emphysematic lung (only 8% reduction 0.25 mm). Alveolar ventilation per unit volume was 30% smaller in emphysema compared to healthy. Destruction of the alveolar walls in emphysema leads to significant differences in flow fields between the healthy and emphysemic lung. Models based on replica geometry provide a useful means to quantify these differences and could ultimately improve our understanding of disease progression.     

A multiscale computational fluid dynamics approach to simulate the micro-fluidic environment within a tissue engineering scaffold with highly irregular pore geometry

Mechanical stimulation can regulate cellular behavior, e.g., differentiation, proliferation, matrix production and mineralization. To apply fluid-induced wall shear stress (WSS) on cells, perfusion bioreactors have been commonly used in tissue engineering experiments. The WSS on cells depends on the nature of the micro-fluidic environment within scaffolds under medium perfusion. Simulating the fluidic environment within scaffolds will be important for gaining a better insight into the actual mechanical stimulation on cells in a tissue engineering experiment. However, biomaterial scaffolds used in tissue engineering experiments typically have highly irregular pore geometries. This complexity in scaffold geometry implies high computational costs for simulating the precise fluidic environment within the scaffolds. In this study, we propose a low-computational cost and feasible technique for quantifying the micro-fluidic environment within the scaffolds, which have highly irregular pore geometries. This technique is based on a multiscale computational fluid dynamics approach. It is demonstrated that this approach can capture the WSS distribution in most regions within the scaffold. Importantly, the central process unit time needed to run the model is considerably low.