This page will look better in a graphical browser that supports web standards, but is accessible to any browser or internet device.

Served by Samwise.


Recirculation of O2 and CO2 between a 2-compartment lung and 2-compartment body.

Model number: 0144

Run Model: 
    Help running a JSim model.
Java runtime required. Web browser must support Java Applets.
(JSim model applet may take 10-20 seconds to load.)


  This oversimplified whole body model captures the uptake of O2 in the
  lung, delivery to tissue, intratissue consumption and CO2 production
  with return of CO2 to the lung where it is lost into the expired air. 
  This lung/body exchange of O2 and CO2, while ridiculously simple, allows 
  one to explore the values for flows, PSs and Volumes, for metabolic rates, 
  and alveolar exchange that give rise to the values in a normally ambient 
  human. The alveolar exchanges are passiveand the tissue consumption is 
  first order Michaelis-Menten. The ratio of CO2 produced to O2 used in the 
  Respiratory Quotient, RQ.  Concentrations of O2 and CO2 in the blood are 
  expressed in the equivalent partial pressures assuming a Henry's Law type 
  of form (concO2 = alphaO2 * pO2) for conversion of partial pressures of 
  O2 and CO2 to the concentrations bound and dissolved in the blood, tissue
  and lung interstitial fluid (isf).    


The transport between and within each of the four regions can be represented by convective transport, e.g. F/V11*(O21-O11), diffusive transport between regions, e.g. PS1O/V11*(Cin-O12), and consumption/production terms, e.g. Vcytox/((Kcytox + O22)*V22)*O22, The equations are as follows:

// Input function
extern real Cin(t) mmHg; //Concentration in inflow stream.

//Initial conditions
when(t=t.min){  O11=0; O12=100; O21=0; O22=0; 
                C11=0; C12=0; C21=0; C22=60;  }

O11:t = F/V11*(O21-O11) + PS1O/V11*(O12-O11);
O12:t = 		  PS1O/V12*(O11-O12) + (PSalv/V12)*(Cin-O12);
O21:t = F/V21*(O11-O21) + PS2O/V21*(O22-O21);
O22:t =                   PS2O/V22*(O21-O22) - Vcytox/((Kcytox + O22)*V22)*O22;
C11:t = F/V11*(C21-C11) + PS1C/V11*(C12-C11);
C12:t =                   PS1C/V12*(C11-C12) - GCair*C12/V12;
C21:t = F/V21*(C11-C21) + PS2C/V21*(C22-C21);
C22:t =                   PS2C/V22*(C21-C22) + RQ*Vcytox/((Kcytox + O22)*V22)*O22;

Where F is blood flowrate, O11-O22 are partial pressures of O2 in each of the regions, C11-C22 are partial pressures of CO2 in each of the regions.
PSalv, PS1O and PS2O are the permeability-surface area product of O2 transport between lung and lung isf, lung isf and blood, and blood and tissue respectively, PS1C and PS2C are the permeability-surface area product of CO2 transport between lung isf and blood, and blood and tissue respectively.

V11-V22 are region volumes, the Vcytox/((Kcytox + O22)*V22)*O22 term is the O2 consumption in tissue, the RQ*Vcytox/((Kcytox + O22)*V22)*O22 is the CO2 production in the tissue and the GCair*C12/V12 term is the exhalation of CO2 from the lung.

Units for concentrations as mmHg: Partial pressures are "activities" equivalent to molar concentrations times an activity coefficient. Just as the p50 for the hemoglobin binding of oxygen can be expressed in mmHg or in mM, so also can the binding of O2 to cytochrome oxidase be expressed in either.

Download JSim model project file



Related Models

Key Terms

Alveolar gas exchange, whole body metabolism, respiratory quotient, oxygen, carbon dioxide, cardiac output, plasma oxygen levels, plasma carbon dioxide, no buffering, no hemoglobin, gas solubilites in blood, air-blood gas exchange

Model Feedback

We welcome comments and feedback for this model. Please use the button below to send comments:

Model History

Get Model history in CVS.

Posted by: BEJ


Please cite in any publication for which this software is used and send an email with the citation and, if possible, a PDF file of the paper to:
Or send a copy to:
The National Simulation Resource, Director J. B. Bassingthwaighte, Department of Bioengineering, University of Washington, Seattle WA 98195-5061.

[This page was last modified 14Mar18, 3:17 pm.]

Model development and archiving support at provided by the following grants: NIH U01HL122199 Analyzing the Cardiac Power Grid, 09/15/2015 - 05/31/2020, NIH/NIBIB BE08407 Software Integration, JSim and SBW 6/1/09-5/31/13; NIH/NHLBI T15 HL88516-01 Modeling for Heart, Lung and Blood: From Cell to Organ, 4/1/07-3/31/11; NSF BES-0506477 Adaptive Multi-Scale Model Simulation, 8/15/05-7/31/08; NIH/NHLBI R01 HL073598 Core 3: 3D Imaging and Computer Modeling of the Respiratory Tract, 9/1/04-8/31/09; as well as prior support from NIH/NCRR P41 RR01243 Simulation Resource in Circulatory Mass Transport and Exchange, 12/1/1980-11/30/01 and NIH/NIBIB R01 EB001973 JSim: A Simulation Analysis Platform, 3/1/02-2/28/07.