Journal of Applied Physiology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Appl Physiol 101: 926-933, 2006. First published June 1, 2006; doi:10.1152/japplphysiol.01548.2005
8750-7587/06 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow All Versions of this Article:
101/3/926    most recent
01548.2005v1
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 ISI Web of Science
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 ISI Web of Science (1)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Lande, B.
Right arrow Articles by Mitzner, W.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Lande, B.
Right arrow Articles by Mitzner, W.

Analysis of lung parenchyma as a parametric porous medium

Boris Lande and Wayne Mitzner

Division of Physiology, Department of Environmental Health Sciences, Johns Hopkins University, Bloomberg School of Public Health, Baltimore, Maryland

Submitted 9 December 2005 ; accepted in final form 28 May 2006


    ABSTRACT
 TOP
 ABSTRACT
 OWEN-LEWIS EQUATIONS OF LUNG...
 INPUT IMPEDANCE OF THE...
 ANALYSIS OF THE LUNG...
 DISCUSSION
 APPENDIX
 REFERENCES
 
The dynamic behavior of the lung in health and disease depends on its viscoelastic properties. To better understand these properties, several mathematical models have been utilized by many investigators. In the present work, we present a new approach that characterizes the dynamics of gas flow into a viscoelastic porous medium that models the lung structure. This problem is considered in terms of the lung input impedance on a macro level and parenchymal tissue impedance on the level of an alveolar wall. We start from a basic theoretical analysis in which macroscopic tissue deformations are represented in accordance with the linearized Navier-Stokes equations. This approach has strong theoretical underpinnings in other situations but has not been applied to analyze the impedance of the inflated lung. Our analysis provides a theoretical basis for analyzing the interaction between flow into the lungs as a biophysical diffusion process and parenchymal viscoelasticity described phenomenologically, within the frameworks of standard viscoelasticity and structural damping. This lung impedance incorporates parameters of porosity, permeability, and viscoelasticity on micro and macro levels of parenchymal tissue. The analysis shows the theoretical basis of the transformation from the impedance of alveolar walls or isolated tissue strips to that of the intact parenchyma. We also show how the loading impedance at the lung boundary may have a significant impact on the dynamic behavior of whole lung viscoelasticity. Our analysis may be useful in directing specific tests of different models and for analyzing experimental measurements of viscoelastic parameters of lung material under normal and pathological conditions.

impedance; admittance; resistance; viscoelasticity; structural damping; tissue resistance


AN IMPORTANT MACROSCOPIC MANIFESTATION of the pathophysiological structural changes in pulmonary disease is the viscoelastic behavior of the lung. Classically, this experimentally derived behavior is studied as a time-dependent measure of pressure and flow (17). The most useful characteristic of these functions is impedance, which is simply the ratio of time-dependent pressure and flow expressed as a complex variable. Impedance then is a function of the frequency, whose magnitude and phase are dependent on the losses and stored energy of the physical substance. In a linear physical system, one can analyze specific models and theoretically calculate the different impedances. To have the most practical utility, the lung impedance should be describable from the intrinsic viscoelastic properties of the parenchymal tissue. The most commonly used model to characterize lung impedance is known as the constant-phase (CP) model (9, 12). This model takes into account the finding that the phase angle of the mechanical impedance is often constant over a wide range of frequencies (9, 12, 23). In addition to its use in lung physiology, the CP model has also been widely used in polymer physics, material mechanics, and geophysics, suggesting a broad utility to characterize material properties (11). The CP model for the lung was originally based on an empirical pressure response to step changes in gas volume as performed with rubber balloons and excised lung tissue (12). The empirical relationships thus derived involved an inverse power law of time with a power exponent less than one, and a Fourier transformation of this functional dependency results in the constant-phase impedance relationship.

The CP model, however, does not directly address the link between tissue properties and lung properties. Although isolated tissue strips and intact lungs both can be modeled with the CP model, the theoretical linkage between them has not been described. This would not be an issue if the lung really were like a rubber balloon, because then the value of the lung input impedance would be identical to the parenchymal tissue impedance (i.e., analogous to the balloon wall in this case, and ignoring airway resistance and gas compression). However, the fact that this is not the case is highlighted in a textbook scanning electron micrograph of a cut lung surface as shown in Fig. 1. The lung is thus far closer in appearance to a sponge than a child's balloon, and this anatomic complexity of a lung structure with such microscopic cavities needs to be considered. A second factor that only exists in the intact lung is the fact that the lung has a finite boundary condition dictated by the visceral pleura. Although this may not be a large constraint at low lung volumes, it may play an increasingly important role at higher inflations.


Figure 1
View larger version (226K):
[in this window]
[in a new window]
 
Fig. 1. Scanning electromicrograph of the lung parenchyma, highlighting its porouslike structure.

 
In the following analysis, we have attempted to theoretically provide the link between the impedance of intact lung parenchyma distal to terminal bronchioles (referred to here as the macro level) and that of the lung tissue at the level of an alveolar wall (micro level). To do this, we need a model that can dynamically handle the spongelike nature of the lung. One such model has considered the lung by modeling it as a porous medium (16). With this model, we have achieved a generalized mathematical description of the link between micromechanical and intact parenchymal impedances. This impedance analysis incorporates the process of flow diffusion in a porous lung medium and the influence of the pleural boundary condition (22) or loading impedance on this process. As will be shown, this theoretical basis can be used to characterize this linkage between several different lung models, including classical models with independent values of elastance and viscosity (2, 7) and models with structural damping that tightly couple elasticity and viscosity, as in the CP model (9, 12).


    OWEN-LEWIS EQUATIONS OF LUNG MECHANICS
 TOP
 ABSTRACT
 OWEN-LEWIS EQUATIONS OF LUNG...
 INPUT IMPEDANCE OF THE...
 ANALYSIS OF THE LUNG...
 DISCUSSION
 APPENDIX
 REFERENCES
 
Owen and Lewis considered the physical system of the lung as a fluid-solid system described by the Navier-Stokes equation, in terms of fluid velocity, pressure, fluid viscosity, Young's modulus (Ep), shear viscosity (µ1p), and density ({rho}s) as characterized by the stress-strain relationships for alveolar tissue (16). Utilizing the homogenization method of analysis (4), they considered the microscopic three-dimensional system to derive macroscopic equations, assuming the lung to be composed of a matrix of collagen, elastin fibers, and water. The micro scale length they used was that of an alveolus, l {cong} 2 · 10–4 m, and the macro scale was related to length of an acinus with L {cong} 10–2 m. Thus the macro scale includes all of the lung tissue distal to terminal bronchioles. Defining {epsilon} = l/L, and noting that {epsilon} << 1, they wrote functional equations for pressure, velocity of fluid, and displacement of tissue as a series of two scale presentations (16). Proceeding with further subsystem analysis and averaging the three-dimensional equations over the unit cell by volume of the alveolus, they were able to develop a fundamental system of equations for lung tissue dynamics. They proceeded to derive a one-dimensional system of equations in terms of the macroscopic physical variables of pressure, fluid velocity, and material displacement in any specific direction. It is worth emphasizing that these equations contain macroscopic viscoelastic parameters that derive from integration in a three-dimensional space. By analogy, this is not unlike the situation relating microscopic random Brownian motion of molecules in three dimensions to create a net pressure that can be measured with a pressure transducer oriented in one specific direction. Although the analysis given by Owen and Lewis was elegant, what was missing from it was an extension of the theory to develop practical equations for lung impedance.

To address this omission, we derive in the following analysis a solution for lung impedance in terms of the physical functions of pressure (p), fluid velocity (v), and solid displacement (u). The time dependency (t{infty} = –{infty}) of all functions can be written as exp (j{omega}t), where {omega} is angular frequency, with time derivatives in the form {partial}F/{partial}t = j{omega}F, where F represents any of functions p, v, or u. In the original analysis, the lung parenchyma was analyzed with a zero boundary condition for the pressure function. Although this situation may sometimes be approximated in a real lung at low stress, in general there is a finite load imposed by the tissue of the boundary. Thus we have extended the analysis to the more generalized situation including parameters of the boundary loading.

The system of equations presented by Owen and Lewis for a one-dimensional dynamic system (16) can be rewritten (ignoring the compressibility and inertance of air) in the following form:

Formula 1(1)

Formula 2(2)

Formula 3(3)
where j = (–1)1/2; the prime symbol indicates differentiation with respect to spatial coordinate x.

The other parameters of these equations are taken from Ref. 16 and shown in Table 1. Although Owen and Lewis described beta as the macroscopic manifestation of the composite microscopic elements (16), they considered it to be negligible for all conditions. However, this variable beta can be considered as a tissue admittance, and our analysis shows that, whereas it can be negligible for p, it cannot be ignored for v'. For this reason, we have added a component proportional to beta in our Eq. 2.


View this table:
[in this window]
[in a new window]
 
Table 1. Parameters of the 1-dimensional porous lung model

 
The system of Eqs. 13 has physical meaning with respect to any diffusion process in a porous medium. With specific reference to the lung, we note here that the original article (16) gave an estimation of the parameters on the basis of summaries of published data, with special emphasis on specific properties of the viscosity in two cases: standard viscoelasticity (SV) and structural damping (SD). Such characterizations were originally considered by Fredberg and Stamenovic (6), on the basis of the assumption that an integral biomechanical element is responsible for both the elastic and viscous stress in the tissue, and as a result the viscosity of the lung tissue is inversely proportional to the frequency. We have further developed and extended this concept in the framework of the parametric model of structural damping, incorporating the frequency dependence of the viscoelastic parameters of a porous lung material, which adds an additional degree of freedom (seeAPPENDIX).


    INPUT IMPEDANCE OF THE LUNG PARENCHYMA
 TOP
 ABSTRACT
 OWEN-LEWIS EQUATIONS OF LUNG...
 INPUT IMPEDANCE OF THE...
 ANALYSIS OF THE LUNG...
 DISCUSSION
 APPENDIX
 REFERENCES
 
Our objective was to transform the system of Eqs. 13 to equations for pressure (p) and flow (f). In the particular case of a one-dimensional system, the value of flow is coupled with fluid velocity as

Formula 4(4)
where S is cross-sectional area penetrated by the flow. Next we substitute Eq. 3 in the right side of the Eq. 2 and then insert this result into a differentiated Eq. 1, to arrive at the equation for the pressure:

Formula 5(5)

This equation results from integration by time of the respective diffusion equations. The coefficient D can be represented in the following form:

Formula 6(6)
where Formula 6= {alpha}2beta (Formula 6 + jFormula 6) = {alpha}2{sigma}({phi} + {alpha}), and with the quantitative values in Table 1 becomes {alpha}Formula 6{approx} 0.35.

Next we substitute the relationship p = D–1 · p'' in Eq. 3 and use this result in place of u' in the differentiated Eq. 1. After integrating this resultant equation (with a zero integration constant), we get an equation describing the functional dependency between p' and v. We obtain the functional dependency between v' and p by simply substituting Eq. 3 into Eq. 2. Finally, using Eq. 4, we arrive at the following two equations to replace Eqs. 13:

Formula 7(7)

Formula 7
where

Formula 8(8)

Formula 8
Here it is important note that zd · yd = D.

This system of Eq. 7 has a well known physical interpretation as equations of a long transmission line (Fig. 2), with specific longitudinal impedance density,

Formula 8
and cross-sectional admittance density,

Formula 8
These parameters give physical meaning to the propagation of the gas through lung parenchyma. This is somewhat analogous to the acoustical propagation of gas pressure perturbations, where the densities of the longitudinal impedance and cross-sectional admittance characterize the compressibility. However, in our situation, we are concerned only with the compressibility (deformation) of lung parenchymal material and not with a gas compressibility, because such effects are relatively minor over a frequency range of less than 20 Hz.


Figure 2
View larger version (7K):
[in this window]
[in a new window]
 
Fig. 2. Equivalent schematic of the diffusion process in a porous lung medium.

 
Figure 2 includes the load impedance, ZL, characterizing the boundary condition with ZL = pL/fL, where pL and fL are pressure and flow at the boundary of the lung, respectively. Note that flow, fL, is proportional to velocity, vL, of the fluid, and that vL exactly equals the velocity of the displacement of the membrane duL/dt at the periphery of the lung. The impedance, ZL, thus characterizes the viscoelastic impedance of the pleural membrane surrounding the lung. The solution of the system of Eq. 7 can be represented as follows:

Formula 9(9)

Formula 9
where

Formula 10(10)
is the propagation constant of the diffusion process and {gamma} and æ are the attenuation and phase constants of the diffusion process, respectively;

Formula 11(11)
is the wave impedance of the diffusion process.

The value of the wave impedance, Z0, plays a very important role in definition of lung impedance; in accordance with Eq. 8, this parameter can be represented as:

Formula 12(12)
where

Formula 13(13)
This formula for B can be simplified after substitution of the parameters from Table 1 into Eq. 12 as following:

Formula 14(14)
Using the values provided in Table 1, we have, in our particular case, B {approx} 2.9.

The next step in this simplification is to separate the value of the transformed parenchymal tissue impedance in formula Eq. 12 as

Formula 15(15)
and hence we have

Formula 16(16)
Equation 16 together with the load impedance ZL can be used to define the impedance of the whole lung. Considering Z as the ratio p0/f0 at x = 0, Eq. 9 yields the following formula for lung impedance,

Formula 17(17)
This formula represents the general relationship between lung input impedance and the material parameters of the lung.


    ANALYSIS OF THE LUNG INPUT IMPEDANCE
 TOP
 ABSTRACT
 OWEN-LEWIS EQUATIONS OF LUNG...
 INPUT IMPEDANCE OF THE...
 ANALYSIS OF THE LUNG...
 DISCUSSION
 APPENDIX
 REFERENCES
 
Equation 17 allows interpretation of the influence of the internal viscoelastic parameters of parenchymal tissue [through parameters Z0 and ({delta}L)] and loading impedance (ZL) on the overall lung impedance. The role of these parameters is dependent on values of the porosity ({phi}), Poisson ratio (Formula 17) of the lung tissue, and the ratio between macroscopic Young's modulus of the lung tissue (Ê) and Young's modulus of the alveolar wall (E) (through parameters {alpha} and {sigma}). The specific relations are shown in Eqs. 6, 8, 15, and 16 and Table 1.

We next estimate significance of the propagation parameter ({delta}L), which can be represented in accordance to Eqs. 6 and 10 as following:

Formula 18(18)
In general, this variable can have a significant effect on the impedance in both normal and pathological conditions. However, in the case when parameters are close to those represented in Table 1, we can use this to simplify the equation. In this case, we have {alpha}{delta}(KFormula 18)–1/2 {approx} 0.018{omega}1/2, and with the approximation of tanh({delta}L) {approx}{delta}L for frequencies less than 10 Hz, Eqs. 16 and 17 give us

Formula 19(19)
In the limiting case when ZL >> Zt we have from Eq. 15:

Formula 20(20)
showing that the lung input impedance equals the transformed parenchymal tissue impedance. We still have to take into account the fact that the transformation coefficient B depends on porosity and the ratio between Young's moduli measured on the macro and micro level. With the quantitative variables in Table 1 taken from Ref. 16, Eq. 20 gives, under standard viscoelasticity (SV) conditions,

Formula 21(21)
and under structural damping (SD) conditions,

Formula 22(22)
Equation 21 allows description of the dynamical behavior of the lung analogous to a classical lumped-parameter (RC) model with normalized parameters of resistance (R = 1) and capacitance (C = 1.6). Although Eq. 22 looks like a characterization of the CP model, this is not strictly true, because if the CP model exponent has a value of –1, this would automatically make the real part of the impedance equal zero. We also note that most experimental measurements find this power exponent of lung impedance to be in the range of 0.69–0.98. It is worth reiterating that Eq. 22 is obtained by substituting the structural damping model of viscosity used in Ref. 16 into Eq. 20. Furthermore, the time or frequency dependence of viscosity (Formula 22) and elasticity (Formula 22) from the structural damping model underlies the frequency dependence of impedance (G and H) in the CP model. Thus, to avoid potential contradictions between Eq. 22 and the CP model with measurements in actual physiological conditions, we offer a phenomenological representation of viscoelastic tissue properties in a framework of one variant of the parametric structural damping model (see APPENDIX). In accordance to this model along with Eqs. A.9A.13 we can rewrite Eq. 15 as follows:

Formula 23(23)
where Rt = BgE{rho}SLS–1 is transformed parenchyma tissue resistance.

Finally, to complete the solution, we still have to include the loading impedance ZL. Initially, we described this parameter as the impedance of the membrane material surrounding the lung parenchyma. Viscoelastic properties of this membrane may depend on several factors but are likely dominated by the mechanical properties of the visceral pleura. There have been a number of studies that have attempted to quantify the dynamic properties of this pleural boundary (3, 12, 20, 22). Simultaneous measurement of lung and pleural impedances showed significant differences between their frequency characteristics (22). This fact supports the use of ZL, in accordance with the parametric model of structural damping, and by analogy with Eq. 23 in the following form:

Formula 24(24)
where RL is the loading resistance.

We can further interpret Eq. 19 using the relationship between the tissue and loading impedances. This formula can be further simplified over the frequency range less than 10 Hz, where terms containing ({delta}L)2 can be neglected. Then, on the basis of Eqs. 19, 23, and 24, we have:

Formula 24
where

Formula 25(25)
Equation 25 demonstrates that, for a normal lung, input impedance can be described as the parallel connection of parenchymal and loading impedances. Certain caveats apply, however, to this consideration, one being that we have not included any discrete upstream airway resistance, and the other that we do not include any heterogeneity of tissue properties. To the extent that airways smaller than terminal bronchioles are included in the porous matrix, their effect will be included, but the more conventional frictional airway resistance of the larger airways is not. It would surely be possible to add a series resistance to an expanded model without much additional complexity, but including the effect of parenchymal heterogeneity would be considerably more cumbersome. Equation 25 also shows the important role of {zeta} on the frequency dependence of lung input impedance. For example, when {zeta} << 1 (i.e., {alpha}1 > {alpha}0), the value of the impedance equals the transformed parenchymal tissue impedance with a frequency dependence proportional to {omega}Formula 25, at least when {omega} is less than ~20 Hz (i.e., minimal gas compressibility). At the other extreme, when {zeta} >> 1 (i.e., {alpha}1 < {alpha}0), the value of the impedance equals the pleural impedance with a frequency dependence proportional to {omega}Formula 25. This analysis provides a theoretical explanation for the experimental observations comparing the dynamic properties of lung parenchymal tissue in vitro with those in the whole intact lung (8). This work showed a greater power of inverse frequency dependence of lung tissue damping compared with that for the inverse frequency dependence of the whole lung.

To summarize, we can now present a formula for lung impedance over an arbitrary frequency range and with arbitrary parameters. Combining Eqs. 16, 17, 23, and 24 yields:

Formula 26(26)
where {zeta}0 = tanh({delta}L)/({delta}L).

This impedance formula contains relationships with viscoelastic parameters of the parenchyma and loading materials on the micro and macro levels, also incorporating parameters of the porosity and permeability. This equation thus characterizes the link between the impedance and the structural and material parameters of the lung.


    DISCUSSION
 TOP
 ABSTRACT
 OWEN-LEWIS EQUATIONS OF LUNG...
 INPUT IMPEDANCE OF THE...
 ANALYSIS OF THE LUNG...
 DISCUSSION
 APPENDIX
 REFERENCES
 
In this paper, we developed an approach of the analysis of dynamic lung biomechanics based on consideration of the Navier-Stokes equations for gas transport, in conjunction with the viscoelastic equations for porous material at the micro and macro level (4, 14, 16, 18, 24). A solution to this problem based on homogenization theory (4) was published a few years ago by Owen and Lewis for the mechanical system of the lung (16), but, unfortunately, this work stopped short of the analysis needed for application to interpreting the lung impedance. The analysis we present now provides the theoretical basis of both the structural damping and standard viscoelasticity models as applied to the lung. We also show how the tissue impedance can be related to the intact lung.

There have been several attempts to connect macroscopic rheology of the parenchymal tissue with the fibrous microstructure of the cell on the biophysical level (1, 15) and also with the molecular theory of polymer viscoelasticity (19). It was shown that either of these approaches can describe viscoelastic deformations of the parenchymal tissue that follow predictions of the CP model. For this purpose, there are at least two ways to characterize the parenchymal tissue impedance: with either a statistical model or a phenomenological description of viscoelastic tissue behavior. In our approach, because we analyzed lung impedance using a porous structure based on deterministic considerations, it was more appropriate to utilize the phenomenological model for the parenchymal tissue impedance. This approach also enabled either of the two common phenomenological models (the CP or RC) to fit with the analysis. In accordance with published descriptions of these models (6, 16), they can be characterized with elements defined either by standard viscoelasticity (SV for the RC model) or structural damping (SD for the CP model). The SV framework uses constant parameters, and SD framework uses time-variable or parametric elements. One parametric model of the viscoelastic deformation was considered within the constraints of quasi-linear viscoelasticity (7) and was shown to have characteristics approximating the CP model. For the purpose of the jointly considering lung input and tissue impedances, we showed how the same mathematical expressions for the parametric model of structural damping could be applied to either case, as arising from the micromechanical deformation. This analysis supports the work of Suki et al. (19) that originally probed the molecular mechanism underlying the CP model. However, more recent analogies of the time variability of parenchymal viscoelasticity to glasses and polymers can also be fit well with our present analytical framework. For example, Cavaille et al. (5) have shown that the compliance of polymers contains four components, among which are two associated with correlated motions of molecules with time-constant parameters as power functions of time. In our model of structural damping we use an identical time variable function for the phenomenological compliance.

One further consideration with regard to our analysis is that we have only considered the whole lung as a homogeneous material, ignoring the potential for parallel regional variations the can lead to heterogeneous behavior. This is an added level of complexity, but if enough information is available on the nature of the heterogeneity, it would surely be possible to use the same analytical model to describe these parallel components. In this case the lung impedance would be considered as an ensemble of the similar impedances with different parameters and in different one-dimensional directions. This analysis, however, is beyond the scope of the present work.

In conclusion, we have applied and extended the theoretical solution describing the dynamics of a porous medium, with parenchyma viscoelasticity presented by Owen and Lewis (16), to derive a general formula for the impedance of the lung. This result utilizing a variant of the parametric model of structural damping provides an explicit relationship between impedance and parameters of viscoelasticity on micro (alveolar wall) and macro (acinar) levels. We have also shown that the dynamic impedance of the lung in the general case is tightly coupled with parenchymal tissue impedance and with loading impedance. In situations in which loading impedance is significantly greater than tissue impedance, the lung input impedance is then proportional to the parenchymal tissue impedance. To describe the CP model of parenchymal viscoelasticity, we used one variant of the parametric model of structural damping. With this approach, our results allow further analysis of lung impedance with viscoelastic tissue parameters at any level of complexity. This analysis may be useful in directing specific tests of different models used to explain experimental results designed to measure viscoelastic parameters of lung structure under normal and pathological conditions.


    APPENDIX
 TOP
 ABSTRACT
 OWEN-LEWIS EQUATIONS OF LUNG...
 INPUT IMPEDANCE OF THE...
 ANALYSIS OF THE LUNG...
 DISCUSSION
 APPENDIX
 REFERENCES
 
One variant of the parametric model of structural damping.   The concept of structural damping was first suggested in the physiological literature by Fredberg and Stamenovic (6). This term is used to describe the phenomenon of frequency independence of the elastic and viscous parameters of biological material (e.g., parameters Ep and µSDp in Table 1). Under conditions of structural damping, the lung impedance can be interpreted within the framework of the CP model. In this model, identical frequency behavior of both parameters can exist with the assumption that one united physical element is responsible for elasticity and viscosity (6). However, their analysis did not describe characteristics of this element that would be consistent with power law of the CP model. It turns out that for this physical element to fit the CP model, one needs to incorporate a slight frequency dependency of the elastic and viscous characteristic of the biological material. On the basis of these considerations and as motivated in the DISCUSSION, we suggest an alternative development of the structural damping concept. This model is analyzed by use of time variable viscosity and elasticity. Although these parameters are proportional to each other on the micro and macro levels (as linked by the constant {sigma} in Table 1), for our purposes it is more convenient to work with the macroparameters of the elasticity Formula 26p and viscosity Formula 26p. If desired, information at the micro level can be easily obtained from the known constant, {sigma}.

To further clarify our concept of the parametric model of structural damping, let us consider the simpler situation in one dimension of a parenchymal strip in vitro with cross-sectional area ST, length LT, and internal force {tau}(t). The linear deformation u(t) of this material will lead to a strain of u(t)/LT and stress of {tau}(t)/ST. Mechanical properties of the material are represented with elasticity Formula 26p and viscosity Formula 26p. To describe the dynamic behavior of standard viscoelasticity in this material, we start from the Maxwell model, describing the relationship between values of the velocity of the linear deformation u(t), i(t) = du/dt, and {tau}(t), the internal force in the material, in the following form:

Formula 26
This formula shows that conventional deformation is described with two separate elements of viscosity and elasticity. With a parametric model of structural damping, we assume the presence of just one elastic element with the following time dependence of capacitance, c(t) = Formula 26p–1(t):

Formula A1(A1)
where Cm–1 is a constant with units of Young's modulus; {tau}{epsilon} ≤ t, k < 1, and {xi} has units of s–1.

We include in this consideration a function of the volume, Q(t):

Formula A1
We then have:

Formula A2(A2)
Equation A2 is a differential equation with time-varying coefficients. Its solution is well known and can be written with help of the following function {varphi}(t,{tau}):

Formula A3(A3)
Using this function, we have the following solution of Eq. A2:

Formula A4(A4)
where the initial tension {tau}0 = 0. Then, from Eqs. A3 and A4, and c({tau}) from Eq. A1, we have:

Formula A5((A5))
We are a seeking specific solution of Eq. A2 under synchronized conditions between the momentum ({tau}) at the onset of the pulse function and the momentum associated with the development of the parametric effect. In this case, the current {tau}c(t) will be considered as a limit of {tau}(t) as {tau}{epsilon} -> {tau}. That is, we have at t0 = 0:

Formula A6(A6)
This result describes a convolution relation between {tau}c(t) and i(t). Taking the Fourier transform of both sides of Eq. A6 gives the following:

Formula A7(A7)
where {omega} is angular frequency; {tau}c({omega}) and i({omega}) are Fourier transforms of {tau}c(t) and i(t), respectively. We use the subscript ts to indicate the special situation of an isolated tissue strip; Zts({omega}) is the Fourier transform of the kernel of Eq. A5:

Formula A8(A8)
where {alpha}0 = 1 – k; {Gamma}({alpha}0) is the gamma function.

In accordance with the theory of the CP model (9, 10), the function Zts({omega}) can be considered as

Formula A9(A9)
where

Formula A9
;

Formula A9
. The function Zts({omega}) is the tissue strip impedance, and the viscoelastic parameters, Formula A9p and Formula A9p, are coupled with this function through the following relationship:

Formula A10(A10)
where

Formula A11(A11)
We can represent Eq. A10 in the following form:

Formula A12(A12)
In the context of our model and in accordance with Eq. 14 in the text, we have the following normalized values in Eq. A12:

Formula A13(A13)

Formula A13
where

Formula A13
The relationships in Eq. A13 are the general equations of the parametric model of structural damping.

Glossary

B
transformation coefficient in the tissue impedance relationship

c(t)
capacitance of the parametric model of the structural damping

Cm
constant quasi-capacitance of the parametric model of the structural damping

D
inversely proportional value to the coefficient diffusion

E
normalized Young's modulus of parenchymal tissue

Ep
physical Young's modulus of parenchymal tissue

Ê
macroscopic Young's modulus

Formula A13
normalized macroscopic Young's modulus of the porous parenchyma material

Formula A13p
physical macroscopic Young's modulus of the porous parenchyma material

Formula A13p(t)
time variable Young's modulus in the parametric model of structural damping

f
fluid flow

fL
fluid flow at the periphery of the lung

gE
normalized equivalent viscosity of the porous lung material

hE
normalized equivalent elasticity of the porous parenchymal material

i(t)
velocity of the linear deformation

j
imagine unit

k
exponent of power law characterizing the parametric capacitance change

K
normalized permeability, proportional to the conductivity of the flow

Kp
physical permeability of the flow

l
length of the alveolar sac

L
length of the acinus

Lt
length of the tissue strip

p
physical function of pressure

pE
equivalent elasticity of the porous lung material

pL
pressure at the periphery of the lung

qE
equivalent viscosity of the porous lung material

Q(t)
function of the quasi-volume in the parametric model of the structural damping

Rt
transformed parenchymal tissue resistance

RL
loading resistance

S
cross-sectional area for flow

ST
cross-sectional area of the tissue strip

t
time coordinate

u
physical function of solid displacement

uL
solid displacement at the lung periphery

v
physical function of fluid velocity

vL
solid velocity at the lung periphery

x
spatial coordinate

yd
cross-sectional admittance density

zd
longitudinal impedance density

z
coefficient of the wave impedance under standard visco-elasticity conditions

zt
normalized lung parenchymal tissue impedance

Z0
wave impedance

ZL
loading impedance

Zt
transformed parenchymal tissue impedance

Zts
mechanical impedance of the tissue strip

{alpha}
transformation coefficient of the elastic parameters of the lung material from micro level to macro level

{alpha}{delta}
modified transformation coefficient

{alpha}0
power of frequency characterizing the parenchymal tissue impedance

{alpha}1
power of frequency characterizing the loading impedance

beta
transformed lung tissue admittance from micro level to macro level

{delta}
propagation constant of the diffusion process

{epsilon}
ratio between the alveolar sac and acinus lengths

{phi}
porosity of the lung structure

{gamma}
attenuation constant of the diffusion process

{Gamma}(...)
gamma function

{eta}
time variable of the integration

æ
wave number of the diffusion process

{varphi}(t,{tau})
multiplier in the Green function of the parametric model

Formula A13
normalized viscosity of the lung porous material

µ1
normalized shear viscosity of the alveolar wall

µ1p
physical shear viscosity of the alveolar wall

µSD
normalized viscosity for structural damping

µSDp
physical viscosity for structural damping

Formula A13
macroscopic Poisson ratio

{rho}s
density of the alveolar tissue

{sigma}
transformation coefficient of the viscoelastic complex from micro level to macro level

{tau}
time variable of integration of the parametric equation

{tau}{epsilon}
time parameter of the quasi-capacitance time function

{tau}(t)
internal force in the material of the parametric model of the structural damping

{omega}
angular frequency

{xi}
normalizing parameter of the quasi-capacitance time function

{zeta}
impedance ratio between parenchymal tissue and pleura materials


    FOOTNOTES
 

Address for reprint requests and other correspondence: W. Mitzner, Division of Physiology, Dept. of Environmental Health Sciences, Johns Hopkins Univ., Bloomberg School of Public Health, 615 North Wolfe St., Baltimore, MD 21205 (e-mail: wmitzner{at}jhsph.edu)

The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.


    REFERENCES
 TOP
 ABSTRACT
 OWEN-LEWIS EQUATIONS OF LUNG...
 INPUT IMPEDANCE OF THE...
 ANALYSIS OF THE LUNG...
 DISCUSSION
 APPENDIX
 REFERENCES
 

  1. Bates JHT. A micromechanical model of lung tissue rheology. Ann Biomed Eng 26: 679–687, 1998.[CrossRef][ISI][Medline]
  2. Bates JHT, Brown KA, and Kochi T. Respiratory mechanics in the normal dog determined by expiratory flow interruption. J Appl Physiol 67: 2276–2285, 1989.[Abstract/Free Full Text]
  3. Bense L, Hök K, Gabriel S, and Wiman LG. Effect of fibrin sealant on the tensile strength of rat visceral pleura. Scand J Thorac Cardiovasc Surg 27: 179–180, 1987.
  4. Burridge R and Keller JB. Poroelasticity equations derived from microstructure. J Acoust Soc Am 70: 1140–1146, 1981.[CrossRef]
  5. Cavaille JY, Perez J, and Johari GP. Molecular theory for the rheology of glasses and polymers. Phys Rev B Condens Matter 39: 2411–2422, 1989.
  6. Fredberg JJ and Stamenovic D. On the imperfect elasticity of lung tissue. J Appl Physiol 67: 2408–2419, 1989.[Abstract/Free Full Text]
  7. Fung YC. Biomechanics: Mechanical Properties of Living Tissues. New York: Springer-Verlag, 1981.
  8. Fust A, Bates JH, and Ludwig MS. Mechanical properties of mouse distal lung: in vivo versus in vitro comparison. Respir Physiol Neurobiol 143: 77–86, 2004.[CrossRef][ISI][Medline]
  9. Hantos Z, Daroczy B, Suki B, Nagy S, and Fredberg JJ. Input impedance and peripheral inhomogeneity of dog lungs. J Appl Physiol 72: 168–178, 1992.[Abstract/Free Full Text]
  10. Hantos Z, Suki B, Czendes T, and Daróczy B. Constant-phase modeling of pulmonary tissue impedance (Abstract). Bull Eur Physiopath Respir 23, Suppl 12: 326s, 1987.
  11. Heymans N. Fractional calculus description of non-linear viscoelastic behavior of polymers. Nonlinear Dynamics 38: 221–231, 2004.[CrossRef]
  12. Hildebrandt J. Comparison of mathematical models for cat lung and viscoelastic balloon derived by Laplace transform methods from pressure-volume data. Bull Math Biophys 31: 651–667, 1969.[ISI][Medline]
  13. Humphrey JD. A possible role of the pleura in lung mechanics. J Biomech 20: 773–777, 1987.[CrossRef][ISI][Medline]
  14. Levy T. Propagation of waves in fluid-saturated porous elastic solid. Int J Eng Sci 17: 1005–1014, 1979.[CrossRef]
  15. Mijailovich SM, Stamenovic D, and Fredberg JJ. Toward a kinetic theory of connective tissue micromechanic. J Appl Physiol 74: 665–681, 1993.[Abstract/Free Full Text]
  16. Owen MR and Lewis MA. The mechanics of lung tissue under high-frequency ventilation. SIAM Appl Math 61: 1731–1761, 2001.[CrossRef]
  17. Peslin R and Fredberg JJ. Oscillation mechanics of the respiratory system. In: Handbook of Physiology. The Respiratory System. Mechanics of Breathing. Bethesda, MD: Am. Physiol. Soc., 1986, sect. 3, vol. III, pt. 1, chapt. 11, p. 145–177.
  18. Singh B. Reflection of plane micropolar viscoelasticity waves at a loosely bonded solid-solid interface. Sadhana Acad P Eng S 27: 493–506, 2002.
  19. Suki B, Barabasi AL, and Lutchen KR. Lung tissue viscoelasticity: a mathematical framework and its molecular basis. J Appl Physiol 76: 2749–2759, 1994.[Abstract/Free Full Text]
  20. Suki B, Daróczy B, and Hantos Z. Mechanical impedance of the canine diaphragm. Part 1. Experimental system and measurements. Med Biol Eng Comput 28: 361–366, 1990.[CrossRef][ISI][Medline]
  21. Suki B and Hantos Z. Viscoelastic properties of the visceral pleura and its contribution to lung impedance. Respir Physiol 90: 271–287, 1992.[CrossRef][ISI][Medline]
  22. Thamrin C, Janosi TZ, Collins RA, Sly PD, and Hantos Z. Sensitivity analysis of respiratory parameter estimates in the constant-phase model. Ann Biomed Eng 32: 815–822, 2004.[CrossRef][ISI][Medline]
  23. Verdier C. Rheological properties of living materials, from cells to tissues. J Theor Med 5: 67–91, 2003.



This article has been cited by other articles:


Home page
J. Appl. Physiol.Home page
W. Mitzner
Counterpoint: Lung Impedance Measurements are not More Useful than Simpler Measurements of Lung Function in Animal Models of Pulmonary Disease
J Appl Physiol, November 1, 2007; 103(5): 1901 - 1903.
[Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow All Versions of this Article:
101/3/926    most recent
01548.2005v1
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 ISI Web of Science
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 ISI Web of Science (1)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Lande, B.
Right arrow Articles by Mitzner, W.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Lande, B.
Right arrow Articles by Mitzner, W.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2006 by the American Physiological Society.