Evaluation AC7-03: Difference between revisions
No edit summary |
m (Mike moved page Lib:Evaluation AC7-03 to Evaluation AC7-03 over redirect) |
||
(14 intermediate revisions by one other user not shown) | |||
Line 4: | Line 4: | ||
}} | }} | ||
__TOC__ | __TOC__ | ||
= | =Flow in a Ventricular Assist Device - Pump Performance & Blood Damage Prediction= | ||
'''Application Challenge AC7-03''' © copyright ERCOFTAC | '''Application Challenge AC7-03''' © copyright ERCOFTAC 2022 | ||
=Evaluation= | =Evaluation= | ||
Line 38: | Line 39: | ||
|} | |} | ||
''' | ''' Effective Stresses <math> \langle \tau_{eff} \rangle </math> in the time-averaged flow field ''' | ||
The computed effective stresses <math> \tau_{eff} </math> (Eq. ( | The computed effective stresses <math> \langle \tau_{eff} \rangle </math> (Eq. (10) in section [[Description AC7-03|Description]]) are presented for both operating points in Fig. 5.2. The stresses from the reference LES are compared to the URANS computations (The instantaneous, effective stresses are dispalyed in Fig. 1.3 in section [[Description AC7-03|Description]]). As can be seen from the LES results, relevant stresses above 9 Pa (threshold for vWF degradation) and 50 Pa (platelet activation) are present within the flow channel of the rotor and the outlet guide vane. In general, the stresses are underpredicted with URANS. Nevertheless, similar hot-spots for significant stresses <math> \geq 9~Pa </math> are observable for the partial load (<math> Q=2.5~l/min </math>). This is in contrast to the nominal operation point (<math> Q=4.5~l/min </math>), where greater deviations in computed stresses are noticable between LES and URANS. Despite the fact that the URANS can reflect the high stresses in the gap vortex and the trailing edge flow regions (red areas marked with (A) and (B)), other relevant regions in the blade channels of the impeller (C) or in the area between the impeller and the outlet guide vane (D) cannot be reflected adequately for the nominal operation point. In these areas of the pump, complex interactions occur between secondary flows are present (explained in Refs. [2] and [27]), which are resolved directly by the LES. In contrast, the URANS turbulence model cannot model the impact of these complex turbulent flow interactions on the effective stress field adequately . | ||
[[Image: | [[Image:AC7-03 LES URANS Effective Stresses.png|620px|center|thumb|Fig. 5.2. Effective stresses <math> \langle \tau_{eff} \rangle </math> of LES and URANS. The figure displays a cylindical cut (explained in Fig. 1.3. in [[Description AC7-03|Description]]) through the rotor and outlet guide vane at a radius of 80% of the outer radius <math> R_2 </math>. Top row: partial load at <math> Q=2.5~l/min </math>. Bottom row: nominal load at <math> Q=4.5~l/min </math>. Note 1: The "hard" change in stresses directly between impeller and outlet guide vane domain is present due to averaging procedures in different frames of references. Note 2: The blades (I,II and 1-3) and the blade channels are mirrored to get an overview of the entire blade channels.]] | ||
''' Hemodynamical Evaluation: Hemolysis Index <math> {MIH} </math> and Volumetric Threshold Analysis <math> I_{\tau_{eff}} </math>''' | ''' Hemodynamical Evaluation: Hemolysis Index <math> {MIH} </math> and Volumetric Threshold Analysis <math> I_{\tau_{eff}} </math>''' | ||
The computed MIH indices (defined by Eq. ( | The computed MIH indices (defined by Eq. (12) in section [[Description AC7-03|Description]]) are shown in Table 5.2. As already displayed in the stress fields in Fig. 5.2., the LES computes higher hemolysis values thans URANS in both operating points. This is partly due to the coarser spatial and temporal resolution of the stresses in URANS (also recognizable from the grid convergence study in Fig. 3.3 in section [[CFD Simulations AC7-03|CFD Simulations]]). But also the already mentioned insufficient resolution of the secondary structures in the blade channels leads to the larger deviations at at <math> Q=4.5~l/min </math>. | ||
{|border="1" cell padding="25" cell spacing="3" align="center" | {|border="1" cell padding="25" cell spacing="3" align="center" | ||
|+align="bottom"|Table 5.2 Modified index of Hemolysis <math> {MIH}~[-] </math> . The constants <math> (C, \alpha, \beta) </math> of Zhang et al. [30] were used in Eq. ( | |+align="bottom"|Table 5.2 Modified index of Hemolysis <math> {MIH}~[-] </math> . The constants <math> (C, \alpha, \beta) </math> of Zhang et al. [30] were used in Eq. (12). The percentage value in brackets indicate the relative deviation between URANS and LES. | ||
! ''<math> MIH </math> by Eq. ( | ! ''<math> MIH </math> by Eq. (12)'' !! colspan="1"| ''partial load <math> Q=2.5~l/min </math>'' !! colspan="1"| ''nominal load <math> Q=4.5~l/min </math>'' | ||
|- | |- | ||
|LES - reference || <math> 83.4 </math>|| <math> 59.1 </math> | |LES - reference || <math> 83.4 </math>|| <math> 59.1 </math> | ||
Line 59: | Line 60: | ||
|} | |} | ||
Tables 5.3. and 5.4. show the percentage of the entire VAD volume in which certain stress thresholds for van-Willebrand-degradation (vWF; <math> >9 Pa </math>), platelet activation (<math> >50 Pa </math>) and hemolysis (<math> >150 Pa </math>) are exceeded. Again, the computed results of URANS are lower than those from LES. Relative deviations of maximal <math> -15.0\% </math> (<math> Q = 4.5~l/min </math>) are observable for the stress thresholds of <math> 9~Pa </math> and <math> 50~Pa </math>. Only for the stress threshold above <math> 150~Pa </math>, larger deviations are present for | Tables 5.3. and 5.4. show the percentage of the entire VAD volume in which certain stress thresholds for van-Willebrand-degradation (vWF; <math> >9 Pa </math>), platelet activation (<math> >50 Pa </math>) and hemolysis (<math> >150 Pa </math>) are exceeded. Again, the computed results of URANS are lower than those from LES. Relative deviations of maximal <math> -15.0\% </math> (<math> Q = 4.5~l/min </math>) are observable for the stress thresholds of <math> 9~Pa </math> and <math> 50~Pa </math>. Only for the stress threshold above <math> 150~Pa </math>, larger deviations are present for both URANS cases, which is due to the coarser near-wall grid density, where the highest stresses are present. These near-wall stresses affect the numerical hemolysis prediction greatly (see Ref. [29] or [31]). | ||
As can be seen from these hemodynamic results, the discrepancies between URANS and LES in the evaluation of blood damage are many times greater than those for the previously evaluated fluid mechanical parameters such as head or efficiency. | As can be seen from these hemodynamic results, the discrepancies between URANS and LES in the evaluation of blood damage are many times greater than those for the previously evaluated fluid mechanical parameters such as head or efficiency. | ||
Line 80: | Line 81: | ||
|- | |- | ||
|URANS || <math> 9.64~(-15.0\%) </math> || <math> 0.84~(-13.4\%) </math> || <math> 0.10~(-23.1\%) </math> | |URANS || <math> 9.64~(-15.0\%) </math> || <math> 0.84~(-13.4\%) </math> || <math> 0.10~(-23.1\%) </math> | ||
|- | |- | ||
|} | |} |
Latest revision as of 10:52, 11 January 2023
Flow in a Ventricular Assist Device - Pump Performance & Blood Damage Prediction
Application Challenge AC7-03 © copyright ERCOFTAC 2022
Evaluation
Experimental Validation of URANS and LES
For simulations in turbopumps, hydraulic characteristics such as the head are among the most important result variables. Also in the field of CFD applications in VADs, it is common to use the head to validate the numerical calculation. In fact, the measurement of the head is the only experimental validation of the simulated flow field in a large number of literature studies (e.g., in [4], [5], [28], [29]). Since the pressure in the pump is coupled via the governing equations to the flow field, the comparison of the heads can be used as the first stage of flow field validation in VADs.
The experimental and numerical results are given in Figure 5.1. For the operating point at , good agreement between numerical and experimental results can be observed with a deviation of for LES and for URANS. For the smaller flow rate the deviations are slightly larger with for LES and for URANS.
In summary, the discrepancy between numerically and experimentally determined head is still within an acceptable range for both LES and URANS. From this, it can be concluded that both numerical models are valid for reproducing the real pressure buildup of the VAD.
Fluid Mechanical & Hemodynamical Evaluation of LES and URANS
Inner Efficiencies
The inner efficiencies according to Eq. (2) of the impeller (index: ) and the whole pump (index: ) are given in Tab. 5.1. The deviation between the URANS and the reference LES case are minor for both operation points with a maximum deviation of . It can be concluded that the turbulence-modelling URANS method can reflect the efficiencies, and hence the global losses, as accurately as the turbulence-resolving LES method.
Flow Rate | Parameter | LES | URANS | Relative deviation to LES [%] |
---|---|---|---|---|
Effective Stresses in the time-averaged flow field
The computed effective stresses (Eq. (10) in section Description) are presented for both operating points in Fig. 5.2. The stresses from the reference LES are compared to the URANS computations (The instantaneous, effective stresses are dispalyed in Fig. 1.3 in section Description). As can be seen from the LES results, relevant stresses above 9 Pa (threshold for vWF degradation) and 50 Pa (platelet activation) are present within the flow channel of the rotor and the outlet guide vane. In general, the stresses are underpredicted with URANS. Nevertheless, similar hot-spots for significant stresses are observable for the partial load (). This is in contrast to the nominal operation point (), where greater deviations in computed stresses are noticable between LES and URANS. Despite the fact that the URANS can reflect the high stresses in the gap vortex and the trailing edge flow regions (red areas marked with (A) and (B)), other relevant regions in the blade channels of the impeller (C) or in the area between the impeller and the outlet guide vane (D) cannot be reflected adequately for the nominal operation point. In these areas of the pump, complex interactions occur between secondary flows are present (explained in Refs. [2] and [27]), which are resolved directly by the LES. In contrast, the URANS turbulence model cannot model the impact of these complex turbulent flow interactions on the effective stress field adequately .
Hemodynamical Evaluation: Hemolysis Index and Volumetric Threshold Analysis
The computed MIH indices (defined by Eq. (12) in section Description) are shown in Table 5.2. As already displayed in the stress fields in Fig. 5.2., the LES computes higher hemolysis values thans URANS in both operating points. This is partly due to the coarser spatial and temporal resolution of the stresses in URANS (also recognizable from the grid convergence study in Fig. 3.3 in section CFD Simulations). But also the already mentioned insufficient resolution of the secondary structures in the blade channels leads to the larger deviations at at .
by Eq. (12) | partial load | nominal load |
---|---|---|
LES - reference | ||
URANS |
Tables 5.3. and 5.4. show the percentage of the entire VAD volume in which certain stress thresholds for van-Willebrand-degradation (vWF; ), platelet activation () and hemolysis () are exceeded. Again, the computed results of URANS are lower than those from LES. Relative deviations of maximal () are observable for the stress thresholds of and . Only for the stress threshold above , larger deviations are present for both URANS cases, which is due to the coarser near-wall grid density, where the highest stresses are present. These near-wall stresses affect the numerical hemolysis prediction greatly (see Ref. [29] or [31]).
As can be seen from these hemodynamic results, the discrepancies between URANS and LES in the evaluation of blood damage are many times greater than those for the previously evaluated fluid mechanical parameters such as head or efficiency.
LES - reference | |||
URANS |
LES - reference | |||
URANS |
Contributed by: B. Torner — University of Rostock, Germany
© copyright ERCOFTAC 2022