A 3D unsteady flow analysis in a doubly constricted arterial vessel

The aim of this study is to examine the interaction between two mild atherosclerotic proliferations spaced apart by a distance S by analyzing their influence on flow structure, pressure drop and stress field in an arterial vessel under pulsatile flow conditions. This has been achieved numerically by employing a time accurate, cell centered finite volume method in solving the Navier-Stokes equations governing the 3D unsteady flow dynamics in a conceptual model of an multiply constricted arterial vessel. In comparison to the pressure drop across a single stenosis, nearly a 50% increase in the late systolic and early diastolic pressure drops has been observed across the two mild constrictions when they are spaced within a distance of S相似文献   

Experimental analysis of the influence of stenotic geometry on steady flow.

We studied the flow behavior under steady flow conditions in four models of cylindrical stenoses at Reynolds numbers from 150 to 920. The flow upstream of the constrictions was always fully developed. The constriction ratios of the rigid tubes (D) to the stenoses (d) were d/D = 0.273; 0.505; 0.548; 0.786. The pressure drop at various locations in the stenotic models was measured with water manometers. The flow was visualized with a photoelasticity apparatus using an aqueous birefringent solution. We also studied the flow behavior at pulsatile flow in a dog aorta with a constriction of 71%. The flow through stenotic geometries depends on the Reynolds number of the flow generated in the tube and the constriction ratio d/D. At low d/D ratios, (with the increased constriction), the flow separation zones (recirculation zones, so-called reattachment length) and flow disturbances increased with larger Reynolds numbers. At lower values, eddies were generated. At high Re, eddies were observed in the pre-stenotic regions. The pressure drop is a function of the length and internal diameter of the stenosis, respective ratio of stenosis to the main vessel and the Reynolds numbers. At low Re-numbers and low d/D, distinct recirculation zones were found close to the stenosis. The flow is laminar in the distal areas. Further experiments under steady and unsteady flow conditions in a dog aorta model with a constriction of 71% showed similar effects. High velocity fluctuations downstream of the stenosis were found in the dog aorta. A videotape demonstrates these results.     

Estimation of valvular resistance of segmented aortic valves using computational fluid dynamics

Aortic valve stenosis is associated with an elevated left ventricular pressure and transaortic pressure drop. Clinicians routinely use Doppler ultrasound to quantify aortic valve stenosis severity by estimating this pressure drop from blood velocity. However, this method approximates the peak pressure drop, and is unable to quantify the partial pressure recovery distal to the valve. As pressure drops are flow dependent, it remains difficult to assess the true significance of a stenosis for low-flow low-gradient patients. Recent advances in segmentation techniques enable patient-specific Computational Fluid Dynamics (CFD) simulations of flow through the aortic valve. In this work a simulation framework is presented and used to analyze data of 18 patients. The ventricle and valve are reconstructed from 4D Computed Tomography imaging data. Ventricular motion is extracted from the medical images and used to model ventricular contraction and corresponding blood flow through the valve. Simplifications of the framework are assessed by introducing two simplified CFD models: a truncated time-dependent and a steady-state model. Model simplifications are justified for cases where the simulated pressure drop is above 10 mmHg. Furthermore, we propose a valve resistance index to quantify stenosis severity from simulation results. This index is compared to established metrics for clinical decision making, i.e. blood velocity and valve area. It is found that velocity measurements alone do not adequately reflect stenosis severity. This work demonstrates that combining 4D imaging data and CFD has the potential to provide a physiologically relevant diagnostic metric to quantify aortic valve stenosis severity.     

Passive movement of human soft palate during respiration: A simulation of 3D fluid/structure interaction

This study reconstructed a three dimensional fluid/structure interaction (FSI) model to investigate the compliance of human soft palate during calm respiration. Magnetic resonance imaging scans of a healthy male subject were obtained for model reconstruction of the upper airway and the soft palate. The fluid domain consists of nasal cavity, nasopharynx and oropharynx. The airflow in upper airway was assumed as laminar and incompressible. The soft palate was assumed as linear elastic. The interface between airway and soft palate was the FSI interface. Sinusoidal variation of velocity magnitude was applied at the oropharynx corresponding to ventilation rate of 7.5L/min. Simulations of fluid model in upper airway, FSI models with palatal Young's modulus of 7539Pa and 3000Pa were carried out for two cycles of respiration. The results showed that the integrated shear forces over the FSI interface were much smaller than integrated pressure forces in all the three directions (axial, coronal and sagittal). The total integrated force in sagittal direction was much smaller than that of coronal and axial directions. The soft palate was almost static during inspiration but moved towards the posterior pharyngeal wall during expiration. In conclusion, the displacement of human soft palate during respiration was mainly driven by air pressure around the surface of the soft palate with minimal contribution of shear stress of the upper airway flow. Despite inspirational negative pressure, expiratory posterior movement of soft palate could be another factor for the induction of airway collapse.     

Mathematical modelling of flow through an irregular arterial stenosis.

A mathematical model of flow through an irregular arterial stenosis is developed. The model is two-dimensional and axi-symmetric with the stenosis outline obtained from a three-dimensional casting of a mildly stenosed artery. Agreement between modelled and experimental pressure drops (obtained from an axi-symmetric machined stenosis with the same profile) is excellent. Results are also obtained for a smooth stenosis model, similar to that used for most mathematical modelling studies. This model overestimates the pressure drop across the stenosis, as well as the wall shear stress and separation Reynolds number. Also, the smooth model predicts one instead of three recirculation zones present in the irregular model. The original stenosis is modified to increase the severity from 48 and 87% areal occlusion, while maintaining the same general shape. This has the effect of increasing the pressure drop by an order of magnitude and decreasing the number of recirculation zones to one, with a lower separation Reynolds number.     

An intergral method for the analysis of flow in arterial stenoses

An approximate solution is presented to the problem of incompressible flow through an axisymmetric constriction. The geometry is intended to simulate an arterial stenosis, and the solution is applicable to both mild and severe stenoses for Reynolds numbers below transition. Theoretical results obtained for specific geometries are given for the velocity distribution, pressure drop, wall shearing stress, and separation phenomena. These results reveal the significant alterations in flow caused by a stenosis. Experiments using model stenoses are described and compared with the theoretical results. Theoretical predictions of pressure drop and separation characteristics are in reasonably good agreement with the experimental observations.     

A fluid--structure interaction finite element analysis of pulsatile blood flow through a compliant stenotic artery.

A new model is used to analyze the fully coupled problem of pulsatile blood flow through a compliant, axisymmetric stenotic artery using the finite element method. The model uses large displacement and large strain theory for the solid, and the full Navier-Stokes equations for the fluid. The effect of increasing area reduction on fluid dynamic and structural stresses is presented. Results show that pressure drop, peak wall shear stress, and maximum principal stress in the lesion all increase dramatically as the area reduction in the stenosis is increased from 51 to 89 percent. Further reductions in stenosis cross-sectional area, however, produce relatively little additional change in these parameters due to a concomitant reduction in flow rate caused by the losses in the constriction. Inner wall hoop stretch amplitude just distal to the stenosis also increases with increasing stenosis severity, as downstream pressures are reduced to a physiological minimum. The contraction of the artery distal to the stenosis generates a significant compressive stress on the downstream shoulder of the lesion. Dynamic narrowing of the stenosis is also seen, further augmenting area constriction at times of peak flow. Pressure drop results are found to compare well to an experimentally based theoretical curve, despite the assumption of laminar flow.     

Noninvasive estimation of pharyngeal airway resistance and compliance in children based on volume-gated dynamic MRI and computational fluid dynamics

Computational fluid dynamics (CFD) analysis was used to model the effect of collapsing airway geometry on internal pressure and velocity in the pharyngeal airway of three sedated children with obstructive sleep apnea syndrome (OSAS) and three control subjects. Model geometry was reconstructed from volume-gated magnetic resonance images during normal tidal breathing at 10 increments of tidal volume through the respiratory cycle. Each geometry was meshed with an unstructured grid and solved using a low-Reynolds number k-ω turbulence model driven by flow data averaged over 12 consecutive breathing cycles. Combining gated imaging with CFD modeling created a dynamic three-dimensional view of airway anatomy and mechanics, including the evolution of airway collapse and flow resistance and estimates of the local effective compliance. The upper airways of subjects with OSAS were generally much more compliant during tidal breathing. Compliance curves (pressure vs. cross-section area), derived for different locations along the airway, quantified local differences along the pharynx and between OSAS subjects. In one subject, the distal oropharynx was more compliant than the nasopharynx (1.028 vs. 0.450 mm(2)/Pa) and had a lower theoretical limiting flow rate, confirming the distal oropharynx as the flow-limiting segment of the airway in this subject. Another subject had a more compliant nasopharynx (0.053 mm(2)/Pa) during inspiration and apparent stiffening of the distal oropharynx (C = 0.0058 mm(2)/Pa), and the theoretical limiting flow rate indicated the nasopharynx as the flow-limiting segment. This new method may help to differentiate anatomical and functional factors in airway collapse.     

Maximal methacholine-induced constriction in rabbit lung: interactions between airways and tissue?

Six anesthetized paralyzed open-chest New Zealand White male rabbits were studied to obtain the maximal or plateau response to the inhalation of methacholine. Tracheal flow, tracheal pressure, and, by use of alveolar capsules, alveolar pressure were measured during tidal mechanical ventilation. We calculated total lung resistance (RL), tissue viscance (Vti), and lung elastance by digital fitting of the equation of motion to changes in tracheal and alveolar pressure. Airways resistance (Raw) was calculated as RL-Vti. Measurements were made under control conditions and after delivery of increasing concentrations of methacholine aerosol (0.5-128 mg/ml). We found that Vti accounted for the major proportion of RL both under control conditions (64.5 +/- 15.9%) and after methacholine-induced constriction (83.6 +/- 11.8%). There was a significant negative correlation between logarithmic percent change in Raw and Vti at the onset of the plateau response (r = 0.973). Furthermore, the slope of the relationship between log change in Vti and log change in Raw during the plateau response was strongly correlated with the degree of tissue response at the onset of the plateau (r = 0.957). Vti was positively correlated with lung elastance both before and during the plateau response (r = 0.946). We propose that the negative correlation between tissue resistance and Raw at the level of the plateau is consistent with a model of a mechanically interdependent lung, where decreases in airway caliber are limited by the constriction of the surrounding parenchyma.     

Computational fluid dynamics modeling of the upper airway of children with obstructive sleep apnea syndrome in steady flow

Computational fluid dynamic (CFD) analysis was used to model the effect of airway geometry on internal pressure in the upper airway of three children with obstructive sleep apnea syndrome (OSAS), and three controls. Model geometry was reconstructed from magnetic resonance images obtained during quiet tidal breathing, meshed with an unstructured grid, and solved at normative peak resting flow. The unsteady Reynolds-averaged Navier-Stokes equations were solved with steady flow boundary conditions in inspiration and expiration, using a two-equation low-Reynolds number turbulence model. Model results were validated using an in-vitro scale model, unsteady flow simulation, and reported nasal resistance measurements in children. Pharynx pressure drop strongly correlated to airway area restriction. Inspiratory pressure drop was primarily proportional to the square of flow, consistent with pressure losses due to convective acceleration caused by area restriction. On inspiration, in OSAS pressure drop occurred primarily between the choanae and the region where the adenoids overlap the tonsils (overlap region) due to airway narrowing, rather than in the nasal passages; in controls the majority of pressure drop was in the nasal passages. On expiration, in OSAS the majority of pressure drop occurred between the oropharynx (posterior to the tongue) and overlap region, and local minimum pressure in the overlap region was near atmospheric due to pressure recovery in the anterior nasopharynx. The results suggest that pharyngeal airway shape in children with OSAS significantly affects internal pressure distribution compared to nasal resistance. The model may also help explain regional dynamic airway narrowing during expiration.     

A threshold lung volume for optimal mechanical effects on upper airway airflow dynamics: studies in an anesthetized rabbit model

Increasing lung volume improves upper airway airflow dynamics via passive mechanisms such as reducing upper airway extraluminal tissue pressures (ETP) and increasing longitudinal tension via tracheal displacement. We hypothesized a threshold lung volume for optimal mechanical effects on upper airway airflow dynamics. Seven supine, anesthetized, spontaneously breathing New Zealand White rabbits were studied. Extrathoracic pressure was altered, and lung volume change, airflow, pharyngeal pressure, ETP laterally (ETPlat) and anteriorly (ETPant), tracheal displacement, and sternohyoid muscle activity (EMG%max) monitored. Airflow dynamics were quantified via peak inspiratory airflow, flow limitation upper airway resistance, and conductance. Every 10-ml lung volume increase resulted in caudal tracheal displacement of 2.1 ± 0.4 mm (mean ± SE), decreased ETPlat by 0.7 ± 0.3 cmH(2)O, increased peak inspiratory airflow of 22.8 ± 2.6% baseline (all P < 0.02), and no significant change in ETPant or EMG%max. Flow limitation was present in most rabbits at baseline, and abolished 15.7 ± 10.5 ml above baseline. Every 10-ml lung volume decrease resulted in cranial tracheal displacement of 2.6 ± 0.4 mm, increased ETPant by 0.9 ± 0.2 cmH(2)O, ETPlat was unchanged, increased EMG%max of 11.1 ± 0.3%, and a reduction in peak inspiratory airflow of 10.8 ± 1.0%baseline (all P < 0.01). Lung volume, resistance, and conductance relationships were described by exponential functions. In conclusion, increasing lung volume displaced the trachea caudally, reduced ETP, abolished flow limitation, but had little effect on resistance or conductance, whereas decreasing lung volume resulted in cranial tracheal displacement, increased ETP and increased resistance, and reduced conductance, and flow limitation persisted despite increased muscle activity. We conclude that there is a threshold for lung volume influences on upper airway airflow dynamics.     

Estimation of inspiratory pressure drop in neonatal and pediatric endotracheal tubes.

Endotracheal tubes (ETTs) constitute a resistive extra load for intubated patients. The ETT pressure drop (DeltaP(ETT)) is usually described by empirical equations that are specific to one ETT only. Our laboratory previously showed that, in adult ETTs, DeltaP(ETT) is given by the Blasius formula (F. Lofaso, B. Louis, L. Brochard, A. Harf, and D. Isabey. Am. Rev. Respir. Dis. 146: 974-979, 1992). Here, we also propose a general formulation for neonatal and pediatric ETTs on the basis of adimensional analysis of the pressure-flow relationship. Pressure and flow were directly measured in seven ETTs (internal diameter: 2.5-7.0 mm). The measured pressure drop was compared with the predicted drop given by general laws for a curved tube. In neonatal ETTs (2.5-3.5 mm) the flow regime is laminar. The DeltaP(ETT) can be estimated by the Ito formula, which replaces Poiseuille's law for curved tubes. For pediatric ETTs (4.0-7.0 mm), DeltaP(ETT) depends on the following flow regime: for laminar flow, it must be calculated by the Ito formula, and for turbulent flow, by the Blasius formula. Both formulas allow for ETT geometry and gas properties.     

Influence of arterial wall compliance on the pressure drop across coronary artery stenoses under hyperemic flow condition

Hemodynamic endpoints such as flow and pressure drop are often measured during angioplasty procedures to determine the functional severity of a coronary artery stenosis. There is a lack of knowledge regarding the influence of compliance of the arterial wall-stenosis on the pressure drop under hyperemic flows across coronary lesions. This study evaluates the influence in flow and pressure drop caused by variation in arterial-stenosis compliance for a wide range of stenosis severities. The flow and pressure drop were evaluated for three different severities of stenosis and tested for limiting scenarios of compliant models. The Mooney-Rivlin model defined the non-linear material properties of the arterial wall and the plaque regions. The non-Newtonian Carreau model was used to model the blood flow viscosity. The fluid (blood)-structure (arterial wall) interaction equations were solved numerically using the finite element method. Irrespective of the stenosis severity, the compliant models produced a lower pressure drop than the rigid artery due to compliance of the plaque region. A wide variation in the pressure drop was observed between different compliant models for significant (90% area occlusion) stenosis with 41.0, 32.1, and 29.8 mmHg for the rigid artery, compliant artery with calcified plaque, and compliant artery with smooth muscle cell proliferation, respectively. When compared with the rigid artery for significant stenosis the pressure drop decreased by 27.7% and 37.6% for the calcified plaque and for the smooth muscle cell proliferation case, respectively. These significant variations in pressure drop for the higher stenosis may lead to misinterpretation and misdiagnosis of the stenosis severity.     

Validation of computational fluid dynamics methodology used for human upper airway flow simulations

An anatomically accurate human upper airway model was constructed from multiple magnetic resonance imaging axial scans. This model was used to conduct detailed Computational Fluid Dynamics (CFD) simulations during expiration, to investigate the fluid flow in the airway regions where obstruction could occur. An identical physical model of the same airway was built using stereo lithography. Pressure and velocity measurements were conducted in the physical model. Both simulations and experiments were performed at a peak expiratory flow rate of 200 L/min. Several different numerical approaches within the FLUENT commercial software framework were used in the simulations; unsteady Large Eddy Simulation (LES), steady Reynolds-Averaged Navier-Stokes (RANS) with two-equation turbulence models (i.e. k?ε, standard k?ω, and k?ω Shear Stress Transport (SST)) and with one-equation Spalart–Allmaras model. The CFD predictions of the average wall static pressures at different locations along the airway wall were favorably compared with the experimental data. Among all the approaches, standard k?ω turbulence model resulted in the best agreement with the static pressure measurements, with an average error of ～20% over all ports. The highest positive pressures were observed in the retroglossal regions below the epiglottis, while the lowest negative pressures were recorded in the retropalatal region. The latter is a result of the airflow acceleration in the narrow retropalatal region. The largest pressure drop was observed at the tip of the soft palate. This location has the smallest cross section of the airway. The good agreement between the computations and the experimental results suggest that CFD simulations can be used to accurately compute aerodynamic flow characteristics of the upper airway.     

Biomechanical Aspects of Compliant Airways due to Mechanical Ventilation

Without proper knowledge of mechanical ventilation effects, physicians can aggravate an existing lung injury. A better understanding of the interaction between airflow and airway tissue during mechanical ventilation will be helpful to physicians so that they can provide appropriate ventilator parameters for intubated patients. In this study, a computational model incorporating the interactions between airflow and airway walls was developed to investigate the effects of airway tissue flexibility on airway pressure and stress. Two flow rates, 30 and 60 l/min, from mechanical ventilation were considered. The transient waveform was active inhalation with a constant flow rate and passive exhalation. Results showed that airway tissue flexibility decreased airway pressure at bifurcation sites by approximately 25.06% and 16.91% for 30 and 60 l/min, respectively, and increased wall shear stress (WSS) by approximately 74.00% and 174.91% for 30 and 60 l/min, respectively. The results from the present study suggested that it is very important to consider the interaction between airflow and airway walls when computational models are developed. Results of this study help to better quantify how the airflow rate used in mechanical ventilation, in conjunction with airway tissue flexibility, affects airway pressure and stresses.     

Noninvasive determination of upper airway resistance and flow limitation.

We have shown that a polynomial equation, FP = AP3 + BP2 + CP + D, where F is flow and P is pressure, can accurately determine the presence of inspiratory flow limitation (IFL). This equation requires the invasive measurement of supraglottic pressure. We hypothesized that a modification of the equation that substitutes time for pressure would be accurate for the detection of IFL and allow for the noninvasive measurement of upper airway resistance. The modified equation is Ft = At3 + Bt2 + Ct + D, where F is flow and t is time from the onset of inspiration. To test our hypotheses, data analysis was performed as follows on 440 randomly chosen breaths from 18 subjects. First, we performed linear regression and determined that there is a linear relationship between pressure and time in the upper airway (R2 0.96 +/- 0.05, slope 0.96 +/- 0.06), indicating that time can be a surrogate for pressure. Second, we performed curve fitting and found that polynomial equation accurately predicts the relationship between flow and time in the upper airway (R2 0.93 +/- 0.12, error fit 0.02 +/- 0.08). Third, we performed a sensitivity-specificity analysis comparing the mathematical determination of IFL to manual determination using a pressure-flow loop. Mathematical determination had both high sensitivity (96%) and specificity (99%). Fourth, we calculated the upper airway resistance using the polynomial equation and compared the measurement to the manually determined upper airway resistance (also from a pressure-flow loop) using Bland-Altman analysis. Mean difference between calculated and measured upper airway resistance was 0.0 cmH2O x l(-1) x s(-1) (95% confidence interval -0.2, 0.2) with upper and lower limits of agreement of 2.8 cmH2O x l(-1) x s(-1) and -2.8 cmH2O x l(-1) x s(-1). We conclude that a polynomial equation can be used to model the flow-time relationship, allowing for the objective and accurate determination of upper airway resistance and the presence of IFL.     

Pressure drops through arterial stenosis models in steady flow condition.

Measures of pressure drops were made in two different plexiglass models of axial-symmetric arterial stenoses. The stenosis models had the same area reduction (86 percent) but were of different length so as to have a different tapering degree. Pressures were measured in steady flow condition at three equidistant points of the stenosis: upstream, in the middle, and downstream. Results indicate that: the upstream-middle pressure drop is independent of tapering degree but is highly influenced by area reduction; moreover it is much greater than the middle-downstream drop. The upstream-middle pressure drop can be accurately predicted by means of a relationship deduced by the momentum equation.     

Upper airway collapsibility and patterns of flow limitation at constant end-expiratory lung volume

The passive pharyngeal critical closing pressure (Pcrit) is measured using a series of pressure drops. However, pressure drops also lower end-expiratory lung volume (EELV), which independently affects Pcrit. We describe a technique to measure Pcrit at a constant EELV. Continuous positive airway pressure (CPAP)-treated obstructive sleep apnea (OSA) patients and controls were instrumented with an epiglottic catheter, magnetometers (to measure change in EELV), and nasal mask/pneumotachograph and slept supine on nasal CPAP. Pcrit was measured in standard fashion and using our novel "biphasic technique" in which expiratory pressure only was lowered for 1 min before the inspiratory pressure was dropped; this allowed EELV to decrease to the drop level before performing the pressure drop. Seven OSA and three controls were studied. The biphasic technique successfully lowered EELV before the inspiratory pressure drop. Pcrit was similar between the standard and biphasic techniques (-0.4 ± 2.6 vs. -0.6 ± 2.3 cmH(2)O, respectively, P = 0.84). Interestingly, we noted three different patterns of flow limitation: 1) classic Starling resistor type: flow fixed and independent of downstream pressure; 2) negative effort dependence within breaths: substantial decrease in flow, sometimes with complete collapse, as downstream pressure decreased; and 3) and negative effort dependence across breaths: progressive reductions in peak flow as respiratory effort on successive breaths increased. Overall, EELV changes do not influence standard passive Pcrit measurements if breaths 3-5 of pressure drops are used. These results also highlight the importance of inspiratory collapse in OSA pathogenesis. The cause of negative effort dependence within and across breaths is not known and requires further study.     

Analysis of the collapsibility of the upper airway in a spectrum of sleep-disordered breathing: a modelling approach

Using a simplified model of the upper airways with two independent collapsible elements (nostrils and hypo-pharynx), we calculated the cross-sectional area of these two elements, taking into account pressure drops. We experimentally measured flow and pressure in the fossa and hypo-pharynx in various syndromes. This allowed us to compare the behaviour of the area supplied by our model with the aerodynamic resistance that is often used to analyse upper airway flow limitation events. We showed that nostril and hypo-pharyngeal areas are better correlated than the resistance values and thus concluded that the pressure divided by the square of the flow is a better parameter for analysing flow limitation in upper airways than resistance. Owing to its simplicity, our model is able to supply the area of the collapsible element in real time, which is impossible with more sophisticated models.     

Effect of lung volume on plateau response of airways and tissue to methacholine in dogs.

We have recently shown in dogs that much of the increase in lung resistance (RL) after induced constriction can be attributed to increases in tissue resistance, the pressure drop in phase with flow across the lung tissues (Rti). Rti is dependent on lung volume (VL) even after induced constriction. As maximal responses in RL to constrictor agonists can also be affected by changes in VL, we questioned whether changes in the plateau response with VL could be attributed in part to changes in the resistive properties of lung tissues. We studied the effect of changes in VL on RL, Rti, airway resistance (Raw), and lung elastance (EL) during maximal methacholine (MCh)-induced constriction in 8 anesthetized, paralyzed, open-chest mongrel dogs. We measured tracheal flow and pressure (Ptr) and alveolar pressure (PA), the latter using alveolar capsules, during tidal ventilation [positive end-expiratory pressure (PEEP) = 5.0 cmH2O, tidal volume = 15 ml/kg, frequency = 0.3 Hz]. Measurements were recorded at baseline and after the aerosolization of increasing concentrations of MCh until a clear plateau response had been achieved. VL was then altered by changing PEEP to 2.5, 7.5, and 10 cmH2O. RL changed only when PEEP was altered from 5 to 10 cmH2O (P < 0.01). EL changed when PEEP was changed from 5 to 7.5 and 5 to 10 cmH2O (P < 0.05). Rti and Raw varied significantly with all three maneuvers (P < 0.05). Our data demonstrate that the effects of VL on the plateau response reflect a complex combination of changes in tissue resistance, airway caliber, and lung recoil.