CFD Rupture Challenge 2013 Phase I — 26グループは同じ形状から同じ答えを出せたか
The Computational Fluid Dynamics Rupture Challenge 2013—Phase I: prediction of rupture status in intracranial aneurysms
同じ形状データを渡したとき、CFD の結果はどれだけ揃うのか。揃わない原因はどこにあるのか。
26グループが同一の2症例を盲検で解析。破裂側の判定は21/26が正解したが、破裂部位の予測は6箇所に散って一致しなかった。ばらつきの主因は解析手法ではなく境界条件で、条件を指定した Phase II ではほぼ同一の流れ場が得られている。
当社の視点
「文献値と自分の結果が合わない」という相談に、最も直接的な根拠を与える一報です。参加グループが自由に選んだ流量は 0.9〜11.36 mL/s と13倍の開きがあり、実測値の 1.76 / 2.5 mL/s から大きく外れたものが含まれていました。入口流速 10 m/s や Reynolds 数 1 という明らかに非現実的な設定も提出されています。つまり結果の散らばりはソルバーの差ではなく、入力条件の差です。Phase II で条件を揃えるとほぼ同一の流れ場になったという事実が、この解釈を裏づけています。
扱う範囲
-
01
チャレンジの設計
破裂・未破裂の中大脳動脈瘤2症例の内腔形状だけを配布し、破裂側と破裂部位を盲検で当てさせた。メッシュ・流量・粘性モデル・境界条件・ソルバーはすべて各グループの自由。26グループ（15か国）と脳外科医43名が参加。
-
02
破裂側の判定
CFD は 21/26（81%）が正解、脳外科医は 39/43（91%）が正解。外科医の根拠はブレブ（56%）と形状の不整（44%）で、形態情報だけでこの水準に達している。
-
03
破裂部位の予測は一致しなかった
予測は瘤内の6箇所に散った。低 WSS と高圧を根拠に小さなドームのブレブを選んだのが10グループ、高 WSS と圧勾配からネックを選んだのが4グループ。脳外科医は43名全員が部位を外した。
-
04
入力条件の散らばり
同じ形状に対して、メッシュ 28.3万〜1797万要素、時間刻み 1×10⁻⁵〜5×10⁻³ 秒、流量 0.9〜11.36 mL/s、Reynolds 数 134〜815 と大きく割れた。メッシュ依存性を検証したと報告したのは26グループ中5グループのみ。
-
05
Phase II が示したこと
境界条件を指定した Phase II では、ほぼすべてのグループが同一の流れ場を再現できた。したがってばらつきの支配要因は CFD の手法そのものではなく、境界条件の選択にある。
同一形状に対して各グループが選んだ条件の幅
| 項目 | 値 | 補足 |
|---|---|---|
| 参加 | 26グループ / 15か国 | 脳外科医43名が独立に参加 |
| 破裂側の判定 | 21/26（81%） | 脳外科医は 39/43（91%） |
| 破裂部位の的中 | CFD は6箇所に分散／外科医 0/43 | 部位予測に合意は成立しなかった |
| メッシュ要素数 | 28.3万〜1797万 | 中央値 210万。依存性を検証したのは5グループのみ |
| 時間刻み | 1×10⁻⁵〜5×10⁻³ 秒 | 1心周期あたり 200〜100,000 ステップ |
| 設定した流量 | 0.9〜11.36 mL/s | 実測は 1.76 mL/s と 2.5 mL/s。中央値 3.5 mL/s |
| 入口流速 | 0.25〜1.51 m/s | 外れ値として 10 m/s を用いたグループがある |
| Reynolds 数 | 134〜815 | 非現実的な Re=1 の報告もあった |
| 血液粘度（ニュートン） | 3.5〜4.01 mPa·s | 18グループ。中央値 4 mPa·s |
| 破裂部位の TAWSS | 0.1〜19 dyn/cm² | 母血管は 50 dyn/cm² 超。OSI は最大 0.42 |
いずれもこの論文が報告した値です。
実務チェック
- 文献値と自分の結果を比べる前に、相手の流量・入口プロファイル・粘性モデル・出口条件を確認する。書かれていなければ比較は成立しない。
- 自分の解析では流量と入口条件の根拠を明記する。実測がないなら、その旨を書く。
- メッシュ依存性を自分で検証する。このチャレンジでは26グループ中5グループしか実施していない。
- 破裂部位や危険部位の特定に単一の指標を使わない。低 WSS だけを根拠にした10グループは部位を外している。
この研究の限界
- 対象は動脈瘤2症例のみ。著者自身が代表性は不十分としています。
- 剛体壁を前提としており、壁厚の分布や壁運動は考慮されていません。
- 破裂部位は自然破裂時ではなく、術中の再破裂を観察したものです。
書誌情報
| 原題 | The Computational Fluid Dynamics Rupture Challenge 2013—Phase I: prediction of rupture status in intracranial aneurysms |
|---|---|
| 著者 | Janiga G, Berg P, Sugiyama S, Kono K, Steinman DA |
| 掲載誌 | AJNR American Journal of Neuroradiology（2015） |
| DOI | 10.3174/ajnr.A4157 |
| PMID | 25500315 |