The novel Mechanical Ventilator Milano for the COVID-19 pandemic.

This paper presents the Mechanical Ventilator Milano (MVM), a novel intensive therapy mechanical ventilator designed for rapid, large-scale, low-cost production for the COVID-19 pandemic. Free of moving mechanical parts and requiring only a source of compressed oxygen and medical air to operate, the MVM is designed to support the long-term invasive ventilation often required for COVID-19 patients and operates in pressure-regulated ventilation modes, which minimize the risk of furthering lung trauma. The MVM was extensively tested against ISO standards in the laboratory using a breathing simulator, with good agreement between input and measured breathing parameters and performing correctly in response to fault conditions and stability tests. The MVM has obtained Emergency Use Authorization by U.S. Food and Drug Administration (FDA) for use in healthcare settings during the COVID-19 pandemic and Health Canada Medical Device Authorization for Importation or Sale, under Interim Order for Use in Relation to COVID-19. Following these certifications, mass production is ongoing and distribution is under way in several countries. The MVM was designed, tested, prepared for certification, and mass produced in the space of a few months by a unique collaboration of respiratory healthcare professionals and experimental physicists, working with industrial partners, and is an excellent ventilator candidate for this pandemic anywhere in the world.

Beat to beat analysis of mechanical heart valves by means of return map.

Three mechanical heart valves (two bileaflet prostheses and a tilting one) were investigated in a basic hardware setup in order to evaluate with a hydrophone their opening and closing action in time and in amplitude of each beat. The recorded signal was then segmented into the series of cycles xi(t) having a temporal duration equal to the working period imposed on the valve. Two return maps were defined, in order to evaluate the degree of dispersion of the resulting scatter plot: (i) the amplitude map xi(t) versus xi+1(t); (ii) the delay map for the closure of the valve within each beat versus the successive ones. To evaluate the results obtained, two indices were proposed based on both the degree of dispersion and the deviation of the regression line of the resulting scatter plot with respect to the bisector of the map plane. The tilting disc valve showed a lower degree of dispersion, both in the amplitude signal and in the closure time delays, with respect to the other two bileaflet heart valves. The methodology proposed here could be regarded as an alternative non-invasive tool to investigate the dynamic behaviour of prosthetic heart valves, especially in the case of their suspected failure.

On the monodimensional approach to the estimation of the highest reynolds shear stress in a turbulent flow.

The measurement of the Reynolds stress tensor, or at least of some of its components, is a necessary step to assess if the turbulence associated with the flow near prosthetic devices can damage blood constituents. Because of the intrinsic three dimensionality of turbulence, in general, a three-component anemometer should be used to measure directly the components of the Reynolds stress tensor. However, this can be practically unfeasible, especially in vivo; therefore, it is interesting to investigate the possibility of characterizing the turbulent flows that may occur in the circulatory system with the monodimensional data that a less complete equipment (e.g., a pulsed ultrasound Doppler) can yield. From the general expression of the Reynolds stress tensor, the highest shear stress can be deduced, as well as the Reynolds normal stress in the main flow direction. The relation between these two quantities, which is an issue already addressed in previous works, can thus be rigorously formulated in terms of some characteristic parameters of the Reynolds stress tensor, the principal normal stresses and the angles that the directions that define them form with the main flow direction. An experimental verification of the ratio of the two above-mentioned quantitites for the flow across bileaflet valves, investigated by means of two-dimensional laser Doppler anemometry, will illustrate the limitations of the monodimensional approach estimating the maximum load on blood constituents.

A discussion on the threshold limit for hemolysis related to Reynolds shear stress.

Turbulence-related damage to blood is a major problem with the use of prosthetic devices, such as mechanical heart valves. An often-cited paper by Sallam and Hwang (1984). Biorheology 21, 783-797) quantified the threshold for hemolysis to be about 400 N m(-2), a value that has hitherto contributed to the evaluation of the potential dangerousness of a medical implantable device. We propose a discussion of the mentioned experiment, based on the application of stress analysis concepts to the original measurements: this is necessary to assess the peak turbulence shear stress value that could have been found in Sallam and Hwangs experiment, with a suitable orientation of the measurement axes. The result of our theoretical discussion is that the threshold value of 400 N m(-2) could probably be considerably underestimated: following this point of view, a 3-D stress analysis shows that the peak turbulence shear stress at the inception of hemolysis should be at least 600 N m(-2). This result, obtained on the basis of the study of RBCs' response to a turbulent environment, indicates that blood particles are probably more resistant to short-time shear stresses than it was thought.

The influence of the leaflets' curvature on the flow field in two bileaflet prosthetic heart valves.

A successful mechanical prosthetic heart valve design is the bileaflet valve, which has been implanted for the first time more than 20 years ago. A key feature of bileaflet valves is the geometry of the two leaflets, which can be very important in determining the flow field. Laser Doppler anemometry (LDA) was used to perform an accurate study of the velocity and turbulence shear stress peak values (TSS(max)) fields at four distances from the valve plane. TSS(max) is a relevant parameter to assess the risk of hemolysis and platelet activation associated to the implantation of a prosthetic device, continuously interacting with blood. Two bileaflet valves were tested: the St. Jude HP and the Sorin Bicarbon, of the same nominal size (19mm). The former has flat leaflets, whereas the latter's leaflets have a cylindrical surface. A high regime (CO: 6l/min) was imposed, in order to test the two valves at maximum Reynolds number and consequent turbulence generation. The flat-leaflet design of the St. Jude generates a TSS field constant with distance; on the contrary, the Bicarbon's shear stress field undergoes an evident development, with an unexpected central peak at a distance comparable to the valve's dimensions (21mm). The two bileaflet valves tested, although very similar in design, behave very differently as for their turbulence properties. In particular, the concept of curved wake leads to conclude that the curvature of the leaflets' surface must be identified as an important parameter, which deserves careful attention in PHV design and development.

Hydraulic functional characterisation of aortic mechanical heart valve prostheses through lumped-parameter modelling.

Lumped-parameter modelling techniques are proposed as a method for studying the hydraulic characteristics of mechanical prosthetic heart valves (PHVs). The global hydraulic behaviour of PHVs in the open position was modelled by taking into account the (nonlinear) resistive and (linear) inertial factors governing the time-dependent relationship between transvalvular pressure drop and fluid flow rate, and neglecting the leaflets' opening and closure transient phenomena. Statistically defined indices associated to the parameters' values attest how properly the model describes PHV hydraulic behaviour. Local fluid dynamics is not modelled with this approach. The proposed method was implemented in a software program and applied to the characterisation of the aortic StJude Medical, StJude Medical Hemodynamic Plus and CarboMedics PHVs, basing on steady- and pulsatile-flow hydraulic-bench experimental data. The results showed that reliable parameters expressing hydraulic resistance can be derived from steady-flow data (R(2)>0.995). Inertance parameters derived from pulsatile-flow experiments are liable to a degree of uncertainty (confidence intervals up to 17%), however, comparing the reconstructed vs. measured pressure drop during systolic time demonstrates that this deficiency is mostly due to the missing description of initial, transient oscillations presumably related to the leaflets' opening (not modelled).

Flow on the symmetry plane of a total cavo-pulmonary connection.

The flow inside a total cavo-pulmonary connection, a bypass operation of the right heart adopted in the presence of congenital malformation, is here studied for a specific geometry which has been recently introduced in clinics. The analysis has been performed by preliminary experimental observation and a novel Navier-Stokes formulation on the symmetry plane. This method, once some basic hypotheses are verified, allows to reproduce the flow on the symmetry plane of a three-dimensional field by using an extension of the two-dimensional approach. The analysis has confirmed the existence of a central vortex showing that it is not a real vortex (i.e. a place with accumulation of vorticity) but, rather, a weakly dissipative recirculating zone. It is surrounded by a shear layer that becomes spontaneously unsteady at moderately high Reynolds number. The topological changes and energy dissipation have been analysed in both cases of unbalanced and of balanced pulmonary artery and caval flows.

Computational model of the fluid dynamics of a cannula inserted in a vessel: incidence of the presence of side holes in blood flow.

Vascular access methods, performed by the insertion of cannulae into vessels, may disturb the physiological flow of blood, giving rise to non-physiological pressure variations and shear stresses. To date, the hydrodynamic behaviour of the cannulae has been evaluated comparing their pressure loss-flow rate relationships, as obtained from in vitro experiments using a monodimensional approach; this methodology neither furnish information about the local fluid dynamics nor the established flow field in specific clinical work conditions. Since the shear stress is a critical factor in the design of artificial circulatory devices, more knowledge should be necessary about the local values assumed by the haemodynamic parameters during cannulation. An alternative way to investigate the fluid dynamic as accurately as possible is given by numeric studies. A 3D model of cannula concentrically placed in a rigid wall vessel is presented, with the finite element methodology used to numerically simulate the steady-state flow field in two different venous cannulation case studies, with two cannulae having a central hole and two or four side holes, respectively, with the same boundary conditions. Lower velocity and shear stress peak values have been computed for the model with four side holes upstream of the central hole, in the region of the cannula where the inlet flows meet and towards cannula's outlet, due to the increased flow symmetry and inlet area with respect to the model with two side holes. Starting from the investigation of different cannula designs, numerically assessing the local fluid dynamics, indications can be drawn to support both the design phase and the device optimal clinical use, in order to limit risks of biomechanical origin. Thus the presence of four side holes implied, as a consequence of the greater inlet area and of the increased symmetry, a less disturbed blood flow, together with reduced shear stress values. Furthermore, results show that the numerical simulations furnished useful informations on the interaction between vessel and cannula, e.g. on the fluid dynamics establishing in the free luminal space left, in the vessel, by the inserted cannula.

19 mm sized bileaflet valve prostheses' flow field investigated by bidimensional laser Doppler anemometry (part I: velocity profiles).

The investigation of the flow field downstream of a cardiac valve prosthesis is a well established task. In particular turbulence generation is of interest if damage to blood constituents is to be assessed. Several prosthetic valve flow studies are available in literature but they generally concern large-sized prostheses. The FDA draft guidance requires the study of the maximum Reynolds number conditions for a cardiac valve model to assess the worst case in turbulence by choosing both the minimum valve diameter and a high cardiac output value as protocol set up. Within the framework of a national research project regarding the characterization of cardiovascular endoprostheses, the Laboratory of Biomedical Engineering is currently conducting an in-depth study of turbulence generated downstream of bileaflet cardiac valves. Four models of 19 mm sized bileaflet valve prostheses, namely St Jude Medical HP Edwards Tekna, Sorin Bicarbon, and CarboMedics, were studied in aortic position. The prostheses were selected for the nominal annulus diameter reported by the manufacturers without any assessment of the valve sizing method. The hemodynamic function was investigated using a bidimensional LDA system. Results concern velocity profiles during the peak flow systolic phase, at high cardiac output regime, highlighting the different flow field features downstream of the four small-sized cardiac valves.

19 mm sized bileaflet valve prostheses' flow field investigated by bidimensional laser Doppler anemometry (part II: maximum turbulent shear stresses)

The investigation of the flow field generated by cardiac valve prostheses is a necessary task to gain knowledge on the possible relationship between turbulence-derived stresses and the hemolytic and thrombogenic complications in patients after valve replacement. The study of turbulence flows downstream of cardiac prostheses, in literature, especially concerns large-sized prostheses with a variable flow regime from very low up to 6 L/min. The Food and Drug Administration draft guidance requires the study of the minimum prosthetic size at a high cardiac output to reach the maximum Reynolds number conditions. Within the framework of a national research project regarding the characterization of cardiovascular endoprostheses, an in-depth study of turbulence generated downstream of bileaflet cardiac valves is currently under way at the Laboratory of Biomedical Engineering of the Istituto Superiore di Sanita. Four models of 19 mm bileaflet valve prostheses were used: St Jude Medical HP, Edwards Tekna, Sorin Bicarbon, and CarboMedics. The prostheses were selected for the nominal Tissue Annulus Diameter as reported by manufacturers without any assessment of valve sizing method, and were mounted in aortic position. The aortic geometry was scaled for 19 mm prostheses using angiographic data. The turbulence-derived shear stresses were investigated very close to the valve (0.35 D0), using a bidimensional Laser Doppler anemometry system and applying the Principal Stress Analysis. Results concern typical turbulence quantities during a 50 ms window at peak flow in the systolic phase. Conclusions are drawn regarding the turbulence associated to valve design features, as well as the possible damage to blood constituents.

The role of usable length in the compliance measurement of vascular prostheses.

In this study we investigated the dependence of the mechanical properties and in particular of the radial compliance of a vascular prosthesis as a function of its usable length. Radial compliance was measured at 60 bpm and in the pressure range 80-120 mmHg. Starting from compliance measurements a simple model was used to calculate the pulse wave velocity and the reflection coefficients between 6 mm and 8 mm grafts (knitted and woven) with iliac and subclavean artery of similar diameter. The results provide an indication of the influence of usable length on the compliance and diameter mismatch at the anasthomosis between graft and host artery.

Morphological analysis of in vivo velocity field in the alteration of the vasomotor tone.

Vessel wall remodeling is involved in atherogenesis and in several important vascular diseases affecting mainly aged and prosthetic implanted patients. This adaptive response to pathological states in arterial hemodynamics strongly suggests that flow-derived stresses act as mechanical stimuli to the release of endothelium-derived vasoactive factors, leading to vascular alterations. As the correlation of intimal hyperplasia (IH) with blood flow alterations in arteries has been shown to be significant, and as it is well-known that clinical procedures carry a substantial risk of development of vascular disease, the relevance of local hemodynamics must be investigated to describe changes in compliance matching in prosthetic applications. The aim of our research is to investigate the use of principal components analysis, together with varimax rotation, in the individuation process of morphological characteristics of real time ultrasound in in vivo recordings of blood flow velocities, as provided by two different carotid perivascular manipulations. This would be of use in the clinical assessment of atherogenesis, hypertension, prosthetic replacement or more in general in all applications in which vascular tone may be impaired. Data recordings refer to previous animal experiments where the Moncada model was investigated by means of an ultrasound profilometer. The present study confirms the feasibility of the proposed analysis to follow vascular pathology evolution, distiguishing between an in progress and a static situation.

Pathological patient in protocol definition for bench testing of mechanical cardiac support system.

Clinical techniques for the restoration of a failing heart are mainly based on the use of mechanical assist devices. In recent years, with the growing need for mechanical circulatory support, these devices have been shown to be a useful therapeutic tool, thanks to their intrinsic capability to unload the failing ventricle, allowing the heart to recover. Mechanical circulatory support systems (MCSS) require an accurate biomechanical characterization of the complex interaction that occurs between the patient and the mechanical support. A protocol for MCSS testing is proposed which takes into account several working conditions, in a modified test mock loop apparatus able to mimic various pathological conditions. Both physiological and pathological conditions can be replicated to show the actual efficacy of a MCSS device in correctly supporting a wide spectrum of ventricular conditions. The test bench is able to simulate the recovery of the pathological condition quite accurately, showing, at the same time, that this set up can be a reliable choice to characterize cardiac support devices. Thus the results of this experimentation can be useful to clinicians in forecasting the response of the heart affected by a cardiac disease and to set appropriate parameters for suitable assistance.

Proposal for a quantitative description of blood spiral flow in medical devices.

The association between specific blood flow patterns and blood behaviour through medical devices suggests that a Lagrangian study may be a useful instrument for the evaluation of the thrombogenic and/or hemolytic potential of certain devices' geometries and biomaterials. In this study a description of blood particle trajectories in terms of their spiral contents is proposed; such a mathematical description for blood spiral flow, computed along several pathlines, is tested for a quantitative determination of the spiralled motion of blood flow into two three-dimensional numerical models, having different design characteristics, of venous cannula inserted in a vessel. As the influence of vortical flow conditions have been observed to have both beneficial and detrimental influence on blood behaviour in terms of blood-device interaction, of the degradation of its components, and of the efficiency of mass-exchange (in red cells oxygenation and plasma filtration, for example), the herein proposed method for the description of spiral laminar motion may be a helpful instrument to build up a tool to investigate, for example, the existence of correlations between level of spiral flow and geometry (as in the present investigated test case), rather than the effects of blood-surface contact. The results obtained in this test case investigation, confirm the effectiveness of the proposed function for a quantitative analysis of spiral flow in medical devices.

Endovascular stents: market vigilance and risk factors.

With the aim of enhancing the safety and reliability level of coronary stents, we analyzed data collected from accident reports drawn from the MAUDE database (Manufacturer and User Facility Device Experience Database) of the FDA from 1996 to 2000. This analysis allowed us to highlight problems related to the use of coronary stents by means of the analysis of these reports at different levels, beginning from the causes that can lead to a certain type of accident up to the possible complication related to that event. Moreover we analyzed the procedure outcomes in terms of stent position inside the patient's body and the possible therapies adopted to solve the problems. The results showed that the most probable event that can lead to an accident is the stent separation from the balloon which, alone, turns up in a number of cases equal to the sum of all the others. This result highlights the importance of the technical skill of the operators accomplished by special training and of the importance of clarity and completeness in the instructions for the use of the device. Another critical point is the reliability of the device which must guarantee an adequate safety level when it is used according to the instructions.

Potential mechanical blood trauma in vascular access devices: a comparison of case studies.

Since vascular access devices may cause disturbances in blood flow, possibly damaging red blood cells (RBCs), the correlated risk of lysis must be assessed. The monodimensional approach for the evaluation of cannulae hydrodynamic behaviour (in vitro measured flow curves) does not furnish information on the local flow field occurring in specific clinical conditions. Researchers consider the prediction of blood trauma, induced by mechanical loading, to optimize the design phase, and to furnish indications on their optimal clinical use. In this study, a model of cannula inserted in a non compliant wall vessel was used as a test bench in a Computational Fluid Dynamics (CFD) problem. By means of CFD the flow field was 3D analysed to achieve information on velocity and shear stress local values, when cannula is used for inflow and outflow cannulation. A prediction of potential blood corpuscle damage, based on a power law, quantified the potential blood damage. Several numerical simulations, with different cannula/vessel flow rate ratios were provided, to investigate the incidence of local sites in the design on blood damaging potential during cannulation. Several regions appeared to be sensitive to the flow rate not only inside the cannula but also in the space between cannula and vessel, suggesting new indications for the assessment of a quality factor based on the evaluation of induced blood cells injury.

Pulsatile flow and atherogenesis: results from in vivo studies.

Compliance mismatch between prosthetic vascular replacement (possibly stented) and native artery is considered to be an important factor in implant failure due, e.g., to vascular remodeling, tissutal growth or intimal hyperplasia (IH). From an in vivo study involving altered vascular mechanics (and, consequently, compliance mismatch), carried out using the Moncada model of atherosclerosis development and smooth muscle cell (SMC) proliferation, the hemodynamic assessment was followed by means of real-time multigated ultrasound profilometry, of collared carotid artery using two different models: non-constrictive and constrictive plastic collars, wrapped around the vessel. The experiments provided the real-time measurement of velocity profiles in vivo and the subsequent estimation of wall shear stresses, locally responsible for the altered hemodynamics. Endothelium modifications were correlated with local hemodynamic alterations by using statistical regression analysis of the development of intimal hyperplasia and the mechanical stimulus applied to the endothelium by means of the two different manipulation models. Different correlations were found between wall shear rate and IH in the two models, showing the importance of the vascular pulsatility in determining SMC proliferation. This result could be useful in minimizing the negative consequences of clinical interventions such as graft and/or stent implantation.

Numerical simulation of a realistic total cavo-pulmonary connection: effect of unbalanced pulmonary resistances on hydrodynamic performance.

Total cavo pulmonary connection (TCPC) is one of the surgical techniques adopted to compensate the failure of the right heart in pediatric patients. The main goal of this procedure is the realization of a configuration for the caval veins and for the pulmonary arteries that can guarantee as low as possible pressure losses and appropriate lung perfusion. Starting from this point of view, a realistic TCPC with extracardiac conduit (TECPC) is investigated by means of Computational Fluid Dynamics (CFD) to evaluate the pressure loss under different pressure conditions, simulating different vessel resistances, on the pulmonary arteries. A total flow of 3 L/min, with a distribution between the inferior vena cava (IVC) and the superior vena cava (SVC) equal to 6/4, was investigated; three different boundary conditions for the pressure were imposed, resulting in three simulations in steady-state conditions, to the right pulmonary artery (RPA) and to the left pulmonary artery (LPA), simulating a balanced (deltaP(LPA-RPA) = 0 mmHg) and two unbalanced pulmonary resistances to blood flow (a pressure difference deltaP(LPA-RPA) = +/- 2 mmHg, respectively). The geometry for the TECPC was realized according to MRI derived physiological values for the vessels and for the configuration adopted for the anastomosis (the extra-cardiac conduit was inclined 22 degrees towards the left pulmonary artery with respect to the IVC axis). The computed power losses agree with previous in vitro Particle Image Velocimetry investigations. The results show that a higher resistance on the LPA causes the greater pressure loss for the TECPC under study, while the minimum pressure loss can be achieved balancing the pulmonary resistances, subsequently obtaining a balanced flow repartition towards the lungs.

A mathematical model of the fetal cardiovascular system based on genetic algorithms as identification technique.

The development of fetal cardiac surgery, considered the ultimate goal in the treatment of congenital cardiac malformations, needs to be supported by detailed knowledge of the blood circulation in the fetal cardiovascular system. The hemodynamic behavior in distal territories is usually inferred from vessel resistance indices, which give limited physiological information. This study presents a mathematical model of the human fetal global cardiovascular system, developed to clarify the relationships and differences existing between upper and lower body circulation. We modelled the heart with two time-varying capacitances, each representing the respective ventricle's pressure-volume relationship. The fetal vascular system was represented using two six-element Windkessel models, for the upper and lower body respectively. We obtained the identification of the set of circuital and elastance function parameters of the model using Genetic Algorithms (GAs), which follow the laws of evolutionary theory. We compared the results of our numerical study on the model identified with data collected from measurements and literature, to validate the proposed global cardiovascular system model of the human fetus. This model is intended as an instrument to investigate the differences in blood distribution between the different vascular districts in the upper and lower fetal body and the role of the aortic isthmus, the small tract of vessel connecting upper and lower fetal vascular beds; it may also represent a useful tool in the assessment of dynamic balance during mechanical assistance of circulation.

Extracorporeal circulation in ewe's foetus: towards a reliable foetal cardiac surgery protocol. A comparison of two cases.

UNLABELLED: Foetal cardiac surgery is the ultimate goal in the treatment of congenital cardiac malformations. The aim of our research is to elucidate some of the features of the necessarily invasive experimental protocol to be used in an animal model of foetal cardiac surgery. In particular, we assessed the foetal placental reactivity to prolonged cardiac bypass in steady-flow conditions. METHODS: Two cases were selected to show the outcome of prolonged (> 30 minutes) extracorporeal circulation (ECC) instituted without oxygenator under steady-flow assistance. Following the instrumentation of the animal (placement of pressure, flow and myocardial fiber length transducers) and the baseline recordings, a 60-minute bypass period was established with an axial turbopump (Hemopump 14 Fr), after systemic heparinisation and artero-venous cannulation. At the end of the circulatory assistance, the cannulae were removed and a 90 minute observation period followed. The cardiac function was assessed by means of indirectly obtained P-V loops. RESULTS: Case A showed a marked reduction in the end-systolic pressure-volume relationship (ESPVR) during ECC, corresponding to a rightward shift of the P-V loop, with a gradual recovery after the assisted circulation. On the contrary, case B was subjected to progressive placental dysfunction, as evidenced by haemogasanalytical data. Consequently, the haemodynamic data also outlined a negative outcome, with high ESPVR values after bypass. CONCLUSIONS: The present study, while confirming the possibility of cardiac intervention in the foetus, underlines the critical role of minimally invasive protocol to limit both foetal stress and placental dysfunction.