A Validated Cfd Model to Predict O 2 and Co 2 Transfer Within Hollow Fiber Membrane Oxygenators

Computational Pulsatile Flow and Efficiency Analysis of Biocompatible Microfluidic Artificial Lungs for Different Fiber Configurations

Average-sized microfluidic artificial lungs consisting of rows and columns of fiber bundles with different column to row aspect ratios (AR) are numerically analyzed for flow characteristics, maximum gas transfer performance, minimum pressure drop, and proper wall shear stress (WSS) values in terms of biocompatibility. The flow is fully laminar and assumed to be incompressible and Newtonian. The transport analysis is performed using a combined convection-diffusion model, and the numerical simulations are carried out with the finite element method. The inlet volumetric flow is modeled as a sinusoidal wave function to simulate the cardiac cycle and its effect on the device performance. The model is first validated with experimental studies in steady-state condition and compared with existing correlations for transient conditions. Then, the validated model is used for a parametric study in both steady and pulsatile flow conditions. The results show that increasing the aspect ratio in fiber configuration leads to converging gas transfer, higher pressure drop, and higher WSS. While determining the optimum configuration, the acceptable shear stress levels play a decisive role to ensure biocompatibility. Also, it is observed that the steady analysis underestimates the gas transfer for higher aspect ratios.

Computational Analysis of the Effects of Fiber Deformation on the Microstructure and Permeability of Blood Oxygenator Bundles

Mechanical loads on the polymeric fibers of oxygenating bundles are commonly present due to bundle press-fitting during device assembly and blood pressure load. However, computational fluid dynamics (CFD) simulations for fiber bundle optimization neglect possible changes in microstructure due to such deformations. The aim of this study is to investigate the impact of fiber deformability on bundle microstructure and fluid dynamics mainly in terms of permeability. Fibers from commercial mats typically used for blood oxygenators were mechanically tested and based on these experimental data, a material model was developed to simulate the structural deformations the fibers undergo under press-fitting and blood pressure loads. Then, CFD simulations were performed on deformed bundle repetitive units to investigate permeability under varying loading conditions. The effects of different bundle geometric parameters on the variation of bundle permeability due to press-fitting were evaluated. Bundle press-fitting results in significant changes in microstructure that are reflected in a bundle permeability more than halved for a 15% press-fitting. This impact on permeability is present in all the simulated fiber bundles and becomes more pronounced as the pitch between fibers and thus bundle porosity decreases. Instead, the analyses on pressurized bundle show only small deformations caused by pressure load, with permeability changes below 1%. While blood pressure effects could be neglected, bundle press-fitting turns out to have a significant impact on bundle microstructure and permeability. Neglecting such microstructure variations during CFD simulations could also lead to incorrect assessment of the local fluid dynamics within the bundle.

Comparison of hemodynamic features and thrombosis risk of membrane oxygenators with different structures: A numerical study

PURPOSE

The geometric structure of the membrane oxygenator can exert an impact on its hemodynamic features, which contribute to the development of thrombosis, thereby affecting the clinical efficacy of ECMO treatment. The purpose of this study is to investigate the impact of varying geometric structures on hemodynamic features and thrombosis risk of membrane oxygenators with different designs.

METHODS

Five oxygenator models with different structures, including different number and location of blood inlet and outlet, as well as variations in blood flow path, were established for investigation. These models are referred to as Model 1 (Quadrox-i Adult Oxygenator), Model 2 (HLS Module Advanced 7.0 Oxygenator), Model 3 (Nautilus ECMO Oxygenator), Model 4 (OxiaACF Oxygenator) and Model 5 (New design oxygenator). The hemodynamic features of these models were numerically analyzed using the Euler method combined with computational fluid dynamics (CFD). The accumulated residence time (ART) and coagulation factor concentrations (C[i], where i represents different coagulation factors) were calculated by solving the convection diffusion equation. The resulting relationships between these factors and the development of thrombosis in the oxygenator were then investigated.

RESULTS

Our results show that the geometric structure of the membrane oxygenator, including the location of the blood inlet and outlet as well as the design of the flow path, has a significant impact on the hemodynamic surroundings within the oxygenator. In comparison to Model 4, which had the inlet and outlet located in the center position, Model 1 and Model 3, which had the inlet and outlet at the edge of the blood flow field, exhibited a more uneven distribution of blood flow within the oxygenator, particularly in areas distant from the inlet and outlet, which was accompanied with lower flow velocity and higher values of ART and C[i], leading to the formation of flow dead zones and an elevated risk of thrombosis. The oxygenator of Model 5 is designed with a structure that features multiple inlets and outlets, which greatly improves the hemodynamic environment inside the oxygenator. This results in a more even distribution of blood flow within the oxygenator, reducing areas with high values of ART and C[i], and ultimately lowering the risk of thrombosis. The oxygenator of Model 3 with circular flow path section shows better hemodynamic performance compared to the oxygenator of Model 1 with square circular flow path. The overall ranking of hemodynamic performance for all five oxygenators is as follows: Model 5 > Model 4 > Model 2 > Model 3 > Model 1, indicating that Model 1 has the highest thrombosis risk while Model 5 has the lowest.

CONCLUSION

The study reveals that the different structures can affect the hemodynamic characteristics inside membrane oxygenators. The design of multiple inlets and outlets can improve the hemodynamic performance and reduce the thrombosis risk in membrane oxygenators. These findings of this study can be used to guide the optimization design of membrane oxygenators for improving hemodynamic surroundings and reducing thrombosis risk.

Three-dimensional membranes for artificial lungs: Comparison of flow-induced hemolysis

BACKGROUND

Membranes based on triply periodic minimal surfaces (TPMS) have proven a superior gas transfer compared to the contemporary hollow fiber membrane (HFM) design in artificial lungs. The improved oxygen transfer is attributed to disrupting the laminar boundary layer adjacent to the membrane surface known as main limiting factor to mass transport. However, it requires experimental proof that this improvement is not at the expense of greater damage to the blood. Hence, the aim of this work is a valid statement regarding the structure-dependent hemolytic behavior of TPMS structures compared to the current HFM design.

METHODS

Hemolysis tests were performed on structure samples of three different kind of TPMS-based designs (Schwarz-P, Schwarz-D and Schoen's Gyroid) in direct comparison to a hollow fiber structure as reference.

RESULTS

The results of this study suggest that the difference in hemolysis between TPMS membranes compared to HFMs is small although slightly increased for the TPMS membranes. There is no significant difference between the TPMS structures and the hollow fiber design. Nevertheless, the ratio between the achieved additional oxygen transfer and the additional hemolysis favors the TPMS-based membrane shapes.

CONCLUSION

TPMS-shaped membranes offer a safe way to improve gas transfer in artificial lungs.

Numerical modeling of pulsatile blood flow through a mini-oxygenator in artificial lungs

While previous in vitro studies showed divergent results concerning the influence of pulsatile blood flow on oxygen advection in oxygenators, no study was done to investigate the uncertainty affected by blood flow dynamics. The aim of this study is to utilize a computational fluid dynamics model to clarify the debate concerning the influence of pulsatile blood flow on the oxygen transport. The computer model is based on a validated 2D finite volume approach that predicts oxygen transfer in pulsatile blood flow passing through a 300-micron hollow-fiber membrane bundle with a length of 254 mm, a building block for an artificial lung device. In this study, the flow parameters include the steady Reynolds number (Re = 2, 5, 10 and 20), Womersley parameter (Wo = 0.29, 0.38 and 0.53) and sinusoidal amplitude (A = 0.25, 0.5 and 0.75). Specifically, the computer model is extended to verify, for the first time, the previously measured O₂ transport that was observed to be hindered by pulsating flow in the Biolung, developed by Michigan Critical Care Consultants. A comprehensive analysis is carried out on computed profiles and fields of oxygen partial pressure (P_O2) and oxygen saturation (S_O2) as a function of Re, Wo and A. Based on the present results, we observe the positive and negative effects of pulsatile flow on P_O2 at different blood flow rates. Besides, the S_O2 variation is not much influenced by the pulsatile flow conditions investigated. While being consistent with a recent experimental study, the computed O₂ volume flow rate is found to be increased at high blood flow rates operated with low frequency and high amplitude. Furthermore, the present study qualitatively explains that divergent outcomes reported in previous in vitro experimental studies could be owing to the different blood flow rates adopted. Finally, the contour analysis reveals how the spatial distributions of P_O2 and S_O2 vary over time.

Suitable CO2 Solubility Models for Determination of the CO2 Removal Performance of Oxygenators

CO₂ removal via membrane oxygenators during lung protective ventilation has become a reliable clinical technique. For further optimization of oxygenators, accurate prediction of the CO₂ removal rate is necessary. It can either be determined by measuring the CO₂ content in the exhaust gas of the oxygenator (sweep flow-based) or using blood gas analyzer data and a CO₂ solubility model (blood-based). In this study, we determined the CO₂ removal rate of a prototype oxygenator utilizing both methods in in vitro trials with bovine and in vivo trials with porcine blood. While the sweep flow-based method is reliably accurate, the blood-based method depends on the accuracy of the solubility model. In this work, we quantified performances of four different solubility models by calculating the deviation of the CO₂ removal rates determined by both methods. Obtained data suggest that the simplest model (Loeppky) performs better than the more complex ones (May, Siggaard-Anderson, and Zierenberg). The models of May, Siggaard-Anderson, and Zierenberg show a significantly better performance for in vitro bovine blood data than for in vivo porcine blood data. Furthermore, the suitability of the Loeppky model parameters for bovine blood (in vitro) and porcine blood (in vivo) is evaluated.

Toward a Long-Term Artificial Lung

Only a very small portion of end-stage organ failures can be treated by transplantation because of the shortage of donor organs. Although artificial long-term organ support such as ventricular assist devices provide therapeutic options serving as a bridge-to-transplantation or destination therapy for end-stage heart failure, suitable long-term artificial lung systems are still at an early stage of development. Although a short-term use of an extracorporeal lung support is feasible today, the currently available technical solutions do not permit the long-term use of lung replacement systems in terms of an implantable artificial lung. This is currently limited by a variety of factors: biocompatibility problems lead to clot formation within the system, especially in areas with unphysiological flow conditions. In addition, proteins, cells, and fibrin are deposited on the membranes, decreasing gas exchange performance and thus, limiting long-term use. Coordinated basic and translational scientific research to solve these problems is therefore necessary to enable the long-term use and implantation of an artificial lung. Strategies for improving the biocompatibility of foreign surfaces, for new anticoagulation regimes, for optimization of gas and blood flow, and for miniaturization of these systems must be found. These strategies must be validated by in vitro and in vivo tests, which remain to be developed. In addition, the influence of long-term support on the pathophysiology must be considered. These challenges require well-connected interdisciplinary teams from the natural and material sciences, engineering, and medicine, which take the necessary steps toward the development of an artificial implantable lung.

Three-dimensional computational model of a blood oxygenator reconstructed from micro-CT scans

Cardiopulmonary bypass procedures are one of the most common operations and blood oxygenators are the centre piece for the heart-lung machines. Blood oxygenators have been tested as entire devices but intricate details on the flow field inside the oxygenators remain unknown. In this study, a novel method is presented to analyse the flow field inside oxygenators based on micro Computed Tomography (μCT) scans. Two Hollow Fibre Membrane (HFM) oxygenator prototypes were scanned and three-dimensional full scale models that capture the device-specific fibre distributions are set up for computational fluid dynamics analysis. The blood flow through the oxygenator is modelled as a non-Newtonian fluid. The results were compared against the flow solution through an ideal fibre distribution and show the importance of a uniform distribution of fibres and that the oxygenators analysed are not susceptible to flow directionality as mass flow versus area remain the same. However the pressure drop across the oxygenator is dependent on flow rate and direction. By comparing residence time of blood against the time frame to fully saturate blood with oxygen we highlight the potential of this method as design optimisation tool. In conclusion, image-based reconstruction is found to be a feasible route to assess oxygenator performance through flow modelling. It offers the possibility to review a product as manufactured rather than as designed, which is a valuable insight as a precursor to the approval processes. Finally, the flow analysis presented may be extended, at computational cost, to include species transport in further studies.

Computational evaluation of the thrombogenic potential of a hollow-fiber oxygenator with integrated heat exchanger during extracorporeal circulation

The onset of thromboembolic phenomena in blood oxygenators, even in the presence of adequate anticoagulant strategies, is a relevant concern during extracorporeal circulation (ECC). For this reason, the evaluation of the thrombogenic potential associated with extracorporeal membrane oxygenators should play a critical role into the preclinical design process of these devices. This study extends the use of computational fluid dynamics simulations to guide the hemodynamic design optimization of oxygenators and evaluate their thrombogenic potential during ECC. The computational analysis accounted for both macro- (i.e., vortex formation) and micro-scale (i.e., flow-induced platelet activation) phenomena affecting the performances of a hollow-fiber membrane oxygenator with integrated heat exchanger. A multiscale Lagrangian approach was adopted to infer the trajectory and loading history experienced by platelet-like particles in the entire device and in a repetitive subunit of the fiber bundles. The loading history was incorporated into a damage accumulation model in order to estimate the platelet activation state (PAS) associated with repeated passes of the blood within the device. Our results highlighted the presence of blood stagnation areas in the inlet section that significantly increased the platelet activation levels in particles remaining trapped in this region. The order of magnitude of PAS in the device was the same as the one calculated for the components of the ECC tubing system, chosen as a term of comparison for their extensive diffusion. Interpolating the mean PAS values with respect to the number of passes, we obtained a straightforward prediction of the thrombogenic potential as a function of the duration of ECC.