朱 橋 輝, 許 少 鋒, 王 子 恒, 陸 俊 杰, 張 學(xué) 昌

(1.浙江大學(xué) 機(jī)械工程學(xué)院, 浙江 杭州 310027;2.浙大寧波理工學(xué)院 機(jī)電與能源工程學(xué)院, 浙江 寧波 315100;3.浙江理工大學(xué) 機(jī)械與自動(dòng)控制學(xué)院, 浙江 杭州 310018 )

0 引 言

人工肝支持系統(tǒng)(artificial liver support system,ALSS)簡稱人工肝,可將人體血液引出至體外機(jī)械裝置幫助凈化人體血液內(nèi)各種有害物質(zhì)、改善血液內(nèi)環(huán)境,暫時(shí)替代和輔助肝臟代謝、解毒等功能[1]．人工肝在臨床上可用于肝病重癥、肝衰竭患者的相關(guān)治療[2]．在對危重癥肺炎患者的治療中[3],人工肝可用于清除炎癥介質(zhì)等有害物質(zhì),補(bǔ)充白蛋白等有益物質(zhì),從而阻斷細(xì)胞因子風(fēng)暴,幫助患者改善呼吸功能,減輕肺炎癥狀,提高危重癥肺炎患者的救治率．中空纖維組件是人工肝的核心裝置,其中有大量的中空纖維管,當(dāng)血液流經(jīng)中空纖維管時(shí),血液內(nèi)的有害物質(zhì)通過中空纖維膜進(jìn)入外腔,從而達(dá)到凈化血液的目的,研究血液在中空纖維管內(nèi)的流動(dòng)行為對研究人工肝具有科學(xué)意義．

Baskurt等[4]的研究表明血液的表觀黏度取決于現(xiàn)有的剪切力,并由紅細(xì)胞比體積、血漿黏度、紅細(xì)胞聚集和紅細(xì)胞的力學(xué)特性決定．Shariatkhah等[5]使用非線性黏彈性模型研究了毛細(xì)血管中周期性發(fā)展的血流,并根據(jù)實(shí)驗(yàn)數(shù)據(jù)與Giesekus 模型推導(dǎo)出了流動(dòng)因子和零剪切速率黏度的最佳值,結(jié)果表明只有非線性黏彈性模型才能準(zhǔn)確描述毛細(xì)血管中血流的實(shí)驗(yàn)數(shù)據(jù)．Apostolidis等[6]研究了動(dòng)脈血流動(dòng)中的非牛頓血液流變學(xué)效應(yīng),結(jié)果顯示采用牛頓模型與非牛頓模型描述血液黏度時(shí)模擬結(jié)果存在顯著差異,這些差異主要?dú)w因于流動(dòng)中低剪切速率區(qū)域和高剪切速率區(qū)域之間存在的耦合．Hund等[7]在Krieger懸浮液模型的基礎(chǔ)上,提出了新的血液黏度模型,解決了先前模型在紅細(xì)胞比體積和剪切速率方面表現(xiàn)出的不連續(xù)性問題．Wu等[8]使用混合理論模擬了微通道的血流,紅細(xì)胞被建模為具有剪切依賴性黏度的非線性流體,并考慮了紅細(xì)胞比體積的影響,結(jié)果顯示在速度場和紅細(xì)胞的體積分?jǐn)?shù)分布方面與實(shí)驗(yàn)觀察結(jié)果非常吻合．Marhefka等[9]對血液在微尺度流動(dòng)進(jìn)行了實(shí)驗(yàn)觀察,結(jié)果顯示血液在微管或小管(直徑小于0.3 mm)內(nèi)流動(dòng)時(shí)在靠近管壁的地方會出現(xiàn)一層薄薄的血漿層,稱為紅細(xì)胞的耗竭層,在耗竭層下游的分支血管中會出現(xiàn)血漿撇渣效應(yīng),導(dǎo)致在較小的血管分支中紅細(xì)胞的減少,可能會損害血液在遠(yuǎn)端毛細(xì)血管中氧氣的運(yùn)輸．其他學(xué)者還發(fā)現(xiàn)了F?hr?us-Lindqvist效應(yīng)、塞子流等現(xiàn)象[10]．

血液的微觀流動(dòng)已經(jīng)得到了大量研究,但對血液中紅細(xì)胞受力的作用機(jī)理和遷移規(guī)律的研究還遠(yuǎn)遠(yuǎn)不夠．為了真實(shí)反映出人工肝中空纖維管內(nèi)的血流動(dòng)力學(xué)和紅細(xì)胞遷移規(guī)律,本文對中空纖維管內(nèi)血液兩相流動(dòng)進(jìn)行三維數(shù)值模擬,構(gòu)建中空纖維管的三維模型,采用歐拉雙流體模型,并使用用戶自定義函數(shù)(UDF)描述紅細(xì)胞升力方程,模擬人工肝中空纖維管內(nèi)血液的流動(dòng)情況,討論壁面升力作用對紅細(xì)胞徑向移動(dòng)的影響,分析中空纖維管內(nèi)血液兩相流體的速度分布、紅細(xì)胞徑向體積分?jǐn)?shù)分布、黏度分布以及壁面條件對紅細(xì)胞體積分?jǐn)?shù)分布的影響,并根據(jù)相關(guān)數(shù)值結(jié)果分析紅細(xì)胞所受升力的作用機(jī)理．研究結(jié)果旨在探討人工肝內(nèi)血液微觀流動(dòng)行為規(guī)律,以期為中空纖維組件的選材與優(yōu)化設(shè)計(jì)提供一定的理論指導(dǎo),提高研發(fā)效率,節(jié)約研發(fā)成本．

1 模型建立

1.1 數(shù)值模型

人體血液是由紅細(xì)胞、白細(xì)胞、血小板、血漿等組成的混合流體,血液中紅細(xì)胞與血漿的比例大于0.99,而白細(xì)胞與血小板所占比例不到0.01,一般忽略不計(jì)．因此血液被看作紅細(xì)胞與血漿組成的兩相系統(tǒng),其中血漿在血液內(nèi)的體積分?jǐn)?shù)εp=0.55,紅細(xì)胞在血液內(nèi)的體積分?jǐn)?shù)εRBC=0.45[11]．對血液流動(dòng)模擬計(jì)算時(shí),血漿被看作液相,紅細(xì)胞被看作顆粒相(固相),由于紅細(xì)胞的體積分?jǐn)?shù)較高,可將顆粒相紅細(xì)胞假設(shè)為擬流體,模擬采用的兩相流模型為歐拉雙流體液-固兩相流模型．

1.1.1 血液兩相流基本方程 采用歐拉雙流體模型,各相均可用統(tǒng)一的質(zhì)量和動(dòng)量守恒方程描述．質(zhì)量守恒方程[12]為

(1)

式中：k為p或RBC,p、RBC分別代表血漿與紅細(xì)胞;ρ為密度;ε為體積分?jǐn)?shù);v為速度矢量．血漿和紅細(xì)胞體積分?jǐn)?shù)之和恒為1,可表示為

(2)

動(dòng)量守恒方程[12]為

(3)

式中：p為流場壓力,τk為應(yīng)力張量,Fk為虛擬質(zhì)量力、浮升力等力源項(xiàng),β為相間動(dòng)量交換系數(shù),g為重力．

1.1.2 本構(gòu)方程 為使上述基本方程組封閉,需給出密度、黏度等參數(shù)對應(yīng)的本構(gòu)方程．混合密度ρmix是由血漿和紅細(xì)胞的密度按其所占的體積分?jǐn)?shù)加權(quán)計(jì)算得到,表達(dá)式為

(4)

1.1.3 血液黏度模型 血液黏度模型是模擬血液動(dòng)力學(xué)流動(dòng)的關(guān)鍵因素．非牛頓剪切稀化模型能夠描述中性稠密浮力懸浮液中高體積分?jǐn)?shù)顆粒的流動(dòng)．采用Carreau-Yasuda模型來模擬血液流動(dòng)的非牛頓特性,血液黏度會隨著紅細(xì)胞體積分?jǐn)?shù)和剪切速率的變化而變化,根據(jù)人體血液的實(shí)驗(yàn)流變學(xué)數(shù)據(jù),計(jì)算得到混合物的具體量綱一黏度表達(dá)式[13-14]為

(5)

(6)

(7)

正常生理?xiàng)l件下,血漿是牛頓流體,血漿黏度[16]為0.001 2 kg/(m·s),而紅細(xì)胞黏度會隨其體積分?jǐn)?shù)的變化而變化,呈現(xiàn)出剪切稀化特性[17]．

液相(血漿)應(yīng)力張量以牛頓形式表示,本構(gòu)方程為

(8)

式中：κp表示血漿的體積黏度,I表示單位應(yīng)力張量．一般條件下,血漿的體積黏度κp設(shè)為0．

固相(紅細(xì)胞)應(yīng)力張量的本構(gòu)方程[18]為

τRBC=-pRBCδ+εRBCμRBC(?v+(?v)T)+

(9)

式中：δ表示Kronecker函數(shù);pRBC表示紅細(xì)胞顆粒間因發(fā)生相互碰撞、排斥作用而產(chǎn)生的壓力,由于紅細(xì)胞的比體積范圍內(nèi)顆粒間產(chǎn)生的壓力微乎其微,可忽略不計(jì);κRBC表示紅細(xì)胞的體積黏度,其數(shù)值為0;紅細(xì)胞黏度μRBC采用Carreau-Yasuda模型計(jì)算得到．

1.1.4 曳力系數(shù) 顆粒相與液相(血漿)之間的曳力采用Schiller-Naumann模型來表示相間動(dòng)量交換系數(shù)β,其表達(dá)式[12]為

(10)

其中曳力系數(shù)Cd表達(dá)式為

(11)

(12)

式中：dRBC表示紅細(xì)胞顆粒的直徑,φ表示形狀因子．紅細(xì)胞顆粒的直徑被假定為8 μm,由于不考慮紅細(xì)胞變形影響,其形狀因子φ為1[15]．

1.1.5 虛擬質(zhì)量力 當(dāng)加速的顆粒相與液相(血漿)發(fā)生碰撞時(shí),由于具有一定慣性作用,顆粒相會受到虛擬質(zhì)量力,其表達(dá)式[19]為

(13)

1.1.6 升力方程 受限剪切流場中,在壁面條件和徑向速度梯度的作用下,考慮到紅細(xì)胞變形作用[20],紅細(xì)胞會受到經(jīng)過其軸心的非慣性升力作用[21-22],該作用力會導(dǎo)致紅細(xì)胞發(fā)生流場內(nèi)的徑向運(yùn)動(dòng),對血液的流場分布、紅細(xì)胞的徑向分布、壁面剪切力產(chǎn)生影響．

兩相流模型中紅細(xì)胞所受的非慣性升力作用,當(dāng)紅細(xì)胞的軸心距離壁面小于5 μm時(shí),升力作用由壁面條件決定;當(dāng)紅細(xì)胞的軸心距離壁面大于5 μm時(shí),升力作用由剪切流場中的速度梯度引起,其表達(dá)式[22-23]為

(14)

式中：U表示量綱一滑移速度,與紅細(xì)胞內(nèi)部、外部流體黏度比有關(guān),在本式中U恒為常數(shù),即U=0.03[19,24];h表示紅細(xì)胞軸心與壁面的距離;f(s)是以s為自變量的量綱一函數(shù)[25-26],s是指與紅細(xì)胞具有相同表面積的球體相比紅細(xì)胞的減少量,假設(shè)紅細(xì)胞在h=5 μm處不再和壁面發(fā)生接觸,由壁面條件所引起的升力[25-26]為0,并假設(shè)f(s)隨著h呈線性變化,當(dāng)h=0 μm時(shí),f(s)=2;當(dāng)h=5 μm時(shí),f(s)=0．

1.2 幾何模型與求解設(shè)置

由于1.1節(jié)中用于血液流動(dòng)模擬的兩相流數(shù)學(xué)模型比較復(fù)雜,在模擬中空纖維管內(nèi)血液流動(dòng)時(shí),為了降低數(shù)值求解過程對計(jì)算機(jī)性能的要求以及相關(guān)計(jì)算參數(shù)的求解,中空纖維管模型采用三維軸對稱結(jié)構(gòu),模型示意圖如圖1所示,中空纖維管截面直徑Ф為100 μm,長度L為5 mm．

圖1 中空纖維管模型示意圖

采用ANSYS ICEM 2020R2對三維結(jié)構(gòu)離散化,模型網(wǎng)格劃分采用六面體網(wǎng)格法,為準(zhǔn)確模擬血液在壁面附近的流動(dòng)情況,通過對壁面邊界層進(jìn)行加密,共設(shè)置了11個(gè)邊界層．其中,模擬時(shí)設(shè)置了3種不同網(wǎng)格數(shù)量的中空纖維管模型進(jìn)行數(shù)值求解,并選用距入口0.5L截面處的紅細(xì)胞徑向體積分?jǐn)?shù)分布結(jié)果進(jìn)行比較,表1給出了不同網(wǎng)格數(shù)量模型下的紅細(xì)胞在壁面鄰近處的最大體積分?jǐn)?shù)比較．由表1可知,不同網(wǎng)格數(shù)量模型下紅細(xì)胞在壁面鄰近處的最大體積分?jǐn)?shù)存在差異,網(wǎng)格數(shù)量為80 735的模型與網(wǎng)格數(shù)量為200 901的模型計(jì)算結(jié)果相差3.404%,而網(wǎng)格數(shù)量為200 901的模型與網(wǎng)格數(shù)量為429 600的模型計(jì)算結(jié)果相差0.488%．因此在保證計(jì)算效率和計(jì)算精度的前提下,本文選用網(wǎng)格數(shù)量為200 901的模型進(jìn)行數(shù)值求解．

表1 不同網(wǎng)格數(shù)量模型下紅細(xì)胞在距入口0.5L截面處的最大體積分?jǐn)?shù)比較

模擬時(shí)血液混合密度為1 059.7 kg/m3．血漿密度為1 030 kg/m3,黏度為0.001 2 kg/(m·s)．紅細(xì)胞顆粒的直徑為8 μm,密度為1 096 kg/m3．紅細(xì)胞的體積分?jǐn)?shù)在入口處恒有εRBC=0.45,液相(血漿)和固相(紅細(xì)胞)在入口處保持相同速度,模擬過程中流速始終保持不變,均為0.2 m/s．出口邊界條件設(shè)置為開放出口,相對壓強(qiáng)為0．根據(jù)中空纖維管直徑與血液流速,計(jì)算得到平均雷諾數(shù)Re為17.66,可看作層流狀態(tài)．使用計(jì)算流體模擬軟件ANSYS FLUENT 2020R2進(jìn)行模擬,迭代計(jì)算采用經(jīng)典相間耦合Simple算法．

2 結(jié)果與討論

2.1 血液兩相流體的速度分布

對中空纖維管內(nèi)血液兩相流體的速度分布特性進(jìn)行分析,在流場作用下,入口邊界條件對血液各相沿徑向速度分布的影響較小,模擬發(fā)現(xiàn)距入口不同位置處,血液兩相流體沿徑向速度分布曲線基本相同．圖2給出了紅細(xì)胞在距入口0.5L截面處的速度分布．由圖2可知,中空纖維管軸心處速度最大,管壁處速度為0,符合層流流動(dòng)特性．根據(jù)圖2的相關(guān)結(jié)果,圖3給出了紅細(xì)胞和血漿在距入口0.5L截面處的速度分布曲線,同時(shí)給出了將血液作為單相牛頓流體的速度分布曲線,橫坐標(biāo)r/r0表示量綱一半徑,r表示紅細(xì)胞或血漿與管內(nèi)中心的距離,r0表示中空纖維管半

圖2 紅細(xì)胞在距入口0.5L截面處的速度分布

圖3 中空纖維管內(nèi)血液兩相流體沿徑向速度

徑．由圖3可知,血液各相的速度呈近似拋物線分布,各相在中空纖維管軸心上的流速大于壁面流速,但血液兩相流模型中血液各相的速度小于血液流動(dòng)(單相牛頓流體)所形成的拋物線速度．血液表現(xiàn)出一定的非牛頓特性,由于顆粒相紅細(xì)胞的存在,顆粒相會拖曳液相(血漿)的流動(dòng),減緩其流動(dòng)速度．其中,徑向各位置處紅細(xì)胞流動(dòng)速度略小于血漿流動(dòng)速度,該結(jié)果與Chandran等[27]采用擴(kuò)散通量模型得到的研究結(jié)論一致．

2.2 紅細(xì)胞徑向體積分?jǐn)?shù)分布

對中空纖維管內(nèi)紅細(xì)胞在距入口不同位置處沿徑向的體積分?jǐn)?shù)分布進(jìn)行統(tǒng)計(jì),共選取了3個(gè)截面,其位置與入口的距離分別為0.2L、0.5L、0.8L,圖4給出了紅細(xì)胞在距入口0.5L截面處沿徑向的體積分?jǐn)?shù)分布．結(jié)果顯示,在徑向上,紅細(xì)胞在壁面處受壁面產(chǎn)生的升力作用較大,紅細(xì)胞遠(yuǎn)離壁面．與入口處紅細(xì)胞體積分?jǐn)?shù)εRBC=0.45相比,壁面處紅細(xì)胞體積分?jǐn)?shù)減小并小于入口處

圖4 在距入口0.5L截面處紅細(xì)胞沿徑向的體積分?jǐn)?shù)分布

紅細(xì)胞體積分?jǐn)?shù),說明紅細(xì)胞向遠(yuǎn)離壁面的方向發(fā)生遷移,與Geislinger等[28]的研究結(jié)論相符．在鄰近壁面處(距離壁面一定距離處),紅細(xì)胞所受的壁面升力作用與流場剪切梯度對紅細(xì)胞產(chǎn)生的升力作用兩者相互平衡[21],此時(shí)紅細(xì)胞體積分?jǐn)?shù)出現(xiàn)最大值并高于入口處紅細(xì)胞體積分?jǐn)?shù)．在管內(nèi)軸心附近處,紅細(xì)胞體積分?jǐn)?shù)略小于入口處紅細(xì)胞體積分?jǐn)?shù)．

根據(jù)不同截面上紅細(xì)胞徑向體積分?jǐn)?shù)分布結(jié)果,圖5給出了紅細(xì)胞在不同截面上沿徑向的體積分?jǐn)?shù)分布曲線．結(jié)果顯示,在管內(nèi)軸心附近,紅細(xì)胞體積分?jǐn)?shù)略小于入口處紅細(xì)胞體積分?jǐn)?shù),而在鄰近壁面處出現(xiàn)最大值,紅細(xì)胞沿徑向體積分?jǐn)?shù)分布呈雙峰狀．從鄰近壁面處至管內(nèi)軸心處,紅細(xì)胞體積分?jǐn)?shù)由最大值逐漸減小至局部最小值;而從鄰近壁面處至壁面處,紅細(xì)胞體積分?jǐn)?shù)急劇減?。酥?對比各截面上紅細(xì)胞在鄰近壁面處的體積分?jǐn)?shù)分布情況,距入口0.2L截面上,紅細(xì)胞體積分?jǐn)?shù)的最大值比其他兩個(gè)截面上的最大值大,隨著各截面與入口距離的增加,各截面上紅細(xì)胞體積分?jǐn)?shù)最大值逐漸減小;而在壁面處,隨著截面與入口距離增加,紅細(xì)胞體積分?jǐn)?shù)逐漸增大,距入口0.8L截面上紅細(xì)胞體積分?jǐn)?shù)最大．

圖5 紅細(xì)胞在不同截面上沿徑向的體積分?jǐn)?shù)分布曲線

2.3 黏度分布

采用血液兩相流模型模擬時(shí),顆粒相紅細(xì)胞被看作擬流體,模擬過程采用Carreau-Yasuda黏度模型描述紅細(xì)胞黏度．圖6給出了紅細(xì)胞在距入口0.5L截面處的黏度分布．紅細(xì)胞在壁面處流動(dòng)速度小,速度梯度大,紅細(xì)胞黏度低;在軸心附近處流動(dòng)速度大,速度梯度小,紅細(xì)胞黏度高．根據(jù)圖6的黏度分布結(jié)果,圖7給出了血液兩相流體沿徑向各位置處的黏度分布曲線．由圖7可知,紅細(xì)胞呈現(xiàn)出非牛頓特性,具有一定的剪切稀化特性．紅細(xì)胞沿徑向黏度變化最大,由于血漿是一種牛頓流體,其黏度不隨體積分?jǐn)?shù)與流速的變化而變化,模擬過程中黏度始終保持不變．血液主要是由紅細(xì)胞與血漿組成的混合溶液,即血液由牛頓流體與非牛頓流體混合而成,整體上表現(xiàn)出一定的非牛頓特性,血液黏度會隨紅細(xì)胞體積分?jǐn)?shù)與剪切速率的變化而變化．

圖6 紅細(xì)胞在距入口0.5L截面處的黏度分布

圖7 血液各相沿徑向的黏度分布曲線

2.4 壁面條件對紅細(xì)胞體積分?jǐn)?shù)分布的影響

由2.2節(jié)可知,紅細(xì)胞在各截面沿徑向的體積分?jǐn)?shù)分布受到壁面升力作用以及血液流動(dòng)對其產(chǎn)生的流體力學(xué)作用的影響．根據(jù)Wu等[8]的相關(guān)研究,壁面對囊泡狀物體的升力會受到壁面剪切速率與囊泡變形作用兩者共同的影響．為了分析紅細(xì)胞在不同壁面條件下所受升力作用對其體積分?jǐn)?shù)分布的影響,本節(jié)共設(shè)置了3組紅細(xì)胞軸心與壁面距離h=0 μm時(shí)的f(s),以模擬血液在中空纖維管不同壁面彈性條件下的流動(dòng)情況．假設(shè)f(s)隨著紅細(xì)胞距壁面的距離h呈線性變化,f(s)越小表征壁面所產(chǎn)生的壁面剪切速率對紅細(xì)胞的變形作用越小時(shí),即當(dāng)紅細(xì)胞距壁面的距離越小,受到來自壁面對紅細(xì)胞產(chǎn)生的非慣性升力越大．圖8給出了不同壁面條件下紅細(xì)胞在距入口0.5L截面處沿徑向體積分?jǐn)?shù)分布曲線．由圖8可知,隨著壁面處的剪切速率與壁面對紅細(xì)胞的變形作用越來越小(即壁面處f(s)越來越小),紅細(xì)胞在鄰近壁面處的體積分?jǐn)?shù)逐漸降低,而在壁面處的體積分?jǐn)?shù)越來越高,表明紅細(xì)胞遠(yuǎn)離壁面的趨勢減弱．

圖8 不同壁面條件下紅細(xì)胞在距入口0.5L截面處沿徑向體積分?jǐn)?shù)分布曲線

為了說明非慣性升力對紅細(xì)胞徑向遷移的影響,圖9給出了紅細(xì)胞在徑向上的遷移機(jī)理示意圖,圖中Fl表示壁面對紅細(xì)胞產(chǎn)生的非慣性升力,Fs表示流場剪切梯度對紅細(xì)胞產(chǎn)生的升力．如圖所示,當(dāng)紅細(xì)胞在壁面處時(shí),所受的升力Fl>Fs,此時(shí)紅細(xì)胞遠(yuǎn)離壁面;當(dāng)紅細(xì)胞在遠(yuǎn)離壁面處于管內(nèi)中心線附近時(shí),主要受流場剪切梯度產(chǎn)生的升力Fs,此時(shí)紅細(xì)胞有沿徑向向壁面遷移的趨勢;當(dāng)紅細(xì)胞在緊鄰壁面處(離壁面距離約為5 μm)時(shí),紅細(xì)胞受到的非慣性升力與剪切梯度產(chǎn)生的升力兩者達(dá)到平衡,即Fl=Fs,此時(shí)紅細(xì)胞會在鄰近壁面處積聚,在該處紅細(xì)胞體積分?jǐn)?shù)出現(xiàn)最大值．當(dāng)壁面處對應(yīng)的f(s)減小時(shí),紅細(xì)胞在壁面處受到的非慣性升力Fl變小,導(dǎo)致徑向上紅細(xì)胞向管內(nèi)中心的遷移趨勢減弱,紅細(xì)胞在鄰近壁面處積聚程度下降,因此紅細(xì)胞在壁面處的體積分?jǐn)?shù)增大．該結(jié)果與Abkarian等[25]模擬的數(shù)值結(jié)果吻合．

圖9 紅細(xì)胞所受非慣性升力與剪切梯度產(chǎn)生的升力作用機(jī)理示意圖

3 結(jié) 論

(1)在層流流場作用下,血液中紅細(xì)胞與血漿的流速呈近似拋物線分布,各相在軸心處流速均大于在壁面處流速,兩相流體速度小于血液單相流動(dòng)形成的拋物線速度,沿徑向各位置處的紅細(xì)胞流速略大于血漿流速．

(2)紅細(xì)胞黏度會隨體積分?jǐn)?shù)與剪切速率變化而變化,壁面處速度梯度大,紅細(xì)胞黏度低;軸心附近處速度梯度小,紅細(xì)胞黏度高,血液整體上表現(xiàn)出非牛頓特性,呈現(xiàn)剪切稀化特性．

(3)紅細(xì)胞在中空纖維管內(nèi)流動(dòng)時(shí),壁面對紅細(xì)胞產(chǎn)生升力作用,使得紅細(xì)胞遠(yuǎn)離壁面,在距離壁面一定位置處,紅細(xì)胞所受壁面升力作用與流場產(chǎn)生的流體力學(xué)作用達(dá)到平衡,紅細(xì)胞體積分?jǐn)?shù)在該位置處出現(xiàn)最大值．隨著截面與入口距離的增加,該截面處紅細(xì)胞體積分?jǐn)?shù)最大值逐漸減小,在壁面處紅細(xì)胞體積分?jǐn)?shù)則逐漸增大．整體上看,紅細(xì)胞沿徑向體積分?jǐn)?shù)分布呈雙峰狀．

(4)壁面對紅細(xì)胞的升力作用會對紅細(xì)胞沿徑向體積分?jǐn)?shù)分布產(chǎn)生影響,隨壁面處剪切速率與壁面對紅細(xì)胞變形作用的減小,紅細(xì)胞在距離壁面一定位置處的體積分?jǐn)?shù)減小,而在壁面處的體積分?jǐn)?shù)增大,紅細(xì)胞由壁面向管內(nèi)中心運(yùn)動(dòng)的趨勢減弱．

本文采用兩相流歐拉方法描述了真實(shí)條件下血液在微觀尺度下的流動(dòng)規(guī)律,準(zhǔn)確設(shè)置了血液兩相流體的黏度模型、曳力模型、紅細(xì)胞升力方程,相關(guān)參數(shù)根據(jù)表達(dá)式使用UDF描述．但由于本文主要采用升力方程來模擬紅細(xì)胞在流動(dòng)時(shí)的受力以反映管壁對紅細(xì)胞遷移的影響,后續(xù)相關(guān)研究還可通過設(shè)置中空纖維管彈性條件,采用流固耦合等方法來模擬管壁對血液流動(dòng)的影響．