Large-eddy simulation of the FDA benchmark blood pump: validation against experiments and implications for turbulent flow mechanisms
8000萬網格大渦模擬解析95%動能，成功破解FDA血泵亂流預測難題。
- 傳統RANS無法解析瞬態渦流，導致擴壓器噴流預測嚴重失真。
- 8000萬網格模型解析超過95%湍流動能，達到高保真模擬標準。
- WALE次網格模型結合瞬態滑動介面，精確重現極端剪切流動。
在流體動力學領域，準確預測心室輔助裝置內的血液流動是一大挑戰。最新研究透過高達 8000 萬網格的大渦模擬結合瞬態技術，成功破解美國 FDA 離心式血泵基準模型的亂流預測難題。數據證實，唯有當解析的湍流動能比例超過 95% 時，計算流體力學才能真正捕捉高度非穩態的血泵流場特徵。
離心式血泵與 FDA 基準模型的亂流預測難題
心室輔助裝置（VAD）是治療心臟衰竭的關鍵技術，近年以磁力懸浮轉子為核心的第三代離心式血泵逐漸成為主流。相較於轉速極高的軸流式血泵，離心式裝置具備較寬闊的流道，能在 2000 至 7000 rpm 的較低轉速下達到每分鐘 3 至 7 公升的流量。此類裝置的雷諾數（Reynolds number）通常介於 10^5 至 10^6 之間，恰好落在流體由層流轉變為亂流的臨界邊緣。由於多種流動型態在泵內同時發展，產生了極度複雜的內部流場。
為了評估計算流體力學（CFD）在醫療器材的預測精準度，美國 FDA 建立了離心泵的基準模型，並由獨立實驗室利用粒子影像測速（PIV，Particle Image Velocimetry）收集實驗數據。多年來，各大研究機構始終無法透過 CFD 穩定且精準地重現擴壓器（diffuser）區域的速度分佈。這項缺失暴露了既有模擬方法的侷限性。準確預測速度場不僅影響水力性能的評估，更是估算高剪切應力是否會造成紅血球破裂等血液損傷指標的核心依據。
突破傳統 RANS 限制的大渦模擬（LES）架構
過去針對 FDA 基準泵的數值模擬，大多仰賴雷諾平均納維-斯托克斯（RANS，Reynolds-Averaged Navier-Stokes）方程。RANS 架構雖然計算成本較低，但它僅能捕捉整體平均流動行為，無法解析瞬態的亂流結構細節。為了解決這項痛點，研究團隊導入了尺度解析的大渦模擬（LES，Large-eddy simulation）技術。LES 透過直接計算大尺度且富含能量的亂流結構，並藉由次網格模型處理無法被網格解析的微小渦流，藉此還原高度非穩態的流場特性。
在處理次網格尺度的湍流應力時，團隊選擇了壁面適應局部渦流黏度（WALE，Wall-Adapting Local Eddy-viscosity）模型。相較於經典的 Smagorinsky 模型，WALE 模型無需依賴近壁阻尼函數，即可在固體邊界附近產生精確的局部耗散，並正確重現渦流黏度的近壁縮放特性。考慮到離心血泵中轉子與外殼之間僅有狹窄的 4 毫米間隙，且在極高轉速下會產生由剪切力主導的強烈流動，WALE 模型的物理特性完美契合了這種極端環境。對照組的 RANS 模擬則採用經典的 k-ω SST 雙方程模型，作為評估 LES 優勢的基準。
解決定子與轉子之間的相對運動，是血泵模擬的另一項核心挑戰。多重參考系（MRF，Multiple Reference Frame）方法習慣將流場視為穩態來估算旋轉效應，但在極窄的間隙中常出現非物理的預測結果。因此，研究採用瞬態滑動介面（SI，Sliding Interface）技術，讓網格隨著轉子進行剛體旋轉，並在交界面上即時交換數據。這種方法對時間步長要求極高，卻能精準呈現葉片掃過時產生的週期性強烈擾動。
8000萬網格與1080核CPU的 LES 解析度驗證
為避免 LES 計算落入網格依賴的陷阱，研究建構了 1000 萬、5100 萬與 8000 萬網格數量的三組模型。這三組模型依序搭配對應於 0.6°、0.3° 與 0.225° 旋轉增量的極小時間步長。網格架構採用混合離散化策略，在壁面附近配置了 14 層邊界層，確保最靠近壁面的 y+ 值小於 1，以直接解析黏性底層內的劇烈速度梯度。在耗能龐大的瞬態求解過程中，最精細的 8000 萬網格模型使用了 1080 個 Intel Xeon 8581 CPU 核心，轉子每旋轉一圈就需要消耗約一天的實際運算時間。
評估 LES 品質的核心指標，在於確認網格是否解析了足夠比例的湍流動能（TKE，Turbulent Kinetic Energy）。根據學界公認的 Pope 標準，一個具備工程參考價值的 LES 至少需要解析 80% 的總湍流動能。測量結果顯示，1000 萬網格的解析度僅落在 80% 左右的邊緣標準，5100 萬網格達到了 91%。而 8000 萬網格則突破了 95% 的高水準，證明幾乎所有主導能量傳遞的渦流運動都已被直接計算，而非依賴次網格模型的經驗公式。
除了檢驗湍流動能比例，研究進一步導入次網格活動參數（Subgrid activity parameter）來量化局部能量傳遞的細節。該參數藉由比較渦流黏度與流體物理黏度的比例，來判定網格模型是否過度依賴次網格的建模耗散。當該參數數值小於 0.3 時，代表絕大部分的能量級聯過程是由直接解析的大尺度網格所主導，建模與離散化誤差皆被控制在 1% 以下。數據分析證實，8000 萬網格的計算結果完美符合了此項嚴苛指標，確保局部流場耗散率的高度保真。
PIV 粒子影像測速與擴壓器流場精準預測
進入實驗數據的驗證階段，團隊將長達 30 圈旋轉的模擬結果進行相位平均，並與 FDA 公布的 PIV 量測數據逐一比對。模擬結果顯示，葉片通道內的流場能迅速收斂，但擴壓器內部的流動異常複雜。特別是在大流量、高雷諾數的操作條件 5 下，噴射流對側的減速區需要極長的採樣週期才能穩定。為忠實還原實驗室的觀測情境，模擬數據嚴格遵照每旋轉 90° 擷取一次樣本的頻率，與 PIV 系統利用相機觸發機制捕捉影像的方式完全一致。
在不同操作條件的嚴格檢驗下，各網格解析度的 LES 展現了明顯的層級差異。在擴壓器區域中，1000 萬與 5100 萬網格高估了噴流的峰值速度，並且噴流核心位置皆出現偏移。唯有最細緻的 8000 萬網格模型，成功鎖定了正確的噴流位置，將大部分速度輪廓線拉回實驗室測得的標準誤差區間內。這項進展解決了過去多項研究皆無法克服的速度高估問題，大幅超越既往的預測水準。
為建立清晰的對照組，研究團隊同步採用結合滑動介面的非穩態 RANS（URANS-SI）以及穩態 RANS（RANS-MRF）進行交叉測試。數據顯示，RANS-MRF 大幅低估了葉片尖端間隙（半徑 26 至 30 毫米處）的流速，並在葉片根部出現不自然的空間振盪。在擴壓器中，無論是 URANS 或穩態 RANS 皆無法正確模擬強烈的間歇性與壁面邊界亂流，導致噴射流預測完全失真。相對而言，結合 SI 耦合的 LES 展現了壓倒性的準確度，證實尺度解析與瞬態模擬在離心式血泵研究中具有無可取代的地位。
大渦模擬精準解析血泵瞬態亂流，為血液相容性評估與設備優化確立了全新科學標竿。