A mathematical model of the controlled plant of the respiratory system

Integrative Computational Models of Lung Structure-Function Interactions

Anatomically based integrative models of the lung and their interaction with other key components of the respiratory system provide unique capabilities for investigating both normal and abnormal lung function. There is substantial regional variability in both structure and function within the normal lung, yet it remains capable of relatively efficient gas exchange by providing close matching of air delivery (ventilation) and blood delivery (perfusion) to regions of gas exchange tissue from the scale of the whole organ to the smallest continuous gas exchange units. This is despite remarkably different mechanisms of air and blood delivery, different fluid properties, and unique scale-dependent anatomical structures through which the blood and air are transported. This inherent heterogeneity can be exacerbated in the presence of disease or when the body is under stress. Current computational power and data availability allow for the construction of sophisticated data-driven integrative models that can mimic respiratory system structure, function, and response to intervention. Computational models do not have the same technical and ethical issues that can limit experimental studies and biomedical imaging, and if they are solidly grounded in physiology and physics they facilitate investigation of the underlying interaction between mechanisms that determine respiratory function and dysfunction, and to estimate otherwise difficult-to-access measures. © 2021 American Physiological Society. Compr Physiol 11:1501-1530, 2021.

Mathematical modeling of extracorporeal CO2 removal therapy : A validation carried out on ten pigs

The extracorporeal CO₂ removal device (ECCO₂RD) is used in clinics to treat patients suffering from respiratory failures like acute respiratory distress syndrome (ARDS) or chronic obstructive pulmonary disease (COPD). The aim of this device is to decarboxylate blood externally with low blood flow. A mathematical model is proposed to describe protective ventilation, ARDS, and an extracorporeal CO₂ removal therapy (ECCO₂RT). The simulations are compared with experimental data carried out on ten pigs. The results show a good agreement between the mathematical simulations and the experimental data, which provides a nice validation of the model. This model is thus able to predict the decrease of PCO₂ during ECCO₂RT for different blood flows across the extracorporeal lung support.

Importance of metabolism variations in a model of extracorporeal carbon dioxide removal

Extracorporeal CO2 Removal device is used in clinics when a patient suffers from a pulmonary insufficiency like Acute Respiratory Distress Syndrome and allows to decarboxylate blood externally. In this work, a model of the respiratory system coupled with such a device is proposed to analyze the decrease of CO2 partial pressure in blood. To validate the model, some parameters are estimated thanks to experimental data. Metabolism is a crucial parameter and we show that its time evolution must be taken into account in order to have correct CO2 partial pressure simulations in arteries and in veins.

Towards a virtual lung: multi-scale, multi-physics modelling of the pulmonary system

The essential function of the lung, gas exchange, is dependent on adequate matching of ventilation and perfusion, where air and blood are delivered through complex branching systems exposed to regionally varying transpulmonary and transmural pressures. Structure and function in the lung are intimately related, yet computational models in pulmonary physiology usually simplify or neglect structure. The geometries of the airway and vascular systems and their interaction with parenchymal tissue have an important bearing on regional distributions of air and blood, and therefore on whole lung gas exchange, but this has not yet been addressed by modelling studies. Models for gas exchange have typically incorporated considerable detail at the level of chemical reactions, with little thought for the influence of structure. To date, relatively little attention has been paid to modelling at the cellular or subcellular level in the lung, or to linking information from the protein structure/interaction and cellular levels to the operation of the whole lung. We review previous work in developing anatomically based models of the lung, airways, parenchyma and pulmonary vasculature, and some functional studies in which these models have been used. Models for gas exchange at several spatial scales are briefly reviewed, and the challenges and benefits from modelling cellular function in the lung are discussed.

Simplified models for gas exchange in the human lungs

This paper presents a hierarchy of models with increasing complexity for gas exchange in the human lungs. The models span from a single compartment, inflexible lung to a single compartment, flexible lung with pulmonary gas exchange. It is shown how the models are related to well-known models in the literature. A long-term purpose of this work is to study nonlinear phenomena seen in the cardio-respiratory system (for example, synchronization between ventilation rate and heart rate, and Cheyne-Stokes respiration). The models developed in this paper can be regarded as the controlled system (plant) and provide a mathematical framework to link between "molecular-level", and "systems-level" models. It is shown how changes in molecular level affect the alveolar partial pressure. Two assumptions that have previously been made are re-examined: (1) the hidden assumption that the air flow through the mouth is equal to the rate of volume change in the lungs, and, (2) the assumption that the process of oxygen binding to hemoglobin is near equilibrium. Conditions under which these assumptions are valid are studied. All the parameters in the models, except two, are physiologically realistic. Numerical results are consistent with published experimental observations.

[Mathematical simulation of the respiratory system (author's transl)]

The respiratory system is described as a feedback control system. The controller consists of the peripheral chemoreceptors and the central chemosensitive structures, the respiratory centre in the medulla oblongata and the thorax-lung pump which they drive. The controlled system is comprised of three compartments (lung, brain and the remaining tissue) connected by the blood circulation. The controlled values are arterial pH and arterial O2 partial pressure and cerebral extracellular pH. Earlier models have been improved by: (1) the dead space description, (2) the thermodynamic formulation of the CO2 dissociation equation and the simple but accurate O2 dissociation equation of the blood, (3) the alteration of the CO2 dissociation equation for the brain and the remaining tissue to accommodate recent results, (4) the application of the one-receptor-theory of central chemosensitivity, (5) the pH dependence of brain circulation, (6) the bicarbonate exchange between blood and extracellular fluid of the brain and (7) the introduction of variable circulation times. Respiratory and metabolic disturbances of the respiratory system are analyzed. The mathematical formulation of the respiratory system is a differential difference equation system. In the steady state the experimental results are reproduced fairly well. A slight discrepancy is found in the simulation of metabolic acidosis. Apparently we have assumed the sensitivity of the peripheral chemoreceptors to be too large so that the respiratory response is not correctly predicted. In the numerical solution there is an overshoot in the on-transient and a damped oscillation in the off-transient of the alveolar CO2 partial pressure during respiratory acidosis. We have varied the parameters to make deviations small. The best agreement seems to result, if the central threshold is near the normal extracellular pH of the brain. A further deviation from experimental findings is that the cerebral CO2 and H+ concentration, the blood circulation of the brain, the alveolar O2 partial tension and the ventilation show a slight oscillation in the off-transient. Except for these discrepancies the experimental results, especially the stability of the extracellular pH of the brain, are reproduced fairly well. During hypoxia there are deviations form the experimental results if the central residual activity is constant and the central threshold deviates from the normal extracellular pH of the brain. But if the central residual activity is pH dependent and if the central threshold is equal to the normal extracellular pH of the brain, then the time course of VE and the other variables agree fairly well with experimental results. There is also a good correspondence between the theoretical and experimental data during hyperoxia. During metabolic acidosis the time constant of the bicarbonate exchange between blood and extracellular fluid of the brain is important. If a time constant of one minute is assumed, then the predicted and the experimental results correspond sufficiently well.