- Open Access
Virtual respiratory system in investigation of CPAP influence on optimal breathing frequency in obstructive lungs disease
© Golczewski and Darowski; licensee BioMed Central Ltd. 2007
Received: 20 June 2007
Accepted: 16 July 2007
Published: 16 July 2007
Continuous Positive Airway Pressure (CPAP) is a commonly accepted method of spontaneous breathing support in obstructive lung disease. Previous work suggested that the cause of the CPAP efficacy in the obstructive lung disease localized in bronchi of middle order (OLDMO) is not as obvious as, for example, in the obstructive sleep apnea. Since CPAP reduces obstruction and the optimal breathing frequency (BF) depends on the obstruction level, it seems to be important to analyze the dependence of the optimal BF on CPAP.
To analyze the support efficacy cause in OLDMO, esp. the relationship between the CPAP value and optimal BF.
Investigations utilized previously built virtual respiratory system. Its most important factors: nonlinear lungs compliance and changeability of nonlinear airway resistance (Raw). Influence of BF and the CPAP value on the tidal volume and minute ventilation was analyzed for four exemplary virtual patients: healthy ("standard") and suffering from moderate, severe, and the very severe OLDMO (the other parameters, esp. respiratory muscles effort, were unchanged). Minute inspiratory work as a criterion of the BF optimization.
CPAP decreased Raw making breathing easier, however, it shifted the working point of the respiratory system towards the smaller lungs compliance making breathing harder. The final result depended on the Raw value: CPAP improved breathing of patients with the serious OLDMO while it worsened healthy person breathing. The optimal CPAP value depended on the Raw value. If a virtual patient suffering from the serious OLDMO was not supported with CPAP, he had to breathe with low frequency because minute ventilation did not rise with BF increase. The optimal BF depended on the CPAP value (the greater the value, the greater the frequency).
The CPAP efficacy depends on the level of OLDMO. CPAP is efficient in the severe OLDMO because it increases the optimal BF, which makes possible less energy-consuming breathing with frequency close to the normal one (greater BF means smaller tidal volume and thus smaller work against lungs compliance).
To avoid ethical and economic problems connected with experiments and investigations on animals or human beings (patients), the Virtual Respiratory System (VRS) has been recently developed in the Institute of Biocybernetics and Biomedical Engineering in Warsaw [1, 2]. Testing and analyses of various artificial ventilation and ventilatory support methods were original applications of VRS [3, 4].
The paper presents analysis of efficacy of spontaneous breathing support by means of the Continuous Positive Airway Pressure (CPAP) in patients suffering from obstructive lung disease. Despite one name, three types of such disease may be distinguished. The type depends on the obstruction localization. If the obstruction concerns the upper airway (e.g. as in the sleep apnea), the cause of the CPAP efficacy seems to be clear: CPAP keeps the upper airway open. Sullivan et al used CPAP in sleep apnea as the first  and nasal CPAP is a standard therapy for obstructive sleep apnea syndrome, now. If the obstruction that concerns smallest bronchi is connected with atelectasis, then the efficacy cause seems to be clear, too: CPAP prevents alveoli from collapse.
However, if the obstruction concerns the bronchi of the middle order (e.g. asthma), the efficacy cause is not so obvious. Although easier inspiration is usually supposed to be the efficacy cause, our previous investigation with VRS  has shown that the expiration period is much more critical in this case. It is because both the obstruction and transmural pressure (pressure inside bronchi minus intrapleural pressure) influence airflow through such bronchi: the smaller the transmural pressure because of increased intrapleural pressure during the expiration, the greater the summarized influence of the obstruction and the pressure on the airflow. Since intrapleural pressure rise depends on the duration of the expiration, and thus on the breathing frequency, the frequency has to be taken into account in an analysis of the CPAP efficacy.
Summarizing: this paper presents detailed analysis of the CPAP efficacy in the case of patients with obstruction localized in the bronchi of middle generations. The analysis concerns dependence of the tidal volume on breathing frequency, obstruction level, and CPAP, when force developed by respiratory muscles is constant.
2.1. Virtual respiratory system
(a) Resistances of the upper airways (the main bronchi, the trachea, etc.).
(b) Resistances that depend on lungs volume (Rvi in Fig. 1) – the resistance of smallest bronchi (i.e. of highest generations), which are component of the lung tissue, and thus their size depends on the current lungs volume.
(c) Resistances that depend on transmural pressure (Rpi in Fig. 1) – the resistance of bronchi of the middle generations, the diameter of which depends on the transmural pressure, i.e. on the difference between pressure inside the bronchi and the intrapleural pressure. Since both the intrapleural pressure and pressure inside the bronchi change during respiration, the bronchi diameter (thus also Rpi) changes, too. In particular, if the transmural pressure is negative during expiration, such bronchi collapse and the airflow does not depend on the pressure gradient. Thus, the simplified idea of "resistance" as the coefficient of proportionality between the flow Q and the pressure gradient Δp (Δp = Q•R) loses the sense. Therefore, the direct dependence between the airflow and the pressures is used in VRS, in fact (the way of the dependence description has been explained in , whereas the final formula is presented in Appendix). Nevertheless, the term "resistance" is used in the paper because of common practice.
In the text below, the "Raw = norm" means the airways of the standard person, i.e. the normal Raw. However, since changeable Rpi and Rvi are the most meaningful components of Raw, it is impossible to characterize this "normal resistance" with simple numeric value(s). Therefore, the "Raw = norm" means that the airflow depends on the pressures as for the standard person. Consequently, an obstruction described below with the "Raw = X·norm" means that the airflows are "X" times smaller than they would be for the standard person under the same instantaneous conditions. Since the analyzed obstruction concerns bronchi of the middle generation, resistance of which is described with the formula A.5 (see Appendix), Raw = X·norm means that the coefficient k in the formula A.5 is X times greater than for the standard person. In this paper, results for exemplary virtual patients, for which X = 1, 4, 8, and 16, have been analyzed. To connect this mathematical description with commonly used clinical classification, spirometry was performed . According to European Respiratory Society and American Thoracic Society criteria, simulated disease with Raw = 4·norm should be classified as a moderate obstruction while with Raw = 8·norm and Raw = 16·norm as severe and very severe obstructions.
2.2. Experimental procedures
Investigation on VRS enables analysis of a problem with the "step by step" method, i.e. it is possible to change a single factor to observe results, then to change another one to observe further results, etc. Such an approach enables easier determination which one(s) from the set of initially postulated factors really influences an analyzed phenomenon. The "step by step" method is usually impossible in the case of living models because of feedback between all factors and parameters.
This paper regards the reason of efficacy of ventilatory support with CPAP in the obstructive disease, when the obstruction concerns the bronchi that depend on the transmural pressure, i.e. when Raw is increased because of increased Rpi. Thus, breathing improvement with CPAP was the analyzed phenomenon. A change of the tidal volume (VT) was treated as an index of the improvement. The values of Raw and/or the breathing frequency were changed to analyze the dependence of VT on CPAP. The other factors remained constant. Especially, the respiratory muscle effort during particular inspirations was unchanged.
Finally, influence of CPAP on minute ventilation for different breathing frequency was analyzed to find the optimal frequency and the CPAP value.
The muscle effort has to be well defined to be used in quantitative analysis, esp. in a computer model. Generally, oxygen consumption is commonly accepted as the best index of the effort of living being. However, such a chemical index cannot be used in simulations of mechanics, which is the matter of this paper. The mechanical work done could be a good index in the case of mechanical (physical) devices. Unfortunately, the physical work is proportional to the movement (e.g. the work is equal to the force times displacement or pressure times volume change) while there is no simple connection between the movement and the organism effort. Isovolumetric or isometric efforts can be examples. Moreover, such a work cannot be a presettable parameter of the simulations because it may be determined with integration of time-related pressures and volumes, i.e. it may be determined only after the simulation.
3.1. CPAP and respiratory system nonlinearity
The value of Raw strongly depends on the bronchi diameter (according to the Poiseuille's law for the ideal tube, the resistance is inversely proportional to the fourth power of the radius). Therefore, since the diameter depends on the transmural pressure (the difference between pressures inside and outside the bronchus), CPAP influences Raw: the greater the CPAP value, the smaller the Raw value. Hence it appears that CPAP may make breathing easier because a smaller part of the force developed by the respiratory muscles has to be used against Raw (which is illustrated by narrow hystereses for the greater CPAP in Fig. 3).
On the other hand, however, CPAP may make breathing harder. Indeed, if the airflow is equal to zero, as at the inspiration or expiration end, then the pressure inside the lungs is equal to the pressure inside the airways. In particular, the pressure inside the lungs at the inspiration beginning is equal to CPAP, which means that the starting lungs volume is greater than would be for the pressure equal to zero. Since the lungs compliance is nonlinear (see the formula A.3), CPAP moves the working point towards the smaller differential compliance (Fig. 2). As the consequence, the same force developed by the respiratory muscles causes smaller VT (the formula 2). Therefore, the greater the CPAP value, the more horizontal the hystereses as shown in Fig. 3.
In the case of the significantly increased Raw, the use of CPAP gives contrary result: CPAP increases VT (Fig. 4 – the results for Raw = 8·norm and Raw = 12·norm), and thus CPAP makes breathing easier. As it will be shown below, a relatively small fall of VT caused by a differential compliance decrease is compensated with a surplus by the fact that the lungs can empty themselves through the Raw decreased by CPAP.
3.2. CPAP and breathing frequency
Summarizing: in the obstructive lung disease, CPAP improves breathing with the normal frequency but it is unprofitable in the case of breathing with low frequency (Fig. 6).
As Fig. 5 shows, the maximal (i.e. end-inspiratory) volume is almost independent from the breathing frequency. Therefore, the frequency influence on VT is mainly caused by its influence on RC. Hence it appears that expiration rather than inspiration is critical in the analyzed obstructive lungs disease. However, although expiration is critical directly, its meaning is connected with its influence on inspiration, esp. its influence on the starting volume (i.e. RC), which determines the working point. For that reason, inspiration, especially the inspiratory work, has to be analyzed, too. In particular, such frequency and CPAP, for which the work is minimal, may be treated as the optimal ones.
The real work of the respiratory system (measured with oxygen consumption) and its mechanical work W (described with the formula dW = p·dV, where: p-pressure, V-volume) are not exactly the same (especially in the case of isovolumetric processes). Nevertheless, changes of the mechanical work may usually estimate changes of the real work.
Let us assume initially, that both Raw and the respiratory system compliance C are constant (the system is linear) and a patient needs ventilation of V liters per the time unit T.
Thus, since , the more shallow and frequent breathing is less energy-consuming. Since the mean airflow is approximately the same in the both cases (it is equal to V/T), the work against Raw is approximately the same, too (it has to be stressed that Raw during inspiration is smaller than during expiration – see the formula A.5). Hence it appears that differences in the work against the elasticity are significant.
The above consideration shows that deep breathing with low frequency is more energy-consuming than more frequent breathing. Since the compliance decreases according to lung volume rise, deep breaths are much more energy-consuming, in fact. Hence it appears that since breathing is one of the most energy-consuming physiological processes, the organism should breathe as frequently and shallowly as it is possible under actual conditions (certainly, since VT has to be greater than the dead space to ventilate alveoli, the breath deepness cannot be too small). To determine the optimal frequency for a particular CPAP value, the inspiratory work against the lungs compliance for constant minute ventilation was approximated for the different frequencies. It was assumed that: (a) the work is proportional to squared value of Ps (as in the formula 3), and (b) the lungs compliance depends on the working point but it does not depend on the breath deepness.
It is commonly known that patients suffering from lung obstructive disease breathe with low frequency. Although it is breathing in the energy-consuming way, they have no choice. As the simulations showed (Fig. 7), if there is no support with CPAP (i.e. CPAP = 0), rise of the frequency (above 8/min for Raw = 8·norm) is pointless because such a rise does not increase the minute ventilation. That seems to be the true reason of the fact that breathing with low frequency is the well-known natural defense mechanism in the case of patients suffering from the obstructive lungs disease. Our explanation, which bases on the simulations and stresses problems during the expiration, is different than the commonly accepted but rather intuitive explanation, which stresses the influence of the work against Raw during the inspiration . If CPAP is used, the above natural defense mechanism is not necessary.
1. It is commonly accepted now that a breathing frequency decrease in obstructive lungs disease is connected with an increase of the inspiratory work against Raw. Our simulations, however, suggest that problems of the expiration rather than the inspiration are the main cause of such frequency fall. The simulations indicated that the minute ventilation rises with the breathing frequency increase only if the frequency is smaller than certain boundary value because of difficulties with lungs emptying during too short expiration. The emptying through the increased Raw needs more time than through the normal Raw, which causes that in the severe obstructive disease the boundary value is smaller than the normal breathing frequency. Therefore, a patient has to breathe slowly since more frequent breathing would not increase the ventilation but it would increase the minute inspiratory work.
2. The simulations suggest that the CPAP efficacy is connected mainly with possibility of effective breathing with the frequency, which is close to the normal one. Therefore, CPAP seems to be helpful only when Raw is increased significantly, i.e. when the frequency that is economic without the CPAP support is smaller than the normal one.
3. CPAP makes breathing easier because it decreases required inspiratory work against the reduced Raw (which is commonly known) as well as the lungs compliance since: (a) the work against the compliance is approximately proportional to the squared tidal volume, (b) the tidal volume is approximately inversely proportional to the breathing frequency, and (c) CPAP enables more frequent breathing.
4. The optimal breathing frequency depends on both the CPAP value and the obstruction level.
Fig. 1 presents the general structure of the computer model that is postulated to be VRS. The main features of the model are:
1) nonlinearity of the meaningful model elements;
2) separation of the lungs and chest;
3) division of the lungs into lobes (i = 1...5 identifies the upper left lobe ....the bottom right lobe, respectively);
4) division of the airway resistance into:
(a) resistances that depend on the lung volume (Rvi in Fig. 1) – resistances of the smallest bronchi that are a component of the lung tissue;
(b) resistances that depend on the transmural pressure (Rpi in Fig. 1) – resistances of the bronchi that may collapse; and
(c) resistances of the large bronchi;
5) many others such as influence of gravity (Gi), air compressibility (Cci), etc.
Data for the model for the "standard" (healthy) human being have been collected on the basis of accessible literature. Some of the model elements in Fig. 1 represent parameters, which are described by single numbers. Such parameters are: Ru, R, RL, Rr, Rti, Rw, and L, LL, Lr – (resistances of mouth, trachea, left and right main bronchi, lung tissue, chest tissue, and inertance of trachea, left and right main bronchi, respectively). R, RL, Rr, L, LL, Lr were calculated from dimensions of trachea or main bronchi. Rti and Rw were estimated on the basis of data from .
Nonlinear model elements are described in the following way:
1. The air compressibility 'compliance' symbolized by Cci in Fig. 1, i.e. the dependence between the alveolar pressure (Patm+Pai), the lobe volume Vai, and the volume (Vai+Vci) that the air from the lobe would occupy at the atmospheric pressure (Patm):
Pai = Patm·Vci/Vai (A1)
where: Ptw – the trans-wall pressure, Vw – the chest volume, coefficients C1, C2, and C3 depend on simulated patient (C1 = 0.3 kPa/l, C2 = 1.722 kPa·l1/2, C3 = 1.5 l for the standard patient). The first part of the formula represents differential linearity of the chest-wall compliance for greater volumes, whereas the second part describes the limit of chest compression.
where: Ptpsi – the static trans-pulmonary (recoil) pressure concerning the lobe of number i, Vai – the volume of this lobe, C4 and C5 depend on simulated patient (C4 = -2.76, C5 = 0.66 l-1 for all i, for the standard homogenous lungs), ui – coefficient determining what part of the whole lung is this particular lobe (u = [0.2, 0.25, 0.15, 0.2, 0.2]). Note that Vai·C5/ui = V·ui·C5/ui = V·C5 for all i (where V- the total lungs volume). It means that it has been assumed initially that each lobe of the homogenous lungs has the same elastic properties as the whole lung.
The formulas (A.2) and (A.3) have been chosen because they well approximate characteristic dependence of the pressures on the volumes. Numeric values of the parameters C1, ..., C5 fit the graphic dependence presented in [9, 10].
4. The resistance Rvi that depends on the lobe volume Vai (aRVi – a proportionality coefficient, ccVai – critical closing volume) is described with the formula:
Rvi = aRVi/(Vai-ccVai) (A4)
Details, esp. the derivation of the formula (as well as explanations of other formulas and model verification) has been presented precisely in .
Coefficients bi determine how quickly the resistance falls with transmural pressure increase (the smaller bi, the quicker the resistance fall). ki describe the resistance for the infinite (theoretically) transmural pressure. Thus, bi are approximately connected with mechanical bronchi wall properties, while ki reflect bronchi internal dimensions. Such mathematical description of Rpi enables to simulate – with one formula – both the airflow auto-limitation employed in spirometry and the commonly known, experimental dependence of the airway resistance on lungs volume.
If a disease is examined, appropriate model parameters are adequately changed to simulate this disease. For example, increase of C4 simulates restrictive lung disease connected with the lung compliance fall, whereas increase of C5 may simulate restrictive lung disease because of the lung (lobe) capacity decrease.
- Darowski M, Kozarski M, Golczewski T: Model studies on respiratory parameters for different lung structures. Biocybernetics and Biomedical Engineering. 2000, 20: 67-77.Google Scholar
- Golczewski T, Kozarski M, Darowski M: The respirator as a user of virtual lungs. Biocybernetics and Biomedical Engineering. 2003, 23 (2): 57-66.Google Scholar
- Golczewski T, Darowski M: Influence of ventilatory mode on respiration parameters – investigation on virtual lungs. Biocybernetics and Biomedical Engineering. 2003, 23 (3): 63-72.Google Scholar
- Darowski M, Gólczewski T, Michnikowski M: Choice of proper lung ventilation method. Biocybernetics and Biomedical Engineering. 2006, 26 (1): 21-37.Google Scholar
- Sullivan CE, Issa FG, Berthon-Jones M, Eves L: Reversal of obstructive sleep apnoea by continuous positive airway pressure applied through the nares. Lancet. 1981, 1: 862-865. 10.1016/S0140-6736(81)92140-1.View ArticleGoogle Scholar
- Golczewski T, Darowski M: Virtual respiratory system for education and research: simulation of expiratory flow limitation for spirometry. Int J Artif Organs. 2006, 29: 961-972.Google Scholar
- Bates JHT, Milic-Emili J: Influence of the viscoelastic properties of the respiratory system on the energetically optimum breathing frequency. Annals of Biomedical Engineering. 1993, 21: 489-499. 10.1007/BF02584331.View ArticleGoogle Scholar
- Tomalak W, Peslin R, Duvivier C: Respiratory tissue properties derived from flow transfer function in healthy humans. J Appl Physiol. 1997, 82 (4): 1098-1106.Google Scholar
- Mead J, Turner JM, Macklem PT, Little JB: Significance of the relationship between lung recoil and maximum expiratory flow. J Appl Physiol. 1967, 22: 95-108.Google Scholar
- Nunn JF: Applied Respiratory Physiology. 1987, Butterworth-Heinemann, London [etc.], 3Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.