Numerical simulation of pulmonary O2 and CO2 exchange

The effects of chemical kinetics on oxygen delivery to tissue

A mathematical model has been formulated to analyze the effect of nonequilibrium kinetics on oxygen delivery to tissue. The model takes into account molecular diffusion, facilitated diffusion in the capillary blood, convection, chemical kinetics of O2 with hemoglobin, and the rate of metabolic consumption. A line iterative technique is described to solve numerically the resulting coupled system of nonlinear partial differential equations with physiologically relevant boundary and entrance conditions. With nonequilibrium kinetics the end-capillary PO2 is found to be lower than that in the venous blood. The effect is more pronounced during hypoxia and anemia. It is found that the tissue PO2 at the lethal corner decreases with the decrease in blood velocity, arterial PO2, hemoglobin concentration, P50, and increase in COHb concentration or metabolic rate, while the difference between end-capillary PO2 and venous PO2 increases, which reflects the effect of nonequilibrium kinetics on the delivery of O2 to tissue. Thus, the consideration of venous PO2 as an indicator of tissue PO2 in clinical and experimental studies may be questionable.

Mathematical model for the computation of alveolar partial pressure of carbon monoxide

A mathematical model is formulated for computing alveolar partial pressure of carbon monoxide (PACO) from that in the atmospheric air. The model takes into account parameters like inspired/expired air flow rates, diffusion capacity of the lung, concentration of CO in the atmospheric air, blood flow rate and the non-linear CO dissociation curve. The effect due to the presence of O2 in the blood on CO dissociation curve is also incorporated. It is shown that for a given atmospheric CO concentration, PACO increases exponentialy with time and attains asymptotic value. Alveolar PCO increases further with the increase in the atmospheric CO concentration. The model can also be used to compute carboxyhaemoglobin levels in the blood as a function of exposure time and the results are comparable with the CFK equation and the values measured experimentally.

A numerical model for the oxygenation of blood in lung capillaries--effect of nth order one-step kinetics of oxygen uptake by haemoglobin

A numerical model is described for the oxygenation of blood in lung capillaries by considering the transport mechanisms of molecular diffusion, convection and the facilitated diffusion due to the presence of haemoglobin. In order to represent the oxygen dissociation curve accurately in the model, the nth order one-step kinetics of oxygen uptake by haemoglobin has been used. The resulting system of coupled, non-linear partial differential equations is solved numerically. It is shown that the blood is required to traverse a larger distance in the capillary before becoming fully oxygenated with nth order one-step kinetics in comparison to first-order one-step kinetics.

A mathematical model for the computation of the oxygen dissociation curve in human blood

The mathematical relations developed by various researchers for the oxygen dissociation curve are reviewed. Using well-known mechanisms of chemical kinetics of various species in the blood, we have developed a mathematical formula to compute the oxygen dissociation curve in the blood showing its dependence on the pH and PCO2. The functional form, proposed here, is much simpler in comparison to those available in the literature for use in the mathematical modelling of O2 transport in the pulmonary and systemic circulations. In the process, the well-known Hill's equation has been generalized showing an explicit dependence on PCO2 and pH. It is shown that the oxygen dissociation curve computed from our comparatively simpler equation, fits in fairly well with the documented data and shows realistic shift with PCO2 and pH.

A numerical study of the nonsteady transport of gases in the pulmonary capillaries

A mathematical model is formulated for simulating the unsteady transport of gases in the blood flowing through the pulmonary capillaries. The formulation takes into account the transport mechanisms of molecular diffusion, convection and facilitated diffusion of the species due to haemoglobin. A time dependent situation is created by allowing to vary suddenly the partial pressures of the gases either in the venous blood or in the alveolar air. A numerical technique is described to solve the resulting time-dependent system of nonlinear coupled partial differential equations with the physiologically relevant boundary, entrance and initial conditions. The time required by the gases to achieve equilibrium is computed. It is shown that the dissolved oxygen takes longest in reaching equilibration whereas the carbon dioxide is the fastest. The various physiologically relevant unsteady situations have been examined.

The process of gas exchange in the pulmonary circulation incorporating the contribution of axial diffusion

A mathematical model is made to describe the process of gas exchange in the pulmonary circulation incorporating the contribution of axial diffusion. The model takes into account the transport mechanisms of molecular diffusion, convection and facilitated diffusion due to the presence of haemoglobin as a carrier of the gases. The mathematical formulation leads to a coupled system of non-linear elliptic partial differential equations. A numerical scheme is described to solve such a system. It is found that the axial diffusion does not have an appreciable effect on the transport of the species in the blood.

A numerical model for blood oxygenation in the pulmonary capillaries--effect of pulmonary membrane resistance

A study of the blood oxygenation in pulmonary capillaries is made by considering the transport mechanisms of molecular diffusion, convection and the facilitated diffusion due to the presence of haemoglobin. The resistance offered by the pulmonary membrane on the transport of gases has been incorporated. The resulting system of coupled, non-linear partial differential equations is solved numerically. It is found that, in the immediate neighbourhood of the entry, the amount of dissolved O2 decreases. This decreases further as the resistance offered by the pulmonary membrane increases. The rate of oxygenation of blood increases as the permeability coefficient for O2 (PO) increases. It is shown that the ideally permeable case for both O2 and CO2 can be approximated by taking PO approximately 10 cm/s. Further, it is shown that the oxygen takes longest and CO2 is the fastest to attain equilibration. The equilibration length increases as the resistance offered by the membrane increases. Finally, some of the pulmonary diseases such as pulmonary oedema and fibrosis have been analyzed.