|Home | About | Journals | Submit | Contact Us | Français|
Global measurements of the lung provided by standard pulmonary function tests do not give insight into the regional basis of lung function and lung disease. Advances in imaging methodologies, computer technologies, and subject-specific simulations are creating new opportunities for studying structure-function relationships in the lung through multi-disciplinary research. The digital Human Lung Atlas is an imaging-based resource compiled from male and female subjects spanning several decades of age. The Atlas comprises both structural and functional measures, and includes computational models derived to match individual subjects for personalized prediction of function. The computational models in the Atlas form part of the Lung Physiome project, which is an international effort to develop integrative models of lung function at all levels of biological organization. The computational models provide mechanistic interpretation of imaging measures; the Atlas provides structural data upon which to base model geometry, and functional data against which to test hypotheses. The example of simulating air flow on a subject-specific basis is considered. Methods for deriving multi-scale models of the airway geometry for individual subjects in the Atlas are outlined, and methods for modeling turbulent flows in the airway are reviewed.
Understanding the complex dynamics of the pulmonary system and the lung’s response to disease, injury, and treatment is fundamental to the advancement of pulmonary medicine; however the global measurements obtained via clinical pulmonary function tests do not adequately capture lung complexity and may only be minimally altered by significant local disease. Quantifiable image-based measures, including evaluation of static and dynamic structure and function, are now recognized as very sensitive indicators of localized subclinical disease and appear to describe complex lung processes much better than global measurements . Dramatic advances in imaging methodologies including multidetector row CT (MDCT), hyperpolarized gas MRI, positron emission tomography (PET), and single photon emission tomography (SPECT), in parallel with advances in computer technologies and subject-specific simulations is creating new opportunities for studying structure-function relationships in the lung through truly multi-disciplinary research. For example, development of a subject-specific model for ventilation of a breathing lung can only be accomplished through multi-disciplinary efforts that require expertise in medical imaging, airway geometric reconstruction, computational techniques, pulmonary physiology and medicine, and fluid mechanics.
Anatomically and biophysically based computational models of human or animal physiology provide a means to integrate vast amounts of data across many spatial and temporal scales into a framework that can be linked to whole-body physiology and clinical medicine. Development of pulmonary models that can be customized to a human or animal subject means that an individual’s structure can be linked to personalized predictions of function [2–5]. These developments form part of the Lung Physiome Project (www.physiome.org.nz), a public domain effort to develop integrative models at all levels of biological organization, from genes to the whole organism. The lung is a unique internal organ in that with each breath it interacts directly with the environment. This can have long term consequences on the structure of the airway tree , or an immediate consequence for the asthmatic lung. Pathologic changes in the lung are slow to develop, heterogeneous, and multiple pathologic processes intermingle; this points to a fundamental need for sensitive, quantitative measures of phenotypic outcomes, and interpretive computational models to differentiate between interacting mechanisms.
A goal of the Lung Physiome is to curate peer-reviewed models and data in a standardized format, facilitating model sharing and minimizing the need for scientists to reimplement existing models. A parallel project developing a digital Human Lung Atlas [7–10] aims to use imaging measures to quantify normal lung structure and pathological ranges of departure from normal, along with standard and imaging-based functional measurements. The Human Lung Atlas provides the primary imaging database for the Lung Physiome, and conversely the Physiome models provide an interpretive tool to link structure and function in the Atlas. This review focuses primarily on the link between the Atlas and the Physiome, considering the derivation of subject-specific models of the airway tree, and their emerging application in studying individualized computation of air flow.
The conducting airways in the mammalian lung comprise an extensive asymmetric branching structure, however all computational studies of function in the conducting airways rely on at least some simplifications of the airway geometry due to constraints on computer memory and computing time. The airway tree can be simplified by assuming that the tree has symmetry  or regular asymmetry of branching  and/or by modelling it as a one-dimensional (1D) structure and integrating the governing equations for function over the circular cross-section of the airways [13–15], or by only simulating within a subset of the tree in three-dimensional (3D) space [16, 17].
Two classic simplifications of the airway tree are Weibel’s symmetric model A , and the Horsfield Δ model with regular asymmetry of branching . Both models are based on measurements made on casts of the airway tree [11, 18]. Both of these models provide a convenient simplification of the complex branching airway tree, and are still widely used for computational studies because they summarize detailed anatomical measurements into idealized structures that are relatively easily used and computationally inexpensive. These ‘simplified’ models are necessary for computing analytic solutions, however, neither model directly includes spatial information, and they cannot be used to predict function on a subject-specific basis.
Weibel and Horsfield relied on measuring airways from casts of the airway tree. Modern imaging provides a high resolution tool with which to measure and define the lung. Early imaging-based studies measured the geometry of the airway tree ; later work used imaging to derive models of the airway tree geometry and connectivity [4, 20–23]. These anatomically based models are specific to an imaged individual rather than seeking to be representative of a population. The airway model geometry can be anatomically-accurate in two ways: 1) it follows the centerlines of the airways to define a 1D tree [4, 20], or 2) it conforms to surface data of the airways segmented from medical imaging [21–23].
The first type of anatomically-accurate, spatially distributed 1D model for the airway tree can be derived by tracing the bifurcating airway paths using imaging of airway casts from an excised lung  or imaging of the in vivo lung . Schmidt et al.  derived a digital reference model of the human bronchial tree over 17 Horsfield orders, using CT scanning of an airway cast. This study provided further insight into airway branching geometry, demonstrating multi-fractal properties in the bronchial airway tree, but it is not a method that can be used to generate geometric models for the in vivo lung. Using an airway cast to define the branching geometry has the advantage of defining connectivity through an extensive airway tree, but casting potentially introduces geometric distortion of the tissue, is representative of a single filling pressure, and it does not extend through the entire conducting airway system. It is also not applicable to deriving individualized models for living subjects.
X-ray CT is currently the imaging modality of choice for comprehensive evaluation of the lung. Scanners are now capable of sub half second data acquisition (330 msec per rotation and faster), allowing for the imaging of not only anatomy  but also ventilation [25–27] and perfusion , providing unprecedented structure-to-function correlations. Conventional x-ray CT is a sensitive predictor of emphysema , smoking-related pathophysiology , and pulmonary hypertension , compared with standards such as pulmonary function tests. An electronic CT-based deformable ‘Human Lung Atlas’ is currently being compiled as a comprehensive tool for image and model based analysis of lung disease [7–10], supported through a Bioengineering Research Partnership through the National Heart, Lung, and Blood Institute of the National Institutes of Health. Medical atlases have traditionally been illustrated texts that fully characterize an anatomical structure. With the advent of modern imaging, atlases can now be based on high resolution volumetric data, including structure and function from a range of imaging modalities. Examples of imaging-based brain atlases can be found at www.brainmaps.org and www.brain-map.org.
The Human Lung Atlas includes male and female subjects spanning several decades of life, integrating both structural and functional information from MDCT imaging and spirometry. The human female lung is smaller in volume relative to body size than the male lung , has smaller mean airway diameters  and a smaller diffusion surface. Functional gender differences include greater expiratory flow limitation and increased work of breathing in the female lung during heavy exercise . One goal of the Atlas is to establish the range of normality of an extensive set of both structural and functional measures of the adult human lung and to catalogue these normative values in a form which can then serve as the basis for detection of early pathology in an individual, as well as a basis for tracking the changes in an individual lung over time. A further goal is to partner imaging-based measures with computational model-based interpretation of physiological mechanisms, for analysis of structure-function relationships in lung disease.
The natural contrast between the low-density lungs and the surrounding high-density chest wall can be used to guide image segmentation . Lung segmentation in the Atlas uses a fully automatic and validated algorithm  to first extract the lung region by gray-level thresholding, then separate the left and right lungs by identifying the anterior and posterior junctions using dynamic programming, and finally smoothing the irregular boundary along the mediastinum to obtain results consistent with manual analysis. The lobar fissures are low contrast surfaces with blurred boundaries in the cross-sectional CT images, and may be incomplete or vary between subjects. Computer-based detection of the fissures is complicated by surrounding vessels or other structures, and noise and artifacts in the images. The Atlas previously used a semi-automatic lobar fissure detection method , but currently utilizes an anatomic pulmonary atlas constructed from 16 volumetric CT scans of 16 subjects to guide fully automatic lobe segmentation . Segmentation and analysis of the intra-thoracic airway trees utilizes the Volumetric Imaging Display and Analysis (VIDA) software . Several iterations of automated algorithms for segmentation of the human airway tree have been developed for the Atlas and incorporated into VIDA. The most recent generations are based on a combination of 3D region growing and two-dimensional (2D) mathematical morphology , and later incorporating fuzzy logic . Using this approach , airway trees can be successfully extracted from most clinical quality standard-dose CT data sets. Segmentation of the human airway tree at close to TLC reliably yields between six to nine generations of branches, and often up to twelve generations depending on image resolution and branch diameter. Many more branch generations can be segmented from a monopodial airway tree – such as that of the ovine lung – due to large-diameter central airways that taper without significantly changing direction .
Extracting the airway tree is only a first step: to perform topology-based analyses, airway tree branches must be identified, either manually  or automatically . The binary airway tree formed in the segmentation step is skeletonized to identify the 3D centerlines of individual branches and branchpoint locations. A sequential 3D thinning algorithm [43, 44] was customized for this application. After thinning, the skeleton is smoothed, pruned, and the complete tree converted into a graph structure. The airway-tree branch-point matching algorithm utilized for the Atlas is able to find corresponding branch-points between two intra-subject scans in-vivo [40, 45]. An automatic method is used for anatomical labeling of the airways. The method is capable of assigning all 33 commonly used anatomical names to a human airway-tree .
Anatomical measurements of the minor and major airway diameter and cross-sectional area are made in the original gray-level MDCT volume . Measurements are made in two-dimensional slices that are resampled perpendicular to the airway centerline. The measurement method has been demonstrated to achieve subvoxel accuracy (typical error ~0.1 mm) on a Plexiglas phantom, independent of airway orientation and X-ray dose .
Tawhai et al.  derive subject-specific airway models for the Atlas subject database by supplementing in vivo imaging-based airway geometry with a volume-filling branching (VFB) method to fill a lung volume with an airway-consistent tree. Combining in vivo imaging with supplemental airways has the advantage of producing a model that is specific to an imaged subject. The airway trees created with the VFB method have been shown to be morphometrically consistent with measurements from cast [11, 12, 48] and imaging  studies. Individual airway trees generated for human and ovine subjects have geometry appropriate to their respective species: the human airways form a relatively symmetric bifurcating tree, whereas the sheep airways branch monopodially and hence are far more asymmetric. The generation algorithm reproduces these features in response to the shape of the lung or lobe boundary . This algorithm is the only method to date that produces airway trees in anatomically-realistic lung shapes, and is the only method that has been used to generate non-human airway trees. The method has also been extended to model the pulmonary vascular trees, including the supernumerary vessels . A model for one Atlas subject is shown in Figure 1.
The resulting model is a 1D tree with branches distributed in 3D space. To use the model for computation, diameters must be defined through the entire structure. For the imaging-based (uppermost) airways the diameters are assigned directly from calculating the diameter based on measurement of airway cross-sectional area in a plane orthogonal to the central axis. The algorithm-based airways can have diameters assigned using Horsfield or Strahler ordering  weighted against a ratio to parent branch diameter.
The model from  is a centreline (1D) tree with associated branch diameters. It is therefore equivalent in basic structure to the Weibel and Horsfield models, but with subject-specific geometry and spatial positioning of the airways. This becomes important when considering interacting function - such as the effect of changing tissue properties on airway tethering and airway collapse – or for comparison of simulated results with spatially distributed experimental results .
The physiologically consistent spatial distribution of airways in this 1D tree structure makes it an adequate and appropriate model for many simulation studies [2, 5, 49]. However, simulation of 3D flow and pressure via computational fluid dynamics (CFD) techniques requires a 3D structure that incorporates the salient features of the airway geometry. Previous models of the airway tree have simplified the computational domain by representing the airways as smooth cylinders, and/or have considered only a few of the up to 30 generations in the conducting airway tree [50–53]. The reason for this is clear: constructing a CFD-appropriate computational mesh of the airways is complex and time consuming, and imaging data for the lower airways is generally not attainable.
3D models of airway geometry will be considered in the next section on CFD simulation.
The lungs are increasingly recognized as a potential route for delivery of systemic drugs such as insulin , because this avoids the hepatic metabolic pathway and may better simulate the physiological response to a food bolus. Penetration and deposition of particles within the airways depends on airway size and branching patterns, which vary between species . Variation in individual airway geometry therefore makes subject-specific models essential for the study of pulmonary air flow and drug delivery. Recent evidence also suggests that early exposure to environmental pollutants has chronic, adverse effects on lung development in children . Furthermore, it has been demonstrated that a strong interaction exists between lung geometry and gas properties , which has major implications in determining gas delivery to and clearance from the lung periphery during ventilation imaging via x-ray CT using xenon gas [25, 27, 56], or MRI using hyperpolarized helium gas [57–60]. The ability to predict air flow and particle deposition on a subject-specific basis is therefore necessary for understanding the correlation between structure and function, for assessing individual differences in vulnerability to airborne pollutants, and to provide insights into inter-subject differences in regional lung function and underlying mechanisms of pathologic developments.
CFD is now becoming a viable tool for simulation of airflow in the human respiratory tract and has been pursued by a number of researchers. The earliest studies of airway CFD used smooth cylindrical tubes merging in a single bifurcation , building to several bifurcations [62–65], and then incorporating imaging-defined airway surface data [21–23, 66–68]. Most of the emphasis of numerical studies has been placed on flow structure and particle deposition for drug delivery. Limitations of airway CFD studies that will be considered here are: 1) the representation of airway geometry, and 2) the treatment of turbulence. Other limitations of CFD studies not considered here but studied by others are the effects of mesh structure and grid convergence , the accuracy of particle tracking algorithms for sedimentation and diffusion in commercial software packages , and the inlet velocity profile . A partial list of earlier CFD studies is given by , summarizing the computational packages used, parameters tested, and success in predicting experimental data. Many studies of airway CFD have been presented in the literature, for example references [16, 17, 21, 50–53, 61–65, 68–74] to name a few; the ones primarily considered here are the most recent studies that are based on medical imaging in an attempt to move towards subject-specific predictions of flow and particle transport.
Several studies have demonstrated the importance of anatomically accurate geometry for valid predictions of flow and particle distribution [16, 75]. Studies based on single and double bifurcation models have been concluded insufficient for analyzing particle transport and deposition in the bronchial airways [52, 71]. The velocity profiles in daughter generations after a bifurcation have been shown to be affected by the angle of rotation, viz. the azimuthal angle, between consecutive bifurcations [17, 75]. Thus, simulations that do not consider the azimuthal angle are limited in their scope. Nowak et al.  demonstrated major differences between Weibel symmetric and imaging-based models derived from a cadaver cast. Using the commercial software FLUENT®, they calculated airflow in generations 0–9 of a Weibel-based airway tree separated into three simulation subunits (generations 0–3, 3–6, and 6–9). Comparison with computation in a more realistic geometry showed that the velocity fields in the more accurate geometry are more complicated than predicted in the symmetric geometry. Furthermore, there was no location-to-location correspondence between the velocity field in the anatomical geometry and the symmetric geometry, and the Weibel symmetric model was found to be insufficient for predicting aerosol deposition in the human lungs.
Recognizing the need for a more accurate description of the airway bifurcation than is provided by symmetric models, Hegedus  presented a mathematical description of a morphologically realistic airway bifurcation, designed for CFD study, and used their method to merge several bifurcations into multiple airway geometries. One important feature of the model was to enforce a smooth transition between the airways, and rounding at the carina. Farkas et al.  used the model from  to piece together a model from generation 1 to generation 5 in the right upper lobe. They found that aerosol deposition is highly dependent on airway geometry, breathing parameters, and particle characteristics.
Lin et al.  demonstrated the process of constructing a realistic geometry of the upper respiratory tract from MDCT images measured from a subject in the Human Lung Atlas. The geometry included a mouth piece, the mouth, the oropharynx, the larynx, and the airways of ~6 generations. Constructing this mesh involved triangulation of the surface geometry in volumetric MDCT imaging and exporting to an STL file, which is a standard data format used in the rapid prototyping industry. The original lung geometry in this format required surface repairing, remeshing, and clipping due to rough surfaces and bad edges. Numerical simulations of airflow in the mouth-to-airway geometry were computed to study the dependence of airflow in the intra-thoracic airways on the complex oropharynx and larynx geometry. Results were compared with simulations in a model comprising only the intra-thoracic airway tree of the same subject. Most prior numerical simulations are based on idealized airway models without including the oropharynx region, therefore requiring some sort of inlet boundary condition to be specified at the trachea. Comparison of results from the two geometries suggested that turbulence induced by the laryngeal jet (and therefore only present in the mouth-to-airway model) can significantly affect airway flow patterns. In the mouth-to-airway model, a curved sheet-like turbulent laryngeal jet was formed, impinging on the rear side of the trachea with high turbulence intensity in the subglottic space and the upper trachea. The velocity distribution in the intra-thoracic airway model was a parabola with its maximum located along the centerline of the trachea and the turbulence intensity essentially zero, thus the flow was laminar. The velocity characteristics in the left and right main-stem bronchi between the two cases were also different.
Several other studies have also used imaging-based central airway geometries for analysis, [22, 23, 66], however the study of airflow in the lung must ultimately be able to span from the mouth to the alveolated airways in the lung periphery. The smaller airways – beyond approximately generations 6 to 9 in the human lung – cannot usually be visualized with current clinical imaging. Lin et al.  presented a method for creating subject-specific 3D and 1D coupled airway mesh structures with seamless transition between the 3D and 1D scales, incorporating the desired level of geometric detail wherever it is needed in the airway tree. Subject-specific 1D models for the entire conducting airway tree  are converted to a high-order (cubic Hermite) 2D surface mesh of the entire domain. The parent and child branches merge at the bifurcations with a smooth, continuous surface. The uppermost airways – for which there is MDCT surface data - are geometry fitted to enforce accurate airway surface geometry. The resulting 2D surface mesh is continuous with the surface of the algorithm-based airways. Any portion of this surface mesh can be converted to a 3D CFD-ready mesh by selecting a region of interest to study in detail. The 3D mesh is created only within those airways, and the remainder of the domain is the original 1D tree. The 3D and 1D trees are a single continuous model that has a different dimension in specified regions. An advantage of this approach is that transport can be studied in - for example – models that include successive generations, to define the conditions under which a 1D model representation will suffice.
Figure 2 illustrates the creation of a 3D CFD mesh for the subject in Figure 1. In Figure 2a a subset of airways are isolated from the full conducting airway tree (in this case every airway in generations 1 to 11), and two selected pathways that extend to the terminal bronchioles (one each in the right upper and right lower lobes). Figure 2b and c illustrate a coarse 3D CFD mesh generated within a single bifurcation structure.
Idealized models of the airway structure generally assume a circular cross-section. The airway cross-section can be quite irregular in shape, particularly in the trachea [78, 79], and would perhaps be better represented by an elliptical cross-section. The impact that this might have on simulation results remains to be studied.
Finlay  estimated the Reynolds number in various generations of simplified cylinder-based airway models, showing that the Reynolds number may be above 2,000 in the larynx - in the turbulent regime - and that the turbulence produced in the oropharynx can be convected to several generations of the tracheobronchial tree.
One approach to avoid computing complex turbulent flows is to simulate flow in small airways that are presumably well into the laminar regime [50–52, 61, 71]. This is not entirely satisfactory for two reasons: 1) in relatively short airways the flow may never be fully developed, and turbulent eddies or flow disturbance may propagate into small airways that would nominally carry laminar flow, and 2) developing an understanding of flow transport phenomena in the airway tree ultimately requires inclusion of all structures from the mouth to the most distal airways. Airway CFD studies that have included turbulence have almost all adopted Reynolds-Averaged Navier-Stokes models (RANS) models. RANS solves the Reynolds-Averaged Navier-Stokes equations, in which only the mean velocity field is resolved and all of the unsteadiness is averaged out. Thus turbulent structures are completely absent; they are simply represented by the Reynolds stress tensor. This tensor can be modeled, for example, by the k-ε model, where k is the turbulent kinetic energy and ε is the turbulent dissipation rate . Experimental and numerical studies have shown that RANS models lack accuracy for this type of study. For example, Heenan et al.  compared measurements of the velocity field in the central plane of the human oropharynx – captured using particle image velocimetry - with numerical results computed using CFX-Tascflow® (AEA Technology Engineering) software, which is based on the RANS model. They found that although the RANS model reproduced the basic features of the flow, it did not capture the increased viscous effects at lower Reynolds numbers, suggesting the need for more development and validation of the RANS CFD methods. Stapleton et al.  also evaluated the applicability of the RANS approach in a model of the oropharynx, and concluded that the RANS k-ε model of turbulence is not accurate for particle tracking in the range of Reynolds numbers found in the oropharynx.
Lin et al.  computed flow in two imaging-based models (mouth-to-airway and intra-thoracic airway) using both Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) to resolve the turbulent structures. DNS directly solves the Navier-Stokes equations to resolve turbulent eddies at all scales, and has previously been applied to simulate low Reynolds number turbulent flows . LES solves the space-filtered Navier-Stokes equations to resolve only large-scale energy-containing turbulent eddies in turbulent flows, but parameterizes small-scale ones that are not resolved with the given grid size. The effect of unresolved eddies is accounted for by the Sub Grid-Scale (SGS) Reynolds stress, which is modeled by a SGS model . In Lin et al. , the comparison of DNS and LES results suggested that LES is sufficient to capture the turbulent laryngeal jet at a Reynolds number of 1,700 in the trachea. The presence of turbulence was found to increase the maximum localized wall shear stress three fold. To understand turbulent coherent spatial structures and their spectral signatures in a deterministic manner, the proper orthogonal decomposition technique  was applied to the velocity fields, showing that the turbulence was primarily attributable to Taylor-Görtler-like vortices. This study demonstrated that turbulence induced by the laryngeal jet could significantly affect airway flow patterns as well as tracheal wall shear stress, which may have pathophysiologic implications. To demonstrate the notion of subject-specific variability and the importance of individualized CFD analysis, Figure 3 shows the characteristics of the turbulent laryngeal jet in two human subjects at inspiration. Due to variations in the tracheal geometrical structure, the jet flow in the subject of Figure 3a impinges on the front side of the tracheal wall, whereas that of Figure 3b  impinges on the rear side. As a result, the location of the maximum wall shear stress is subject-dependent. It is also noted that the jet speed in the subject of Figure 3a is much higher than that of Figure 3b because the former has a much narrower glottis than the latter.
To account for the multi-scale effects of pulmonary airflow and to predict realistic regional ventilation, a 1D–3D coupled CFD simulation has been demonstrated by Lin et al. . For this type of simulation - assuming use of a 3D CT-based airway tree model that contains about 6 generations of airways and 5 million tetrahedral elements - the computational time needed for a breathing cycle of 4.8 seconds is approximately 90 hours running on 128 parallel CPU cores.
Airway CFD is a multi-scale problem, computed in what is in reality a deforming body that interacts with the surrounding tissue as the lung expands and recoils during ventilation. Predicting realistic distribution of inspired air and particulates therefore requires consideration of airway-tissue coupling through simulation of fluid structure interaction (FSI)  or by linking to continuum models for the tissue stress and strain . This essentially provides an organ-to-tissue link. While this is an extremely important issue for realistic simulation of air flow, to date it has not been addressed in any published CFD modeling studies. To model airway deformation, Lin et al.  proposed two approaches. The first is achieved by mapping CT images at various lung volumes via numerical optimization, which may be most suitable for the study of structure-function relationships at the organ scale; the second approach is based upon an FSI solver, which may be most appropriate for study of interplay between gas flow and lung tissue at the local scale. The 1D–3D modeling approach proposed by Lin et al.  holds promise for efficient coupling with other phenomena in the lung. For example, coupling with models of lung soft tissue mechanics, or distributed models of cell function. These types of coupled model lend themselves to studying how airflow distribution is altered by changes to the regional tissue properties in the lung: an important analysis for understanding delivery of aerosolized medicine to subjects with different pathologies. As noted by Lin et al. , maximum localized airway wall shear stress is increased three-fold by turbulence. This may have pathophysiologic implications for reactive inflammation and fibrosis. Linking down another scale to – interactions with the epithelial cell – is therefore an essential next step. Warren et al. [85, 86] have developed a model of airway epithelial cell regulation of the periciliary liquid. The epithelial cell is known to respond to shear, so interaction between the organ and tissue level models (airway CFD, computational tissue mechanics, and FSI) and the model of the epithelial cell (a system of ordinary differential equations) will provide a whole lung prediction of how regional alteration to ventilation distribution impacts on cellular regulation of the surface liquid, and therefore determines regional clearance of mucus.
A further example of a Physiomic model of the pulmonary system that was partly developed from the Human Lung Atlas and that spans organ to cell is the work by Burrowes et al. [2, 87–89] studying perfusion of the lung. In this series of papers blood is modeled as a two-phase fluid in a morphometrically based segmented (tube) model of the alveolo-capillary network, and as a Newtonian fluid in the arteries and veins. Flow in the elastic large vessels is estimated using 1D Navier-Stokes equations, and includes a gravitational term. The transit time for blood cells through the microcirculation is calculated for different regional pressure boundary conditions. A current limitation of the work is that the microcirculatory and arterial/venous models remain as separate entities; coupling these systems is necessary for studying regional gas exchange .
The Human Lung Atlas and the Lung Physiome are two major synergistic projects that together and individually are providing new insights into lung structure and function. The Atlas provides a comprehensive database of subjects with different pathologies from which computational models can be developed to predict disease-specific function and to interrogate hypotheses such as the principles that govern the regional distribution of ventilation and perfusion . Anatomical customized models of different individuals or different species can be used to probe the sensitivity of predictions to differences in geometry, whether subtle (normal temporal change in airway caliber of an individual) or major (between species). An important feature of an anatomically-based model is that it should be able to be reduced to a more simplified form, to facilitate understanding of the more complex model via analytic solutions or established mathematical analysis . The anatomical models reviewed here are not the ‘best’ tool for all computational studies, but they are necessary for airway CFD and other situations where heterogeneity in the tree structure , or spatially-dependent interaction with the surrounding tissue , is important.
This work was supported in part by NIH Grants R01-HL-064368 and R01-EB-005823.
Merryn H Tawhai, The University of Auckland.
Eric A Hoffman, The University of Iowa.
Ching-Long Lin, The University of Iowa.