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


     


J Appl Physiol 86: 273-284, 1999;
8750-7587/99 $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 Mates vanLöbensels, E.
Right arrow Articles by Hlastala, M. P.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Mates vanLöbensels, E.
Right arrow Articles by Hlastala, M. P.
Vol. 86, Issue 1, 273-284, January 1999

Modeling diffusion limitation of gas exchange in lungs containing perfluorocarbon

Elisabeth Mates vanLöbensels1, Joseph C. Anderson2, Jacob Hildebrandt1, and Michael P. Hlastala1

Departments of 1 Physiology and Biophysics and 2 Chemical Engineering, University of Washington, Seattle, Washington 98195

    ABSTRACT
Top
Abstract
Introduction
References

We reported changes in alveolar-arterial PO2 gradient, ventilation-perfusion heterogeneity, and arterial-alveolar PCO2 gradient during partial liquid ventilation (PLV) in healthy piglets (E. A. Mates, P. Tarczy-Hornoch, J. Hildebrandt, J. C. Jackson, and M. P. Hlastala. In: Oxygen Transport to Tissue XVII, edited by C. Ince. New York: Plenum, 1996, vol. 388, p. 585-597). Here we develop two mathematical models to predict transient and steady-state (SS) gas exchange conditions during PLV and to estimate the contribution of diffusion limitation to SS arterial-alveolar differences. In the simplest model, perfluorocarbon is represented as a uniform flat stirred layer and, in a more complex model, as an unstirred spherical layer in a ventilated terminal alveolar sac. Time-dependent solutions of both models show that SS is established for various inert and respiratory gases within 5-150 s. In fluid-filled unventilated terminal units, all times to SS increased sometimes by hours, e.g., SF6 exceeded 4 h. SS solutions for the ventilated spherical model predicted minor end-capillary disequilibrium of inert gases and significant disequilibrium of respiratory gases, which could explain a large portion of the arterial-alveolar PCO2 gradient measured during PLV (14). We conclude that, during PLV, diffusion gradients for gases are generally small, except for CO2.

liquid breathing; perfluorocarbon liquids; mathematical model; gas exchange

    INTRODUCTION
Top
Abstract
Introduction
References

PARTIAL LIQUID VENTILATION (PLV) is a technique of ventilatory support in which the air spaces of the lung are partially replaced with liquid perfluorocarbon (PFC) and then periodically insufflated with O2-enriched gas with use of a conventional mechanical ventilator. PLV was first described by Fuhrman et al. (3) and has been shown to improve oxygenation and lung mechanics in animal models and in humans with acute respiratory distress syndrome (2, 4, 8, 10, 11). We previously showed that PLV in healthy piglets causes mild increases in arterial-alveolar PO2 and PCO2 gradients [(A-a)DO2 and (a-A)DCO2] (12-14). Compared with conventional gas ventilation with 100% O2, there was a 50% increase in ventilation-perfusion (VA/Q) heterogeneity and a 50% increase in O2 shunt, both of which can contribute to the alveolar-arterial difference. We hypothesized, but were unable to verify experimentally, that a diffusion barrier exists across the PFC in the lung periphery and that it is responsible for a significant portion of measured increases in alveolar-arterial differences in healthy animals during PLV. To test the feasibility of this hypothesis, we developed two mathematical models of gaseous diffusion in partially PFC-filled lung subunits.

In our experimental studies we used the multiple inert gas elimination technique (MIGET) to measure VA/Q heterogeneity in healthy piglets during PLV (13, 14, 20). The use of this method raised the question of whether inert and respiratory gas exchange reaches steady state during PLV within a time frame similar to conventional gas ventilation. Steady state refers to the condition in which, given a constant source of a gas infused into mixed venous blood, the ratio of input to output partial pressures across the lung (i.e., Pa/P<SUB><OVL>v</OVL></SUB> and PE/P<SUB><OVL>v</OVL></SUB>, where Pa, P<SUB><OVL>v</OVL></SUB>, and PE are arterial, mixed venous, and expired pressures, respectively) does not change with time and there is no further storage or net loss of mass within the lung over time. Using a very simple model, we showed previously that the time to steady state for SF6 (a gas used in MIGET to estimate shunt) was prohibitively long because of its high solubility in PFC vs. blood (13). This required us to modify MIGET by eliminating SF6 from the analysis, inasmuch as it did not satisfy the underlying assumption that steady-state conditions exist (13). With the more sophisticated models described here, we were able to refine and verify these original predictions and further explore the effects of PFC on attainment of steady-state gas exchange for the remaining five inert gases as well as O2 and CO2. We are also able to explore the effect of PFC dose on diffusion-limited gas transport in the alveolus.

In recent publications, PLV has been shown to improve gas exchange in humans with acute lung injury (4, 8, 11). We have focused our efforts on studying the effects of PLV in healthy animals to shed light on the fundamental differences in gas exchange between gas- and liquid-filled lungs. Many of the equations in traditional gas exchange theory are based on the assumptions that steady-state mass flux exists and that there is a negligible diffusion barrier in the alveolus (e.g., Berggren shunt and Bohr dead space). These assumptions need to be critically evaluated in the novel situation of a fluid-filled lung. Despite mild increases in (A-a)DO2 and (a-A)DCO2 during PLV in healthy animals, oxygenation and ventilation can be achieved surprisingly well through a liquid-filled lung. The success of PLV in a clinical setting may depend on altering our thinking about shunt and dead space when we add a high-solubility fluid with diffusion resistance to the air space of the lung. The models described here have been helpful in exploring these ideas.

    MATHEMATICAL MODELS

In prior publications we presented two different models of gas exchange during PLV: 1) a two-compartment well-mixed model used to estimate times to steady state (13) and 2) a spherical gas exchange model used to estimate steady state arterial-alveolar differences across a PFC diffusion barrier (14). Here we expand on both models, adding a gas compartment to the one-dimensional well-mixed model, providing time- and space-dependent numerical solutions to the spherical model, and providing a full discussion of the underlying assumptions and model behavior. We explore solutions to the time rate of change of partial pressures of O2, CO2, and six MIGET gases in PFC after a step change in input partial pressures.

A comparison of two separate model configurations is particularly enlightening, since the in vivo PFC-filled alveolus probably includes some features of both. The well-stirred compartment model reflects a PFC layer with complete convective mixing and no diffusion limitation within the gas exchange unit, whereas the spherical shell model imitates a perfectly still diffusion barrier interposed between gas and blood. The true nature of gas exchange in PFC lies somewhere between these models. With each breath, PFC probably moves in and out of some alveoli and small airways and exists as small stagnant puddles in others.

Glossary

Solubility of a tracer gas in a solvent (ml gas · 100 ml solvent-1 · Torr-1)

C

Concentration of a tracer gas in a solvent (ml gas/ml solvent)

D

Molecular diffusion coefficient (cm2/s)

M

Mass of tracer gas in a solvent (ml gas)

MIGET

Multiple inert gas elimination technique

n

Number of gas exchange units in a piglet lung

P

Partial pressure of a tracer gas (Torr)

PFC

Perfluorochemical

PLV

Partial liquid ventilation

Q

Blood flow (ml/s)

RR

Respiratory rate (min-1)

r

Radial distance from center of gas compartment (cm)

rc

Radius of gas exchange unit at the capillary boundary (cm)

rg

Radius of gas compartment (cm)

t

Time (s)

tau

Time constant (s)

T

Temperature (°C or K)

T98

Time to 98% of steady state (s)

VD

Dead space (ml/breath)

VT

Tidal volume (ml/breath)

a

Arterial

A

Alveolar

b

Blood

c

Capillary

g

Alveolar gas

gi

Inspired gas

pfc

Perfluorocarbon

<OVL>v</OVL>

Mixed venous

Model assumptions. In both models we assumed that the blood and gas compartments on either side of the PFC are well mixed. The models also assume that diffusion barriers at the capillary membrane and the PFC-gas interface are negligible. Because the presence of the PFC in the alveolar space does not affect gas exchange properties of the alveolar-capillary membrane, the assumption of complete equilibrium across the membrane is as valid as in the gas-filled lung. Blood flow and ventilation are assumed continuous and nonpulsatile (i.e., Q and VA are constant).

P<SUB><OVL>v</OVL></SUB> of tracer gases was assumed to be constant, and variation in body tissue partial pressures was assumed to be negligible. In the experimental situation, P<SUB><OVL>v</OVL></SUB> of the inert gases will vary slightly with time as the body comes to a new steady state after a perturbation in gas exchange. We believed that this variation was small, inasmuch as the body tissues were previously equilibrated with inert gas and the recirculated component is a small fraction of the total P<SUB><OVL>v</OVL></SUB>. The error introduced by this assumption will lead to a slight underestimation of the true time to reach steady state.

Time-dependent gas exchange in a well-stirred three-compartment model. Figure 1 schematically describes this model, in which blood is delivered to the capillary compartment at a flow rate Q (ml/s) and ventilation through the gas compartment occurs at a rate VA (ml/s). A tracer gas may enter the gas exchange unit dissolved in blood at partial pressure P<SUB><OVL>v</OVL></SUB> or via ventilation at partial pressure Pgi. Mass balance for the tracer in three compartments is given by Eq. 1 with the assumption that the PFC layer is well mixed. Thus the gas partial pressures in all compartments are equal (i.e., Pc = Ppfc = Pg) at time t
<FR><NU>d(<IT>M</IT><SUB>c</SUB> + <IT>M</IT><SUB>pfc</SUB> + <IT>M</IT><SUB>g</SUB>)</NU><DE>d<IT>t</IT></DE></FR> = <A><AC>Q</AC><AC>˙</AC></A>⋅(C<SUB><OVL>v</OVL></SUB> − C<SUB>c</SUB>) − <A><AC>V</AC><AC>˙</AC></A><SC>a</SC>⋅(C<SUB>g</SUB> − C<SUB>gi</SUB>) (1)
Converting mass (M = C · V) and concentration (C = beta · P) to partial pressures (P), applying the assumption Pc = Ppfc = Pg (i.e., well-mixed with no diffusion gradients), and rearranging into the standard form for a first-order differential equation tau  · P + P = K (where tau  is the time constant and K is the steady-state asymptotic value of P)
<FR><NU>V<SUB>c</SUB>⋅&bgr;<SUB>b</SUB> + V<SUB>pfc</SUB>⋅&bgr;<SUB>pfc</SUB> + V<SUB>g</SUB>⋅&bgr;<SUB>g</SUB></NU><DE><A><AC>V</AC><AC>˙</AC></A><SC>a</SC>⋅&bgr;<SUB>g</SUB>+ <A><AC>Q</AC><AC>˙</AC></A>⋅&bgr;<SUB>b</SUB></DE></FR> ⋅ <FR><NU>dP<SUB>pfc</SUB></NU><DE>d<IT>t</IT></DE></FR>
 + P<SUB>pfc</SUB> = <FR><NU><A><AC>V</AC><AC>˙</AC></A><SC>a</SC>⋅&bgr;<SUB>g</SUB>⋅P<SUB>gi</SUB> + <A><AC>Q</AC><AC>˙</AC></A>⋅&bgr;<SUB>b</SUB>⋅P<SUB><OVL>v</OVL></SUB></NU><DE><A><AC>V</AC><AC>˙</AC></A><SC>a</SC>⋅&bgr;<SUB>g</SUB> + <A><AC>Q</AC><AC>˙</AC></A>⋅&bgr;<SUB>b</SUB></DE></FR> (2)
tau  can be expressed as follows
&tgr; = <FR><NU>V<SUB>c</SUB> + V<SUB>pfc</SUB>⋅<FR><NU>&bgr;<SUB>pfc</SUB></NU><DE>&bgr;<SUB>b</SUB></DE></FR> + V<SUB>g</SUB>⋅<FR><NU>&bgr;<SUB>g</SUB></NU><DE>&bgr;<SUB>b</SUB></DE></FR></NU><DE><A><AC>V</AC><AC>˙</AC></A><SC>a</SC>⋅<FR><NU>&bgr;<SUB>g</SUB></NU><DE>&bgr;<SUB>b</SUB></DE></FR> + <A><AC>Q</AC><AC>˙</AC></A></DE></FR> (3)
The standard solution to Eq. 2 takes the form
P = <IT>K</IT>⋅(1 − <IT>e</IT><SUP>−<IT>t</IT>/&tgr;</SUP>) (3b)
where K = [VA · Pgi · (beta g/beta b) + Q · P<SUB><OVL>v</OVL></SUB>]/[VA · (beta g/beta b) + Q].


View larger version (20K):
[in this window]
[in a new window]
 
Fig. 1.   Schematic of 3-compartment model with well-mixed blood, PFC, and alveolar gas compartments. Capillary compartment is perfused at a rate Q at partial pressure P<SUB><OVL>v</OVL></SUB>. Capillary, PFC, and gas compartment partial pressures vary with time but are uniformly mixed and in spatial equilibrium at each time t. Partial pressure vs. distance along capillary is illustrated for a point in time. See Glossary for definition of abbreviations.

The rate at which Ppfc approaches steady-state equilibrium is determined by tau , the time for the exponential term to decrease by 63%. At 4tau , steady-state equilibrium is >98% complete. The standard MIGET theory assumes that Pc = Pg = constant; i.e., after a change in the infusate, the time at which gas exchange measurements are taken is much longer than tau , so the exponential term in Eq. 3 becomes negligible.

Equation 3 shows that when PFC is present in the alveolus and beta pfc > beta b, tau  is prolonged, especially if beta pfc is greater than both beta b and beta g. For gases in which this holds true, larger volumes of PFC result in longer times to equilibrium. For O2, tau  is actually prolonged in the absence of PFC, because beta g > beta pfc. It is also prolonged as VA approaches zero (i.e., shunt), because PFC must equilibrate to a higher final value, i.e., input partial pressures P<SUB><OVL>v</OVL></SUB> or Pgi. When VA is nonzero, the steady-state partial pressure (K) is less than input partial pressure and tau  is accordingly shorter. Increasing VA or Q shortens tau  for all gases.

Time-dependent gas exchange in a spherical shell with radial diffusion. To simulate gas exchange in a functional subunit of lung (Fig. 2), we chose a spherically shaped structure with an outer layer of capillary blood surrounding a layer of PFC that, in turn, surrounds a gas-filled center. The branching, space-filling nature of lung architecture is too complex for small-scale mathematical modeling. We chose to model gas exchange at the level of the terminal alveolar duct and represented them as smooth spheres. If the anatomic subunit is larger than this, the surface area of a smooth sphere would greatly underestimate the surface-to-volume ratio. On the other hand, representing a structure as small as an alveolus by a closed sphere would overestimate the ratio, since alveoli are roughly hexagonal cups. We therefore compromised on a structure the size of a single terminal alveolar sac to be portrayed by a sphere with dimensions derived accordingly.


View larger version (27K):
[in this window]
[in a new window]
 
Fig. 2.   Schematic of spherical shell model representing a terminal sac partially filled with PFC. Size of unit approximates a terminal alveolar sac. Gases enter model via mixed venous blood at pressure P<SUB><OVL>v</OVL></SUB> in capillary compartment or through inspired gas at pressure Pgi. Gas diffuses radially (r) through PFC and is removed from system by ventilation at pressure Pg or by blood flow at pressure Pc. PFC-gas and capillary-PFC interfaces are located at r = rg and r = rc, respectively. Only radial gradients in gas concentration occur. Differential volume element for analysis is a spherical shell of thickness dr. See Glossary for definition of abbreviations.

We assumed that the capillary and alveolar gas compartments were individually well mixed and that uniform radial diffusion occurred in the PFC. Mass exchange between the compartments is dependent on the interfacial area bounding two adjacent regions. The area of the capillary-PFC interface is fixed at 4pi r2c. The area of the inner gas space (4pi r2g) depends on the volume of PFC administered and on total lung volume. PFC is assumed to distribute uniformly as a spherical shell with the ventilated gas "hole" in the center. As the hole radius approaches zero, the unit becomes "flooded" with PFC. As rg approaches rc the model represents a gas-filled lung with no diffusion gradient (see PARAMETER ESTIMATES for description of actual dimensions used).

We use three coupled differential equations to describe mass flux between blood, PFC, and gas. Equation 4 represents the rate of change of mass (beta  · V · P) of a dissolved gas in the capillary blood compartment. It is equal to the rate of gas delivery to the capillary space via blood flow, the rate of gas removal via blood flowing out of the capillary, and the rate of diffusive gas flux across the alveolar capillary membrane into the PFC. Equation 5 describes radial diffusion in the PFC shell, which has spherical symmetry (1). Equation 6 represents the rate of change of mass in the central air space determined by addition of gas via inspiration, subtraction of gas removed by expiration, and subtraction of gas diffusing across the air-liquid interface from the PFC layer adjacent to the compartment
&bgr;<SUB>b</SUB> ⋅ V<SUB>c</SUB> ⋅ <FR><NU>dP<SUB>c</SUB></NU><DE>d<IT>t</IT></DE></FR> = <A><AC>Q</AC><AC>˙</AC></A> ⋅ &bgr;<SUB>b</SUB> ⋅ (P<SUB><OVL>v</OVL></SUB> − P<SUB>c</SUB>)
 − (<IT>D</IT><SUB>pfc</SUB> ⋅ &bgr;<SUB>pfc</SUB> ⋅ 4 ⋅ &pgr; ⋅ <IT>r</IT><SUP>2</SUP><SUB>c</SUB>) ⋅ <FENCE><FR><NU>∂P<SUB>pfc</SUB></NU><DE>∂<IT>r</IT></DE></FR></FENCE><SUB><IT>r</IT>=<IT>r</IT><SUB>c</SUB></SUB> (4)
&bgr;<SUB>pfc</SUB> ⋅ V<SUB>pfc</SUB> ⋅ <FR><NU>∂P<SUB>pfc</SUB>(<IT>r</IT>, <IT>t</IT>)</NU><DE>∂<IT>t</IT></DE></FR> = <IT>D</IT><SUB>pfc</SUB> ⋅ &bgr;<SUB>pfc</SUB> ⋅ V<SUB>pfc</SUB>
 ⋅ <FENCE><FR><NU>∂<SUP>2</SUP>P<SUB>pfc</SUB></NU><DE>∂<IT>r</IT><SUP>2</SUP></DE></FR> + <FR><NU>2</NU><DE><IT>r</IT></DE></FR> ⋅ <FR><NU>∂P<SUB>pfc</SUB></NU><DE>∂<IT>r</IT></DE></FR></FENCE> (5)
&bgr;<SUB>g</SUB> ⋅ V<SUB>g</SUB> ⋅ <FR><NU>dP<SUB>g</SUB></NU><DE>d<IT>t</IT></DE></FR> = <A><AC>V</AC><AC>˙</AC></A><SC>a</SC> ⋅ &bgr;<SUB>g</SUB> ⋅ (P<SUB>gi</SUB>− P<SUB>g</SUB>)
 + (<IT>D</IT><SUB>pfc</SUB> ⋅ &bgr;<SUB>pfc</SUB> ⋅ 4 ⋅ &pgr; ⋅ <IT>r</IT><SUP>2</SUP><SUB>g</SUB>) ⋅ <FENCE><FR><NU>∂P<SUB>pfc</SUB></NU><DE>∂<IT>r</IT></DE></FR></FENCE><SUB><IT>r</IT>=<IT>r</IT><SUB>g</SUB></SUB> (6)

The system of three partial differential equations was solved numerically to determine the partial pressure profiles in the PFC layer from the capillary-PFC interface to the PFC-gas interface. Spatial derivatives were determined by finite difference, and time derivatives were solved using LSODE, a time-integrating algorithm developed by Hindmarsh (7). The executable program was submitted as a batch job in which each simulation was solved numerically using an IBM model RS6000 computer running Unix version 4.2. Pc and Pg are equal to Ppfc(r) at the rc and rg boundaries. The time to steady-state equilibrium (T98) was defined as the time for the numerical solutions to converge to 98% of the analytically determined Pc and Pg for a steady-state gas diffusion in a spherical shell, as defined by Crank (1) (see Eqs. 7-11). The two calculated mass flow rates across the capillary-PFC and PFC-gas boundaries were nearly equal at "steady state" by use of these criteria.

Steady-state gas exchange in a spherical shell with radial diffusion. Under steady-state conditions, the time rate of change of compartmental partial pressures is zero and mass flow is equal across all boundaries. We used Crank's (1) steady-state solution to Eq. 5 describing the concentration profile as a function of radial position [C(r)] in a spherical shell to simplify the above system of equations and to analytically calculate blood-gas partial pressure differences
C(<IT>r</IT>) = <FR><NU><IT>r</IT><SUB>g</SUB> ⋅ &bgr;<SUB>pfc</SUB> ⋅ P<SUB>g</SUB> ⋅ (<IT>r</IT><SUB>c</SUB> − <IT>r</IT>) + <IT>r</IT><SUB>c</SUB> ⋅ &bgr;<SUB>pfc</SUB> ⋅ P<SUB>c</SUB> ⋅ (<IT>r</IT> − <IT>r</IT><SUB>g</SUB>)</NU><DE><IT>r</IT> ⋅ (<IT>r</IT><SUB>c</SUB> − <IT>r</IT><SUB>g</SUB>)</DE></FR>
 = &bgr;<SUB>pfc</SUB> ⋅ P<SUB>pfc</SUB>(<IT>r</IT>) (7)
Differentiating Eq. 7 with respect to r, evaluating beta pfc · dPpfc/dr at r = rc and also at r = rg, and then substituting into Eqs. 4 and 6 gives
<A><AC>Q</AC><AC>˙</AC></A> ⋅ &bgr;<SUB>b</SUB> ⋅ P<SUB><OVL>v</OVL></SUB> = <FENCE><A><AC>Q</AC><AC>˙</AC></A> ⋅ &bgr;<SUB>b</SUB> + <IT>D</IT><SUB>pfc</SUB> ⋅ &bgr;<SUB>pfc</SUB> ⋅ 4&pgr; <FR><NU><IT>r</IT><SUB>c</SUB> ⋅ <IT>r</IT><SUB>g</SUB></NU><DE><IT>r</IT><SUB>c</SUB> − <IT>r</IT><SUB>g</SUB></DE></FR></FENCE> ⋅ P<SUB>c</SUB>
− <FENCE><IT>D</IT><SUB>pfc</SUB> ⋅ &bgr;<SUB>pfc</SUB> ⋅ 4&pgr; <FR><NU><IT>r</IT><SUB>c</SUB> ⋅ <IT>r</IT><SUB>g</SUB></NU><DE><IT>r</IT><SUB>c</SUB> − <IT>r</IT><SUB>g</SUB></DE></FR></FENCE> ⋅ P<SUB>g</SUB> (8)
P<SUB>gi</SUB> ⋅ &bgr;<SUB>g</SUB> ⋅ <A><AC>V</AC><AC>˙</AC></A><SC>a</SC> = − <FENCE><IT>D</IT><SUB>pfc</SUB> ⋅ &bgr;<SUB>pfc</SUB> ⋅ 4&pgr; <FR><NU><IT>r</IT><SUB>c</SUB> ⋅ <IT>r</IT><SUB>g</SUB></NU><DE><IT>r</IT><SUB>c</SUB> − <IT>r</IT><SUB>g</SUB></DE></FR></FENCE> ⋅ P<SUB>c</SUB>
+ <FENCE><A><AC>V</AC><AC>˙</AC></A><SC>a</SC> ⋅ &bgr;<SUB>g</SUB> + <IT>D</IT><SUB>pfc</SUB> ⋅ &bgr;<SUB>pfc</SUB> ⋅ 4&pgr; <FR><NU><IT>r</IT><SUB>c</SUB> ⋅ <IT>r</IT><SUB>g</SUB></NU><DE><IT>r</IT><SUB>c</SUB> − <IT>r</IT><SUB>g</SUB></DE></FR></FENCE> ⋅ P<SUB>g</SUB> (9)
Equations 8 and 9 constitute simultaneous equations in two unknowns (Pc and Pg). Substituting K1 = Dpfc · 4pi · (rc · rg)/(rc - rg) · beta pfc/beta g and solving for Pc and Pg
P<SUB>c</SUB> = <FR><NU>P<SUB><OVL>v</OVL></SUB> <FENCE>1 + <FR><NU><A><AC>V</AC><AC>˙</AC></A><SC>a</SC></NU><DE><IT>K</IT><SUB>1</SUB></DE></FR> + <FR><NU><A><AC>V</AC><AC>˙</AC></A><SC>a</SC>⋅&bgr;<SUB>g</SUB></NU><DE><A><AC>Q</AC><AC>˙</AC></A>⋅&bgr;<SUB>b</SUB></DE></FR>⋅<FR><NU>P<SUB>gi</SUB></NU><DE>P<SUB><OVL>v</OVL></SUB></DE></FR></FENCE></NU><DE>1 + <A><AC>V</AC><AC>˙</AC></A><SC>a</SC> <FENCE><FR><NU>1</NU><DE><IT>K</IT><SUB>1</SUB></DE></FR> + <FR><NU>&bgr;<SUB>g</SUB></NU><DE><A><AC>Q</AC><AC>˙</AC></A>⋅&bgr;<SUB>b</SUB></DE></FR></FENCE></DE></FR> (10)
P<SUB>g</SUB> = <FR><NU>P<SUB><OVL>v</OVL></SUB> <FENCE>1 + <FR><NU><A><AC>V</AC><AC>˙</AC></A><SC>a</SC></NU><DE><IT>K</IT><SUB>1</SUB></DE></FR>⋅<FR><NU>P<SUB>gi</SUB></NU><DE>P<SUB><OVL>v</OVL></SUB></DE></FR> + <FR><NU><A><AC>V</AC><AC>˙</AC></A><SC>a</SC>⋅&bgr;<SUB>g</SUB></NU><DE><A><AC>Q</AC><AC>˙</AC></A>⋅&bgr;<SUB>b</SUB></DE></FR>⋅<FR><NU>P<SUB>gi</SUB></NU><DE>P<SUB><OVL>v</OVL></SUB></DE></FR></FENCE></NU><DE>1 + <A><AC>V</AC><AC>˙</AC></A><SC>a</SC> <FENCE><FR><NU>1</NU><DE><IT>K</IT><SUB>1</SUB></DE></FR> + <FR><NU>&bgr;<SUB>g</SUB></NU><DE><A><AC>Q</AC><AC>˙</AC></A>⋅&bgr;<SUB>b</SUB></DE></FR></FENCE></DE></FR> (11)
and
P<SUB>c</SUB> − P<SUB>g</SUB> = <FR><NU>P<SUB><OVL>v</OVL></SUB> − P<SUB>gi</SUB></NU><DE>1 + <FR><NU><IT>K</IT><SUB>1</SUB></NU><DE><A><AC>V</AC><AC>˙</AC></A><SC>a</SC></DE></FR> + <FR><NU><IT>K</IT><SUB>1</SUB></NU><DE><A><AC>Q</AC><AC>˙</AC></A></DE></FR> <FR><NU>&bgr;<SUB>g</SUB></NU><DE>&bgr;<SUB>b</SUB></DE></FR></DE></FR> (12)
At the extremes of no PFC (rc = rg) and Pgi = 0, Eqs. 10 and 11 reduce to the MIGET equations for retention (R) and excretion (E): R = E = lambda b/(lambda b + VA/Q), where lambda b = beta b/beta g. Notice that the capillary-to-gas partial pressure (Pc - Pg) difference is dependent on the absolute values of VA and Q.

    PARAMETER ESTIMATES

Parameter values were chosen to correspond to the dimensions of lung structure and function of healthy piglets weighing 2-4 kg. Piglets this size typically have a functional residual capacity of 30 ml/kg and respiratory rate (RR) of 20 breaths/min. For calculation purposes, an average weight of 2.5 kg was used. As discussed above, our gas exchange unit represents a terminal sac in the lung of a piglet. Haefeli-Bleuer and Weibel (5) measured the outer diameter of human terminal sacs (an alveolar duct plus 2 alveoli in total width) to be 656 ± 127 µm. Tenney and Remmers (18) showed that species variation in alveolar diameter was correlated to metabolic rate per unit body weight, with adult pig alveolar diameter ~91% of the diameter of human alveoli (656 × 0.91 = 597 µm). On the basis of these data we chose an end-inspiratory rc of 300 µm. Surface area and volume of a single spherical unit are therefore 0.0113 cm2 and 0.000113 cm3, respectively. The number (n) of terminal sacs or gas exchange units in a piglet lung was then determined by the ratio of end-inspiratory lung volume [(48 ml/kg) × (2.5 kg) = 120 ml at rc = 300 µm] to gas exchange unit volume (1.13 × 10-5 ml): 1,062,000 units/lung, which we rounded to 1 × 106. End-inspiratory lung volume was determined as the sum of functional residual capacity lung volume (30 ml/kg), tidal volume (VT, 15 ml/kg), and 3 ml/kg associated with positive end-expiratory pressure of 5 cmH2O used in all our experimental work (12). If there are 20 alveoli per terminal gas exchange unit, there would be ~20 × 106 alveoli/piglet. Lung volume is obviously not constant throughout the respiratory cycle. We evaluated the steady-state model (Eqs. 10 and 11) for several lung volumes in the range of tidal breathing, i.e., rc of 270 and 300 µm, to illustrate the impact of lung volume on (A-a)DO2 and (a-A)DCO2. We did not simulate tidal breathing in the sense of second-to-second variation in Vg.

Ventilation per gas exchange unit (VA) was determined using our typical experimental VT of 15 ml/kg (12), estimated dead space (VD) of 4.5 ml/kg, RR of 20 min-1, M of 2.5 kg, and n as follows: VA = (VT - VD) · RR · M/n = 8.74 × 10-6 ml/s. Blood flow per gas exchange unit (Q) was derived from average piglet cardiac output of 500 ml/min (12) divided by n: 8.33 × 10-6 ml/s. The capillary blood volume was derived on the basis of anatomic data that show pulmonary capillaries to cover 75% of the alveolar surface (i.e., capillary surface area per spherical model unit = 0.75 × 0.0113 cm2) and have a thickness equivalent to the red cell diameter (5 µm), giving a Vc per unit of 4.24 × 10-6 ml. The volume of PFC per gas exchange unit (Vpfc) was determined from the total dose of PFC divided by n. For example, a dose of 30 ml/kg in a 2.5-kg piglet results in a total dose of 75 ml, or Vpfc of 7.5 × 10-5 ml/unit. PFC layer thickness is dependent on rc and the volume of PFC present, with the assumption that PFC is distributed as a spherical shell with a gas hole in the center (Fig. 2). Normal parameter values for the spherical model under matched VA/Q and Vpfc of 30 ml/kg are summarized in Table 1.

                              
View this table:
[in this window]
[in a new window]
 
Table 1.   Normal model parameters for VA/Q = 1 and Vpfc = 30 ml/kg

Values of beta b and beta pfc for inert gases were obtained from experimental measurements of gas solubility in pig blood and in the PFC perflubron (C8F17Br, LiquiVent, Alliance Pharmaceutical, San Diego, CA) (12). The "solubility" of a tracer gas in the gas phase (beta g) is defined in the classic paper by Piiper et al. (16) as 0.00132 Torr-1 (=1/760 at sea level).

The solubility of O2 and CO2 in blood was determined by the slope of the curve of gas content vs. partial pressure. This relationship is nonlinear over the physiological range of partial pressures of these gases because of chemical binding in the blood. O2 combines with Hb, resulting in an S-shaped concentration vs. pressure curve in the partial pressure range 0-150 Torr. For PO2 >150 Torr, the concentration vs. partial pressure curve is linear, because Hb is saturated, and for O2, beta b is the same as in plasma: 0.003 ml · 100 ml solvent-1 · Torr-1. For PO2 <150 Torr, beta b for O2 is much higher; e.g., at PO2 of 40 Torr it is 0.06 ml · 100 ml solvent-1 · Torr-1 as determined by the slope of the O2 content (CO2, ml O2/100 ml blood) vs. PO2 (Torr) curve generated by the subroutines of Olszowka and Farhi (15). For the steady-state partial pressure differences calculated using Eqs. 10 and 11, we used only beta b for O2 of 0.003, because for all the experimental data against which we are comparing model results arterial PO2 (PaO2) was >150 Torr (12). The solubility of CO2 in blood is a function of dissolved CO2 as well as CO2 converted to HCO-3. The content (CCO2) vs. PCO2 curve is approximately linear within 40-80 Torr PCO2. With use of the blood-gas routines of Olszowka and Farhi, beta b for CO2 was determined from the slope of CCO2 vs. PCO2 over this range and was found to be 0.779 ml · 100 ml blood-1 · Torr-1. O2 and CO2 solubilities in PFC were provided by Alliance Pharmaceutical (Table 2).

                              
View this table:
[in this window]
[in a new window]
 
Table 2.   Solubilities of 6 inert gases, O2, and CO2 in blood and PFC

Few molecular diffusion coefficients (Dpfc) of dissolved gases in PFC are precisely known. Tham et al. (19) measured Dpfc of O2 and CO2 in three perfluorochemicals (Caroxin-D, Caroxin-F, and FC-80), finding the average diffusion coefficient for O2 in PFC to be 5.61 × 10-5 cm2/s at 37°C with a range of 5.57-5.65 × 10-5 cm2/s and for CO2 in PFC at 37°C to be 4.36 × 10-5 cm2/s with a range of 4.21-4.48 × 10-5 cm2/s. The diffusion coefficients of O2 and CO2 in H2O at 37°C are 3.3 × 10-5 and 2.6 × 10-5 cm2/s, respectively (6).

We used the average value of the CO2 diffusion coefficient as measured by Tham et al. (19) to estimate Dpfc of each respiratory gas in perflubron, the PFC used in our experiments. There are no experimental data available measuring diffusivity in PFC of the six inert gases used in MIGET (9, 20). Their diffusivities in H2O at 37°C are 1.63 × 10-5 cm2/s for SF6, 1.96 × 10-5 cm2/s for ethane, 1.84 × 10-5 cm2/s for cyclopropane, 1.28 × 10-5 cm2/s for halothane, 0.85 × 10-5 cm2/s for ether, and 1.62 × 10-5 cm2/s for acetone (17, 21). Because their diffusivities in H2O are only slightly less than those of CO2 in H2O, we chose the value of Dpfc for CO2 in PFC to represent the diffusivity of the six inert gases in the absence of experimental data.

    RESULTS

Solutions for both of the models were well behaved with no instances of negative results or mass imbalance. Partial pressures at the boundaries between compartments were continuous. The numerically integrated time- and space-dependent solutions for the spherical model converged on the analytic steady-state solutions. For each of the eight gases simulated, the time to steady-state equilibrium was estimated by two independent models, and the times generated by both models were within 30% of each other and usually within 10%.

Time to reach steady-state equilibrium. Figure 3 illustrates the time rate of change of partial pressure of the eight gases in the simpler well-mixed three-compartment model with Vpfc of 30 ml/kg after a step change in the input partial pressure of each gas. For O2 this involved setting Pgi at 650 Torr and P<SUB><OVL>v</OVL></SUB> at 40 Torr and for the remaining 7 gases Pgi at 0 Torr and P<SUB><OVL>v</OVL></SUB> at 1 Torr at time 0. Figure 3A illustrates the application of Eq. 3 for normal conditions of matched VA and Q (VA/Q = 1). Figure 3B illustrates the same for near-zero ventilation (shunt conditions). Because O2 is delivered by ventilation, PO2 was not simulated for shunt conditions. Whenever VA is negligible, Eq. 3 shows that the final value is always P<SUB><OVL>v</OVL></SUB>, and the time constants are lengthened. Both features are apparent in Fig. 3B. Gases with the lowest beta b/beta pfc ratio (i.e., SF6) take the longest to equilibrate, because PFC acts as a large capacitor that fills slowly when there is great disparity in solubilities.


View larger version (25K):
[in this window]
[in a new window]
 
Fig. 3.   Normalized partial pressure (P/P<SUB><OVL>v</OVL></SUB>) vs. time of 6 inert gases, O2, and CO2 in well-mixed 3-compartment model with Vpfc = 30 ml/kg. P represents Pc = Ppfc = Pg in well-mixed model, beta  values are as described in PARAMETER ESTIMATES. Pgi = 0 and P<SUB><OVL>v</OVL></SUB> = 1 for 6 inert gases and CO2 (smooth lines); for O2 Pgi = 650 and P<SUB><OVL>v</OVL></SUB> = 0. A: solutions for normal conditions (VA = 8.74 × 10-6 ml/s and Q = 8.33 × 10-6 ml/s) after a step change in P<SUB><OVL>v</OVL></SUB> or Pgi at time 0. B: solutions for shuntlike conditions (VA = 0 and Q = 8.33 × 10-6 ml/s) after a step change in input partial pressures. T98 values are given in Table 3.

Figure 4 demonstrates the time and space rate of change in the spherical gas exchange unit with 30 ml/kg PFC and matched VA and Q (as described in PARAMETER ESTIMATES). Figure 4A shows successive time traces of PCO2 vs. radial distance from the capillary through PFC to the central gas region. After a step change in P<SUB><OVL>v</OVL></SUB> from 0 to 40 Torr, PCO2 increases in the gas exchange unit until it converges on the steady-state value. Figure 4B shows similar successive time traces of PO2 vs. radial distance through the PFC after a step change in Pgi from 0 to 650 Torr.


View larger version (23K):
[in this window]
[in a new window]
 
Fig. 4.   Examples of time- and space-dependent changes in PCO2 (A) and PO2 (B) in spherical shell model. Multiple tracings represent P in PFC layer vs. r at 2-s intervals after a step change in P<SUB><OVL>v</OVL></SUB> of 40 Torr PCO2 (A) or in Pgi of 650 Torr PO2 with P<SUB><OVL>v</OVL></SUB> = 40 (B). Solutions converge on steady-state values. Pc = Ppfc at r = rc and Pg = Ppfc at r = rg. Model parameters for these solutions are as follows: Vpfc = 30 ml/kg, rc = 300 µm, rg = 210 µm, VA = 8.74 × 10-6 ml/s, and Q = 8.33 × 10-6 ml/s. Vc, D, and beta  values are defined in PARAMETER ESTIMATES.

Table 3 reports the T98 for eight gases in each of the two models with Vpfc of 30 ml/kg. T98 values were defined slightly differently for the two models. In the well-mixed model T98 was defined as 4tau in Eq. 3; for the spherical model it was the time at which the time-dependent solutions (Eqs. 4-6) converged to 98% of the analytic steady-state solutions (Eqs. 10 and 11). We evaluated the model for three conditions to illustrate the range of T98 likely to be encountered in the lung partially filled with PFC: matched VA and Q, VA approximately zero with Q normal (shunt), and Q near zero with VA normal (dead space).

                              
View this table:
[in this window]
[in a new window]
 
Table 3.   T98 after a step change in input partial pressure (P<SUB><OVL>v</OVL></SUB> or Pgi) with Vpfc = 30 ml/kg and rc = 300 µm

For VA and Q well-matched (VA/Q = 1), all times to steady state were <3 min. The gas with the longest time to steady state was cyclopropane followed by ether, SF6, halothane, O2, CO2, and acetone. Under shunt conditions all times to steady state were prolonged (except for acetone, which is insensitive to shunt), with SF6 having the longest times at ~5 h. The time to steady state for O2 was also markedly prolonged at ~27 min, whereas that for CO2 remained short at 15-20 s. Under dead space conditions the times were intermediate, with the longest being for acetone at ~26 min. CO2 equilibration times were mildly prolonged under these conditions, ~95 s.

Steady-state gas exchange in a PFC-filled spherical shell. Steady-state partial pressure differences of inert and respiratory gases were calculated from Eqs. 10 and 11. SF6 was left out of the following discussion, since it was not included in our experimental MIGET analysis (12, 13) because of its prohibitively long time to reach steady state under shunt conditions. Figure 5 shows Pc-Pg differences of seven gases normalized by input partial pressure (P<SUB><OVL>v</OVL></SUB> for the 5 inert gases and CO2 and Pgi for O2) vs. rg for VA of 8.74 × 10-6 ml/s and Q of 8.33 × 10-6 ml/s. An rg of 0 corresponds to a flooded terminal sac with no gas compartment, and rg of 300 µm corresponds to a gas exchange unit with no PFC. Values of rg equal to 210, 250, and 280 µm correspond to the three doses of PFC used in our experimental work: 30, 20, and 10 ml/kg, respectively (12). CO2 shows the largest difference at all values of rg, with the Pc-Pg difference nearly 10% of the input pressure at rg of 210 µm. The partial pressure gradient of O2 is very low at the same dose (<1% of Pgi), rising only when rg becomes very small as the gas exchange unit becomes flooded with PFC. The inert gases also show a negligible Pc-Pg difference for rg of 210 µm, with halothane having the largest Pc-Pg difference at 3% of P<SUB><OVL>v</OVL></SUB> followed by cyclopropane, ethane, acetone, and ether. As the PFC layer increases in thickness, the Pc-Pg difference rises exponentially, approaching P<SUB><OVL>v</OVL></SUB> values for CO2 and the inert gases and Pgi for O2. For these simulations, P<SUB><OVL>v</OVL></SUB> for O2 was set to zero and Pgi to 1 for the sake of comparison.


View larger version (17K):
[in this window]
[in a new window]
 
Fig. 5.   Spherical shell model. Steady-state normalized and capillary-to-gas partial pressure (Pc-Pg) differences are shown vs. PFC thickness for 5 inert gases plus O2 and CO2. Radial thickness of PFC is varied from 0 (no PFC present, rg = rc) to 300 µm (all PFC, rg = 0). Arrows under abscissa indicate rg corresponding to 30, 20, and 10 ml/kg doses of PFC. Vc, D, and beta  are defined in PARAMETER ESTIMATES; Vpfc and Vg vary with rg. Pgi = 0 and P<SUB><OVL>v</OVL></SUB> = 1 for all inert gases. VA = 8.74 × 10-6 ml/s and Q = 8.33 × 10-6 ml/s. Top curve is for CO2, bottom curve is for ether, and other curves are for gases in sequence shown. O2 and ethane are indistinguishable. Capillary-to-gas disequilibrium is, by far, the largest for CO2.

Figure 6 illustrates the effect of gas exchange unit volume ("lung volume") on partial pressure difference of O2 and CO2. Although our model does not incorporate features of tidal breathing, we explored the effect of varying the gas exchange unit volume between the extremes of end inspiration (rc = 300 µm) and end expiration (rc = 270 µm). This might be equivalent to breath-holding maneuvers at the extremes of cyclic breathing. For both gases, the Pc-Pg difference increased at the lower lung volume for all PFC doses. The percent increase in the Pc-Pg difference was greater with larger doses of PFC. The Pc-Pg difference for CO2 with a P<SUB><OVL>v</OVL></SUB> of 40 Torr and 30 ml/kg PFC in the lung varied from 3.7 Torr at the large lung volume to 9.6 Torr at the lower lung volume. At the small PFC dose of 10 ml/kg, the Pc-Pg difference for CO2 varied from 0.8 to 1.2 with the change in lung volume. We previously showed (a-A)DCO2 in healthy animals with 30 ml/kg PFC in the lungs to be 12 Torr (12). The difference for Pc-Pg difference for O2 varied in a similar manner with an increase from 6 to 18 Torr as the gas exchange unit volume decreased with 30 ml/kg PFC in the lung and Pgi of 650 Torr.


View larger version (13K):
[in this window]
[in a new window]
 
Fig. 6.   Pc-Pg difference for O2 and CO2 vs. PFC dose at a high gas-exchange unit volume (bullet ) and a low gas-exchange unit volume () for VA = 8.74 × 10-6 ml/s and Q = 8.33 × 10-6 ml/s. D and beta  values are defined in PARAMETER ESTIMATES. A: Pg-Pc difference for O2 for P<SUB><OVL>v</OVL></SUB> = 40 Torr and Pgi = 650 Torr rises with dose and with smaller volumes. B: Pc-Pg difference for CO2 with P<SUB><OVL>v</OVL></SUB> = 40 Torr and Pgi = 0 Torr shows the same trend as in A.

We examined the impact of varying VA and Q independently on the Pc-Pg difference for inert gases, O2, and CO2. At PFC thicknesses up to 100 µm (PFC dose ~30 ml/kg), varying VA and Q had a small impact on MIGET gas Pc-Pg differences. At PFC thicknesses >100 µm, the gradients increased exponentially, as in the case of matched VA and Q (Fig. 5). The Pc-Pg difference for the inert gases never exceeded 10% of P<SUB><OVL>v</OVL></SUB> over this range of VA and Q. Each gas was affected to a different degree depending on their relative solubilities. Figure 7 illustrates the effect of varying VA and Q on O2 and CO2. Figure 7A shows minimal effect on the Pg-Pc difference for O2 with varying VA over a range from 0.1 to 10 times the average ventilation of a terminal alveolar sac (8.74 × 10-6 ml/s) with Q fixed (8.33 × 10-6 ml/s). Figure 7B shows an 8-fold increase in the Pg-Pc difference for O2 with a 10-fold increase in Q. The gradient drops to near zero as Q decreases to 0.1 its average value. Changes in partial pressure differences of CO2 with varying VA and Q are shown in Fig. 7, C and D. There is a 3-fold increase in the Pc-Pg difference with VA 10 times its average value, and the gradient drops to near zero with VA at 0.1 its average value. The Pc-Pg difference for CO2 drops in half with a decrease in Q but is essentially unchanged with a 10-fold increase in Q. Comparison of the solutions 10 × VA and 0.1 × Q in Fig. 7, A and B, as well as 7, C and D, illustrates that the Pc-Pg gradient is different for each condition, despite equivalent VA/Q ratios.


View larger version (19K):
[in this window]
[in a new window]
 
Fig. 7.   Pg-Pc difference for O2 (A and B) and Pc-Pg differences for CO2 (C and D) vs. PFC thickness, with VA or Q varied while the other is fixed at a mean value (VA = 8.74 × 10-6 ml/s and Q = 8.33 × 10-6 ml/s). Radial thickness of PFC ranges from 0 (no PFC present) to 100 µm; arrows under abscissa indicate rg corresponding to Vpfc = 10, 20, and 30 ml/kg with rc = 300 µm. Vc, D, and beta  values are defined in PARAMETER ESTIMATES. For CO2, P<SUB><OVL>v</OVL></SUB> = 40 Torr and Pgi = 0 Torr; for O2, P<SUB><OVL>v</OVL></SUB> = 40 Torr and Pgi = 650 Torr. A: Pg-Pc difference for O2 is relatively insensitive to variation in VA from 0.1 to 10 times its average value. B: Pg-Pc for O2 rises 8-fold with a 10-fold increase in Q and drops to near zero at 0.1 its average value. C and D: sensitivity of Pc-Pg difference for CO2 to variations in VA and Q from 0.1 to 10 times their average value.

    DISCUSSION

Evaluation of model assumptions. We had two specific questions in mind when developing these models of gas exchange in a terminal sac filled with PFC: 1) Do gases that are exchanged in a PFC-containing alveolus reach steady state at usual respiratory rates? 2) How large are the alveolar-arterial differences as a result of diffusion across PFC barriers? Two different models were developed in an attempt to answer these questions. The well-mixed three-compartment model provided a simple approach to estimating time to steady state. Its major assumptions are that neither diffusion times in the PFC nor the geometry of a gas-exchanging subunit significantly affect the solutions. By contrast, our spherical model explicitly incorporated the diffusion gradients and more realistic geometry but, despite major mathematical differences, the results showed very close agreement with the well-mixed compartment model predictions of time to steady state (Table 3).

Both models depict gas exchange in a single terminal alveolar sac. Parameters such as Q, VA, and Vpfc were arrived at by partitioning an equal amount of Q, VT, and Vpfc to all terminal sacs in the lung. The lung is not homogeneous in its distribution of any of these parameters, and application of model results to interpretation of experimental data must be done with this in mind. In reality, there will be a heterogeneous distribution of gas exchange units ranging from completely PFC filled to partially PFC filled to completely gas filled that are ventilated and perfused in some heterogeneous distribution. Measured arterial and expired gas partial pressures are weighted averages of gas exchange subunits. Model predictions of gas exchange in a single terminal sac help us explore the range of possible alveolar Pc-Pg differences due to diffusion limitation and provide a gross approximation to overall lung arterial-alveolar differences.

Additionally, the choice of a spherical shape of our gas exchange unit to approximate the terminal alveolar sac likely overestimates the diffusion barrier somewhat. A terminal sac is not a smooth sphere but, rather, a cluster of cup-shaped alveoli opening up to a common duct. There are sheets of perfused alveolar-capillary membrane extending inward toward the duct that increase the surface area for exchange and bring those parts of the membrane close to the PFC-gas interface. This would be equivalent to "thinning" the PFC spherical shell in our model and decreasing the Pc-Pg gradient for ventilated units. On the other hand, regions of shunt where VA is zero probably behave similarly to the model as the PFC pool equilibrates with mixed venous blood and geometry becomes irrelevant.

We feel justified in our choice of inert gas diffusion coefficients on the basis of the fact that the inert gases and CO2 had similar diffusion coefficients in H2O and that all should have increased diffusivity in perflubron because it is a nonpolar solvent. The rate of diffusion of a molecule through a fluid medium depends on the "effective radius" of the molecule, a function of molecular size and van der Waals interactions with neighboring molecules. Increased Dpfc for O2 and CO2 in PFC compared with H2O suggests that the molecules have smaller effective volumes in PFC because of reduced van der Waals interactions. Although there are certain to be discrepancies between the true diffusion coefficients of these gases in PFC and our approximated Dpfc, model results show little dependence of our time- or space-dependent solutions on diffusive resistance. As we demonstrate, the disparity in Pc-Pg gradients for different gases with the same Dpfc (i.e., CO2 vs. ether) supports the conclusion that minor variations in Dpfc will not significantly affect our model results.

Time to reach steady state. Of particular interest to us was whether the inert gases used in MIGET would reach steady state during PLV within the time period of our experimental measurements (12, 13). In using MIGET to assess VA/Q heterogeneity in healthy piglets during PLV, we modified the standard protocol (9, 20) to incorporate a 60-min equilibration period between experimental conditions (15 min is more common). Results from our two models suggest that all gases come to equilibrium well within this time period with the exception of SF6 under shunt conditions (>4 h). We showed this previously and eliminated SF6 from MIGET analyses during PLV (12, 13). The remaining five inert gases reach steady state within the 1-h time frame. The next longest equilibration time was for acetone, which took ~26 min to come to steady state under "dead space" conditions (Q = 0). Except for one case, O2 and CO2 reached steady state in <2 min for the range of possible VA and Q that might exist during PLV. O2 took 26.6 min to reach steady state under shunt conditions because of the slow delivery rate.

We did not incorporate the periodic nature of VA and Q in our model, but this would be a useful extension. It would be interesting to see if O2 and CO2 reach steady state, despite breath-to-breath variations in Pg and pulsatile changes in Pc that occur over a 2- to 4-s time period. Intuition leads us to think that a gas exchange unit would reach steady state about an average value of Pc and Pg, filtering out second-to-second fluctuations.

Partial pressure differences at steady state. PFC acts as a mild diffusion barrier for all gases in the steady state, creating a Pc-Pg difference that increases with volume of liquid in the alveolar space (Fig. 5). Less intuitive is the fact that the partial pressure gradient for each gas is different on the basis of the relative solubility of the gas in blood, PFC, and the gas phase. The presence of a partial pressure gradient in alveoli during steady-state gas exchange has several important consequences. Gas exchange efficiency is reduced with overall arterial-alveolar partial pressure gradients increased compared