Computational and Mathematical Methods in Medicine

Volume 2013 (2013), Article ID 679308, 12 pages

http://dx.doi.org/10.1155/2013/679308

## On the Stability of Lung Parenchymal Lesions with Applications to Early Pneumothorax Diagnosis

^{1}Thomas Jefferson High School for Science and Technology, Alexandria, VA 22312, USA^{2}Mathematical Sciences, George Mason University, Fairfax, VA 22030, USA

Received 1 December 2012; Revised 26 March 2013; Accepted 30 March 2013

Academic Editor: Thierry Busso

Copyright © 2013 Archis R. Bhandarkar et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Spontaneous pneumothorax, a prevalent medical challenge in most trauma cases, is a form of sudden lung collapse closely associated with risk factors such as lung cancer and emphysema. Our work seeks to explore and quantify the currently unknown pathological factors underlying lesion rupture in pneumothorax through biomechanical modeling. We hypothesized that lesion instability is closely associated with elastodynamic strain of the pleural membrane from pulsatile air flow and collagen-elastin dynamics. Based on the principles of continuum mechanics and fluid-structure interaction, our proposed model coupled isotropic tissue deformation with pressure from pulsatile air motion and the pleural fluid. Next, we derived mathematical instability criteria for our ordinary differential equation system and then translated these mathematical instabilities to physically relevant structural instabilities via the incorporation of a finite energy limiter. The introduction of novel biomechanical descriptions for collagen-elastin dynamics allowed us to demonstrate that changes in the protein structure can lead to a transition from stable to unstable domains in the material parameter space for a general lesion. This result allowed us to create a novel streamlined algorithm for detecting material instabilities in transient lung CT scan data via analyzing deformations in a local tissue boundary.

#### 1. Introduction

Spontaneous pneumothorax is a form of sudden lung collapse closely associated with risk factors such as lung cancer and emphysema. The disorder represents a pressing clinical challenge, affecting between 20 and 40% of patients with major trauma [1]. The onset of pneumothorax has been attributed to the rupture of quasi spherical parenchymal lesions or air blebs of the visceral pleural membrane [2]. The rupture of these lesions in turn releases air into the pleural cavity, and the consequent pressure buildup leads to the collapse of the lung [3]. Because pneumothorax has a relatively large recurrence rate of approximately 54% in the first four years after surgery, wholly understanding and quantifying the mechanisms behind lesion rupture are important [4].

The etiology and trigger mechanisms behind lesion rupture, however, still remain elusive [2, 5]. We hypothesize that lesion rupture and instability are closely associated with elastodynamic strain of the visceral pleural membrane from pulsatile air flow and changes in the constitutive protein composition of tissue. In order to assess the validity of this hypothesis, we constructed a biomechanical model based on the principles of finite strain continuum mechanics and fluid-structure interaction. Our proposed model closely aligns with the behavior of true biological soft tissue through the coupling of membrane energy limitations and growth and remodeling driven by collagen-elastin dynamics. Through the use of an exhaustive ordinary differential equation stability analysis, we isolated several instability regions in the material parameter space of a general lesion.

Based on these mathematical results, we developed an algorithm to rapidly assay clinical lung CT scan data for structural instabilities. The algorithm is based on cataloging and processing local deformations in a tissue boundary isolated via Sobel Edge Detection. The local deformation of the boundary is then fit to our biomechanical model and registered as either stable or unstable based on whether or not the interpolated material parameters lie in the instability region. Figure 1 shows a general outline of our novel early diagnostic method. In the future, the algorithm may be real time coupled with either ultrasound or CT scan data to quantify patient risk of pneumothorax earlier and improve the current diagnostic benchmark of only 75% lesion sensibility [6]. Unlike competitive diagnostic procedures for pneumothorax such as infrared thermography [1], our proposed method relies on detecting pathological hallmarks in tissue behavior before acute lung collapse actually occurs.

#### 2. Methodology

Our study focused on the construction of a theoretically based biomechanical model for lung parenchymal lesion rupture and then the implementation of this model towards a streamlined algorithm for early pneumothorax diagnosis. We began with the theoretical development of a model for simple isotropic tissue deformation and then conducted several initial numerical investigations. Next, we derived mathematical instabilities for our system and then translated these mathematical instabilities to physically relevant structural instabilities via the incorporation of an energy limiter. We then incorporated collagen-elastin dynamics (dynamic alteration of the constitutive proteins in the tissue) to demonstrate how changes in tissue structure may lead to lesion rupture. Lastly, we developed an algorithm for determining lesion stability from CT scan data by fitting our biomechanical model onto tissue deformation data.

##### 2.1. Mathematical Formulation

We modeled the parenchymal lesions with a spherical membrane geometry with thickness significantly smaller than the radius. Several parameters were needed in order to fully characterize the lesion geometry. The nondimensional stretch ratio was defined as , where represents the deformed radius varying with time and is the original undeformed radius. Assuming membrane incompressibility, the constant volume condition implies that the deformed thickness is , where represents the original undeformed thickness. Figure 2 depicts the geometry of these quasi-spherical lesions, including and .

###### 2.1.1. Constitutive Parenchymal Wall Model

Based on the findings of Tai and Lee [7], who concluded that lung tissue exhibits less than 10% anisotropy of the mean deformation, we pursued an isotropic model for the lung parenchymal wall. We had based a collagen-only isotropic pseudostrain-energy function for the wall of the parenchymal lesion on the work of Denny and Schroter [8]. Previous works by Denny and Schroter [9, 10] had also considered collagen and elastin using separate models [9], which were enhanced with viscoelastic contributions [10]. We adapted the Denny-Schroter constitutive relationship for our 1D formulation through multiplication by the deformed thickness :
where and are constants with dimensions N·m^{−2}, is dimensionless, and (where is the Green-Lagrangian strain tensor). Note that for simplicity, our model does not account for large stiffnesses for large strains. However, this is something that we hope to incorporate in future work. In order to facilitate the stability analysis and further calculations, this pseudostrain-energy function was approximated by a Taylor series centered about , yielding
If the 2D deformation gradient is diag, then the Cauchy stress resultant tensor for inner membrane stress can be represented by the following equation [11]:
Thus, the isotropic stress resultant () is given by the following:

###### 2.1.2. Model for Breathing Motion

We used a Fourier series developed by Jakuszkin [12] to model the pulsatile nature of breathing. The driven flow value of air through the lungs is Table 1 contains the constants , , and [12].

Note that while we have employed a sinusoidal model for simplicity, we can extend the model to account for erratic breathing behavior, which is not considered here. The driven flow value is related to the pulsatile pressure through the relation , where was the respiratory resistance of the bronchiole network. Because is 1.0 mmHg·s·L^{−1} [13] and the original pressure was atmospheric pressure (760 mmHg), the final pressure series was as follows:

###### 2.1.3. Dynamic Equation System

We obtained the dynamic equation for the wall through linear momentum balance to yield where is the density of the membrane and is the inner pressure minus the outer pressure (transmural pressure). The Navier-Stokes equations in spherical coordinates were used to model the pressure that the pleural fluid exerted on the lung parenchyma. The resultant pressure was as follows [14]: where is a constant value representing the pressure at a significantly large distance from the lung and is the density of the fluid. Note the absence of a viscous term as well as the same convention for the pressure at infinity as in the published work by Shah and Humphrey [14]. Combining the equations for the dynamic wall, the outer pressure, the dynamic radial stress which contributes the factor of , and viscoelasticity, we arrived at the following differential equation for the air-tissue-pleural fluid system: where and are the dynamic viscosities of the pleural fluid and membrane, respectively.

###### 2.1.4. Nondimensionalization

Nondimensionalization is a process that helps to accommodate multiscale variations in the data. This technique is a common engineering practice which essentially simplifies multiphysics problems with different measured units involved. Because our final governing equations couple three different physical aspects (the air, the membrane, and the pleural fluid), this method was necessitated. Moreover, nondimensionalization recovers the characteristic properties of the stretch ratio and stretch rate. We followed this standard which is commonly used in published works on biomechanics such as Shah and Humphrey [14]. The following nondimensionalized quantities were required: The resulting nondimensionalized equation system was as follows:

###### 2.1.5. Parameter Table

Table 2 is a complete collection of the relevant parameter values for the model.

##### 2.2. Initial Numerical Investigations

We conducted several numerical investigations in order to quantify both system dynamics and stability. The previous dynamical system was numerically evolved in MATLAB through a fourth-order explicit Runge-Kutta scheme. Plots of the stretch ratio versus time reveal sustained stable oscillations of the membrane for normal lung parenchymal tissue parameters of the Denny-Schroter model: , , and . Variations in material stiffness proved an effective tool in analyzing the effects of parameter variation on system dynamics. Figure 3 shows that as material stiffness increased in the model (for test values of ), the amplitude of the oscillations decreased as one would expect intuitively. Figure 4 shows that for the general range of , amplitude can be seen to monotonically decrease and begin to level off near zero. For the same window of stiffness constants, oscillation frequency reaches a peak near and then decreases monotonically from that value.

Figure 5 demonstrates that the system is stable even with minor perturbations in initial conditions for that class of material parameters. This suggests that there may be conditions where the bleb will not rupture. Figure 6 is explained in the preceding paragraph as a relative force proportion graph, which illustrates the dynamic behavior of the different forces including internal pressure, fluid structure, internal membrane, radial stress, and viscoelasticity. The graph suggests that friction is not a significant factor leading to pneumothorax, and this aligns with other published works. The numerical phase plane analysis in Figure 5 of versus further suggests the existence of an orbit about a stable equilibria. The spiral nature of the superimposed vector field, representing the stretch rate versus the derivative of the stretch rate , suggests an asymptotically stable node at for normal lung parenchymal tissue parameters. A relative force proportion graph in Figure 6 reveals the forces of internal pressure and inner membrane stress as the major driving forces behind this oscillation, whereas the viscoelastic and fluid forces effects are negligible. Furthermore, several unstable regions had been isolated through variation of the material parameters such as for , , and , but a rigorous mathematical analysis had been called for to fully determine whether the instabilities were truly physical or just due to numerical effects.

##### 2.3. Stability Results

###### 2.3.1. Determining System Equilibria

In order to simplify the following equilibria identification and stability analysis, we made the system autonomous by approximating . The system is at an equilibrium for any point such that and . This implies that for any equilibrium point. The equilibrium stretch ratio may be found by solving as follows: where , , and . For physically meaningful values of , the positive equilibrium point is the one of significance. For the normal lung parenchymal tissue parameters, as suggested by the initial numerical investigation.

###### 2.3.2. Characterization of Unstable Nodes

We conducted a general ordinary differential equation (ODE) stability analysis on the physical system in order to evaluate a lung parenchymal lesion with arbitrary material parameters for physical instabilities. Given a system with Jacobian matrix , ODE instability is prescribed for . For the previous dynamical system in and , the Jacobian matrix about the positive equilibrium point reduces to

The constraint of physically relevant material parameters imposes the inequalities , , and , whereas the constraint of defined real elements of the Jacobian matrix imposes , , and . Further note that from the definition of , is always greater than zero for physically relevant material constants .

Theorem 1. *For any values of and such that , , and , the system is unstable because .*

The expression is less than zero for all such that . The sign of this expression is in turn governed by the following three subexpressions. The subexpression of is consciously excluded from the sign analysis because it is invariantly positive.

Lemma 2. *The sign of and cannot both be positive for physically relevant parameters. **The lemma is proven by contradiction. By definition, implies . Isolation of reduces the inequality to . Likewise, reduces to and places bounds on . However, the inequality reduces to , which cannot be true from the restrictions of and . *

Lemma 3. *The expression
**
is greater than zero for any value of such that either and or and . Similarly, the expression is less than zero for any value of such that either and or and .**The inequality reduces to . Because both sides of the inequality are positive, the inequality may be squared and the terms brought to one side to give . Because only values of , and such that are considered, the expression is true . Thus, for either and or and . Likewise, for and or and . *

Lemma 4. *For physically relevant restriction of , the expression
**
is greater than zero for any value of such that any one of the following three conditions are met: and ; and ; , , and .**For either the conditions of and or and , the expression is always greater than zero. Hence, because is invariantly positive, the entire expression (15) will be greater than zero for these conditions. Given such that and , the inequality reduces to or . This in turn reduces to or .*

Lemma 5. *For physically relevant restriction of , the expression (15) is less than zero for any value of such that , , and . **Either the conditions of and or and make the statement unconditionally false, leaving only the case of , . The fact that is proven likewise as in Lemma 4.*

Lemma 6. *The expression
**
is greater than zero for all and less than zero for all . By definition, and , and thus the result follows. **The expression is true such that the product of expressions (14), (15), and (16) is positive. Upon coupling of the previous inequality systems, the only scenario that does not lead to contradictions is for (14), (15), and (16) all positive. This case implies that for , , and the system is unstable because . *

Theorem 7. *For any values of and such that , , , and , the system is unstable because .*

The expression is less than zero for all such that . The sign of this expression is in turn governed by subexpressions (15) and (17).

Lemma 8. *For physically relevant parameters (),
**
is less than zero for all or for , , and .**For all , expression (17) is invariantly negative, and thus, the statements are unconditionally true. This leaves only the case of and . Manipulation of the expression yields
**
or . This in turn yields , which upon isolation of gives the expression .*

The expression is true such that the product of expressions (15) and (17) is positive. For the signs of (15) and (17) to both be negative, Lemma 5 and Lemma 8 give , , , and . The case of expression (15) and (17) both greater than zero is consciously omitted from the proof because the analysis leads to an instability region already enclosed by Theorem 1. Figure 7 shows 2D portraits of the regions obtained from these instability criteria for a certain selection of material parameters , , and .

##### 2.4. Introducing Energy Limitations of Tissue

In order to accurately capture the biomechanics of lung parenchymal failure, we translated ODE instabilities to physical mechanical failure via the introduction of a finite energy limiter . Similar to Volokh and Vorp [19], we built a finite energy constitutive model based on the framework of the hyperelastic model which was constructed as follows: Note that as , a finite value of energy before rupture in contrast to the hyperelastic model where may grow unbounded as . Similarly, in the limiting case , as can be seen via Taylor expansion. Reformulating the Cauchy stress resultant tensor given by [14] with the new strain-energy function yields Thus, the new isotropic stress resultant () is given by the following:

As demonstrated in Figure 8, numerical simulations with the new isotropic stress resultant demonstrate rupture midexpansion as opposed to expansion to infinity in the original constitutive model.

##### 2.5. Collagen-Elastin Dynamics

The final component of the proposed model for parenchymal lesion growth and rupture is a means by which stable mechanical configurations may enter unstable domains. The introduction of collagen-elastin dynamics (a dynamic protein framework of the tissue) elucidates mechanisms for how a stable lesion may eventually rupture. We introduced a new constitutive model through the creation of a linear combination of Denny and Schroter’s [8] collagen-only strain energy function and Humphrey and Yin’s [20] elastin-only strain energy function:
where and are dimensionless quantities representing active collagen and elastin numbers. Similarly as before, we introduced a finite energy limiter to this constitutive model, and we recalculated the new isotropic strain tensor. Unlike the previous developments of the paper, collagen-elastin dynamics called for a return to the initial momentum balance due to the changing mass of the protein matrix. Given the mass of collagen and elastin defined by and (where and are equal to 5.10 kg·m^{−3} [21]), the new momentum balance was as follows:
where is the velocity at the boundary of an expanding biological membrane given by . With the nondimensionalized term , the first two differential equations of the dynamical system became
Based on histochemical evidence that had determined that collagen levels increase with increased stress and that elastin levels decrease because of elastolytic processes near rupture [22], we proposed the following differential equations for collagen-elastin dynamics:
Because and are positive constants, we devised the differential equations such that little change occurs in the protein matrix in scenarios far from rupture and larger changes occur as the membrane nears rupture. We numerically solved the series of four differential equations as a whole in MATLAB, revealing escapes from stable orbits for certain rate kinetics of the elastin-collagen dynamics that are depicted in Figure 9.

##### 2.6. Parameter Estimation from *In Vivo* Data

In order to build a database for parameter estimation on the proposed model, we collected *in vivo* lung data in the form of CT scans from the VIA/I-ELCAP Public Lung Image Database, which were about 250 mm × 150 mm in size. The Sobel Edge Detection algorithm provided a means to quantify the deformation of the lung boundary and to extract material parameters from *in vivo* clinical data for the strain-energy function. This algorithm works by estimating the horizontal and vertical spatial gradients at the interior pixels, and , through application of a matrix mask in order to detect edges. Pixels that meet a certain threshold for are counted as edges (Aybar, “Sobel Edge Detection”). This method was preferred over methods such as the Canny algorithm due to the coupling of speed and required accuracy for this problem. In order to facilitate the detection of the edges, we recolored the lung images to a 50% black and white scheme with 50% softening. After recoloration and softening, we ran the images through the Sobel algorithm implemented in Python, which returned a list of edge pixels for each image as coordinate pairs. We first transformed these pairs into a new coordinate system with the center of the image as the origin then categorized the pairs into subsections corresponding to sectors defined by angles. We then calculated values for the stretch factor, kinetic energy per square meter , and radial velocity in each radial direction.

Next, we applied nonlinear regression using the Levenberg-Marquadt algorithm on each subsection in order to extract the three material parameters , , and as well as the corresponding coefficients of determination ( values). We fit the velocity and kinetic energy data to the Denny-Schroter second-order Taylor series approximation. After we determined these parameters, they were classified as either stable or unstable corresponding to the previously derived stability criteria. Lastly, our program recolored the boundary of the lung parenchyma based on the stability of the corresponding subsections, with green corresponding to stable subsections and red corresponding to unstable subsections. The general schema for our computational algorithm is shown in Figure 10. Thus, our algorithm is streamlined towards interpolating material parameters from deformations in lung CT scan data and identifying specific areas of risk in the lung parenchyma.

Motion artifacts, while important and are often observed in CT data acquisition, are not considered in this paper. The images are represented using intensity values at various pixel locations. This paper used only one dataset as a proof of concept. Our algorithm looks at small sections of the lung boundary, allowing for the isotropic assumption to be valid. However, it must be noted that this assumption must be removed for more general studies.

#### 3. Conclusion and Future Directions

Our proposed biomechanical model for lung parenchymal lesions demonstrates that elastodynamic strain of the visceral pleural membrane from pulsatile air flow and changes in the constitutive protein composition of tissue are key pathological factors in pneumothorax. Through proving general instability results, our study demonstrated that these factors biject to physically unstable domains in the material parameter space for a general lesion. Specifically, our discovery that certain rate kinetics of collagen-elastin dynamics may lead to unstable mechanical configurations aligns well with clinical studies, which have found that certain connective tissue diseases lead to spontaneous pneumothorax [23, 24]. Future work on the theoretical biomechanics aspect of our model will incorporate COMSOL multiphysics simulation software integrated with finite element mesh processing to solve the problems of nonaxisymmetry and anisotropic lesion growth. Also, it would be interesting to study how variations in the dynamic viscosities affect the behavior of parenchymal lesions. Clinical studies have indicated that gravity is a significant factor on the air accumulated in the pleural cavity. We hope to consider these effects in future work. Quantification of material properties will require efficient inverse parameter estimation studies, which we hope to incorporate in the future. Such studies will help us to develop efficient simulation tools [25, 26]. Other methods for lung deformation registration such as the dictionary learning approach used by Zhang et al. [27] will be considered.

Our algorithm for determining lesion stability from CT scan data relies on fitting our biomechanical model to deformations in a local tissue boundary. Using Sobel Edge Detection and CT scans from a lung image database, we could isolate unstable regions of local lung parenchymal tissue. The algorithm is a robust method for quantifying patient risk of pneumothorax and may be streamlined to work in the future with in vivo lung data collected from either ultrasound or CT scans. Overall, our computational solution for lung parenchymal lesion detection and patient-specific structural instability profiling is a feasible alternative diagnostic strategy and has the potential to surpass current diagnostic benchmarks for pneumothorax.

#### References

- P. B. Rich, G. R. Dulabon, C. D. Douillet et al., “Infrared thermography: a rapid, portable, and accurate technique to detect experimental pneumothorax,”
*Journal of Surgical Research*, vol. 120, no. 2, pp. 163–170, 2004. View at Publisher · View at Google Scholar · View at Scopus - W. C. Yu, Y. C. Yeung, Y. Chang et al., “Use of endobronchial one-way valves reveals questions on etiology of spontaneous pneumothorax: report of three cases,”
*Journal of Cardiothoracic Surgery*, vol. 4, no. 1, article 1749, p. 63, 2009. View at Publisher · View at Google Scholar · View at Scopus - G. P. Currie, R. Alluri, G. L. Christie, and J. S. Legge, “Pneumothorax: an update,”
*Postgraduate Medical Journal*, vol. 83, no. 981, pp. 461–465, 2007. View at Publisher · View at Google Scholar · View at Scopus - C. A. Laituri, P. A. Valusek, D. C. Rivard et al., “The utility of computed tomography in the management of patients with spontaneous pneumothorax,”
*Journal of Pediatric Surgery*, vol. 46, no. 8, pp. 1523–1525, 2011. View at Google Scholar - C. G. Ball, S. M. Hameed, D. Evans, J. B. Kortbeek, and A. W. Kirkpatrick, “Occult pneumothorax in the mechanically ventilated trauma patient,”
*Canadian Journal of Surgery*, vol. 46, no. 5, pp. 373–379, 2003. View at Google Scholar · View at Scopus - N. Nathan, J. Guilbert, M. Larroquet, M. Lenoir, A. Clement, and R. Epaud, “Efficacy of blebs detection for preventive surgery in children's idiopathic spontaneous pneumothorax,”
*World Journal of Surgery*, vol. 34, no. 1, pp. 185–189, 2010. View at Publisher · View at Google Scholar · View at Scopus - R. C. Tai and G. C. Lee, “Isotropy and homogeneity of lung tissue deformation,”
*Journal of Biomechanics*, vol. 14, no. 4, pp. 243–252, 1981. View at Google Scholar · View at Scopus - E. Denny and R. C. Schroter, “A model of non-uniform lung parenchyma distortion,”
*Journal of Biomechanics*, vol. 39, no. 4, pp. 652–663, 2006. View at Publisher · View at Google Scholar · View at Scopus - E. Denny and R. C. Schroter, “Relationships between alveolar size and fibre distribution in a mammalian lung alveolar duct model,”
*ASME Journal of Biomechanical Engineering*, vol. 119, no. 3, pp. 289–297, 1997. View at Google Scholar · View at Scopus - E. Denny and R. C. Schroter, “Viscoelastic behavior a lung alveolar duct model,”
*ASME Journal of Biomechanical Engineering*, vol. 122, no. 2, pp. 143–151, 2000. View at Publisher · View at Google Scholar · View at Scopus - G. David and J. D. Humphrey, “Further evidence for the dynamic stability of intracranial saccular aneurysms,”
*Journal of Biomechanics*, vol. 36, no. 8, pp. 1143–1150, 2003. View at Publisher · View at Google Scholar · View at Scopus - K. Jakuszkin, “Analysis of properties of a lung mechanical model during artificial ventilation using measurement station,”
*Metrology and Measurement Systems*, vol. 17, no. 3, pp. 427–438, 2010. View at Google Scholar · View at Scopus - A. Ben-Tal, “Simplified models for gas exchange in the human lungs,”
*Journal of Theoretical Biology*, vol. 238, no. 2, pp. 474–495, 2006. View at Publisher · View at Google Scholar · View at Scopus - A. D. Shah and J. D. Humphrey, “Finite strain elastodynamics of intracranial saccular aneurysms,”
*Journal of Biomechanics*, vol. 32, no. 6, pp. 593–599, 1999. View at Publisher · View at Google Scholar · View at Scopus - K. Amjadi, G. G. Alvarez, E. Vanderhelst, B. Velkeniers, M. Lam, and M. Noppen, “The prevalence of blebs or bullae among young healthy adults: a thoracoscopic investigation,”
*Chest*, vol. 132, no. 4, pp. 1140–1145, 2007. View at Publisher · View at Google Scholar · View at Scopus - O. Yetkin, I. Tek, F. Yetkin, and N. Numanoglu, “Role of pleural viscosity in the differential diagnosis of exudative pleural effusion,”
*Respirology*, vol. 12, no. 2, pp. 267–271, 2007. View at Publisher · View at Google Scholar · View at Scopus - J. Rubins and H. Manning, “Pleural effusion workup,” 2012, http://emedicine.medscape.com/article/299959-workup.
- S. Girnyk, A. Barannik, E. Barannik, V. Tovstiak, A. Marusenko, and V. Volokhov, “The estimation of elasticity and viscosity of soft tissues in vitro using the data of remote acoustic palpation,”
*Ultrasound in Medicine and Biology*, vol. 32, no. 2, pp. 211–219, 2006. View at Publisher · View at Google Scholar · View at Scopus - K. Y. Volokh and D. A. Vorp, “A model of growth and rupture of abdominal aortic aneurysm,”
*Journal of Biomechanics*, vol. 41, no. 5, pp. 1015–1021, 2008. View at Publisher · View at Google Scholar · View at Scopus - J. D. Humphrey and F. C. Yin, “A new constitutive formulation for characterizing the mechanical behavior of soft tissues,”
*Biophysical Journal*, vol. 52, no. 4, pp. 563–570, 1987. View at Google Scholar · View at Scopus - Y. Takubo, T. Hirai, S. Muro, K. Kogishi, M. Hosokawa, and M. Mishima, “Age-associated changes in elastin and collagen content and the proportion of types I and III collagen in the lungs of mice,”
*Experimental Gerontology*, vol. 34, no. 3, pp. 353–364, 1999. View at Publisher · View at Google Scholar · View at Scopus - D. S. Faffe and W. A. Zin, “Lung parenchymal mechanics in health and disease,”
*Physiological Reviews*, vol. 89, no. 3, pp. 759–775, 2009. View at Publisher · View at Google Scholar · View at Scopus - G. R. Junior, A. S. Rodrigues, and L. C. Meguins, “Thoracoscopic surgery for simultaneous bilateral spontaneous pneumothorax,”
*The Internet Journal of Thoracic and Cardiovascular Surgery*, vol. 12, no. 1, 2008. View at Google Scholar - H. J. Metivier, R. Masse, F. Vai, and R. G. Pariente, “Increased content of neutrosoluble collagen and dialyzable hydroxyproline in panacinar emphysema and spontaneous pneumothorax,”
*Experimental Lung Research*, vol. 3, no. 2, pp. 129–136, 1982. View at Google Scholar · View at Scopus - A. Santhanam, C. Fidopiastis, K. Langen, S. Meeks, P. Davenport, and J. P. Rolland, “Real-time visualization of subject specific lung dynamics,”
*IEEE Computer-based Medical Systems*, 2006. View at Google Scholar - W. P. Segars, T. S. Lee, and B. M. W. Tsui, “Simulation of motion defects in the 4D NCAT cardiac model,”
*Journal of Nuclear Medicine*, vol. 44, p. 142P, 2003. View at Google Scholar - S. Zhang, Y. Zhan, and D. N. Metaxas, “Deformable segmentation via sparse representation and dictionary learning,”
*Medical Image Analysis*, vol. 16, no. 7, pp. 1385–1396, 2012. View at Google Scholar