Sevoflurane consumption pattern by individual anaesthesiologists varies widely despite using the same high-end workstations in the same hospital.

Volatile anaesthetics are potent greenhouse gasses but contemporary workstations enable considerable savings while improving patient safety. Institutions may provide this technology to reduce the ecological footprint but proper training and motivation is required to maximize their ecologic and financial benefit. This study aims to compare the sevoflurane consumption of 22 anaesthesiologists in a medium sized hospital 4 years after flow-i workstations (Getinge, Sweden) entered into service, in three airway approaches: intubated patients, laryngeal mask ventilation, and mask anaesthesia. Typical sevoflurane consumption for each anaesthesiologist was defined as the mean cumulative consumption in the chronologically first 50 cases meeting the inclusion criteria for each airway group in 2019. The potential savings, if everyone were to adopt the approach of the more economical anaesthesiologists (15th percentile), was calculated. The CO2 equivalent emissions were calculated using a GWP20 of 702 and a GWP100 of 195. The median [range] consumption after 45 min was 10.9 [7.5-18.4] ml in intubated patients and 9.0 [7.4-15.3] ml in patients with laryngeal mask, and 9.9 [3.4-20.9] ml after 8 min with mask ventilation. This corresponds to a double to six fold consumption between the least and most wasteful approach. The typical CO2 equivalent emissions (GWP20) per anaesthesiologist varied between 8.0 and 19.6 kg/45 min in intubated airways, between 7.9 and 16.3 kg/45 min in LMA, and between 3.6 and 22.3 kg/8 min in mask ventilation. Despite using the same workstations in the same hospital, the typical sevoflurane consumption differed dramatically between 22 anaesthesiologists. In addition to providing advanced workstations, proper education is required to achieve the behavior change needed to reduce the pollution and financial waste associated with volatile anaesthetics.

Tomographic PIV in a model of the left ventricle: 3D flow past biological and mechanical heart valves.

Left ventricular flow is intrinsically complex, three-dimensional and unsteady. Its features are susceptible to cardiovascular pathology and treatment, in particular to surgical interventions involving the valves (mitral valve replacement). To improve our understanding of intraventricular fluid mechanics and the impact of various types of prosthetic valves thereon, we have developed a custom-designed versatile left ventricular phantom with anatomically realistic moving left ventricular membrane. A biological, a tilting disc and a bileaflet valve (in two different orientations) were mounted in the mitral position and tested under the same settings. To investigate 3D flow within the phantom, a four-view tomographic particle image velocimetry setup has been implemented. The results compare side-by-side the evolution of the 3D flow topology, vortical structures and kinetic energy in the left ventricle domain during the cardiac cycle. Except for the tilting disc valve, all tested prosthetic valves induced a crossed flow path, where the outflow crosses the inflow path, passing under the mitral valve. The biological valve shows a strong jet with a peak velocity about twice as high compared to all mechanical heart valves, which makes it easier to penetrate deeply into the cavity. Accordingly, the peak kinetic energy in the left ventricle in case of the biological valve is about four times higher than the mechanical heart valves. We conclude that the tomographic particle imaging velocimetry setup provides a useful ground truth measurement of flow features and allows a comparison of the effects of different valve types on left ventricular flow patterns.

Effective Orifice Area during Exercise in Bileaflet Mechanical Valve Prostheses.

BACKGROUND: The aims of this study were to investigate the evolution of the transprosthetic pressure gradient and effective orifice area (EOA) during dynamic bicycle exercise in bileaflet mechanical heart valves and to explore the relationship with exercise capacity. METHODS: Patients with bileaflet aortic valve replacement (n = 23) and mitral valve replacement (MVR; n = 16) prospectively underwent symptom-limited supine bicycle exercise testing with Doppler echocardiography and respiratory gas analysis. Transprosthetic flow rate, peak and mean transprosthetic gradient, EOA, and systolic pulmonary artery pressure were assessed at different stages of exercise. RESULTS: EOA at rest, midexercise, and peak exercise was 1.66 ± 0.23, 1.56 ± 0.30, and 1.61 ± 0.28 cm2, respectively (P = .004), in aortic valve replacement patients and 1.40 ± 0.21, 1.46 ± 0.27, and 1.48 ± 0.25 cm2, respectively (P = .160), in MVR patients. During exercise, the mean transprosthetic gradient and the square of transprosthetic flow rate were strongly correlated (r = 0.65 [P < .001] and r = 0.84 [P < .001] for aortic valve replacement and MVR, respectively), conforming to fundamental hydraulic principles for fixed orifices. Indexed EOA at rest was correlated with exercise capacity in MVR patients only (Spearman ρ = 0.68, P = .004). In the latter group, systolic pulmonary artery pressures during exercise were strongly correlated with the peak transmitral gradient (ρ = 0.72, P < .001). CONCLUSIONS: In bileaflet mechanical valve prostheses, there is no clinically relevant increase in EOA during dynamic exercise. Transprosthetic gradients during exercise closely adhere to the fundamental pressure-flow relationship. Indexed EOA at rest is a strong predictor of exercise capacity in MVR patients. This should be taken into account in therapeutic decision making and prosthesis selection in young and dynamic patients.

Automatic detection of oesophageal intubation based on ventilation pressure waveforms shows high sensitivity and specificity in patients with pulmonary disease.

BACKGROUND: Unrecognised endotracheal tube misplacement in emergency intubations has a reported incidence of up to 17%. Current detection methods have many limitations restricting their reliability and availability in these circumstances. There is therefore a clinical need for a device that is small enough to be practical in emergency situations and that can detect oesophageal intubation within seconds. In a first reported evaluation, we demonstrated an algorithm based on pressure waveform analysis, able to determine tube location with high reliability in healthy patients. The aim of this study was to validate the specificity of the algorithm in patients with abnormal pulmonary compliance, and to demonstrate the reliability of a newly developed small device that incorporates the technology. MATERIALS AND METHODS: Intubated patients with mild to moderate lung injury, admitted to intensive care were included in the study. The device was connected to the endotracheal tube, and three test ventilations were performed in each patient. All diagnostic data were recorded on PC for subsequent specificity/sensitivity analysis. RESULTS AND DISCUSSION: A total of 105 ventilations in 35 patients with lung injury were analysed. With the threshold D-value of 0.1, the system showed a 100% sensitivity and specificity to diagnose tube location. CONCLUSION: The algorithm retained its specificity in patients with decreased pulmonary compliance. We also demonstrated the feasibility to integrate sensors and diagnostic hardware in a small, portable hand-held device for convenient use in emergency situations.

Comparison of symmetric hemodialysis catheters using computational fluid dynamics.

PURPOSE: Symmetric-tip dialysis catheters have become alternative devices because of low access recirculation and ease of tip positioning. Flow characteristics of three symmetric catheters were compared based on computational fluid dynamics (CFD) as they relate to catheter function. MATERIALS AND METHODS: In Palindrome, GlidePath, and VectorFlow catheters, a computational fluid dynamics-based approach was used to assess (i) regions of flow separation, which are prone to thrombus development; (ii) shear-induced platelet activation potency; (iii) recirculation; and (iv) venous outflow deflection. A steady-state, laminar flow model simulated catheter tip position within the superior vena cava. Catheter performance was investigated at high hemodialysis flow rate (400 mL/min). Blood was assumed as a Newtonian fluid. RESULTS: Wide regions of flow separation downstream of the Palindrome side slot and close to the distal tip were observed in forward and reversed line configurations. Geometric asymmetry of the distal guide wire aperture of the GlidePath catheter produced the highest levels of inverted velocity flow when run in reversed configuration. The lowest mean shear-induced platelet activation was exhibited by GlidePath and VectorFlow catheters; the Palindrome catheter exhibited 152% higher overall platelet activation potency. All catheters were associated with a recirculation close to zero; the helically contoured lumens of the VectorFlow catheter produced the greatest amount of deflection of venous flow away from the arterial lumen. CONCLUSIONS: The VectorFlow catheter produced less shear-induced platelet activation than the Palindrome catheter and less flow separation than the Palindrome and GlidePath catheters irrespective of line configuration. These findings have potential implications for differences in thrombogenic risk during clinical performance of these catheters.

Functional respiratory imaging as a tool to personalize respiratory treatment in patients with unilateral diaphragmatic paralysis.

A completely different treatment approach was chosen for 2 patients with unilateral diaphragmatic paralysis and complaints of dyspnea despite similar anatomic and physiologic abnormalities. These decisions were supported by results obtained by functional respiratory imaging (FRI). FRI generated functional information on lobar ventilation and local drug deposition. In the first patient, some lobes were poorly ventilated, and drug deposition simulation showed that some regions were undertreated. This patient underwent diaphragmatic plication to restore ventilation. In the second patient, all lobes were still ventilated. A conservative approach with regular follow-ups was chosen to wait for spontaneous recovery of the diaphragmatic function. Both patients improved subjectively and objectively. These cases demonstrate how novel medical imaging techniques such as FRI can be used to personalize respiratory treatment in patients with unilateral diaphragmatic paralysis.

Change in upper airway geometry between upright and supine position during tidal nasal breathing.

BACKGROUND: As the upper airway is the most important limiting factor for the deposition of inhalation medication in the lower airways, it is interesting to assess how its morphology varies between different postures. The goal of this study is to compare the upper airway morphology and functionality of healthy volunteers in the upright and supine positions during tidal nasal breathing and to search for baseline indicators for these changes. This is done by performing three-dimensional measurements on computed tomography (CT) and cone beam computed tomography (CBCT) scans. METHODS: This prospective study was approved by all relevant institutional review boards. All patients gave their signed informed consent. In this study, 20 healthy volunteers (mean age, 62 years; age range, 37-78 years; mean body mass index, 29.26; body mass index range, 21.63-42.17; 16 men, 4 women) underwent a supine low-dose CT scan and an upright CBCT scan of the upper airway. The (local) average (Savg) and minimal (Smin) cross-sectional area, the position of the latter, the concavity, and the airway resistance were examined to determine if they changed from the upright to the supine position. If changes were found, baseline parameters were sought that were indicators for these differences. RESULTS: There were five dropouts due to movement artifacts in the CBCT scans. Savg and Smin were 9.76% and 26.90% larger, respectively, in the CBCT scan than in the CT scan, whereas the resistance decreased by 26.15% in the upright position. The Savg of the region between the hard palate and the bottom of the uvula increased the most (49.85%). In people with a high body mass index, this value changed the least. The airway resistance in men decreased more than in women. CONCLUSIONS: This study demonstrated that there are differences in upper airway morphology and functionality between the supine and upright positions and that there are baseline indicators for these differences.

Feasibility of a priori numerical assessment of plaque scaffolding after carotid artery stenting in clinical routine: proof of concept.

PURPOSE: Carotid artery stenting (CAS) is an alternative procedure for the treatment of severely stenosed carotid artery lesions in high-risk patients. Appropriate patient selection and stent design are paramount to achieve a low stroke and death rate in these complex high-risk procedures. This study introduces and evaluates a novel virtual, patient-specific, pre-operative environment to quantify scaffolding parameters based on routine imaging techniques. METHODS: Two patients who underwent CAS with two different sizes of the Acculink stent (Abbott Vascular, Santa Clara, CA, USA) were studied. Pre-operative data were used to build the numerical models for the virtual procedure. Numerical results were validated with post-operative angiography. Using novel virtual geometrical tools, incomplete stent apposition, free cell area and largest fitting sphere in the stent cell were evaluated in situ as quantitative measures of successful stent placement and to assess potential risk factors for CAS complications. RESULTS: A quantitative validation of the numerical outcome with post-operative images noted differences in lumen diameter of 5.31 ± 8.05% and 4.12 ± 9.84%, demonstrating the reliability of the proposed methodology. The quantitative measurements of the scaffolding parameters on the virtually deployed stent geometry highlight the variability of the device behavior in relation to the target lesion. The free cell area depends on the target diameter and oversizing, while the largest fitting spheres and apposition values are influenced by the local concavity and convexity of the vessel. CONCLUSIONS: The proposed virtual environment may be an additional tool for endovascular specialists especially in complex anatomical cases where stent design and positioning may have a higher impact on procedural success and outcome.

A computational study of the hemodynamic impact of open- versus closed-cell stent design in carotid artery stenting.

The aim of this study is to analyze the shape and flow changes of a patient-specific carotid artery after carotid artery stenting (CAS) performed using an open-cell (stent-O) or a closed-cell (stent-C) stent design. First, a stent reconstructed from micro-computed tomography (microCT) is virtually implanted in a left carotid artery reconstructed from CT angiography. Second, an objective analysis of the stent-to-vessel apposition is used to quantify the lumen cross-sectional area and the incomplete stent apposition (ISA). Third, the carotid artery lumen is virtually perfused in order to quantify its resistance to flow and its exposure to atherogenic or thrombogenic hemodynamic conditions. After CAS, the minimum cross-sectional area of the internal carotid artery (ICA) (external carotid artery [ECA]) changes by +54% (-12%) with stent-O and +78% (-17%) with stent-C; the resistance to flow of the ICA (ECA) changes by -21% (+13%) with stent-O and -26% (+18%) with stent-C. Both stent designs suffer from ISA but the malapposed stent area is larger with stent-O than stent-C (29.5 vs. 14.8 mm(2) ). The untreated vessel is not exposed to atherogenic flow conditions whereas an area of 67.6 mm(2) (104.9) occurs with stent-O (stent-C). The area of the stent surface exposed to thrombogenic risk is 5.42 mm(2) (7.7) with stent-O (stent-C). The computer simulations of stenting in a patient's carotid artery reveal a trade-off between cross-sectional size and flow resistance of the ICA (enlarged and circularized) and the ECA (narrowed and ovalized). Such a trade-off, together with malapposition, atherogenic risk, and thrombogenic risk is stent-design dependent.

Evaluation of alternatives for dysfunctional double lumen central venous catheters using a two-compartmental mathematical model for different solutes.

Double lumen (DL) central venous catheters (CVC) often suffer from thrombosis, fibrin sheet formation, and/or suction towards the vessel wall, resulting in insufficient blood flow during hemodialysis. Reversing the catheter connection often restores blood flows, but will lead to higher recirculation. Single lumen (SL) CVCs have often fewer flow problems, but they inherently have some degree of recirculation. To assist bedside clinical decision making on optimal catheter application, we investigated mathematically the differences in dialysis adequacy using different modes of access with CVCs.
A mathematical model was developed to calculate reduction ratio (RR) and total solute removal (TSR) of urea, methylguanidine (MG), beta-2-microglobulin (ß2M), and phosphate (P) during different dialysis scenarios: 4-h dialysis with a well-functioning DL CVC (DL-normal, blood flow QB 350 ml/min), dysfunctional DL CVC (DL-low flow, QB 250), reversed DL CVC (DL-reversed, QB 350, recirculation R = 10%) and 12 Fr SL CVC (effective QB 273). 
With DL-normal as reference, urea RR was decreased by 3.5% (DL-reversed), 13.0% (SL), and 15.6% (DL-low flow), while urea TSR was decreased by 3.3% (DL-reversed), 13.2% (SL), and 13.5% (DL-low flow). The same trend was found for MG and P. However, ß2M RR decreased only 1.5% with SL CVC although TSR decrease was 17.2%, while RR decreased 21.1% with DL-low flow although TSR decrease was only 4.9%.
In the case of dysfunctional DL CVCs, reversing the catheter connection and restoring the blood flow did not impair TSR, with 10% recirculation. The SL CVC showed suboptimal TSR results that were similar to those of the dysfunctional DL CVC.

The accuracy of ultrasound volume flow measurements in the complex flow setting of a forearm vascular access.

PURPOSE: Maturation of an arterio-venous fistula (AVF) frequently fails, with low post-operative fistula flow as a prognostic marker for this event. As pulsed wave Doppler (PWD) is commonly used to assess volume flow, we studied the accuracy of this measurement in the setting of a radio-cephalic AVF. METHODS: As in-vivo validation of fistula flow measurements is cumbersome, we performed simulations, integrating computational fluid dynamics with an ultrasound (US) simulator. Flow in the arm was calculated, based on a patient-specific model of the arm vasculature pre and post AVF creation. Raw ultrasound signals were subsequently simulated, from which Doppler spectra were calculated in both a proximal and a distal location. RESULTS: The velocity component in the direction of the PWD-US beam (vPWD), in a centered, small, sample volume, can be captured accurately using PWD spectrum mean-tracking (maximum bias [mB] 8.1%). However, when deriving flow rate from these measurements, a high degree of inaccuracy occurs. First, the angle-correction of vPWD towards the velocity along the axis of the vessel is largely influenced by the radial velocity components in the complex flow field (mB=16.3%). Second, the largest error is introduced when transferring the centerline velocity to the cross-sectional mean velocity without any knowledge of the flow profile (mB=97.7%). CONCLUSIONS: In the setting of a forearm AVF, flow estimates based on PWD are hampered by the complex flow patterns. Overall, flow estimation based on centerline measurement, analyzed by mean-tracking of the RF-spectral estimates, under the assumption of a parabolic flow profile, appeared to provide the most reasonable values.

Modeling hemodynamics in vascular networks using a geometrical multiscale approach: numerical aspects.

On the one hand the heterogeneity of the circulatory system requires the use of different models in its different compartments, featuring different assumptions on the spatial degrees of freedom. On the other hand, the mutual interactions between its compartments imply that these models should preferably not be considered separately. These requirements have led to the concept of geometrical multiscale modeling, where the main idea is to couple 3D models with reduced 1D and/or 0D models. As such detailed information on the flow field in a specific region of interest can be obtained while accounting for the global circulation. However, the combination of models with different mathematical features gives rise to many difficulties such as the assignment of boundary conditions at the interface between two models and the development of robust coupling algorithms, as the subproblems are usually solved in a partitioned way. This review aims to give an overview of the most important aspects concerning 3D-1D-0D coupled models. In addition, some applications are presented in order to illustrate the potentialities of these coupled models.

Slow extended nocturnal home hemodialysis shows superior adequacy compared to in-center dialysis: a mathematical analysis.

Extended nocturnal home hemodialysis has gained renewed interest. However, no removal data for single/double needle (lumen) (SL and DL, respectively) or for low/high blood flow in extended dialysis are available. Therefore, we studied dialysis adequacy in different nocturnal home hemodialysis strategies. Coupling a kinetic with a dialyzer model, we calculated a reduction ratio from pre- to post-dialysis (RR) and total solute removal (TSR) of urea, methylguanidine (MG), ß2-microglobulin, and phosphate. Simulations were done for dialysis with blood flow Qb350 ml/min (DL-4h), extended DL high flow with Qb350 (DL-HF-8h) and low flow with Qb175 (DL-LF-8h), and SL with Qb273 (SL-8h). Compared to DL-4h, TSR was 28-59% larger for DL-HF-8h. TSR was most increased for ß2-microglobulin (18%) with DL-LF-8h, and for MG (35%) with SL-8h. Furthermore, RRs were equal (DL-LF-8h), higher (SL-8h), and even more increased (DL-HF-8h) for all studied solutes. In the home setting, DL-LF-8h and SL-8h are safe and promising strategies.

Assessment of variables affecting late left ventricular flow propagation velocity.

UNLABELLED: Early colour M-mode flow propagation velocity (Vpe) in the left ventricle is a well-known non-invasive index for assessing left ventricular relaxation. However, the utility and determinants of late colour M-mode flow propagation (Vpa) have received little attention to date. Vpa as a representation of the left ventricular vortex travelling velocity during late filling could have a distinct role in differentiating potential subgroups in diastolic failure. The aim of the present study was to establish the normal values of late flow propagation in a healthy population of various ages (18-79 years), and to examine the general and echocardiographic variables that affect Vpa. METHODS: We studied 75 apparently healthy subjects (age range, 18-79 years; 38 women, 37 men) as part of an outpatient clinic check-up screening. General parameters were recorded, including age, gender, height, weight, blood pressure, and heart rate. In addition, conventional grey-scale M-mode, 2D, as well as colour M-mode, 2D, and pulsed wave (tissue) Doppler echocardiographic parameters were obtained in a single centre and using a single operator setting. Backward linear regression analysis (dependent variable: Vpa) was performed to find the optimal model, taking into account multicollinearity and maximum coefficient of determination (R2). Due to the heteroscedasticity of the collected data, a logarithmic transformation was used. In addition, separate linear backward regression analysis was performed for the male and female subgroups. RESULTS: Vpa values were 26-179 cm/s. The optimal regression model after elimination included the following variables: age (beta = 0.684, P < 0.001), height (beta = 0.521, P < 0.001), gender (beta = 0.343, P < 0.05), left ventricular Vpe (beta = 0.299, P < 0.01), left ventricular posterior systolic (M-mode) wall thickness (beta = 0.288, P < 0.01), interventricular septum thickness diastole (beta = 0.346, P < 0.005), transmitral Doppler E-wave deceleration time apical 4-chamber (beta = -0.297, P < 0.05), and tissue Doppler peak E-wave mitral annulus (beta = 0.459, P < 0.005). The total coefficient of determination (R2) for this model was 0.540 (P < 0.001); 0.673 (P < 0.001) for men and 0.645 (P < 0.001) for women. CONCLUSION: Vpa, representing left ventricular vortex travelling velocity during late filling, shows a large range of values in normal healthy subjects. It is mainly depending on age, gender and left ventricular mass. Moreover, substantially different determinants are found between men and women. Further study is required to explore these findings.

A fast strong coupling algorithm for the partitioned fluid-structure interaction simulation of BMHVs.

The numerical simulation of Bileaflet Mechanical Heart Valves (BMHVs) has gained strong interest in the last years, as a design and optimisation tool. In this paper, a strong coupling algorithm for the partitioned fluid-structure interaction simulation of a BMHV is presented. The convergence of the coupling iterations between the flow solver and the leaflet motion solver is accelerated by using the Jacobian with the derivatives of the pressure and viscous moments acting on the leaflets with respect to the leaflet accelerations. This Jacobian is numerically calculated from the coupling iterations. An error analysis is done to derive a criterion for the selection of useable coupling iterations. The algorithm is successfully tested for two 3D cases of a BMHV and a comparison is made with existing coupling schemes. It is observed that the developed coupling scheme outperforms these existing schemes in needed coupling iterations per time step and CPU time.

Computational flow dynamics and preclinical assessment of a novel hemodialysis catheter.

We examined flow characteristics and recirculation of a novel dialysis catheter with helically contoured lumens, which optimize the vectors of blood entering and leaving the device. Computational flow dynamics were used to measure shear stress, residence time (RT), Platelet Lysis Index (PLI), and recirculation of the study catheter, termed the VectorFlow catheter, compared with the Palindrome catheter. Bench and animal models were used to measure recirculation compared with other dialysis catheter designs. The VectorFlow catheter was associated with an 18% reduction in mean shear stress compared with the Palindrome catheter, at 13.6 and 16.6 Pa, respectively. Low PLI was seen with the VectorFlow (PLI = 0.019) and Palindrome (PLI = 0.015) catheters. The VectorFlow catheter was associated with a % RT >0.030 seconds of 1.2%, compared to 0.9% for the Palindrome. Calculated recirculation of the VectorFlow, Palindrome, and step-tip catheter designs in reversed position were 1.2 × 10(-4) %, 0.039%, and 13.6%, respectively. In vitro, the VectorFlow catheter had no detectable recirculation (0%) compared to 7.3-9.5% with the Palindrome at flow rates of 400-600 ml/minute. In vivo, the VectorFlow catheter had no detectable recirculation, compared to 6.7-12% for the Palindrome. The VectorFlow catheter was associated with an absence of detectable recirculation with favorable flow dynamics.

The upstream boundary condition influences the leaflet opening dynamics in the numerical FSI simulation of an aortic BMHV.

In this paper, the influence of the upstream boundary condition in the numerical simulation of an aortic bileaflet mechanical heart valve (BMHV) is studied. Three three-dimensional cases with different upstream boundary conditions are compared. The first case consists of a rigid straight tube with a velocity profile at its inlet. In the second case, the upstream geometry is a contracting left ventricle (LV), positioned symmetrically with respect to the valve. In the last case, the LV is positioned asymmetrical with respect to the valve. The cases are used to simulate the same three-dimensional BMHV. The change in time of the LV volume is calculated such that the flow rate through the valve is identical in each case. The opening dynamics of the BMHV are modelled using fluid-structure interaction. The simulations show that differences occur in the leaflet movement of the three cases. In particular, with the asymmetric LV, one of the leaflets impacts the blocking mechanism at its open position with a 34% higher velocity than when using the velocity profile, and with an 88% higher velocity than in the symmetric LV case. Therefore, when simulating such an impact, the upstream boundary condition needs to be chosen carefully.