MATLAB-Based Algorithm and Software for Analysis of Wavy Collagen Fibers.

Knowledge of soft tissue fiber structure is necessary for accurate characterization and modeling of their mechanical response. Fiber configuration and structure informs both our understanding of healthy tissue physiology and of pathological processes resulting from diseased states. This study develops an automatic algorithm to simultaneously estimate fiber global orientation, abundance, and waviness in an investigated image. To our best knowledge, this is the first validated algorithm which can reliably separate fiber waviness from its global orientation for considerably wavy fibers. This is much needed feature for biological tissue characterization. The algorithm is based on incremental movement of local regions of interest (ROI) and analyzes two-dimensional images. Pixels belonging to the fiber are identified in the ROI, and ROI movement is determined according to local orientation of fiber within the ROI. The algorithm is validated with artificial images and ten images of porcine trachea containing wavy fibers. In each image, 80-120 fibers were tracked manually to serve as verification. The coefficient of determination R2 between curve lengths and histograms documenting the fiber waviness and global orientation were used as metrics for analysis. Verification-confirmed results were independent of image rotation and degree of fiber waviness, with curve length accuracy demonstrated to be below 1% of fiber curved length. Validation-confirmed median and interquartile range of R2, respectively, were 0.90 and 0.05 for curved length, 0.92 and 0.07 for waviness, and 0.96 and 0.04 for global orientation histograms. Software constructed from the proposed algorithm was able to track one fiber in about 1.1 s using a typical office computer. The proposed algorithm can reliably and accurately estimate fiber waviness, curve length, and global orientation simultaneously, moving beyond the limitations of prior methods.

Pressure-volume mechanics of inflating and deflating intact whole organ porcine lungs.

Pressure-volume curves of the lung are classical measurements of lung function and are impacted by changes in lung structure due to disease or shifts in air-delivery volume or cycling rate. Diseased and preterm infant lungs have been found to show heterogeneous behavior which is highly frequency dependent. This breathing rate dependency has motivated the exploration of multi-frequency oscillatory ventilators to deliver volume oscillation with optimal frequencies for various portions of the lung to provide more uniform air distribution. The design of these advanced ventilators requires the examination of lung function and mechanics, and an improved understanding of the pressure-volume response of the lung. Therefore, to comprehensively analyze whole lung organ mechanics, we investigate six combinations of varying applied volumes and frequencies using ex-vivo porcine specimens and our custom-designed electromechanical breathing apparatus. Lung responses were evaluated through measurements of inflation and deflation slopes, static compliance, peak pressure and volume, as well as hysteresis, energy loss, and pressure relaxation. Generally, we observed that the lungs were stiffer when subjected to faster breathing rates and lower inflation volumes. The lungs exhibited greater inflation volume dependencies compared to frequency dependencies. This study's reported response of the lung to variations of inflation volume and breathing rate can help the optimization of conventional mechanical ventilators and inform the design of advanced ventilators. Although frequency dependency is found to be minimal in normal porcine lungs, this preliminary study lays a foundation for comparison with pathological lungs, which are known to demonstrate marked rate dependency.

Effects of tissue degradation by collagenase and elastase on the biaxial mechanics of porcine airways.

BACKGROUND: Common respiratory illnesses, such as emphysema and chronic obstructive pulmonary disease, are characterized by connective tissue damage and remodeling. Two major fibers govern the mechanics of airway tissue: elastin enables stretch and permits airway recoil, while collagen prevents overextension with stiffer properties. Collagenase and elastase degradation treatments are common avenues for contrasting the role of collagen and elastin in healthy and diseased states; while previous lung studies of collagen and elastin have analyzed parenchymal strips in animal and human specimens, none have focused on the airways to date. METHODS: Specimens were extracted from the proximal and distal airways, namely the trachea, large bronchi, and small bronchi to facilitate evaluations of material heterogeneity, and subjected to biaxial planar loading in the circumferential and axial directions to assess airway anisotropy. Next, samples were subjected to collagenase and elastase enzymatic treatment and tensile tests were repeated. Airway tissue mechanical properties pre- and post-treatment were comprehensively characterized via measures of initial and ultimate moduli, strain transitions, maximum stress, hysteresis, energy loss, and viscoelasticity to gain insights regarding the specialized role of individual connective tissue fibers and network interactions. RESULTS: Enzymatic treatment demonstrated an increase in airway tissue compliance throughout loading and resulted in at least a 50% decrease in maximum stress overall. Strain transition values led to significant anisotropic manifestation post-treatment, where circumferential tissues transitioned at higher strains compared to axial counterparts. Hysteresis values and energy loss decreased after enzymatic treatment, where hysteresis reduced by almost half of the untreated value. Anisotropic ratios exhibited axially led stiffness at low strains which transitioned to circumferentially led stiffness when subjected to higher strains. Viscoelastic stress relaxation was found to be greater in the circumferential direction for bronchial airway regions compared to axial counterparts. CONCLUSION: Targeted fiber treatment resulted in mechanical alterations across the loading range and interactions between elastin and collagen connective tissue networks was observed. Providing novel mechanical characterization of elastase and collagenase treated airways aids our understanding of individual and interconnected fiber roles, ultimately helping to establish a foundation for constructing constitutive models to represent various states and progressions of pulmonary disease.

Biaxial mechanical properties of the bronchial tree: Characterization of elasticity, extensibility, and energetics, including the effect of strain rate and preconditioning.

Distal airways commonly obstruct in lung disease and despite their importance, their mechanical properties are vastly underexplored. The lack of bronchial experiments restricts current airway models to either assume rigid structures, or extrapolate the material properties of the trachea to represent the small airways. Furthermore, past works are exclusively limited to uniaxial testing; investigating the multidirectional tensile loads of both the proximal and distal pulmonary airways is long overdue. Here we present comprehensive mechanical and viscoelastic properties of the porcine airway tree, including the trachea, trachealis muscle, large bronchi, and small bronchi, via measures of elasticity, extensibility, and energetics to explore regional and directional dependencies, cross-examining strain rate and preconditioning effects using planar equibiaxial tensile tests for the first time. We find bronchial regions are notably heterogeneous, where the trachea exhibits greater stiffness, energy loss, and preconditioning sensitivity than the smaller airways. Interestingly, the trachealis muscle is similar to the distal bronchi, despite being anatomically located adjacent to the proximal ring. Tissues are anisotropic and axially stiffer under initial loading, losing more energy with greater stress relaxation circumferentially. Strain rate dependency is also noted, where tissues are more energetically efficient at the faster strain rate, likely attributable to the microstructure. Findings highlight assumptions of homogeneity and isotropy are inadequate, and enable the improvement of aerosol flow and dynamic airway deformation computational predictive models. These results provide much needed fundamental material properties for future explorations contrasting healthy versus diseased pulmonary airway mechanics to better understand the relationship between structure and lung function. STATEMENT OF SIGNIFICANCE: We present comprehensive multiaxial mechanical tensile experiments of the proximal and distal airways via measures of maximum stress, initial and ultimate moduli, strain and stress transitions, hysteresis, energy loss, and stress relaxation, and further assess preconditioning and strain rate dependencies to examine the relationship between lung function and structure. The mechanical response of the bronchial tree demonstrates significant anisotropy and heterogeneity, even within the tracheal ring, and emphasizes that contrary to past studies, the behavior of the proximal airways cannot be extended to distal bronchial tree analyses. Establishing these material properties is critical to advancing our understanding of airway function and in developing accurate computational simulations to help diagnose and monitor pulmonary diseases.

Positive- and Negative-Pressure Ventilation Characterized by Local and Global Pulmonary Mechanics.

Rationale: There is continued debate regarding the equivalency of positive-pressure ventilation (PPV) and negative-pressure ventilation (NPV). Resolving this question is important because of the different practical ramifications of the two paradigms. Objectives: We sought to investigate the parallel between PPV and NPV and determine whether or not these two paradigms cause identical ventilation profiles by analyzing the local strain mechanics when the global tidal volume (Vt) and inflation pressure was matched. Methods: A custom-designed electromechanical apparatus was used to impose equal global loads and displacements on the same ex vivo healthy porcine lung using PPV and NPV. High-speed high-resolution cameras recorded local lung surface deformations and strains in real time, and differences between PPV and NPV global energetics, viscoelasticity, as well as local tissue distortion were assessed. Measurements and Main Results: During initial inflation, NPV exhibited significantly more bulk pressure-volume compliance than PPV, suggestive of earlier lung recruitment. NPV settings also showed reduced relaxation, hysteresis, and energy loss compared with PPV. Local strain trends were also decreased in NPV, with reduced tissue distortion trends compared with PPV, as revealed through analysis of tissue anisotropy. Conclusions: Apparently, contradictory previous studies are not mutually exclusive. Equivalent changes in transpulmonary pressures in PPV and NPV lead to the same changes in lung volume and pressures, yet local tissue strains differ between PPV and NPV. Although limited to healthy specimens and ex vivo experiments in the absence of a chest cavity, these results may explain previous reports of better oxygenation and less lung injury in NPV.

Associating local strains to global pressure-volume mouse lung mechanics using digital image correlation.

Pulmonary diseases alter lung mechanical properties, can cause loss of function, and necessitate use of mechanical ventilation, which can be detrimental. Investigations of lung tissue (local) scale mechanical properties are sparse compared to that of the whole organ (global) level, despite connections between regional strain injury and ventilation. We examine ex vivo mouse lung mechanics by investigating strain values, local compliance, tissue surface heterogeneity, and strain evolutionary behavior for various inflation rates and volumes. A custom electromechanical, pressure-volume ventilator is coupled with digital image correlation to measure regional lung strains and associate local to global mechanics by analyzing novel pressure-strain evolutionary measures. Mean strains at 5 breaths per minute (BPM) for applied volumes of 0.3, 0.5, and 0.7 ml are 5.0, 7.8, and 11.3%, respectively, and 4.7, 8.8, and 12.2% for 20 BPM. Similarly, maximum strains among all rate and volume combinations range 10.7%-22.4%. Strain values (mean, range, mode, and maximum) at peak inflation often exhibit significant volume dependencies. Additionally, select evolutionary behavior (e.g., local lung compliance quantification) and tissue heterogeneity show significant volume dependence. Rate dependencies are generally found to be insignificant; however, strain values and surface lobe heterogeneity tend to increase with increasing rates. By quantifying strain evolutionary behavior in relation to pressure-volume measures, we associate time-continuous local to global mouse lung mechanics for the first time and further examine the role of volume and rate dependency. The interplay of multiscale deformations evaluated in this work can offer insights for clinical applications, such as ventilator-induced lung injury.

Developing a Lung Model in the Age of COVID-19: A Digital Image Correlation and Inverse Finite Element Analysis Framework.

Pulmonary diseases, driven by pollution, industrial farming, vaping, and the infamous COVID-19 pandemic, lead morbidity and mortality rates worldwide. Computational biomechanical models can enhance predictive capabilities to understand fundamental lung physiology; however, such investigations are hindered by the lung's complex and hierarchical structure, and the lack of mechanical experiments linking the load-bearing organ-level response to local behaviors. In this study we address these impedances by introducing a novel reduced-order surface model of the lung, combining the response of the intricate bronchial network, parenchymal tissue, and visceral pleura. The inverse finite element analysis (IFEA) framework is developed using 3-D digital image correlation (DIC) from experimentally measured non-contact strains and displacements from an ex-vivo porcine lung specimen for the first time. A custom-designed inflation device is employed to uniquely correlate the multiscale classical pressure-volume bulk breathing measures to local-level deformation topologies and principal expansion directions. Optimal material parameters are found by minimizing the error between experimental and simulation-based lung surface displacement values, using both classes of gradient-based and gradient-free optimization algorithms and by developing an adjoint formulation for efficiency. The heterogeneous and anisotropic characteristics of pulmonary breathing are represented using various hyperelastic continuum formulations to divulge compound material parameters and evaluate the best performing model. While accounting for tissue anisotropy with fibers assumed along medial-lateral direction did not benefit model calibration, allowing for regional material heterogeneity enabled accurate reconstruction of lung deformations when compared to the homogeneous model. The proof-of-concept framework established here can be readily applied to investigate the impact of assorted organ-level ventilation strategies on local pulmonary force and strain distributions, and to further explore how diseased states may alter the load-bearing material behavior of the lung. In the age of a respiratory pandemic, advancing our understanding of lung biomechanics is more pressing than ever before.

Local and global growth and remodeling in calcific aortic valve disease and aging.

Aging and calcific aortic valve disease (CAVD) are the main factors leading to aortic stenosis. Both processes are accompanied by growth and remodeling pathways that play a crucial role in aortic valve pathophysiology. Herein, a computational growth and remodeling (G&R) framework was developed to investigate the effects of aging and calcification on aortic valve dynamics. Particularly, an algorithm was developed to couple the global growth and stiffening of the aortic valve due to aging and the local growth and stiffening due to calcification with the aortic valve transient dynamics. The aortic valve dynamics during baseline were validated with available data in the literature. Subsequently, the changes in aortic valve dynamic patterns during aging and CAVD progression were studied. The results revealed the patterns in geometric orifice area reduction and an increase in the valve stress during local and global growth and remodeling of the aortic valve. The proposed algorithm provides a framework to couple mechanobiology models of disease growth with tissue-scale transient structural mechanics models to study the biomechanical changes during cardiovascular disease growth and aging.

Novel Mechanical Strain Characterization of Ventilated ex vivo Porcine and Murine Lung using Digital Image Correlation.

Respiratory illnesses, such as bronchitis, emphysema, asthma, and COVID-19, substantially remodel lung tissue, deteriorate function, and culminate in a compromised breathing ability. Yet, the structural mechanics of the lung is significantly understudied. Classical pressure-volume air or saline inflation studies of the lung have attempted to characterize the organ's elasticity and compliance, measuring deviatory responses in diseased states; however, these investigations are exclusively limited to the bulk composite or global response of the entire lung and disregard local expansion and stretch phenomena within the lung lobes, overlooking potentially valuable physiological insights, as particularly related to mechanical ventilation. Here, we present a method to collect the first non-contact, full-field deformation measures of ex vivo porcine and murine lungs and interface with a pressure-volume ventilation system to investigate lung behavior in real time. We share preliminary observations of heterogeneous and anisotropic strain distributions of the parenchymal surface, associative pressure-volume-strain loading dependencies during continuous loading, and consider the influence of inflation rate and maximum volume. This study serves as a crucial basis for future works to comprehensively characterize the regional response of the lung across various species, link local strains to global lung mechanics, examine the effect of breathing frequencies and volumes, investigate deformation gradients and evolutionary behaviors during breathing, and contrast healthy and pathological states. Measurements collected in this framework ultimately aim to inform predictive computational models and enable the effective development of ventilators and early diagnostic strategies.

Introducing a Custom-Designed Volume-Pressure Machine for Novel Measurements of Whole Lung Organ Viscoelasticity and Direct Comparisons Between Positive- and Negative-Pressure Ventilation.

Asthma, emphysema, COVID-19 and other lung-impacting diseases cause the remodeling of tissue structural properties and can lead to changes in conducting pulmonary volume, viscoelasticity, and air flow distribution. Whole organ experimental inflation tests are commonly used to understand the impact of these modifications on lung mechanics. Here we introduce a novel, automated, custom-designed device for measuring the volume and pressure response of lungs, surpassing the capabilities of traditional machines and built to range size-scales to accommodate both murine and porcine tests. The software-controlled system is capable of constructing standardized continuous volume-pressure curves, while accounting for air compressibility, yielding consistent and reproducible measures while eliminating the need for pulmonary degassing. This device uses volume-control to enable viscoelastic whole lung macromechanical insights from rate dependencies and pressure-time curves. Moreover, the conceptual design of this device facilitates studies relating the phenomenon of diaphragm breathing and artificial ventilation induced by pushing air inside the lungs. System capabilities are demonstrated and validated via a comparative study between ex vivo murine lungs and elastic balloons, using various testing protocols. Volume-pressure curve comparisons with previous pressure-controlled systems yield good agreement, confirming accuracy. This work expands the capabilities of current lung experiments, improving scientific investigations of healthy and diseased pulmonary biomechanics. Ultimately, the methodologies demonstrated in the manufacturing of this system enable future studies centered on investigating viscoelasticity as a potential biomarker and improvements to patient ventilators based on direct assessment and comparisons of positive- and negative-pressure mechanics.

Characterizing the viscoelasticity of extra- and intra-parenchymal lung bronchi.

Pulmonary disease is known to cause remodeling of tissue structure, resulting in altered viscoelastic properties; yet the foundation for understanding this phenomenon is still nascent and will enable scientific insights regarding lung functionality. In order to characterize the viscoelastic response of pulmonary airways, uniaxial tensile experiments are conducted on porcine extra- and intra-parenchymal bronchial regions, measuring both axially and circumferentially oriented tissue. Anisotropic and heterogeneous effects on preconditioning and hysteresis are substantial, linking to energy dissipation expectancies. Stress relaxation is rheologically modeled using several classical configurations of discrete spring and dashpot elements; among them, Standard Linear Solid (SLS) and Maxwell-Weichart exhibit better fit performance. Enhanced fractional order derivative SLS (FSLS) model is also evaluated through use of a hybrid spring-pot of order α. FSLS outperforms the conventional models, demonstrating superior representation of the stress-relaxation curve's initial value and non-linear asymptotic decent. FSLS parameters exhibit notable orientation- and region-specific values, trending with observed tissue structural constituents, such as glycosaminoglycan and collagen. To the best of our knowledge, this work is the first to characterize proximal and distal bronchial energy efficiency and contextualize tissue biochemical composition in view of experimental measures and viscoelastic trends. Results provide a foundation for future investigations, particularly for understanding the role of viscoelasticity in diseased states.

Mechanics of pulmonary airways: Linking structure to function through constitutive modeling, biochemistry, and histology.

Breathing involves fluid-solid interactions in the lung; however, the lack of experimental data inhibits combining the mechanics of air flow to airway deformation, challenging the understanding of how biomaterial constituents contribute to tissue response. As such, lung mechanics research is increasingly focused on exploring the relationship between structure and function. To address these needs, we characterize mechanical properties of porcine airways using uniaxial tensile experiments, accounting for bronchial orientation- and location- dependency. Structurally-reinforced constitutive models are developed to incorporate the role of collagen and elastin fibers embedded within the extrafibrillar matrix. The strain-energy function combines a matrix description (evaluating six models: compressible NeoHookean, unconstrained Ogden, uncoupled Mooney-Rivlin, incompressible Ogden, incompressible Demiray and incompressible NeoHookean), superimposed with non-linear fibers (evaluating two models: exponential and polynomial). The best constitutive formulation representative of all bronchial regions is determined based on curve-fit results to experimental data, accounting for uniqueness and sensitivity. Glycosaminoglycan and collagen composition, alongside tissue architecture, indicate fiber form to be primarily responsible for observed airway anisotropy and heterogeneous mechanical behavior. To the authors' best knowledge, this study is the first to formulate a structurally-motivated constitutive model, augmented with biochemical analysis and microstructural observations, to investigate the mechanical function of proximal and distal bronchi. Our systematic pulmonary tissue characterization provides a necessary foundation for understanding pulmonary mechanics; furthermore, these results enable clinical translation through simulations of airway obstruction in disease, fluid-structure interaction insights during breathing, and potentially, predictive capabilities for medical interventions. STATEMENT OF SIGNIFICANCE: The advancement of pulmonary research relies on investigating the biomechanical response of the bronchial tree. Experiments demonstrating the non-linear, heterogeneous, and anisotropic material behavior of porcine airways are used to develop a structural constitutive model representative of proximal and distal bronchial behavior. Calibrated material parameters exhibit regional variation in biomaterial properties, initially hypothesized to originate from tissue constituents. Further exploration through biochemical and histological analysis indicates mechanical function is primarily governed by microstructural form. The results of this study can be directly used in finite element and fluid-structure interaction models to enable physiologically relevant and more accurate computational simulations aimed to help diagnose and monitor pulmonary disease.

Mechanical properties of the airway tree: heterogeneous and anisotropic pseudoelastic and viscoelastic tissue responses.

Airway obstruction and pulmonary mechanics remain understudied despite lung disease being the third cause of death in the United States. Lack of relevant data has led computational pulmonary models to infer mechanical properties from available material data for the trachea. Additionally, the time-dependent, viscoelastic behaviors of airways have been largely overlooked, despite their potential physiological relevance and utility as metrics of tissue remodeling and disease progression. Here, we address the clear need for airway-specific material characterization to inform biophysical studies of the bronchial tree. Specimens from three airway levels (trachea, large bronchi, and small bronchi) and two orientations (axial and circumferential) were prepared from five fresh pig lungs. Uniaxial tensile tests revealed substantial heterogeneity and anisotropy. Overall, the linear pseudoelastic modulus was significantly higher axially than circumferentially (30.5 ± 3.1 vs. 8.4 ± 1.1 kPa) and significantly higher among circumferential samples for small bronchi than for the trachea and large bronchi (12.5 ± 1.9 vs. 6.0 ± 0.6 and 6.6 ± 0.9 kPa). Circumferential samples exhibited greater percent stress relaxation over 300 s than their axial counterparts (38.0 ± 1.4 vs. 23.1 ± 1.5%). Axial and circumferential trachea samples displayed greater percent stress relaxation (26.4 ± 1.6 and 42.5 ± 1.7%) than corresponding large and small bronchi. This ex vivo pseudoelastic and viscoelastic characterization reveals novel anisotropic and heterogeneous behaviors and equips us to construct airway-specific constitutive relations. Our results establish necessary fundamentals for airway mechanics, laying the groundwork for future studies to extend to clinical questions surrounding lung injury, and further directly enables computational tools for lung disease obstruction predictions. NEW & NOTEWORTHY Understanding the mechanics of the lung is necessary for investigating disease progression. Trachea mechanics comprises the vast majority of ex vivo airway tissue characterization despite distal airways being the site of disease manifestation and occlusion. Furthermore, viscoelastic studies are scarce, whereas time-dependent behaviors could be potential physiological metrics of tissue remodeling. In this study, the critical need for airway-specific material properties is addressed, reporting bronchial tree anisotropic and heterogeneous material properties.

Elastosis during airway wall remodeling explains multiple co-existing instability patterns.

Living structures can undergo morphological changes in response to growth and alterations in microstructural properties in response to remodeling. From a biological perspective, airway wall inflammation and airway elastosis are classical hallmarks of growth and remodeling during chronic lung disease. From a mechanical point of view, growth and remodeling trigger mechanical instabilities that result in inward folding and airway obstruction. While previous analytical and computational studies have focused on identifying the critical parameters at the onset of folding, few have considered the post-buckling behavior. All prior studies assume constant microstructural properties during the folding process; yet, clinical studies now reveal progressive airway elastosis, the degeneration of elastic fibers associated with a gradual stiffening of the inner layer. Here, we explore the influence of temporally evolving material properties on the post-bifurcation behavior of the airway wall. We show that a growing and stiffening inner layer triggers an additional subsequent bifurcation after the first instability occurs. Evolving material stiffnesses provoke failure modes with multiple co-existing wavelengths, associated with the superposition of larger folds evolving on top of the initial smaller folds. This phenomenon is exclusive to material stiffening and conceptually different from the phenomenon of period doubling observed in constant-stiffness growth. Our study suggests that the clinically observed multiple wavelengths in diseased airways are a result of gradual airway wall stiffening. While our evolving material properties are inspired by the clinical phenomenon of airway elastosis, the underlying concept is broadly applicable to other types of remodeling including aneurysm formation or brain folding.

Systems biology and mechanics of growth.

In contrast to inert systems, living biological systems have the advantage to adapt to their environment through growth and evolution. This transfiguration is evident during embryonic development, when the predisposed need to grow allows form to follow function. Alterations in the equilibrium state of biological systems breed disease and mutation in response to environmental triggers. The need to characterize the growth of biological systems to better understand these phenomena has motivated the continuum theory of growth and stimulated the development of computational tools in systems biology. Biological growth in development and disease is increasingly studied using the framework of morphoelasticity. Here, we demonstrate the potential for morphoelastic simulations through examples of volume, area, and length growth, inspired by tumor expansion, chronic bronchitis, brain development, intestine formation, plant shape, and myopia. We review the systems biology of living systems in light of biochemical and optical stimuli and classify different types of growth to facilitate the design of growth models for various biological systems within this generic framework. Exploring the systems biology of growth introduces a new venue to control and manipulate embryonic development, disease progression, and clinical intervention.

Patient-Specific Airway Wall Remodeling in Chronic Lung Disease.

Chronic lung disease affects more than a quarter of the adult population; yet, the mechanics of the airways are poorly understood. The pathophysiology of chronic lung disease is commonly characterized by mucosal growth and smooth muscle contraction of the airways, which initiate an inward folding of the mucosal layer and progressive airflow obstruction. Since the degree of obstruction is closely correlated with the number of folds, mucosal folding has been extensively studied in idealized circular cross sections. However, airflow obstruction has never been studied in real airway geometries; the behavior of imperfect, non-cylindrical, continuously branching airways remains unknown. Here we model the effects of chronic lung disease using the nonlinear field theories of mechanics supplemented by the theory of finite growth. We perform finite element analysis of patient-specific Y-branch segments created from magnetic resonance images. We demonstrate that the mucosal folding pattern is insensitive to the specific airway geometry, but that it critically depends on the mucosal and submucosal stiffness, thickness, and loading mechanism. Our results suggests that patient-specific airway models with inherent geometric imperfections are more sensitive to obstruction than idealized circular models. Our models help to explain the pathophysiology of airway obstruction in chronic lung disease and hold promise to improve the diagnostics and treatment of asthma, bronchitis, chronic obstructive pulmonary disease, and respiratory failure.

On the Role of Mechanics in Chronic Lung Disease.

Progressive airflow obstruction is a classical hallmark of chronic lung disease, affecting more than one fourth of the adult population. As the disease progresses, the inner layer of the airway wall grows, folds inwards, and narrows the lumen. The critical failure conditions for airway folding have been studied intensely for idealized circular cross-sections. However, the role of airway branching during this process is unknown. Here, we show that the geometry of the bronchial tree plays a crucial role in chronic airway obstruction and that critical failure conditions vary significantly along a branching airway segment. We perform systematic parametric studies for varying airway cross-sections using a computational model for mucosal thickening based on the theory of finite growth. Our simulations indicate that smaller airways are at a higher risk of narrowing than larger airways and that regions away from a branch narrow more drastically than regions close to a branch. These results agree with clinical observations and could help explain the underlying mechanisms of progressive airway obstruction. Understanding growth-induced instabilities in constrained geometries has immediate biomedical applications beyond asthma and chronic bronchitis in the diagnostics and treatment of chronic gastritis, obstructive sleep apnea and breast cancer.