計(jì)算流體力學(xué)(CFD)在心血管系統(tǒng)中的應(yīng)用是其在醫(yī)學(xué)領(lǐng)域中最成熟、最深入的方向之一。它通過(guò)數(shù)值模擬血液流動(dòng)的動(dòng)力學(xué)特征如,流速、壓力、剪切應(yīng)力等,結(jié)合患者的影像數(shù)據(jù)重建個(gè)性化血管模型,為理解心血管疾病的發(fā)病機(jī)制、優(yōu)化診療方案提供了量化依據(jù)。以下是其核心應(yīng)用場(chǎng)景及解決的問(wèn)題:
1、動(dòng)脈粥樣硬化的機(jī)制研究與風(fēng)險(xiǎn)預(yù)測(cè)
動(dòng)脈粥樣硬化是冠心病、腦梗死等疾病的核心病因,其斑塊易在血管分叉、彎曲處,如頸動(dòng)脈分叉、冠狀動(dòng)脈左前降支的形成,而血流動(dòng)力學(xué)異常是關(guān)鍵誘因。
應(yīng)用場(chǎng)景:
可解決的問(wèn)題:
2、動(dòng)脈瘤的生長(zhǎng)與破裂風(fēng)險(xiǎn)評(píng)估
動(dòng)脈瘤,如顱內(nèi)動(dòng)脈瘤、腹主動(dòng)脈瘤,是血管壁局部膨出形成的“薄弱點(diǎn)”,其破裂可能導(dǎo)致致命性出血,而血流對(duì)瘤壁的力學(xué)作用是決定其生長(zhǎng)和破裂的核心因素。
應(yīng)用場(chǎng)景:
可解決的問(wèn)題:
量化破裂風(fēng)險(xiǎn):對(duì)未破裂動(dòng)脈瘤,通過(guò)CFD參數(shù)區(qū)分“高風(fēng)險(xiǎn)”與“低風(fēng)險(xiǎn)”動(dòng)脈瘤,避免過(guò)度治療。
優(yōu)化治療方案:評(píng)估手術(shù)干預(yù)后的血流動(dòng)力學(xué)變化,預(yù)測(cè)治療效果,如是否有效降低瘤壁壓力。
3、心臟瓣膜疾病的評(píng)估與人工瓣膜優(yōu)化
心臟瓣膜的狹窄或反流會(huì)導(dǎo)致血流紊亂,而人工瓣膜的設(shè)計(jì)需平衡血流動(dòng)力學(xué)性能與并發(fā)癥風(fēng)險(xiǎn)。
應(yīng)用場(chǎng)景:
4、先天性心臟病的血流動(dòng)力學(xué)分析
先天性心臟病,如室間隔缺損、法洛四聯(lián)癥、主動(dòng)脈縮窄等,會(huì)導(dǎo)致血流路徑異常,引發(fā)心功能損害,而CFD可量化異常血流對(duì)心臟和血管的影響。
應(yīng)用場(chǎng)景:
可解決的問(wèn)題:
5、支架植入與血管重建的優(yōu)化
冠心病、外周動(dòng)脈疾病等需通過(guò)支架植入擴(kuò)張狹窄血管,但術(shù)后可能出現(xiàn)再狹窄或支架血栓,血流動(dòng)力學(xué)是關(guān)鍵影響因素。
應(yīng)用場(chǎng)景:
可解決的問(wèn)題:
6、心功能評(píng)估與心力衰竭機(jī)制研究
心力衰竭的核心是心臟泵血功能下降,而心室血流動(dòng)力學(xué)與心肌收縮功能密切相關(guān)。
應(yīng)用場(chǎng)景:
可解決的問(wèn)題:
門靜脈是肝臟的“物流總閘”,其位置如下圖(a)所示。血流量占心輸出量的25 %,卻沒(méi)有任何瓣膜;一旦血流動(dòng)力學(xué)失衡,就可能引發(fā)門脈高壓、血栓乃至肝性腦病。受制于倫理與成像極限,活體測(cè)量門靜脈全域流速、壁面剪切應(yīng)力(WSS)幾乎不可能。
于是我們把血管“搬進(jìn)”計(jì)算機(jī),截取一段如下圖(b)所示,用CFD讓在血管內(nèi)的行為變得可觀測(cè)、可量化、可干預(yù)。
在此案例中,采用上海積鼎信息科技有限公司的自研CFD軟件VirtualFlow 2025的非結(jié)構(gòu)求解器進(jìn)行仿真。該求解器支持任意的非結(jié)構(gòu)網(wǎng)格,仿真精度不亞于主流上軟,同時(shí)支持原生GPU計(jì)算,可大幅提升仿真速度。
a、CPU求解結(jié)果與主流商軟對(duì)比
采用VirtualFlow進(jìn)行網(wǎng)格劃生成軟件生成四面體網(wǎng)格,分別通過(guò)國(guó)外F軟件和自研求解器用該網(wǎng)格進(jìn)行仿真。從仿真的結(jié)果對(duì)比可以看出血管壁面上的壓力分布基本一致,最大高壓絕對(duì)值幾乎一致,最低低壓絕對(duì)值相差10%左右。
b、GPU求解結(jié)果與主流商軟對(duì)比
積鼎科技VirtualFlow 2025版已支持GPU求解能力,其特色表現(xiàn)為:
1. 同時(shí)支持CPU和GPU計(jì)算
2. 支持任意多面體網(wǎng)格,適用于具有復(fù)雜幾何外形的快速網(wǎng)格生成
3. 采用穩(wěn)定高效的數(shù)值求解方法,可以求解多種流動(dòng)特別是多相流問(wèn)題
4. 支持定常和非定常求解
5. 支持多種湍流模型,包括但不限于k-epsilon、sst k-omega等
6. 具有全面的前后處理功能
7. 具有腳本功能,能夠支持自動(dòng)化運(yùn)行
8. AI友好,可與AI模型連接
仿真結(jié)果對(duì)比
從上圖的仿真結(jié)果比較來(lái)看,VirtualFlowGPU版與商軟國(guó)外S軟件的計(jì)算結(jié)果幾乎一致,可見(jiàn)VirtualFlow GPU版本求解器數(shù)學(xué)模型的精度是可以信賴的。
仿真性能對(duì)標(biāo)
國(guó)外S軟件核并行時(shí)間
自研GPU版1卡并行時(shí)間
國(guó)外S軟件在16核配置下仿真500步,耗時(shí)378s。相同的網(wǎng)格與模型配置下,采用VirtualFlow非結(jié)構(gòu)求解器求解,在1塊GPU顯卡的情況下,同樣計(jì)算500步,需要51s??梢?jiàn)自研非結(jié)構(gòu)求解器的GPU版本計(jì)算效率是商軟CPU計(jì)算的7~8倍。
c、與離散元軟件DEMms耦合
在多物理場(chǎng)耦合數(shù)值方法及 CFD-DEM 單向耦合解析的基礎(chǔ)上,針對(duì)紅細(xì)胞與血流場(chǎng)的相互作用,雖正常生理狀態(tài)下紅細(xì)胞體積分?jǐn)?shù)為 35%~50%,屬于高濃度顆粒體系,二者相互作用具有強(qiáng)耦合特性,但為演示單向耦合,特設(shè)計(jì)如下測(cè)試算例。
該測(cè)試算例的核心是對(duì)紅細(xì)胞相關(guān)參數(shù)進(jìn)行特殊設(shè)定,以滿足單向耦合的適用條件。在算例中,我們將紅細(xì)胞的體積分?jǐn)?shù)大幅降低至5% 以下,使其處于稀薄顆粒體系范疇。此時(shí),可近似認(rèn)為紅細(xì)胞對(duì)血流場(chǎng)無(wú)顯著擾動(dòng),僅需考慮血流場(chǎng)對(duì)紅細(xì)胞的作用力,符合單向耦合 “上游計(jì)算結(jié)果驅(qū)動(dòng)下游仿真,下游結(jié)果不反作用于上游” 的特征。
基于積鼎的顆粒流求解器DEMms,與非結(jié)構(gòu)求解器耦合,實(shí)現(xiàn)紅細(xì)胞在血管中流動(dòng)的耦合仿真。通過(guò)下圖耦合仿真結(jié)果可見(jiàn),紅細(xì)胞在中間交匯處容易聚集。通過(guò)耦合仿真,可以對(duì)血栓形成的研究體統(tǒng)理論依據(jù)。
1. VirtualFlow2025版中的求解器可分結(jié)構(gòu)模塊求解器和非結(jié)構(gòu)模塊求解器。比于傳統(tǒng)國(guó)外F軟件等商軟,只能提供一種模式,軟件的仿真場(chǎng)景有更大適用性。
2. 同時(shí)支持CPU和GPU計(jì)算
a. GPU原生支持,性能不亞于國(guó)際主流GPU求解器
b. CPU支持多線程和多進(jìn)程混合并行模式
c.同時(shí)適配國(guó)際主流和國(guó)內(nèi)多種CPU/GPU
3. 支持任意多面體網(wǎng)格,適用于具有復(fù)雜幾何外形的快速網(wǎng)格生成
4. 采用穩(wěn)定高效的數(shù)值求解方法,可以求解多種流動(dòng)特別是多相流問(wèn)題
5. 支持定常和非定常求解
6. 支持多種湍流模型,包括但不限于k-epsilon、sst k-omega等
具有全面的前后處理功能
a. 高度自動(dòng)化網(wǎng)格劃分能力
b.高性能仿真結(jié)果后處理
7. 具有腳本功能,能夠支持自動(dòng)化運(yùn)行
8. AI友好,可與AI模型連接
四、總結(jié)
CFD在心血管系統(tǒng)中的應(yīng)用核心是將“形態(tài)學(xué)影像”轉(zhuǎn)化為“力學(xué)量化指標(biāo)”,解決了傳統(tǒng)影像學(xué)無(wú)法回答的“力學(xué)機(jī)制”問(wèn)題:
隨著CFD技術(shù)與醫(yī)學(xué)研究的結(jié)合,其在臨床中的轉(zhuǎn)化應(yīng)用將更加廣泛,推動(dòng)心血管疾病診療從“定性描述”向“定量預(yù)測(cè)”發(fā)展。
官方微信
掃碼關(guān)注