基于數值模擬的混流式血泵結構改進

許斌，羅基平，黃典貴

（1.上海理工大學能源與動力工程學院，上海市 200093；2.上海市動力工程多相流動與傳熱重點實驗室，上海市 200093）

心室輔助設備（ventricular assist device，VAD）可以輔助或代替心臟作功輸送血液，維持人體的正常生理運轉，對心衰患者的治愈及生存起到重要作用。VAD 最為核心的部件就是血泵，血泵性能的好壞直接決定VAD 的整體性能優劣[1]。

血泵的性能評估主要包括提供泵血功能的水力性能以及血液相容性[2]。水力性能是指通過血泵轉子葉輪做功傳遞給血液的能量以及能量轉換程度的大小，一般由揚程和水力效率來評估，對于血泵實際應用的需求來說，則是滿足心力衰竭患者體內血液正常生理循環所需的流量及血壓，這是血泵設計符合預期目標的前提條件。血液相容性是指人體內部血液對于外源材料或物質相互作用的適應程度，目前旋轉類血泵對于人體血液的危害主要有溶血及血栓。

混流式血泵能有效緩解軸流式血泵轉速過高加劇血液損傷、離心式血泵容易產生流動死區產生血栓等問題，從而能夠減小對血液的損傷程度，改善血液的流動狀況[3-5]。Wu等[6]通過水力性能實驗探究的方法對一款微型混流式血泵進行流體動力學分析，將螺旋葉片數為3，4 的葉輪結構進行對比，探究混流式血泵流動損失及流動狀況，并利用油滴法在葉輪輪轂上繪制條紋線來探究葉輪輪轂面上的流動情況。研究發現：葉輪的二次流及葉頂間隙泄露與血泵的工作狀況及結構設計有著重要關聯，血泵在低流量下更容易產生流動分離等現象。Shu等[7]利用粒子圖像測速實驗（particle image velocimetry，PIV）對同軸混流式血泵內部泰勒渦進行探究，研究發現：泰勒渦能夠增加環狀間隙流動的混合，縮短紅細胞在該高剪切速率區域的停留時間，能夠對血泵內部的溶血和血栓等現象產生起到抑制作用。

在數值模擬方面，藺嫦燕等[8]利用數值模擬及體外溶血實驗探究的方法對螺旋混流式血泵的流體力學性能及溶血性能進行分析，研究發現：數值模擬結果與體外實驗結果較為吻合，螺旋混流式血泵對血液的損傷程度遠低于對照組的軸流式血泵。Arvand等[9]利用計算流體力學方法（computational fluid dynamics，CFD）及體外性能測試對3 款不同葉片結構的閉式葉輪混流式血泵進行性能分析，研究發現：CFD 模擬的揚程-流量曲線與實驗結果吻合，其所研究的范圍內，葉片形狀為直葉片時的閉式葉輪結構能夠取得更好的溶血性能。羅基平等[10]設計并采用數值模擬方法研究了帶導葉的閉式混流血泵，通過分析該類型血泵的流場特性及壓力分布情況，探討其水力性能以及可能對紅細胞造成的損傷程度。與半開式葉輪結構混流式血泵相比，所研究的血泵具有更好的水力性能及避免血液損傷的能力。Olia等[11]指出血細胞長期暴露于高剪切力、與異物接觸碰撞、空化等因素會對血液造成損傷。

血液流經血泵時，由于內部流動情況復雜，會產生溶血和血栓現象，威脅人體生理健康，如何降低溶血程度及避免不均勻流動產生血栓現象，在血泵的設計及結構改進中尤為重要。血泵的結構參數不僅決定著整體性能的好壞，而且對泵體內部流場分布存在重要影響。通過數值模擬方法探究不同葉片數及葉片厚度對血泵水力性能與血液相容性的影響，得到具有水力性能良好、血液相容性較高的血泵結構參數。

1 數值模擬方法及驗證

1.1 閉式葉輪混流式血泵幾何模型

血泵模型為一款帶閉式葉輪的導葉式混流血泵，如圖1 所示[10]。血泵內部通流部分結構主要為轉子葉輪及后導葉兩部分。轉子葉輪將旋轉得到的動能轉換成流經血泵內部血液的壓能，實現輔助心臟為人體正常生理活動供壓；后導葉起到導流作用，將血液從血泵內部均勻排出，改善血液流動環境。血泵基本結構參數如表1 所示。

表1 混流式血泵的特征參數Tab.1 Characteristic parameters of mixed flow blood pump

圖1 混流式血泵模型Fig.1 Mixed flow blood pump model

1.2 數值模擬方法

假設血液為牛頓流體，血液粘度為3.5×10-3Pa·s，密度為1 055 kg/m3[12]。設定進口流量為5 L/min，出口壓力為20 000 Pa，壁面設置為無滑移條件。采用壓力耦合方程組的半隱式方法對壓力與速度耦合求解，壓力求解用二階中心離散差分法。采用RANS 方法對混流式血泵三維流場進行求解，湍流模型采用標準k-ε模型。

血泵的流場計算域如圖2 所示，計算域分為進口、出口、轉子葉輪、后導葉；進出口段使用結構化網格進行劃分，幾何形狀復雜的轉子葉輪和導葉采用非結構網格，對葉片前、尾緣區進行局部加密，圖3 展示了葉輪與導葉網格。此外，為保證計算的準確性，進行了網格獨立性驗證，結果見表2。當網格量達到3 015 000 時，血泵的揚程和效率隨著網格量的增加變化很小，故選用3 015 000 這套網格對血泵的水力性能與溶血性能進行探討。

圖2 混流式血泵流場計算域Fig.2 Simulation region of the mixed flow blood pump

圖3 網格劃分示意圖Fig.3 Schematic diagram of meshing

表2 網格無關性驗證Tab.2 Grid-independent verification

1.3 數值模擬方法的實驗驗證

根據美國食品監測局（food and drug administration，FDA）的基準泵（benchmark pump）幾何數據[13]，建立計算模型，并采用CFD 方法在文獻中給出的不同流量、轉速的6 個工況條件下進行數值計算，其中，在工況1，4 時，泵轉速為2500 r/min；在工況2，3，5，6 時，泵轉速為3 500 r/min。將本文CFD 模擬結果與FDA 提供的試驗結果、24 組CFD 計算結果及其平均值進行對比，結果如圖4所示。本研究得到的CFD 模擬結果與實驗數據及其他學者模擬的平均值吻合良好，因此，本文采用的數值模擬方法是準確的。

圖4 試驗結果對比Fig.4 Comparison of CFD results

2 不同結構參數對混流式血泵性能的影響

2.1 模型建立

血泵葉片厚度受到尺寸大小的限制，通常在血泵研究當中的葉片厚度為0.5～1 mm 之間[14-15]，本文取血泵厚度c分別為0.6，0.8，1 mm，葉片數Z為3，4，5，6。由于葉片厚度與葉片數的變化均會改變流道的通流面積，故需對葉片厚度與葉片數的影響進行綜合考慮，制定不同葉片數及葉片厚度的血泵模型進行探究。

2.2 水力性能

對于血泵來說，滿足血液的供壓需求是非常重要的。人體的正常供壓在80～120 mmHg 之間，本文中的血泵設計目標是在5 L/min 的標況下達到100 mmHg 的供血需求，在滿足此條件下，揚程接近設計目標，則定義為更符合預期需求[16]。

圖5 和圖6 分別為5 L/min 工況下不同結構參數混流式血泵的揚程與效率曲線，可以看出，在標準工況5 L/min 下，隨著葉片數的增加，揚程均呈現先增大后減小的趨勢，葉片厚度較大時，供壓不足，無法滿足人體正常供血需求。隨葉片數的增加，所有葉片厚度下，血泵工作效率均呈先增大后減小的趨勢；當葉片數較少時，葉片厚度大則效率越高；當葉片數較多時，葉片厚度大則效率低。當Z=4、c=0.6，Z=5、c=0.6及c=0.8，Z=6、c=0.6 能夠滿足100 mmHg 供壓需求，同時具有較高的效率。

圖5 不同結構參數的混流式血泵揚程變化Fig.5 Lift of mixed flow blood pump with different structural parameters

圖6 不同結構參數的混流式血泵效率變化Fig.6 Efficiency of mixed flow blood pump with different structural parameters

2.3 溶血性能

血液紅細胞在血泵中流動時，流場中剪切應力的作用會破壞其結構，破壞程度劇烈則會嚴重威脅人體的正常生理健康。故對于血泵來說，是否具有良好的血液相容性是設計者需要重點考慮的問題。溶血現象主要是由于血液中多種原因導致的紅細胞破裂，溶血程度主要與細胞所受的剪切應力及其作用時間有關，血液損傷程度由紅細胞所受剪切應力與作用時間共同決定，量化標準是溶血指數[17-18]，其表達式為

式中：dHb為游離血紅蛋白量；Hb為血紅蛋白總量；τ為剪切應力標量；t為暴露時間。通過拉格朗日法進行跡線追蹤，在流場中提取血液紅細胞粒子在運動過程中受到的剪切應力標量與時間段進行計算，則可以得到紅細胞在該時間段落上的損傷值

其中剪切應力的標量可由流場計算得到的流場當中的應力張量進行轉換得到，參照Bludszuweit[19]提出的計算方式

式中：τii和τjj為正應力分量；τij為切應力分量。而紅細胞在血泵內部經過時，沿著紅細胞運動軌跡某一時刻受到的累計溶血指數預測值HIp,i為

這樣可以采用若干條紅細胞跡線溶血值相加求平均得到血泵內部的溶血指數值

將粒子直徑設為正常紅細胞的平均直徑7 um進行追蹤，提取不同粒子跡線中的切應力和運動時間數據，計算得到跡線的溶血指數值。圖7 所示為追蹤的一個紅細胞粒子在血泵內部的剪切應力隨運動軌跡及運動時間而變化的示意圖，可以看到該粒子在經過運動情況最為復雜的葉輪區域時所受剪切應力較高。

圖7 紅細胞粒子剪切應力變化情況Fig.7 Shear stress changes of erythrocyte particles

首先對粒子跡線數目對溶血預測結果的無關性進行驗證，如表3 所示，在跡線數目為100 時，隨著跡線數目的增加，溶血指數變化幅度較小，故選用100 條跡線即可滿足跡線數目無關性的要求。

表3 跡線數目無關性驗證Tab.3 Verification of trace number independence

不同葉輪結構的混流式血泵溶血指數如圖8所示。可以看到，相比于原模型，其他結構模型的溶血指數均存在一定的下降，這說明原模型的溶血程度較高，存在改進的空間。為了降低混流式血泵的溶血程度，改善血泵的血液相容性，對原模型進行結構改進是很有必要的。在相同的葉片厚度下，溶血指數隨著葉片數的增加呈現先減小后增大的趨勢，而當葉片數越靠近4 時，血液紅細胞所受到的損傷程度越小。

圖8 不同結構參數的溶血指數分布Fig.8 Distribution of hemolysis index with different structural parameters

2.4 流場分布

流場中存在的不規則流動以及緩流、回流區域容易造成血液的血栓，從而造成血管的堵塞，血栓的形成對血液的流動以及人體的正常生理健康有很大的隱患。圖9 分別顯示了12 種不同結構模型在標準工況5 L/min 下的葉輪流場跡線分布。可以看出，在葉片數為3，4 時，不同葉片厚度下的葉片前緣靠近吸力面一側均產生了一個明顯的流動低速區，速度分布紊亂。當葉片數較少時，葉片進口速度方向改變，液流撞擊到葉片前緣靠近壓力面一側，在前緣區產生脫流現象，引起血液在葉輪背面流道的不穩定流動，由于葉片前緣低速區的存在，容易致使該部分區域血液凝結形成血栓等現象危害人體健康。而葉片數為5、葉片厚度為0.8 mm、1 mm 及葉片數為6 時，流線分布較為合理，沒有明顯的回流及流動停滯現象，更符合血泵穩定均勻流動的需求，能有效地防止血栓現象的產生。

圖9 不同葉片數的葉輪流線圖Fig.9 Impeller streamline diagram with different numbers of blades

2.5 模型改進

綜合考慮揚程、效率、溶血指數以及葉輪內部流場分布情況，本文確定最優葉輪葉片數為5、葉片厚度為0.8 mm。此時，水力性能及溶血性能雖不如葉片數為4、葉片厚度為0.6 mm 的葉輪結構，但是其葉片前緣附近未產生低速回流區，流場分布更加均勻穩定，能有效防止血栓發生。

表4 為混流式血泵結構改進前后的性能對比情況。可以看出，改進后模型揚程為101.57，雖然較原模型有所下降，但與混流式血泵的設計目標更為接近，符合混流式血泵設計預期的血液供壓需求，且改進后水力效率有所提高，有助于血泵的安全高效運行；改進后的血泵溶血指數由原來的4.99×10-6減少到4.26×10-6，比原模型降低14.65%，能夠有效降低對血液紅細胞的損傷，達到預期優化的目標。

表4 結構改進前后血泵性能對比Tab.4 Performance comparison of mixed flow blood pump

圖10 為結構模型改進前后在不同葉高處壓力分布云圖。可以看到，在葉片進口處存在局部低壓區，這會導致空化現象的產生。當模型改進后，低壓區的面積減小，這在一定程度上能有效緩解空化現象的產生，降低葉輪葉片的損壞，增加血泵裝置的使用壽命。

圖10 葉柵壓力分布Fig.10 Pressure distribution in cascade

3 結論

本文以閉式葉輪混流式血泵為研究對象，采用數值模擬方法研究不同葉片數及葉片厚度對混流式血泵水力性能及溶血性能的影響，在所研究的范圍內得到了以下結論：

a.相同的葉片厚度下，溶血指數隨著葉片數的增加呈現減小后增加的趨勢，而當葉片數數目越靠近4 時，血液紅細胞所收到的損傷程度更小；葉片數為5 和6 時，流線分布較為合理，更適合內部的流動需求。

b.當葉片數為5、葉片厚度為0.8 mm 時，在流動均勻不變的情況下，揚程較原模型更符合預期設計目標，更易調節；溶血指數下降了14.65%，可有效降低對血液的損傷；葉片進口低壓區域減小，有效抑制了空化現象產生；內部流動均勻穩定，未出現回流、流動死區等問題，可有效防止血栓發生。

c.較原模型相比，改進后葉片數更少，葉片厚度更大，這在一定程度上能夠降低制造加工的難度及成本。