Journal of Applied Physiology Fuel your research with LabChart
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Appl Physiol 82: 233-239, 1997;
8750-7587/97 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Hill, M. J.
Right arrow Articles by Lambert, R. K.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Hill, M. J.
Right arrow Articles by Lambert, R. K.

Journal of Applied Physiology
Vol. 82, No. 1, pp. 233-239, January 1997
GAS EXCHANGE, MECHANICS, AND AIRWAYS

Effects of surface tension and intraluminal fluid on mechanics of small airways

Mark J. Hill1, Theodore A. Wilson1, and Rodney K. Lambert2

1 Department of Aerospace Engineering and Mechanics, University of Minnesota, Minneapolis, Minnesota 55455; and 2 Department of Physics and Biophysics, Massey University, Palmerston North, New Zealand

ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
APPENDIX
ACKNOWLEDGEMENTS
FOOTNOTES
REFERENCES


ABSTRACT

Hill, Mark J., Theodore A. Wilson, and Rodney K. Lambert. Effects of surface tension and intraluminal fluid on the mechanics of small airways. J. Appl. Physiol. 82(1): 233-239, 1997.---Airway constriction is accompanied by folding of the mucosa to form ridges that run axially along the inner surface of the airways. The muscosa has been modeled (R. K. Lambert. J. Appl. Physiol. 71: 666-673, 1991 [Medline] ) as a thin elastic layer with a finite bending stiffness, and the contribution of its bending stiffness to airway elastance has been computed. In this study, we extend that work by including surface tension and intraluminal fluid in the model. With surface tension, the pressure on the inner surface of the elastic mucosa is modified by the pressure difference across the air-liquid interface. As folds form in the mucosa, intraluminal fluid collects in pools in the depressions formed by the folds, and the curvature of the air-liquid interface becomes nonuniform. If the amount of intraluminal fluid is small, <2% of luminal volume, the pools of intraluminal fluid are small, the air-liquid interface nearly coincides with the surface of the mucosa, and the area of the air-liquid interface remains constant as airway cross-sectional area decreases. In that case, surface energy is independent of airway area, and surface tension has no effect on airway mechanics. If the amount of intraluminal fluid is >2%, the area of the air-liquid interface decreases as airway cross-sectional area decreases, and surface tension contributes to airway compression. The model predicts that surface tension plus intraluminal fluid can cause an instability in the area-pressure curve of small airways. This instability provides a mechanism for abrupt airway closure and abrupt reopening at a higher opening pressure.

mathematical model; airway compliance; airway closure; airway opening pressure


INTRODUCTION

REPORTS THAT DESCRIBE ridges that run axially along the inner surface of airways can be traced back to the last century. Until recently, this feature of airway geometry has been largely ignored, and the airways have been modeled as smooth-bored tubes that expand or contract uniformly around their circumference. For example, in their theoretical studies of surface tension and intraluminal fluid, Halpern and Grotberg (5, 6), Johnson et al. (10), and Otis et al. (17) modeled the airway lumen as circular. The work of James et al. (8, 9) brought new attention to the folding of the airway epithelium. They showed that the length of the perimeter of the basement membrane remains constant as airway area decreases, even for the large area decreases that are caused by smooth muscle contraction, and they used basement membrane length to characterize the intrinsic size of airways that were fixed during bronchoconstriction. The constant perimeter is accommodated inside a smaller outer circumference by folding of the mucosa to form ridges that penetrate into the airway lumen. Yager et al. (20) reported observations of airways in which the interstices between the folds were filled with fluid, and this observation added to the interest in mucosal folding. Lambert (12) then suggested that the basement membrane, which contains elastin and collagen, might resist bending and that its bending stiffness might be an important component of airway elastance. Lambert et al. (13) used this model to interpret data that showed correlations between number of folds, airway area, and wall thickness in sheep.

Lambert (12) modeled the airway as a composite structure consisting of an elastic inner layer, which represents the mucosa, an intermediate fluid layer, which represents the submucosa, and an outer layer of smooth muscle, which imposes a compressive stress on the inner layers when it contracts. Surface tension was not included in the model. He computed luminal area as a function of transmucosal pressure for this model. Yager et al. (20) discussed the effect of intraluminal fluid on luminal area and estimated the effect of surface tension on the airway area-transmural pressure curve in the case where the interstices of the folds were completely filled and the air-liquid interface formed a cylindrical surface inside the folds. In this study, we describe a model in which surface tension and intraluminal fluid are added to the Lambert model. The effect of surface tension is analyzed for the case where fluid partially fills the folds. We find that surface tension has no effect on airway area if the amount of intraluminal fluid is small and the fluid lining layer is thin. In that case, the air-liquid interface coincides with the epithelial surface, and the length of the air-liquid interface is constant because the perimeter of the mucosa is constant. Therefore, surface energy is independent of luminal area, and surface tension does not affect airway mechanics. Thus one of the principal consequences of mucosal folding is to remove the compressive effect of surface tension if the airway is relatively dry. However, with fluid volumes on the order of a few percent of luminal volume, fluid pools form in the folds, and the length of the air-liquid interface changes as airway area changes and the shape of the folds changes. In that case, the model predicts that surface tension can cause area-pressure instabilities and fluid-balance instabilities.


METHODS

Model. The model for the airway is the same as that of Lambert (12). The mucosa is modeled as a thin elastic cylinder with Young's modulus (E), thickness (t), and undeformed radius (R) that is loaded by the pressures that are applied to its inner and outer surfaces. The submucosa is modeled as a fluid in which the pressure (Psub) is assumed to be uniform. In Lambert's model, the air pressure in the lumen (Pair) was applied to the inner surface of the mucosa, and transmucosal pressure, Pair - Psub, was uniform. For negative values of Pair - Psub, the mucosa is compressed, and at a critical negative transmucosal pressure, Pcr = - (D/R3)(n2 - 1), where D is the bending stiffness of the mucosa and n is the number of folds in the bending mode, the mucosa buckles. He analyzed the buckling of the mucosa for Pair - Psub < Pcr.

We have modified Lambert's model by including a lining layer of intraluminal fluid with surface tension (gamma ) at the interface between the fluid and the air in the lumen. If the cylinder is circular and the fluid layer is thin, the curvature of the air-liquid interface is ~1/R and Pfl = Pair - gamma /R, where Pfl is the pressure applied to the inner surface of the mucosa. Therefore, with intraluminal fluid and surface tension, buckling occurs at a less negative value of Pair - Psub. If the fluid layer is thin, the folds that form when the cylinder buckles quickly penetrate the fluid layer. When that occurs, the curvature of the air-liquid interface is no longer uniform around the circumference of the wall, and the pressure applied to the inner surface is no longer uniform. We have analyzed the buckling of the cylinder under this nonuniform loading.

Governing equations. A segment of the folded elastic wall is shown in Fig. 1. The elastic wall carries a shear force (V), a tensile force (T), and a bending moment (M) that vary with position (s) measured along the perimeter of the cylinder. The inclination of the tangent to the wall is denoted as theta . The equations of equilibrium for the wall are
d<IT>T</IT>/d<IT>s</IT> + <IT>V</IT>(d&thgr;/d<IT>s</IT>) = 0 (1)
d<IT>V</IT>/d<IT>s</IT> − <IT>T</IT>(d&thgr;/d<IT>s</IT>) = P<SUB>sub</SUB> − P<SUB>fl</SUB> (2)
d<IT>M</IT>/d<IT>s</IT> − <IT>V</IT> = 0 (3)
In addition, M is given by the product of D and the difference between the curvature of the wall (dtheta /ds) and the undeformed curvature (1/R)
<IT>M</IT> = <IT>D</IT>(d&thgr;/d<IT>s</IT> − 1/<IT>R</IT>) (4)

Fig. 1. Schematic diagram of model for mucosa and intraluminal fluid. Mucosa is modeled as a thin elastic sheet that carries shear and tensile forces (V, T) and bending moment (M). Coordinates x and y lie normal to and along plane of symmetry of a fold. s, Distance along mucosa from midplane; s*, boundary of pool; theta , angle between tangent to sheet and horizontal; a, radius of curvature of surface of pool; beta , angle subtended by arc that forms surface of pool; Pair, Pfl, and Psub: pressures in air in lumen of airway, in fluid layer, and in submucosa that surround mucosa, respectively. A pool of fluid covers mucosa to point s = s*.
[View Larger Version of this Image (23K GIF file)]

Pfl that appears in Eq. 2 differs from the uniform Pair by the pressure difference across the air-liquid interface. The pressure difference across the air-liquid interface is given by the product of surface tension and interface curvature. Surface tension is assumed to be uniform, but the curvature of the air-liquid interface varies around the circumference. The curvature of the interface above the fluid pool must be uniform if Pfl is uniform within the pool. The radius of curvature of this surface is denoted as a. Beyond the boundary of the pool (s*), the air-liquid interface coincides with the inner surface of the wall, and the interface curvature is the same as wall curvature (dtheta /ds). Therefore, Pfl is given by
<AR><R><C>P<SUB>fl</SUB> = </C></R><R><C> </C></R></AR><AR><R><C>P<SUB>air</SUB> − &ggr;/<IT>a</IT>,</C></R><R><C>P<SUB>air</SUB> − &ggr;(d&thgr;/d<IT>s</IT>),</C></R></AR>  <AR><R><C><IT>s</IT> < <IT>s</IT>*</C></R><R><C><IT>s</IT> > <IT>s</IT>*</C></R></AR> (5)
Solutions for n identical symmetric folds around the circumference can be generated from solutions for one interval of length (pi R/n), where R is the radius of the undeformed cylinder. The point s = 0 is taken to be the point of maximum wall curvature at the depth of the fold, and the point s = pi R/n is the point of minimum curvature. The boundary conditions for V and theta  at these points are
<IT>V</IT>(0) = <IT>V</IT>(&pgr;<IT>R</IT>/<IT>n</IT>) = 0     (6)
&thgr;(0) = 0;  &thgr;(&pgr;<IT>R</IT>/<IT>n</IT>) = &pgr;/<IT>n</IT> (7)
The interface above the pool must join the interface on the wall smoothly; that is, the circular arc of radius a and subtended angle beta  must meet the wall and be tangent to the wall at s = s*. This requires that the solution satisfy the following conditions
<IT>a</IT> = <IT>x</IT>(<IT>s</IT>*)/sin &thgr;(<IT>s</IT>*) (8)
&bgr;(<IT>s</IT>*) = &thgr;(<IT>s</IT>*) (9)
Finally, the volume of intraluminal fluid (Vfl) is specified. This is given by
 <IT>V</IT><SUB>fl</SUB>/2<IT>n</IT> = <LIM><OP>∫</OP><LL>0</LL><UL><IT>s</IT>*</UL></LIM> <IT>x</IT> sin &thgr;d<IT>s</IT> − <FR><NU>1</NU><DE>2</DE></FR> &bgr;(<IT>s</IT>*)<IT>a</IT><SUP>2</SUP> + <FR><NU>1</NU><DE>2</DE></FR> <IT>a</IT><SUP>2</SUP> cos &thgr;(<IT>s</IT>*) sin &thgr;(<IT>s</IT>*) (10)
where x is the perpendicular distance between the midplane and the wall.

Equations 1-4, for the loading described by Eq. 5 and given values of Pair - Psub, Vfl, and n, were integrated numerically from s = 0 to s = pi R/n. Because the values of T and dtheta /ds at s = 0 are not imposed by the boundary conditions and the values of a and s* are not known a priori, the solution method required a search through values of T(0), dtheta /ds(0), a, beta , and s* to find solutions that satisfied the boundary conditions at s = pi  R/n, given by Eqs. 6 and 7, and the conditions at s* given by Eqs. 8-10. Solutions were obtained for a series of values of Pair - Psub up to the value of Pair - Psub at which the sides of the folds touched at the midplane of the fold.

Particular values of the parameters were used in the calculations. On the basis of the data of Codd et al. (2), Young's modulus E was taken to be 200 cmH2O, the thickness-to-radius ratio (t/R) was taken as 0.04, and the bending stiffness D/R3 was taken to be 0.0015 cmH2O. The value of R was chosen to be 0.5 mm. Solutions were obtained for two values of gamma , 10 and 20 dyn/cm, respectively; for a series of values of Vfl, ranging from 1 to 4% of the area of the undeformed cylinder, pi R2; but for only one value of n, n = 16. This value of n was chosen because it lies in the middle of the range of n observed by Lambert et al. (13).


RESULTS

The area-pressure curves that were computed from the model are shown in Fig. 2. In these curves, airway area (A), normalized by the undeformed area (Ao), is plotted against the pressure difference across the liquid lining and mucosa, Pair - Psub. For pressures greater than the buckling pressure, the cylinder remains circular. The specific elastance, Ao[d(Pa - Psub)/dA] of the circular cylinder is Et/2R. Therefore, for Pair - Psub greater than the critical buckling pressure, the curves shown in Fig. 2 have slope 2R/Et. At the buckling pressure, area decreases sharply with decreasing pressure. The curves for fluid volumes <1% nearly coincide with Lambert's curve for n = 16 (12). With intraluminal fluid, the buckling pressure is increased by gamma /R. However, if the amount of fluid is small, the fluid pools are small, the air-liquid interface nearly coincides with the wall surface, and the area-pressure curve is nearly the same as the area-pressure curve for the dry airway. If the amount of intraluminal fluid is greater, the area-pressure curves depart from those for dry airways. The shapes of a half-fold at various points along the area-pressure curves are also shown. The pressure in the fluid pool (Pfl) relative to gas pressure (Pair) is plotted vs. Vfl at fixed values of Pair - Psub in Fig. 3.


Fig. 2. Area-pressure curves for 2 values of surface tension (gamma ), namely, gamma  = 10 (A) and 20 dyn/cm (B), respectively. Ratio of luminal area (A) to undeformed area (Ao; A/Ao) is shown vs. pressure difference across mucosa (Pair - Psub) for different volumes of intraluminal fluid, expressed as %undeformed luminal volume. Shapes of mucosal fold and fluid pool are shown for 3 cases (a-c).
[View Larger Version of this Image (16K GIF file)]


Fig. 3. Pfl relative to Pair in lumen for Pair - Psub = -0.3 cmH2O vs. volume of fluid (Vfl) in airway lumen as %luminal volume.
[View Larger Version of this Image (10K GIF file)]


DISCUSSION

Modeling assumptions. The model for the mucosa that is analyzed here is an idealized and simplified model. First, the geometry is idealized. It is assumed that the folding is uniform along the axis of the airway so that the deformation is two-dimensional. It is also assumed that the 16 folds around the circumference are identical, whereas the folding seen in airways is irregular. The material properties of the the airway are also modeled simply. The mucosa is modeled as a linear elastic material. Although geometric nonlinearities are included in the analysis of the folding, material nonlinearities are not considered. Finally, the submucosa has been modeled as a fluid in which the pressure is uniform. The material properties of the submucosa are unknown. It appears that the submucosa must be highly deformable because the deformations that occur during mucosal folding are large, but it is not known that the submucosa is completely inelastic. Also, the submucosa is thin compared with the depth of the folds. Contact forces between the mucosa and the surrounding smooth muscle layer or blood vessels in the submucosa may impose a nonuniform stress on the submucosal side of the mucosa.

Particular values for the parameters of the model have been chosen. The values of some of these parameters are uncertain, and the values of some are known to vary over a wide range. The model equations can be nondimensionalized (12), and the nondimensionalization yields scaling laws for the dependence of the variables on the parameter values. These laws allow the results reported here to be scaled for other parameter values. For example, the bending stiffness (D) of the mucosa is not well established, and the size of airways, which is described by the parameter R, varies over a wide range. Scaling laws show that the pressure differences shown in Fig. 2 are proportional to D/R3. Because D is proportional to t3, D/R3 would be independent of airway size if the material properties of the mucosa were the same for all airways and mucosal thickness were scaled with R. Observed values of n also vary considerably. Roughly speaking, Pair - Psub for a given deformation scales with n2. The value of gamma  in the airways is also uncertain. Alveolar surface tension has been measured and is known to change as lung volume changes, but the distribution of gamma  along the airways is unknown. The effect of gamma  scales with gamma /R. Thus, for a given surface tension, the effects of intraluminal fluid and surface tension are greater in smaller airways than in larger airways. If the value of R had been chosen to be 1 mm instead of 0.5 mm, the curves for gamma  = 20 dyn/cm for that airway would be the same as the curves for gamma  = 10 dyn/cm that are shown. We chose to analyze the deformation of a small airway to emphasize the effect of surface tension. Finally, we show results for intraluminal fluid volumes of 1-4% of undeformed luminal area. Yager et al. (21) found intraluminal fluid volumes of ~2% of luminal area at both high and lower lung volumes. For fluid volumes greater than ~6%, the fluid simply floods the space between the folds, and the air-liquid interface is circular. In that case, transmucosal pressure is uniform, and the effect of surface tension on mucosal folding can be obtained from Lambert's results (12).

We assumed that the mucosa folded in a 16-lobe mode. The buckling pressure for this mode is larger than the buckling pressure for lower values of n, and one would expect that buckling with lower values of n would be favored. Nonuniformities in mucosal properties around the circumference of the airway could affect the buckling mode. Lambert (12) suggested that the higher value of n is imposed by the geometric condition that the folded mucosa fit within the surrounding ring of smooth muscle, and Mitzner and Wagner (15) suggested that folds form over blood vessels that impose an initial deformation on the mucosa. Thus, although the mechanism is unknown, multifold buckling occurs, and we have chosen a value of n that lies in the middle of the observed range. However, the computed curves do not take account of the mechanisms that produce multifold buckling, and the sharp change in the slope of the curves shown in Fig. 2 seems unrealistic. It seems likely that the transitions between the segment of the curve for A/Ao >1 to the segments for A/Ao <1 are smoother than those shown in Fig. 2.

Smooth muscle contraction simultaneously imposes both a compressive pressure difference across the mucosa and maintains a circular outer boundary that confines the mucosa. If there were no smooth muscle tone and transmural pressure were negative, the entire airway wall, mucosa, submucosa, and smooth muscle could buckle in the lower energy two-lobe mode. Excised airways with no muscle tone collapse in this two-lobe mode. However, in intact airways, the surrounding parenchyma also tends to maintain a circular outer boundary for the airway. Certainly, reports of ridges in the surface of airways are not restricted to preparations in which smooth muscle has been activated by bronchoconstrictive agents. In the Appendix, we analyze the effect of the surrounding parenchyma on the critical buckling pressure and conclude that, in the absence of smooth muscle tension, the critical buckling pressure for 16-lobe buckling of the mucosa within a circular outer boundary is smaller than the critical pressure for two-lobe buckling of the entire airway wall if transpulmonary pressure is >1 cmH2O. Because passive or active smooth muscle tension and parenchymal stiffness both contribute to maintaining a circular outer boundary, the outer boundary of intact airways may remain circular, and the model described here may be an appropriate model for elastance of the mucosa of normal airways.

Finally, our analysis describes the configuration of the mucosa and intraluminal fluid in static equilibrium. The equilibrium configuration depends on Pair - Psub, and if Pair - Psub changes, the system must relax from the old to the new configuration. The relaxation time for redistribution of intraluminal fluid can be estimated by using the equations that describe two-dimensional viscous fluid flow. The volume flow rate (Q) in the fluid layer is given by Q = (h3/12µ)(Delta p/d) where h is the thickness of the layer, µ is fluid viscosity, and Delta p is the difference between the pressures at two points separated by a distance (d). To change the shape of the fluid layer, a volume of the order of dh must collect, and the time required to transport this volume of fluid is dh/Q. For a value of µ equal to the viscosity of water, Delta p equal to the difference between pressure in the pool and pressure in the thin layer, and h and d equal to pool dimensions, the relaxation time is 0.1 s. Thus, for changes of configuration that occur over time scales that are >0.1 s, our assumption that pressure is uniform in the pool seems justified. However, Q is proportional to h3. An initially uniform layer of fluid would redistribute on a time scale of 0.1 s, but as the fluid layer on the ridge thins, flow becomes smaller and the rate of thinning decreases. In addition, our solution contains a discontinuity in the curvature of the air-liquid interface and hence, a discontinuity in fluid pressure, at the boundary between the pool and the thin layer. Our model does not address the complex mechanics of the fluid layer that covers the ridge; we have simply assumed that this layer is relatively thin.

Airway mechanics. The primary purpose of this work is to investigate the effects of surface tension and intraluminal fluid on the folding of the mucosa that occurs when airways are compressed. A primary, and perhaps surprising, result is that surface tension has no effect on airway compression if the amount of intraluminal fluid is small. This result can be understood immediately by considering the effect of the folds on surface energy. By folding, the perimeter of the airway remains nearly constant as airway area decreases. If the amount of fluid is small so that the pools of intraluminal fluid are small, the air-liquid interface lies near the inner surface of the airway, and the area of the air-liquid interface and hence the surface energy also remain constant. If surface energy is independent of airway area, surface tension cannot affect the pressure-area curve. This result can also be obtained from Eq. 2. If the air-liquid interface coincides with the airway wall, the value of Pfl is Pair - gamma (dtheta /ds). Another term in Eq. 2, T(dtheta /ds), has the same form as the term gamma (dtheta /ds). Therefore, a change in gamma  can be balanced by an equal and opposite change in T, leaving all other equations unchanged. Thus surface tension can be balanced by wall compression with no change in the area-pressure curve, and Lambert's original model (12), in which surface tension was neglected, is valid for any value of surface tension if the airway is relatively dry. Conversely, changes in luminal area would have little effect on surface tension in that case. In the parenchyma, surface tension changes because alveolar surface area changes with lung volume. In the airways, if the amount of intraluminal fluid is small, airway perimeter is nearly independent of airway luminal area, and airway wall area changes only by the change in airway length as lung volume changes. Therefore, the direct effect of lung volume on surface tension in the airways is weak.

The combination of intraluminal fluid and surface tension affect the area-transmucosal pressure curves as shown in Fig. 2. Before buckling of the mucosa, the fluid forms a thin layer around the circumference of the lumen, and the pressure in the fluid layer is lower than luminal pressure because of the pressure drop across the air-liquid interface. Therefore, the compressive pressure applied to the mucosa is >Pair - Psub by gamma /R or 0.2 and 0.4 cmH2O for R of 0.5 mm and gamma  of 10 and 20 dyn/cm, respectively. Thus, with intraluminal fluid, buckling occurs at values of Pair - Psub that are greater than the value for the dry tube. For 1% fluid volume, the fronts of the folds quickly push through the thin fluid lining, and fluid collects in the depths of the folds. However, the curvature of the surface of the pools nearly equals the curvature of the mucosa, and the region covered by the pools is small so that the area-pressure curves are displaced only slightly to the right along the pressure axis. For 2% fluid, the area-pressure curve lies a bit farther to the right down to A/Ao of ~0.2. At that value of A/Ao, fluid fills most of the space within the fold and the air-liquid interface is near the neck between adjacent folds. The neck is narrow and the radius of curvature of the air-liquid interface is small as it passes throught the neck. As a result, pressure in the fluid pool is much more negative than luminal pressure, and the region covered by the pool is subjected to a large compressive stress. This buckled configuration is maintained by the negative pressure in the pool, and Pair - Psub is more positive. In fact, this configuration occurs at positive values of Pair - Psub for both gamma  = 10 and gamma  = 20 dyn/cm. After the air-liquid interface passes through the neck, its radius of curvature increases, pressure in the fluid pool is less negative, and more negative values of Pair - Psub are required to decrease area further. The same sequence occurs as A/Ao decreases with 3% fluid. However, the air-liquid interface passes through the neck at a slightly higher value of A/Ao, the minimum radius of curvature as the interface passes through the neck is larger, and the maximum value of Pair - Psub is smaller than for 2% fluid. The sequence is illustrated by diagrams that show three configurations, all of which occur at Pair - Psub = -0.5.

The combination of surface tension and intraluminal fluid can produce convoluted area-pressure curves. The segments of the curves with positive slope are locally stable branches. That is, if Pair - Psub is held constant, more work is required to deform the tube along the equilibrium curve than would be provided by the given Pair - Psub. The segments of the curve with negative slope are locally unstable. Thus configuration b shown in Fig. 2 is unstable. If Pair - Psub were held constant, the tube would jump to either configuration a or c. The combination of stable and unstable segments leads to hysteresis in airway closing and opening. For example, in the instance where gamma  = 10 and there is 3% fluid volume, area would decrese smoothly along the upper segment of the area-pressure curve down to A/Ao = 0.25 as Pair - Psub decreased to -0.6. For a further decrease in Pair - Psub, A/Ao would drop to the value A/Ao = 0.1 on the lower segment of the curve. If Pair - Psub were then increased, area would increase along the lower segment until Pair - Psub reached -0.3. For any further increase, A/Ao would jump from 0.15 to 0.8. For some values of gamma  and fluid volume, namely, gamma  = 10 dyn/cm and 2% fluid or gamma = 20 dyn/cm and 2-4% fluid, positive values of Pair - Psub are required to reopen the airway.

Recent studies (3, 16, 18) have shown that a positive transmural pressure is required to open airways that have closed. The airways appear to open abruptly at a critical positive transmural pressure. In some modeling studies, the positive opening pressure has been ascribed to the pressure required to drive a miniscus axially into the closed segment of the airway (3). Our model describes a different mechanism that could maintain the airway nearly closed at small positive transmural pressures and produce a sudden opening at a critical pressure. The mechanism does not require an axial nonuniformity between open and closed segments of the airway; it is an instability in the opening of an axially uniform tube with fluid-filled folds.

Fluid instabilities. In addition to the mechanical instability that can occur for certain values of gamma  and fluid volume, a fluid instability can also occur. The mechanism is described by Fig. 3, in which pressure in the fluid pool, relative to luminal pressure, is plotted against fluid volume at a fixed value of Pair - Psub. This figure shows that at Pair - Psub = -0.3, fluid pressure decreases as fluid volume increases because the radius of curvature of the air-liquid interface decreases as the volume of the pool increases. In that situation, a uniform distribution of fluid along the axis of the airway is unstable. If fluid volume were slightly greater at one point along the axis, pressure would be lower at that point and fluid would flow toward it. Pressure would fall further at the point where volume was high and would rise at other points, driving further flow. For gamma  = 20 dyn/cm, fluid accumulation would be terminated by a mechanical instability that occurs at a fluid volume of 3%.

Both the mechanical and the fluid instabity mechanisms are driven by the decrease in radius of curvature as the air-liquid interface approaches the neck between adjacent folds. Folds that are one-half or more filled with fluid are susceptible to these instabilities. Therefore, it would be unlikely to find airways in which folds are one-half filled with fluid; it would be more likely to observe airways that are nearly dry or airways in which the folds are filled past the neck.

Comparison with airway elastance. Our modeling of the effects of surface tension and intraluminal fluid describes potential mechanisms for instabilities that may play a role in airway closure and reopening, and we have argued that this model of mucosal bending is pertinent both to airway compression driven by smooth muscle contraction and compression by transmural pressure differences due to flow. However, there is a fundamental question about the significance of these results; it is not clear that the magnitudes of mucosal bending stiffness and surface tension are large enough for these mechanisms to be significant. Because surface tension affects the area curves shown in Fig. 2, surface tension is significant, relative to mucosal bending stiffness, in 1-mm airways. However, it is not clear that the magnitudes of both bending stiffness and surface tension in small airways are significant. Figure 2 shows the relationship between area and transmucosal pressures. Transmucosal pressure is important if it is significant compared with the pressure generated by smooth muscle tension or if it is a significant fraction of transmural pressure. The range of pressures shown in Fig. 2 is small, on the order of 1 cmH2O. Maximum smooth muscle contraction generates compressive pressures of 30-40 cmH2O (4). Therefore, mucosal bending stiffness offers little resistance to maximum smooth muscle contraction. However, the bending stiffness of the mucosa may have been underestimated. Also, airway wall thickening occurs in asthma, and bending stiffness is proportional to the cube of thickness. A doubling of mucosal thickness would produce an increase of nearly a factor of ten in bending stiffness.

Although it appears that the mechanisms described here would not be significant during maximum muscle contraction, they may be significant at lower levels of smooth muscle tension. To determine whether they contribute significantly to normal airway elastance, the theoretical curves shown in Fig. 2 should be compared with experimental curves of airway area vs. transmural pressure. A precise comparison is not possible at present for a number of reasons. First, the pressures shown in Fig. 2 are the difference between gas pressure in the lumen and pressure in the submucosa. The pressure difference across the muscle and other elastic components of the airway wall must be added to the pressure differences shown in Fig. 2 to obtain transmural pressure. Therefore, plots of airway area vs. transmural pressure would be expected to be shifted to the right along the pressure axis, and the slopes of the curves would be reduced because of the elastance of the smooth muscle and other components of the airway wall, but the magnitudes of these shifts are unknown. Also, the relationship between Ao, the area of the undeformed mucosa, and maximum airway area for large transmural pressures would be needed to align the ordinate in Fig. 2 with the ordinate of experimental area-pressure curves. At present, the transmural pressure at which mucosal folding begins is not known.

Although a precise comparison between the theoretical curves shown in Fig. 2 and experimental data is not possible, an approximate comparison can be made between the magnitude of the elastance of the model mucosa and the magnitude of airway elastance. The value of specific elastance, Ao d(Pair - Psub)/dA, for the model curves is ~1 cmH2O. Gunst and Stropp (4) report area-transmural pressure curves of fluid-filled 4-mm canine airways. To obtain a value of specific elastance that is comparable to the model elastance, an approximate Ao value for the canine airways is needed. The report of Yager et al. (21) shows that the epithelium of airways of guinea pig lungs fixed at a transpulmonary pressure of 5 cmH2O is fairly smooth. The mucosa of sheep airways fixed at zero transpulmonary pressure (13) have deep folds. We will assume that folding begins at a transmural pressure of ~2 cmH2O and take this value of A as the value of Ao. Then, the data of Gunst and Stropp (4) yield a value of specific elastance of 2-3 cmH2O. The model mucosa would provide a significant fraction of the elastance of these 4-mm airways, and because airway elastance decreases with decreasing airway size (7, 14), mucosal elastance would provide a larger component of the elastance of smaller airways.


ACKNOWLEDGEMENTS

We thank Wayne Mitzner for pointing out the early literature on mucosal folding.


FOOTNOTES

   This work was supported by National Heart, Lung, and Blood Institute Grant HL-49788.

Address for reprint requests: T. A. Wilson, 107 Akerman Hall, 110 Union St. SE, Minneapolis, MN 55455.

Received 13 March 1996; accepted in final form 9 September 1996.


APPENDIX

The effect of the surrounding parenchyma on the buckling of an airway can be modeled by representing the surrounding parenchyma as an elastic continuum and adding this feature to the model for a thin-walled elastic tube. This model can be analyzed by combining two results from classic elasticity theory. The first is the analysis of buckling of a thin-walled tube supported by a Winkler foundation (1). A Winkler foundation provides a restoring stress that is proportional to the local displacement. The critical buckling pressure (Pcr W) for a thin-walled tube supported by a Winkler foundation is given by the equation
P<SUB>cr W</SUB> = − (<IT>n</IT><SUP>2</SUP> − 1)(<IT>D</IT>w/<IT>R</IT><SUP>3</SUP>) − <IT>kR</IT>/(<IT>n</IT><SUP>2</SUP> − 1) (A1)
where n is the number of folds in the buckling mode, Dw is the bending stiffness of the tube wall, R is the tube radius, and k is the spring constant of the Winkler foundation. The first term on the right side of Eq. A1 describes the buckling pressure of an isolated thin-walled elastic tube. The second term describes the additional compressive pressure required to overcome the restoring forces provided by the foundation. In the buckled state, the normal stress applied by the foundation and the radial displacement of the tube are both proportional to cos (npsi ) where psi  is the azimuthal angle around the axis of the tube. The ratio of the two is k. The effective k of an elastic continuum can be obtained from a second result from classic elasticity theory (19). The stress and displacements in an elastic body loaded by normal stresses acting at a circular inner boundary can be expressed as series that reduce to Fourier series in the azimuthal angle at the boundary. The ratio of normal stress to radial displacement for each term in the series provides k for that order. The ratio depends on the order, or buckling mode (n), and is given by
<IT>k</IT> = <IT>E</IT><SUB>par</SUB>(<IT>n</IT><SUP>2</SUP> − 1)/<IT>R</IT>(2<IT>n</IT> + 1 − <IT>v</IT>) (A2)
In Eq. A2, Epar and v are the Young's modulus and Poisson's ratio of the parenchyma. Substituting this value for k into Eq. A1 yields the following expression for Pcr W
P<SUB>cr W</SUB> = − (<IT>n</IT><SUP>2</SUP> − 1)(<IT>D</IT>w/<IT>R</IT><SUP>3</SUP>) − <IT>E</IT><SUB>par</SUB> /(2<IT>n</IT> + 1 − <IT>v</IT>) (A3)

The objective here is to compare the buckling pressures for two different buckling geometries, buckling of the entire wall in the n = 2-lobe mode and buckling of the mucosa in the n = 16-lobe mode, with the outer boundary of the airway remaining circular. If Pcr W, the buckling pressure given by Eq. A3 for n = 2, is more negative than -0.4, the critical pressure for mucosal buckling shown in Fig. 2, mucosal buckling will occur. Thus we seek the value of Epar for which Pcr W <-0.4 cmH2O
3(<IT>D</IT>w/<IT>R</IT><SUP>3</SUP>) + <IT>E</IT><SUB>par</SUB> /(5 − <IT>v</IT>) > 0.4 (A4)
The value of the bending stiffness of the airway wall (Dw) is unknown. It is clear that Dw must be greater than the bending stiffness of the mucosa (D) because the airway wall includes the mucosa and the smooth muscle layer and is much thicker than the mucosa. However, a conservative estimate of the value of Epar for which Eq. A4 is satisfied is obtained by assuming that Dw is <102 D and that the first term in Eq. A4 is negligible. The value of v for the parenchyma is ~0.4. Thus a conservative estimate of the value of Epar that will prevent buckling of the entire airway in the n = 2-lobe mode is Epar >2 cmH2O. Young's modulus has not been measured at very low transpulmonary pressures, but an extrapolation of the data of Lai-Fook et al. (11) to lower pressures yields the estimate Epar >2 cmH20 for transpulmonary pressure >1 cmH2O.


REFERENCES

1. Brush, D. O., and B. O. Almroth. Buckling of Bars, Plates, and Shells. New York: McGraw-Hill, 1975, p. 133-134.
2. Codd, S. L., R. K. Lambert, M. R. Alley, and R. J. Pack. Tensile stiffness of ovine tracheal wall. J. Appl. Physiol. 76: 2626-2635, 1994.
3. Gaver, D. P., R. W. Samsel, and J. Solway. Effects of surface tension and viscosity on airway reopening. J. Appl. Physiol. 69: 74-85, 1990. [Abstract/Free Full Text]
4. Gunst, S. J., and J. Q. Stropp. Pressure-volume and length-stress relationships in canine bronchi in vitro. J. Appl. Physiol. 64: 2522-2531, 1988. [Abstract/Free Full Text]
5. Halpern, D., and J. B. Grotberg. Fluid-elastic instability of a thin film coating a flexible tube. J. Fluid Mech. 194: 619-632, 1992.
6. Halpern, D., and J. B. Grotberg. Surfactant effects on fluid-elastic instabilities of liquid-lined flexible tubes: a model of airway closure. J. Biomech. Eng. 119: 271-277, 1993.
7. Hyatt, R. E., R. Sittipong, S. Olafsson, and W. A Porter. Some factors determining pulmonary pressure-flow behavior at high rates of airflow. In: Airway Dynamics, edited by A. Bouhuys. Springfield, IL: Thomas, 1970, p. 43-83.
8. James, A. L., J. C. Hogg, L. A. Dunn, and P. D. Paré. The use of the internal perimeter to compare airway size and to calculate smooth muscle shortening. Am. Rev. Respir. Dis. 138: 136-139, 1988. [Medline]
9. James, A. L., P. D. Paré, and J. C. Hogg. Effects of lung volume, bronchoconstriction, and cigarette smoke on morphometric airway dimensions. J. Appl. Physiol. 64: 913-919, 1988. [Abstract/Free Full Text]
10. Johnson, M., R. D. Kamm, L. Lo, A. H. Shapiro, and T. J. Pedley. The nonlinear growth of surface tension-driven instabilities of a thin annular film. J. Fluid Mech. 233: 191-196, 1991.
11. Lai-Fook, S. J., T. A. Wilson, R. E. Hyatt, and J. R. Rodarte. Elastic constants of inflated dog lungs. J. Appl. Physiol. 40: 508-513, 1976. [Abstract/Free Full Text]
12. Lambert, R. K. Role of bronchial basement membrane in airway collapse. J. Appl. Physiol. 71: 666-673, 1991. [Abstract/Free Full Text]
13. Lambert, R. K., S. L. Codd, M. R. Alley, and R. J. Pack. Physical determinants of bronchial mucosal folding. J. Appl. Physiol. 77: 1206-1219, 1994. [Abstract/Free Full Text]
14. Martin, H. B., and D. F. Proctor. Pressure-volume measurements on dog bronchi. J. Appl. Physiol. 13: 337-343, 1958. [Abstract/Free Full Text]
15. Mitzner, W., and E. Wagner. Bronchial vessels in the submucosa precipitate buckling of the airway epithelium (Abstract). Am. Rev. Respir. Dis. 197: A507, 1993.
16. Naureckas, E. T., C. A. Dawson, B. S. Gerber, D. P. Gaver, H. L. Gerber, J. H. Linehan, J. Solway, and R. W. Samsel. Airway reopening pressure in isolated rat lungs. J. Appl. Physiol. 76: 1372-1377, 1994. [Abstract/Free Full Text]
17. Otis, D. R., M. Johnson, T. J. Pedley, and R. D. Kamm. The role of pulmonary surfactant in airway closure: a computational study. J. Appl. Physiol. 75: 1323-1333, 1993. [Abstract/Free Full Text]
18. Suki, B., A. Barbarasi, Z. Hantos, F. Petak, and H. E. Stanley. Avalanches and power-law behavior in lung inflation. Nature Lond. 368: 615-618, 1994. [Medline]
19. Timoshenko, S., and J. N. Goodier. Theory of Elasticity. New York: McGraw-Hill, 1951, p. 116-121.
20. Yager, D., J. P. Butler, J. Bastacky, E. Israel, G. Smith, and J. M. Drazen. Amplification of airway constriction due to liquid filling of airway interstices. J. Appl. Physiol. 66: 2873-2884, 1989. [Abstract/Free Full Text]
21. Yager, D., T. Cloutier, H. Feldman, J. Bastacky, J. M. Drazen, and R. D. Kamm. Airway surface liquid thickness as a function of lung volume in small airways of the guinea pig. J. Appl. Physiol. 77: 2333-2340, 1994. [Abstract/Free Full Text]

0161-7567/97 $5.00 Copyright © 1997 the American Physiological Society



This article has been cited by other articles:


Home page
J. Appl. Physiol.Home page
R. L. Sorkness, E. R. Bleecker, W. W. Busse, W. J. Calhoun, M. Castro, K. F. Chung, D. Curran-Everett, S. C. Erzurum, B. M. Gaston, E. Israel, et al.
Lung function in adults with stable but severe asthma: air trapping and incomplete reversal of obstruction with bronchodilation
J Appl Physiol, February 1, 2008; 104(2): 394 - 403.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
R. L. Sorkness and A. Tuffaha
Contribution of airway closure to chronic postbronchiolitis airway dysfunction in rats
J Appl Physiol, March 1, 2004; 96(3): 904 - 910.
[Abstract] [Full Text] [PDF]


Home page
Adv. Physiol. Educ.Home page
H. D. Prange
LAPLACE'S LAW AND THE ALVEOLUS: A MISCONCEPTION OF ANATOMY AND A MISAPPLICATION OF PHYSICS
Advan Physiol Educ, March 1, 2003; 27(1): 34 - 40.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
S. N. Ghadiali, J. Banks, and J. D. Swarts
Effect of surface tension and surfactant administration on Eustachian tube mechanics
J Appl Physiol, September 1, 2002; 93(3): 1007 - 1014.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
M. Okazawa, P. D. Pare, and R. K. Lambert
Compliance of peripheral airways deduced from morphometry
J Appl Physiol, December 1, 2000; 89(6): 2373 - 2381.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
R. K. Lambert and P. D. Pare
Lung parenchymal shear modulus, airway wall remodeling, and bronchial hyperresponsiveness
J Appl Physiol, July 1, 1997; 83(1): 140 - 147.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Hill, M. J.
Right arrow Articles by Lambert, R. K.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Hill, M. J.
Right arrow Articles by Lambert, R. K.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online