Study on the optimal elastic modulus of flexible blades for right heart assist device supporting patients with single-ventricle physiologies

Background: Patients with single-ventricle physiologies continue to experience insufficient circulatory power after undergoing palliative surgeries. This paper proposed a right heart assist device equipped with flexible blades to provide circulatory assistance for these patients. The optimal elastic modulus of the flexible blades was investigated through numerical simulation.

Methods: A one-way fluid-structure interaction (FSI) simulation was employed to study the deformation of flexible blades during rotation and its impact on device performance. The process began with a computational fluid dynamics (CFD) simulation to calculate the blood pressure rise and the pressure on the blades’ surface. Subsequently, these pressure data were exported for finite element analysis (FEA) to compute the deformation of the blades. The fluid domain was then recreated based on the deformed blades’ shape. Iterative CFD and FEA simulations were performed until both the blood pressure rise and the blades’ shape stabilized. The blood pressure rise, hemolysis risk, and thrombosis risk corresponding to blades with different elastic moduli were exhaustively evaluated to determine the optimal elastic modulus.

Results: Except for the case at 8,000 rpm with a blade elastic modulus of 40 MPa, the pressure rise associated with flexible blades within the studied range (rotational speeds of 4,000 rpm and 8,000 rpm, elastic modulus between 10 MPa and 200 MPa) was lower than that of rigid blades. It was observed that the pressure rise corresponding to flexible blades increased as the elastic modulus increased. Additionally, no significant difference was found in the hemolysis risk and thrombus risk between flexible blades of various elastic moduli and rigid blades.

Conclusion: Except for one specific case, deformation of the flexible blades within the studied range led to a decrease in the impeller’s functionality. Notably, rotational speed had a more significant impact on hemolysis risk and thrombus risk compared to blade deformation. After a comprehensive analysis of blade compressibility, blood pressure rise, hemolysis risk, and thrombus risk, the optimal elastic modulus for the flexible blades was determined to be between 40 MPa and 50 MPa.

1 Introduction

Congenital heart defects are the most prevalent type of birth defect, occurring in approximately 9 out of every 1,000 live births (1). Among these, 9%–12% are classified as single ventricle physiologies (2). Patients with single-ventricle physiologies often have atrial and ventricular septal defects, leading to the mixing of arterial and venous blood within the heart chambers. A portion of the venous blood enters the systemic circulation without being oxygenated, causing organ hypoxia. Typically, patients undergo three-stage palliative surgeries: the Norwood, Glenn, and Fontan procedures (3). After the Fontan procedure, patient's superior and inferior vena cava (SVC and IVC, respectively) are connected to the pulmonary artery (PA), forming a total cavopulmonary connection (TCPC) structure. This setup forces the single ventricle to pump blood through both systemic and pulmonary circulations in the absence of a functional right ventricle. Over time, this can result in high central venous pressure and low pulmonary artery pressure (known as Fontan failure) (4). Approximately 40% of patients eventually suffer from single ventricle pumping failure (5).

Assistive devices play a critical role in augmenting circulatory power for these patients. In 2010, Throckmorton et al. (6) introduced an axial blood pump into the IVC to aid blood flow to the lungs. However, as the outlet of the blood pump is oriented towards the SVC, the blood flow pressurized by the pump collides with the blood flow within the SVC. This collision leads to an increase in blood pressure inside the SVC and hinders blood return. Then in 2013, Throckmorton et al. (7) implanted axial blood pumps in both SVC and IVC to provide dual-lumen assistance. In 2014, Wang et al. (8) used a paired umbrella double-lumen cannula (Maquet Cardiovascular, Wayne, NJ) and extracorporeal blood pump CentriMag (Thoratec, Pleasanton, CA) for total cavopulmonary assist. Animal studies revealed that this setup effectively improved hemodynamics in failing Fontan sheep. Subsequently, in 2019, Wang et al. (9) combined a double-lumen cannula with a two-valve extracardiac conduit and CentriMag to achieve effective support. In 2016, Gandolfo et al. (10) designed a T-shaped mechanically assisted TCPC model using the Jarvik Child 2,000 axial pump (Jarvik Heart, New York, NY, USA), which showed promising results in vivo. Although these techniques can provide total cavopulmonary assist, they have drawbacks such as the need for multiple devices, limitations on patient activity due to extracorporeal blood pumps, and the necessity for vascular reconfiguration.

To address the issues mentioned above, this study proposed a novel right heart assist device, as shown in Figure 1. The main components of this device are made of flexible materials. During minimally invasive surgery, the device is in a contracted state with a smaller diameter. Once implanted in the human body, the device expands and begins to function. It is anticipated that the deformation of flexible components during operation will impact the device's performance. Consequently, ensuring rigidity in its expanded state while maintaining flexibility during folding poses the principal challenge in developing this device. The material selection for each component is crucial. This research focuses on a key component—the blades, investigating their optimal elastic modulus. A one-way FSI simulation was employed to study blades’ deformation during rotation, and their impact on blood pressure increase, hemolysis risk, and thrombosis risk. The optimal elastic modulus is determined through a comprehensive analysis of these indicators.

Figure 1

Figure 1. Contraction state and post-implantation state of the flexible right heart assist device. (1. bearing seat; 2. shaft; 3. protective stent; 4. impeller; 5. constrictor ring; 6. motor; 7. cable; 8. inferior vena cava; 9. superior vena cava; 10. pulmonary artery).

2 Materials and methods

2.1 Flexible right heart assist device and computational models

Figure 1 shows the proposed right heart assist device, consisting of a protective stent, a bearing seat, a bearing, a shaft, an impeller, a constrictor ring, a motor, and cables. The protective stent is made of Nitinol. One end of the stent is fixed to the bearing seat, while the other end connects to the constrictor ring and can move axially. The impeller is fabricated from flexible materials such as polyurethane through an integrated molding process. It consists of two parts: the hub and the blades, with the base of the blades affixed to the surface of the hub. The hub's upper end is mounted on the shaft, and its lower end can move axially. The blades can adopt two configurations: curled and unfolded. When the assist device contracts, the lower end of the hub moves downward, and the blades curl. As the device expands, both the hub and the blades revert to their initial shape. The motor is powered by cables, and drives the rotation of the shaft and impeller. During surgery, the entire device is compressed within a delivery catheter and implanted into the TCPC via the femoral vein. Once released, the protective stent expands to tightly fit against the inner walls of the vena cava, stabilizing the device. The rotating impeller aids in directing blood within the SVC and IVC into the PA. Additionally, the hub is uniquely designed in a spindle shape to prevent blood from colliding. The design goal of this assist device is to provide a pressure rise of 2–25 mmHg for the total systemic blood flow (3l /min 11) in Fontan patients when operating at 4,000–8,000 rpm (12).

This paper focuses on flexible blades, therefore adopts simplified computational models. The TCPC model is an ideal circular tube model derived from patient's data (13). According to the principle of equal area, the diameters of the SVC and IVC are set at 13 mm, while those of the LPA and RPA are set at 9 mm. In CFD simulation, only the impeller is retained. The base of the blades is fixed to the hub's surface to ensure the model's consistency with reality. The minimum and maximum diameters of the impeller are 0.4 mm and 10.9 mm respectively, the number of blades is 12, and the thickness of the blades is 0.2 mm. In FEA simulation, only the blades are retained. The deformation of the hub is ignored. The impeller and blade models presented in this paper are the optimized versions after structural optimization. Figure 2 illustrates the final computational models used in this study.

Figure 2

Figure 2. CFD and FEA simulation model. (A) Ideal TCPC structure; (B) blades of right heart assist device.

2.2 One-way FSI methodology, boundary conditions and solving settings

In this study, one-way FSI simulation was used to analyze the flexible blades’ deformation and its effect on the device performance, as shown in Figure 3. First, ANSYS Fluent 2020R2 (Canonsburg, USA) was used to model the fluid dynamics in TCPC and obtain the blade surface pressure. Subsequently, these pressure data were utilized in ANSYS Structural 2020R2 (Canonsburg, USA) to compute blades’ deformation. The deformation of the blades was quantitatively described by the displacement of nodes on the blades. After the deformed shape of the blades was exported, the fluid domain was remeshed. Subsequently, the CFD simulation was performed again. The cycle was ended when the difference in blood pressure rise between two consecutive CFD simulations was less than 5%. Finally, the last CFD result was used to calculate the hemolysis risk and thrombus risk.

Figure 3

Figure 3. One-way FSI simulation process.

In CFD simulation, flow boundary conditions were defined at the SVC and IVC inlets, with flow rates of 1.2l /min and 1.8l /min respectively (11). The LPA and RPA outlets were set as pressure boundary conditions, with a pressure value of 1,750 Pa, corresponding to the average pulmonary artery pressure during a patient's cardiac cycle (13). Blood was modeled as a Newtonian fluid with a density of 1,060 kg/m³ and a dynamic viscosity of 3.5 cP. This article analyzes the operation of the assist device at its lowest and highest rotational speeds. The flow domain is divided into rotating and stationary regions. The rotational speed of the rotating region is set at 4,000 rpm and 8,000 rpm. Furthermore, rotational boundary conditions are applied to the surfaces of the hub and blades, with the rotational speed settings matching those of the rotating region. The remaining surfaces are treated as stationary boundaries. Additionally, the SST k-ω model was used to simulate turbulent conditions, as this model can accurately resolve the flow fields both at the wall and in the bulk flow regime (12, 14, 15). The fluid domain was meshed as a polyhedral grid using Fluent Meshing 2020R2 (Canonsburg, USA). After grid independence study, the fluid mesh consisting of approximately 1.35 million elements was used for CFD simulation. Inflation layers were also used in the fluid meshes to achieve the SST k-ω turbulence model requirement (y + less than 2). A steady-state solver was employed, with the convergence criteria was set to a maximum residual of below the 1 × 10⁻⁴ threshold.

In FEA simulation, the blade material was set based on the polyurethane material properties. There is a wide variety of linear elastic polyurethanes, with elastic moduli ranging from 6.2 MPa to 27,600 MPa (16). A pre-simulation (a single round of one-way FSI) was conducted to narrow down the range of elastic moduli for the study, as shown in Figure 4. At 4,000 rpm and 8,000 rpm, the maximum deformation of the blades showed a consistent pattern. When the elastic modulus is within 50 MPa, the maximum deformation decreases rapidly as the elastic modulus increases. Between 50 MPa and 200 MPa, the decrease in maximum deformation is more gradual. Above 200 MPa, both the absolute value and the variation in maximum deformation are extremely small. Based on these results, the elastic moduli selected for the study were 10 MPa, 20 MPa, 30 MPa, 40 MPa, 50 MPa, 100 MPa, and 200 MPa. In the actual configuration, the base of the blades is fixed to the hub's surface, leading to the application of a fixed constraint at the base of the blades. The structural mesh of the blade, generated using HyperMesh (Altair, USA), is a hexahedral mesh and consists of two layers along the thickness direction. After the grid independence study, the structural mesh consisting of 36,096 elements was used for FEA simulation. A transient solver was employed, with a total computational duration of 1.1 s and a time step of 0.001 s. The convergence criteria was set as changes of the maximum deformation, the minimum deformation, and the average deformation are all less than 1%.

Figure 4

Figure 4. Maximum deformation of blades with different elastic moduli at rotational speeds of 4,000 rpm and 8,000 rpm. (These results were obtained from one round of one-way FSI simulation).

2.3 Hemolysis risk assessment model

During the high-speed rotation of the impeller, red blood cells are destroyed due to shear forces, resulting in hemolysis. The hemolysis risk can be quantitatively evaluated using the Normalized index of hemolysis (NIH). Numerical simulations and hemolysis experiments employ different calculation methods to derive this index. In numerical simulation, NIH is commonly derived through a power law model, as detailed in Equations 1–3 (17–19).

$\sigma ={\left({\displaystyle \frac{1}{6}\sum {({\sigma }_{\mathrm{ii}}-{\sigma }_{\mathrm{jj}})}^2+\sum {\sigma }_{\mathrm{ij}}^2}\right)}^{\frac{1}{2}}(1)$

$D={\sum }_{\mathrm{inlet}}^{\mathrm{outlet}}C\cdot {\sigma }^{\alpha }\cdot \mathrm{\Delta }{T}^{\beta }(2)$

$NIH=0.00015\mathbf{\ast }D(3)$

where σ is the scalar stress, $\mathrm{\Delta }T$ is the exposure time, and $C,\alpha ,\beta$ are constants obtained from hemolysis experiments. Mitamura et al. (20) derived a set of parameters based on the experimental data from axial blood pumps: $C=3.62\times {10}^{-5},\alpha =2.416,\beta =0.785$. Heuser et al. (21) used experimental data from a Couette viscometer to fit another set of parameters: $C=1.8\times {10}^{-6},\alpha =1.991,\beta =0.765$. The working principle of the impeller in this paper differs considerably from that of the axial blood pump and Couette viscometer, and is closer to that of the centrifugal blood pump. Therefore, this paper uses experimental data of centrifugal blood pumps, provided by the U.S. Food and Drug Administration (FDA), to determine the parameters of the hemolysis risk assessment model (22).

In addition, Blackshear et al. (23) conducted a hemolysis experiment and proposed the relationship between shear stress experienced by red blood cells and their exposure time:

$({\tau }_{\mathrm{rbc}}^2)(t)=C_1(4)$

In the above equation $C_1$ is a constant of proportionality. The results of Blackshear et al. indicate that the hemolysis risk is proportional to the exposure time and the square of the shear stress. Combined with Equations 1–4, the final hemolysis risk assessment model used in this paper is as follows:

$\sigma ={\left({\displaystyle \frac{1}{6}\sum {({\sigma }_{\mathrm{ii}}-{\sigma }_{\mathrm{jj}})}^2+\sum {\sigma }_{\mathrm{ij}}^2}\right)}^{\frac{1}{2}}(5)$

$NIH={\sum }_{\mathrm{inlet}}^{\mathrm{outlet}}A\cdot {\sigma }^{2B}\cdot \mathrm{\Delta }{T}^B(6)$

In Equation 6, A represents the overall scaling factor, set at 2.7 × 10⁻¹⁰. This value was derived by multiplying the scaling factor C, provided by Heuser (21), with the scaling factor of 0.00015 from Equation 3. The parameter B, set at 0.55, was determined from experimental data, as detailed in the Supplementary Document.

2.4 Thrombus risk assessment model

Thrombus formation is another common complication associated with rotary assist devices (24). It involves the generation, transport, and cascade reactions of numerous components. The adhesion and activation of platelets are critical steps in thrombus formation. Cell adhesion only occurs in areas where the wall shear strain rate $(\dot{\gamma })$ is lower than 100 /s (25). Therefore, the shear strain rate can be used to identify areas prone to thrombus growth (26). Its calculation method is shown in Equation 7. An accelerated thrombus formation model is utilized to quantitatively assess the thrombus risk (27). This model includes three convection-diffusion equations (Equations 8–10), describing the generation and transport of non-activated platelets (NP), activated platelets (AP), and adenosine diphosphate (ADP). The concentration of AP is ultimately used to represent the thrombus risk.

$\dot{\gamma }={\left[2\left({\left({\displaystyle \frac{\mathrm{\partial }u}{\mathrm{\partial }x}}\right)}^2+{\left({\displaystyle \frac{\mathrm{\partial }v}{\mathrm{\partial }y}}\right)}^2+{\left({\displaystyle \frac{\mathrm{\partial }w}{\mathrm{\partial }z}}\right)}^2\right)+{\left({\displaystyle \frac{\mathrm{\partial }u}{\mathrm{\partial }y}+{\displaystyle \frac{\mathrm{\partial }v}{\mathrm{\partial }x}}}\right)}^2+{\left({\displaystyle \frac{\mathrm{\partial }u}{\mathrm{\partial }z}+{\displaystyle \frac{\mathrm{\partial }w}{\mathrm{\partial }x}}}\right)}^2+{\left({\displaystyle \frac{\mathrm{\partial }v}{\mathrm{\partial }z}+{\displaystyle \frac{\mathrm{\partial }w}{\mathrm{\partial }y}}}\right)}^2\right]}^{\frac{1}{2}}(7)$

${\displaystyle \frac{\mathrm{\partial }({\varphi }_a)}{\mathrm{\partial }t}+(u\cdot \mathrm{\nabla }){\varphi }_a=D_a{\mathrm{\nabla }}^2{\varphi }_a+\{[A_C(ADP)]{\varphi }_n+[A_M({\varphi }_f,\tau )]({\varphi }_a+{\varphi }_n)\}}(8)$

${\displaystyle \frac{\mathrm{\partial }({\varphi }_n)}{\mathrm{\partial }t}+(u\cdot \mathrm{\nabla }){\varphi }_n=D_n{\mathrm{\nabla }}^2{\varphi }_n-\{[A_C(ADP)]{\varphi }_n+[A_M({\varphi }_f,\tau )]({\varphi }_a+{\varphi }_n)\}}(9)$

${\displaystyle \frac{\mathrm{\partial }(ADP)}{\mathrm{\partial }t}+(u\cdot \mathrm{\nabla })ADP=D_{ADP}{\mathrm{\nabla }}^{2}ADP+R_{ADP}\{[A_C(ADP)]{\varphi }_n+[A_M({\varphi }_f,\tau )]({\varphi }_a+{\varphi }_n)\}}(10)$

class="mb15">In the above equations, D_x is the respective diffusivity coefficient, A_M and A_C denotes mechanical activation and chemical activation, and R_ADP is the amount of ADP contained in a platelet. Platelet activation (change from NP to AP) occurs either with mechanical cues (Equation 11) or chemical cues (Equation 12).

$A_M({\varphi }_f,\tau )=(1-{\varphi }_f)C_1^{\frac{1}{{\beta }_1}}{\beta }_1{\varphi }_f^{\frac{{\beta }_1-1}{{\beta }_1}}{\tau }^{\frac{{\alpha }_1}{{\beta }_1}}(11)$

$A_C(\mathrm{ADP})=\{\begin{array}{c}{\displaystyle \frac{\mathrm{ADP}}{\mathrm{AD}{\mathrm{P}}_t\cdot t_{\mathrm{ADP}}},\mathrm{ADP}\ge \mathrm{AD}{\mathrm{P}}_t}\\ 0,\mathrm{ADP}<\mathrm{AD}{\mathrm{P}}_t\end{array}(12)$

where ${\varphi }_f$ is the ratio of the number of activated platelets to total platelets count, ${\alpha }_1$, ${\beta }_1$ and $C_1$ are the power law model parameter, $\mathrm{AD}{\mathrm{P}}_t$ is the threshold for chemical platelet activation, $t_{\mathrm{ADP}}$ is the characteristic time of platelet activation, and τ is the shear stress calculated from Equations 13–15.

$\tau ={\displaystyle \frac{1}{\sqrt{3}}\sqrt{{\sigma }_{\mathrm{xx}}^2+{\sigma }_{\mathrm{yy}}^2+{\sigma }_{\mathrm{zz}}^2-{\sigma }_{\mathrm{xx}}{\sigma }_{\mathrm{yy}}-{\sigma }_{\mathrm{yy}}{\sigma }_{\mathrm{zz}}-{\sigma }_{\mathrm{xx}}{\sigma }_{\mathrm{zz}}+3({\sigma }_{\mathrm{xy}}^2+{\sigma }_{\mathrm{yz}}^2+{\sigma }_{\mathrm{xz}}^2)}}(13)$

${\sigma }_{\mathrm{ij}}=\mu \left({\displaystyle \frac{\mathrm{\partial }{u}_i}{\mathrm{\partial }{x}_j}+{\displaystyle \frac{\mathrm{\partial }{u}_j}{\mathrm{\partial }{x}_i}}}\right)-\rho \overline{u_i^{\mathrm{\prime }}{u}_j^{\mathrm{\prime }}}(14)$

$-\rho \overline{u_i^{\mathrm{\prime }}{u}_j^{\mathrm{\prime }}}={\mu }_t\left({\displaystyle \frac{\mathrm{\partial }{u}_i}{\mathrm{\partial }{x}_j}+{\displaystyle \frac{\mathrm{\partial }{u}_j}{\mathrm{\partial }{x}_i}}}\right)-{\displaystyle \frac{2}{3}\rho k{\delta }_{\mathrm{ij}}}(15)$

In the above equations μ is the dynamic viscosity, ${\mu }_t$ is the eddy viscosity, ρ is the fluid density, k is the turbulent kinetic energy, ${\delta }_{\mathrm{ij}}$ is the Kronecker delta, $u_i$ is the velocity component and $u_i$ is the spatial component.

3 Results

3.1 Variation of blades’ deformation and pressure rise during FSI simulation iterations

Figure 5; Table 1 illustrate the changes in blades’ deformation and blood pressure rise during FSI simulation iterations. Figures 5A,C,D sequentially show the blades’ deformation during three rounds of FSI simulation. It is evident that as the iterations progress, the solid body gradually aligns with the dashed frame. Additionally, the maximum deformation of the blades gradually decreases, indicating that the blade shape is stabilizing. According to Table 1, the pressure rises in both the SVC and IVC also tend to stabilize during the iteration process. The pressure rises of the last two iterations differing by less than 5%. These results validate that the one-way FSI simulation used in this paper can achieve a stable blade shape and the corresponding pressure rise.

Figure 5

Figure 5. Blades' deformation during the one-way FSI simulation iterations. (A) Blades' deformation obtained from the first round of FSI; (B) Blades' deformation under stress; (C) Blade deformation obtained from the second round of FSI; (D) Blade deformation obtained from the third round of FSI. (Dotted frame represents the blade shape before deformation, solid body represents the blade shape after deformation).

Table 1

Table 1. Maximum blade deformation and blade pressure rise during three rounds of one-way FSI simulation.

Moreover, stress exists within the blade in the second and subsequent rounds of FSI simulation. Figures 5A,B were compared with each other to verify the accuracy of applied stress. Figure 5A shows the blade bending under surface pressure. Figure 5B depicts the blade returning to vertical shape after the removal of surface pressure, while stress is retained. The maximum deformations in Figures 5A,B are 1.55 mm and 1.505 mm, respectively. The difference is only 3%, indicating that the stress is sufficiently precise.

3.2 Blood pressure rise corresponding to blades with different elastic moduli at two rotational speeds

Figure 6 displays the blood pressure rise corresponding to flexible blades and rigid blades. DP-Total represents the total pressure rise, obtained by the weighted sum of DP-SVC and DP-IVC, based on the flow rates in SVC and IVC (DP-Total = DP-SVC × 0.4 + DP-IVC × 0.6). At a speed of 8,000 rpm, blades with an elastic modulus of 10 MPa deform excessively and lose their functionality entirely. Consequently, data related to this condition are not included in the analysis.

Figure 6

Figure 6. Blood pressure rise corresponding to blades with different elastic moduli at two rotational speeds. (A) At 4,000 rpm; (B) at 8,000 rpm.

From Figure 6A, it can be observed that at 4,000 rpm, the pressure rises corresponding to flexible blades show an increasing trend. As the elastic modulus increases, DP-SVC, DP-IVC and DP-Total all gradually increase. The pressure rises corresponding to flexible blades with elastic moduli of 50 MPa, 100 MPa, and 200 MPa are relatively close to those of rigid blades. Figure 6B reveals that at 8,000 rpm, the pressure rises corresponding to flexible blades also trends upward. However, when the elastic modulus of the blades increases from 40 MPa to 50 MPa, the pressure rise decreases instead. The pressure rises corresponding to flexible blades with elastic moduli of 40 MPa, 50 MPa, 100 MPa, and 200 MPa are relatively close to those of rigid blades. Notably, the pressure rise corresponding to the 40 MPa elastic modulus flexible blade is higher than that of the rigid blades. Comprehensive analysis indicates that in most cases, blade deformation leads to a decrease in functionality. However, at the highest rotational speed, there are instances where blade deformation results in an increase in functionality.

3.3 Hemolysis risk corresponding to blades with different elastic moduli at two rotational speeds

Figure 7 shows the hemolysis risk associated with flexible blades and rigid blades. After completing the streamline number independence test, 1,100 streamlines were released at both the SVC and IVC inlets. The NIH values for each streamline were calculated and average values were used to represent the hemolysis risk in the blood flow of the SVC and IVC (NIH-SVC and NIH-IVC). The overall hemolysis risk (NIH-Total) was derived by the weighted sum of NIH-SVC and NIH-IVC (NIH-Total = NIH-SVC × 0.4 + NIH-IVC × 0.6). The results indicate that at both rotational speeds, the hemolysis risks associated with flexible and rigid blades are relatively similar. At 4,000 rpm, the NIH-Total for flexible blades ranges between 7.26 × 10⁻⁴ and 7.67 × 10⁻⁴ g/100 L, with a maximum difference of 5.6%. The NIH-Total for rigid blades is 7.03 × 10⁻⁴ g/100 L. Therefore, the difference in NIH-Total between flexible and rigid blades ranges from 3.3% to 9.1%. At 8,000 rpm, the NIH-Total for flexible blades ranges between 15.01 × 10⁻⁴ and 16.65 × 10⁻⁴ g/100 L, with a maximum difference of 10.9%. The NIH-Total for rigid blades is 15.48 × 10⁻⁴ g/100 L.As a result, the difference in NIH-Total between flexible and rigid blades ranges from 0.2% to 7.6%.

Figure 7

Figure 7. Hemolysis risk corresponding to blades with different elastic moduli at two rotational speeds. (A) At 4,000 rpm; (B) at 8,000 rpm.

Under the same rotational speed, flexible blades with different elastic moduli have varying shapes, yet the corresponding hemolysis risks are similar. This suggests that the impact of rotational speed on hemolysis risk is greater than that of blade deformation. Moreover, literature data indicate that the NIH values for HeartMateII (Thoratec Inc., Pleasanto, California, USA) and CH-VAD (CH Biomedical, Inc., Suzhou, China) at medium speeds are 5.83 × 10⁻³ g/100 L and 1.35 × 10⁻³ g/100 L, respectively (28). The highest NIH-Total in this study is 1.665 × 10⁻³ g/100 L, falling between the NIH values of HeartMateII and CH-VAD. This result demonstrates that the assist device proposed in this paper is reliable in terms of hemolysis.

3.4 Thrombus risk corresponding to blades with different elastic moduli at two rotational speeds

Figure 8 presents the results of thrombus formation simulation, corresponding to a rotational speed of 4,000 rpm and a blade elastic modulus of 10 MPa. As mentioned earlier, the shear strain rate (SSR) is utilized to determine areas prone to thrombus growth. The concentration of AP is used for quantitatively assessing the thrombus risk. To enhance readability, the AP concentration is scaled and represented as the SAP index. This index has a minimum value of 0, corresponding to the initial AP concentration (25 × 10¹²/m³). Each unit increase indicates an increase of one millionth of the initial AP concentration. Figure 8A shows the distribution of wall SSR. It is observed that regions with SSR less than 100 /s are only present on the walls of the SVC, LPA and RPA. There are no low SSR areas on the blades’ surfaces. Figure 8B, based on Figure 8A, displays grids where the SAP exceeds 23 in a dark red color. Detailed examination through close-up and cross-sectional views reveals that all four vessels exhibit wall grids with higher SAP. However, IVC lacks regions with low SSR, and the high SAP at the LPA and RPA outlets is due to the transport of AP. Only the wall of the SVC displays areas with both low SSR and high SAP, indicating that thrombus is more likely to form on the SVC wall. Figure 8C shows streamlines released from the IVC. Due to the larger flow volume in the IVC, a portion of the blood enters the SVC, leading to a collision within the SVC. This collision area has a low SSR, which facilitates the deposition of AP and thus the formation of thrombus. At both rotational speeds, the thrombus formation scenarios for different blades are similar to those depicted in Figure 8B. Thrombus are more likely to growth on the wall of SVC, hence separate illustrations are not provided.

Figure 8

Figure 8. Simulation results of thrombus formation. (Corresponding to rotational speed of 4,000 rpm, blade elastic modulus of 10 MPa. (A) Distribution of wall shear strain rate; (B) thrombus distribution; (C) IVC streamline diagram).

A statistical analysis of the mean SAP values in the SVC for different blades provides the following insights: At 4,000 rpm, the mean SAP value for flexible blades with various elastic moduli ranges between 2.60 and 2.68, with a maximum difference of 3%. The mean SAP value for rigid blades is 2.60. As a result, the maximum difference in the thrombus risk between flexible and rigid blades is 3%. At 8,000 rpm, the mean SAP value for flexible blades with different elastic moduli ranges from 14.92 to 15.88, with a maximum difference of 6%. The mean SAP value for rigid blades is 15.24. Therefore, the maximum difference in thrombus risk between flexible and rigid blades is 4.2%. These results indicate that, in terms of thrombus formation, the impact of rotational speed is also much greater than that of blade deformation.

4 Discussion

The Fontan surgery successfully diverts venous blood into the lungs by constructing the TCPC structure. The issue of mixing between arterial and venous blood has been resolved. However, long-term reliance on a single ventricle to power both systemic and pulmonary circulations can lead to various complications, such as elevated central venous pressure, increased pre-load on the single ventricle, and insufficient pulmonary perfusion. Patients are at high risk of developing single ventricle heart failure. Numerous scholars have proposed a variety of solutions to provide additional power for the Fontan circulation. Based on the location of the assistive device, these solutions can be categorized into extracorporeal and intracorporeal assistance. Trusty et al. (29) used pulsatile blood pumps, conduits, and baffles to explore efficient extracorporeal assistance methods in an in vitro loop. Their findings revealed that without partitions, blood pumped out of the pump's outlet tube could be reabsorbed by the inlet tube, causing recirculation. Effective power assistance was only possible by completely separating the inlet and outlet tubes. Wang et al. (8, 9) used dual-lumen catheters and extracorporeal blood pumps for Fontan support. The successful assistance achieved was also due to the use of valves or umbrella-like membranes that separated the inlet and outlet of the dual-lumen catheters. Therefore, providing extracorporeal assistance to patients requires the addition of specific obstructions to the existing TCPC structure, increasing the number of devices and the complexity of surgery. Intracorporeal assistance methods are more convenient as they do not alter the existing TCPC structure. Among these, collapsible assist devices have a smaller size during implantation, leading to reduced invasiveness. After implantation, these devices expand, increasing in diameter. Compared to rigid devices of similar implantation size, collapsible devices operate at lower speeds to achieve the same pressure rise. Therefore, the temperature increase, mechanical wear, and blood damage during the operation of collapsible pumps are all lower (30).

In the field of ventricular assist devices, the folding mechanisms of flexible blades mainly include mechanical folding, material folding, and hybrid folding. Mechanical folding typically employs hinged structures (31–33), requiring locking components to ensure the blades’ operational shape. However, the use of complex blood-immersed joint significantly increases thrombus risks. Material folding relies on the material’s properties for folding, offering simpler structures (34–36). This approach requires the blades to exhibit two distinct states: one for elastic deformation into the stored position, and another to resist hydrodynamic forces during operation. The hybrid folding mechanism combines the advantages of the previous two methods. It involves designing unique structures for flexible blades to achieve structural anisotropy (37, 38). This mechanism allows the device to be both easily collapsible and capable of maintaining its operational shape. In the latter two folding mechanisms, the blade deformation and its impact on device performance are particularly crucial for product development. FSI simulation can rapidly explore relevant patterns, offering guidance for material selection and structural design.

For high-speed rotating devices, conducting two-way FSI simulation can yield results that are closer to reality. However, this method requires a significant amount of resources. At the early stage of research, employing one-way FSI simulation can also reveal certain patterns while substantially reducing resource consumption. The results demonstrate that the one-way FSI simulation described in this paper can accurately model the stable shape of flexible blades. This allows for the investigation of the impact of blade deformation on device performance, aiding in the determination of the optimal blade elastic modulus. It is important to note that one-way FSI simulation neglects the impact of blade deformation velocities and accelerations on the fluid, failing to capture transient flow fluctuations. This leads to some deviation in the calculations of blood pressure rise, hemolysis risk, and thrombus risk. However, these neglected fluctuations are primarily near the impeller, and their impact area is limited. As presented in the paper, blood pressure rise is largely determined by the rotational speed and blade morphology. Furthermore, the hemolysis risk and thrombus risk are mainly influenced by the rotational speed. Compared to these factors, the impact of flow fluctuations is relatively minor. Although one-way FSI simulation overlooks the influence of structural dynamics on the fluid, it still accurately depicts the impact of the main factors on performance indicators. This method can be further applied to the entire impeller to explore the effect of hub deformation. For flexible left ventricular assist device, such as impella ECP (Abiomed, Danvers, MA), one-way FSI simulation is equally applicable.

This study has the following limitations: (1) It did not analyze the entire device, thus ignoring the impact of components other than the impeller on blood pressure rise, hemolysis risk, and thrombus risk; (2) Fluctuations in rotational speed caused by the operation of the brushless motor were not considered; (3) The one-way FSI simulation can't capture blades’ status at each moment; (4) No experiments were conducted to verify the accuracy of the simulation results; (5) This study employs steady-state boundary conditions and an ideal TCPC model for research, neglecting the effects of blood pulsatility and real vascular geometry. The future research plans include: (1) Modeling and simulating the entire device; (2) Constructing an in-vitro experimental system, using a high-speed camera to capture the shape of flexible blades during rotation, and comparing it with the simulation results; (3) Conducting hydraulic experiments, particle image velocimetry (PIV) experiments, hemolysis experiments, and thrombosis experiments to verify the simulation results of pressure rise, velocity field, hemolysis risk, and thrombus risk; (4) Collecting data on rotational speeds during device operation and applying these varying speeds in FLUENT using a user-defined function (UDF); (5) Exploring two-way FSI simulation methods and using transient simulation to study the impact of blade flutter on the flow field; (6) Performing patient-specific studies using measured boundary conditions and real TCPC models.

5 Conclusion

This article proposed a right heart assist device equipped with flexible blades. A one-way FSI simulation was used to study the deformation of the flexible blades and its impact on the device's performance. The simulation results indicate that within the studied range (rotational speeds of 4,000 rpm and 8,000 rpm, elastic modulus between 10 MPa and 200 MPa), deformation of the flexible blades leads to a decrease in functionality. However, within a certain range of elastic modulus, the pressure rise associated with flexible blades is close to that of rigid blades. Additionally, compared to blade deformation, rotational speed has a more significant impact on the hemolysis risk and thrombus risk. After a comprehensive analysis of blade compressibility, blood pressure rise, hemolysis risk, and thrombus risk, the optimal elastic modulus for the flexible blades is determined to be between 40 MPa and 50 MPa.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Ethics Statement

Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent from the patients/participants was not required to participate in this study in accordance with the national legislation and the institutional requirements.

Author contributions

TC: Conceptualization, Formal Analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – original draft, Writing – review & editing. XC: Conceptualization, Formal Analysis, Investigation, Methodology, Software, Validation, Writing – review & editing. XL: Methodology, Software, Writing – review & editing. HZ: Data curation, Resources, Writing – review & editing. SW: Conceptualization, Funding acquisition, Project administration, Supervision, Writing – review & editing, Resources.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article.

This research was funded by National Natural Science Foundation of China, grant number 32071310.

Conflict of interest

XL is employed by Shanghai Key Laboratory of Interventional Medical Devices and Equipment, Shanghai MicroPort Medical Group Co., Ltd, Shanghai, China.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2024.1377765/full#supplementary-material

References

1. Botto LD, Correa A, Erickson JD. Racial and temporal variations in the prevalence of heart defects. Pediatrics. (2001) 107(3):e32–32. doi: 10.1542/peds.107.3.e32

PubMed Abstract | Crossref Full Text | Google Scholar

2. Pike NA, Evangelista LS, Doering LV, Griffin DK, Lewis AB, Child JS. Clinical profile of the adolescent/adult fontan survivor. Congenit Heart Dis. (2011) 6(1):9–17. doi: 10.1111/j.1747-0803.2010.00475.x

PubMed Abstract | Crossref Full Text | Google Scholar

3. Trusty PM, Slesnick TC, Wei ZA, Rossignac J, Kanter KR, Fogel MA, et al. Fontan surgical planning: previous accomplishments, current challenges, and future directions. J Cardiovasc Transl. (2018) 11:133–44. doi: 10.1007/s12265-018-9786-0

Crossref Full Text | Google Scholar

4. Khairy P, Fernandes SM, Mayer JE, Triedman JK, Walsh EP, Lock JE, et al. Long-term survival, modes of death, and predictors of mortality in patients with fontan surgery. Circulation. (2008) 117:85–92. doi: 10.1161/CIRCULATIONAHA.107.738559

PubMed Abstract | Crossref Full Text | Google Scholar

5. Rodefeld MD, Frankel SH, Giridharan GA. Cavopulmonary assist: (em)powering the univentricular fontan circulation. Semin Throac Cardiov. (2011) 14:45–54. doi: 10.1053/j.pcsu.2011.01.015

Crossref Full Text | Google Scholar

6. Bhavsar SS, Moskowitz WB, Throckmorton AL. Interaction of an idealized cavopulmonary circulation with mechanical circulatory assist using an intravascular rotary blood pump. Artif Organs. (2010) 34:816–27. doi: 10.1111/j.1525-1594.2010.01080.x

PubMed Abstract | Crossref Full Text | Google Scholar

7. Throckmorton AL, Lopez Isaza S, Moskowitz W. Dual-pump support in the inferior and superior vena cavae of a patient-specific fontan physiology. Artif Organs. (2013) 37:513–22. doi: 10.1111/aor.12039

PubMed Abstract | Crossref Full Text | Google Scholar

8. Wang D, Gao G, Plunkett M, Zhao G, Topaz S, Ballard-Croft C, et al. A paired membrane umbrella double-lumen cannula ensures consistent cavopulmonary assistance in a fontan sheep model. J Thorac Cardiovasc Surg. (2014) 148:1041–7. doi: 10.1016/j.jtcvs.2014.04.051

PubMed Abstract | Crossref Full Text | Google Scholar

9. Zhou C, Wang D, Condemi F, Zhao G, Topaz S, Ballard-Croft C, et al. Avalonelite double lumen cannula for total cavopulmonary assist in failing fontan sheep model with valved extracardiac conduit. Asaio J. (2019) 65:361–6. doi: 10.1097/MAT.0000000000000817

PubMed Abstract | Crossref Full Text | Google Scholar

10. Gandolfo F, Brancaccio G, Donatiello S, Filippelli S, Perri G, Iannace E, et al. Mechanically assisted total cavopulmonary connection with an axial flow pump: computational and in vivo study. Artif Organs. (2016) 40:43–9. doi: 10.1111/aor.12641

PubMed Abstract | Crossref Full Text | Google Scholar

11. Zhao ZX, Chen T, Liu XD, Wang SZ, Lu HY. Numerical investigation of an idealized total cavopulmonary connection physiology assisted by the axial blood pump with and without diffuser. Comput Model Eng Sci. (2020) 125(3):1173–84. doi: 10.32604/cmes.2020.013702

Crossref Full Text | Google Scholar

12. Hirschhorn M, Bisirri E, Stevens R, Throckmorton AL. Fluid-structure interaction analysis of a collapsible axial flow blood pump impeller and protective cage for fontan patients. Artif Organs. (2020) 44:E337–47. doi: 10.1111/aor.13685

PubMed Abstract | Crossref Full Text | Google Scholar

13. Chen T, Liu X, Si B, Feng Y, Zhang H, Jia B, et al. Comparison between single-phase flow simulation and multiphase flow simulation of patient-specific total cavopulmonary connection structures assisted by a rotationally symmetric blood pump. Symmetry (Basel). (2021) 13:912. doi: 10.3390/sym13050912

Crossref Full Text | Google Scholar

14. Jagani J, Mack E, Gong J, Untaroiu A. Effect of stent design parameters on hemodynamics and blood damage in a percutaneous cavopulmonary assist device. ASME. (2018) 51562:V002T11A007. doi: 10.1115/FEDSM2018-83210

Crossref Full Text | Google Scholar

15. Pei J, Meng F, Li Y, Yuan S, Chen J. Fluid–structure coupling analysis of deformation and stress in impeller of an axial-flow pump with two-way passage. Adv Mech Eng. (2016) 8:168781401664626. doi: 10.1177/1687814016646266

Crossref Full Text | Google Scholar

16. TPU. Unspecified (2024). Available online at: https://www.plasdata.com/ (accessed January 25, 2024).

17. Bludszuweit C. Three-dimensional numerical prediction of stress loading of blood particles in a centrifugal pump. Artif Organs. (1995) 19:590–6. doi: 10.1111/j.1525-1594.1995.tb02386.x

PubMed Abstract | Crossref Full Text | Google Scholar

18. Giersiepen M, Wurzinger LJ, Opitz R, Reul H. Estimation of shear stress-related blood damage in heart valve prostheses–in vitro comparison of 25 aortic valves. Int J Artif Organs. (1990) 13:300–6. doi: 10.1177/039139889001300507

PubMed Abstract | Crossref Full Text | Google Scholar

19. Throckmorton AL, Kishore RA. Design of a protective cage for an intravascular axial flow blood pump to mechanically assist the failing fontan. Artif Organs. (2009) 33:611–21. doi: 10.1111/j.1525-1594.2009.00779.x

PubMed Abstract | Crossref Full Text | Google Scholar

20. Mitamura Y, Nakamura H, Sekine K, Kim DW, Youzu R, Kawada S, et al. Prediction of hemolysis in rotary blood pumps with computational fluid dynamics analysis. J CongestHeart Fail Circ Supp. (2001) 1(4):331–6. doi: 10.1201/B14731-35

Crossref Full Text | Google Scholar

21. Heuser G, Opitz R. A couette viscometer for short time shearing in blood. Biorheology. (1980) 17(1-2):17–24. doi: 10.3233/BIR-1980-171-205

PubMed Abstract | Crossref Full Text | Google Scholar

22. Malinauskas RA, Hariharan P, Day SW, Herbertson LH, Buesen M, Steinseifer U, et al. FDA benchmark medical device flow models for CFD validation. Asaio J. (2017) 63:150–60. doi: 10.1097/MAT.0000000000000499

PubMed Abstract | Crossref Full Text | Google Scholar

23. Blackshear PL, Dorman FD, Steinbach JH. Some mechanical effects that influence hemolysis. ASAIO J. (1965) 11(1):112–7. doi: 10.1097/00002480-196504000-00022

Crossref Full Text | Google Scholar

24. Copeland JG, Copeland H, Gustafson M, Mineburg N, Covington D, Smith RG, et al. Experience with more than 100 total artificial heart implants. J Thorac Cardiovasc Surg. (2012) 143:727–34. doi: 10.1016/j.jtcvs.2011.12.002

PubMed Abstract | Crossref Full Text | Google Scholar

25. Goel MS, Diamond SL. Adhesion of normal erythrocytes at depressed venous shear rates to activated neutrophils, activated platelets, and fibrin polymerized from plasma. Blood. (2002) 100:3797–803. doi: 10.1182/blood-2002-03-0712

PubMed Abstract | Crossref Full Text | Google Scholar

26. Menichini C, Xu XY. Mathematical modeling of thrombus formation in idealized models of aortic dissection: initial findings and potential applications. J Math Biol. (2016) 73:1205–26. doi: 10.1007/s00285-016-0986-4

PubMed Abstract | Crossref Full Text | Google Scholar

27. Blum C, Groß-Hardt S, Steinseifer U, Neidlin M. An accelerated thrombosis model for computational fluid dynamics simulations in rotary blood pumps. Cardiovasc Eng Techn. (2022) 13:638–49. doi: 10.1007/s13239-021-00606-y

Crossref Full Text | Google Scholar

28. Berk ZBK, Zhang J, Chen Z, Tran D, Griffith BP, Wu Z. Evaluation of in vitro hemolysis and platelet activation of a newly developed maglev LVAD and two clinically used LVADs with human blood. Artif Organs. (2019) 43:870–9. doi: 10.1111/aor.13471

PubMed Abstract | Crossref Full Text | Google Scholar

29. Trusty PM, Tree M, Vincent D, Naber JP, Maher K, Yoganathan AP, et al. In vitro examination of the VentriFlo true pulse pump for failing fontan support. Artif Organs. (2019) 43:181–8. doi: 10.1111/aor.13301

PubMed Abstract | Crossref Full Text | Google Scholar

30. Hsu PL, McIntyre M, Kuetting M, Jack P, Christina E, Rudiger A, et al. Review of recent patents on foldable ventricular assist device. Recent Patents Biomed Eng (Discontinued). (2012) 5(3):208–22. doi: 10.2174/1874764711205030208

Crossref Full Text | Google Scholar

31. Reitan Ö, Epple K. Catheter Pump for Circulatory Support. United States Patent US 20110034874 (2011).

Google Scholar

32. Siess T. Foldable Intravascularly Inserted Blood Pump. United States Patent US 20080103591 (2008).

Google Scholar

33. Röhn D. Radially Compressible and Expandable Rotor for a Fluid Pump. Would Intellectual Property Organization Patent WO 2011076439 (2011).

Google Scholar

34. McBride MW, Mallison TM, Dillon GP, Campbell RL, Boger DA, Hambrie SA, et al. Expandable Impeller Pump. United States Patent US 7927068 (2011).

Google Scholar

35. Schmitz-Rode T, Günther RW. Self-deploying Axial-flow Pump Introduced Intravascularly for Temporary Cardiac Support. United States Patent US 6533716 (2003).

Google Scholar

36. Pfeffer JG, Schmitz-Rode T, Günther RW. Catheter Device. United States Patent US 20090093764 (2009).

Google Scholar

37. Töllner T. Compressible and Expandable Blade for a Fluid Pump. Would Intellectual Property Organization Patent WO 2010149393 (2010).

Google Scholar

38. Töllner T, Scheckel M. Radially Compressible and Expandable Rotor for a Pump Having an Impeller Blade. Would Intellectual Property Organization patent WO 2012007140 (2012).

Google Scholar