マイクロ流体力学
マイクロ流体力学の理論基礎
マイクロ流体力学の基礎
先生、マイクロ流体力学ってマクロの流体力学と何が違うんですか？ スケールが小さいだけ？
スケールの変化は物理の「支配者」を変える。マクロスケールでは慣性力が支配的だが、マイクロスケール（代表長さ $L \sim 1\text{--}100\,\mu\text{m}$）では粘性力と表面張力が支配的になる。
具体的に数値で比較しよう。水流速 $U = 1\,\text{mm/s}$、チャネル幅 $D = 100\,\mu\text{m}$ の場合:
Re << 1 なので慣性力は無視でき、Stokes流れ（クリーピングフロー）が支配する。
Re = 0.1 だと完全に粘性支配ですね。乱流の心配はゼロ。
そうだ。さらに表面張力の重要度を示す毛管数（Capillary number）は:
Ca << 1 なので、液滴・気泡の形状は表面張力が決定する。
Stokes方程式
Re << 1 の極限では、Navier-Stokes方程式の慣性項が消えてStokes方程式になる。
Stokes方程式には以下の重要な性質がある:
- 線形: 解の重ね合わせが可能
- 時間可逆: 外力を反転すれば流体は元の状態に戻る（混合が困難な理由）
- 瞬時応答: 慣性がないので圧力変化が即座に全体に伝搬
時間可逆ということは、マイクロスケールでは撹拌しても混ざらないんですか？
通常の方法では混ざらない。これがマイクロミキサーの設計が難しい理由だ。拡散に頼るか、カオス的対流（流路のジグザグ構造など）を利用するしかない。
マイクロ流体特有の物理
マイクロスケールで重要になる物理現象を整理しよう。
| 現象 | 支配パラメータ | マクロでの影響 | マイクロでの影響 |
|---|---|---|---|
| 表面張力 | Ca, We, Bo | 通常は無視 | 支配的 |
| 電気浸透流 (EOF) | $\zeta$ 電位、Debye長 | 無視 | 駆動力として利用 |
| スリップ流 | Knudsen数 Kn | no-slip | Kn > 0.01 でスリップ |
| 拡散混合 | Peclet数 Pe | 対流支配 | 拡散支配 |
| 接触角ヒステリシス | 前進/後退接触角 | 重要でない | デバイス動作を支配 |
Bond数 $\text{Bo} = \rho g L^2 / \sigma$ もマイクロスケールでは非常に小さくなる（$L = 100\,\mu\text{m}$ で $\text{Bo} \sim 10^{-6}$）。つまり重力は完全に無視でき、宇宙でも地上でも同じ振る舞いをする。
重力が効かないのは直感に反しますが、スケールの効果なんですね。
マイクロ流体工学の誕生——1990年代のμTAS革命とラブオンチップ
マイクロ流体工学（Microfluidics）の黎明期は1990年代初頭、スイスのETHでManz & Widmer(1990)が「μTAS（Micro Total Analysis System）」概念を提唱したことに始まる。半導体製造技術（フォトリソグラフィ）でガラス基板上に幅数十〜数百μmの流路を作製し、試薬の混合・分離・検出をチップ上で完結させる「ラブオンチップ（Lab-on-Chip）」の概念だ。この流路スケールではReynolds数が1未満のStokes流れが支配的で、乱流混合が使えないためT字型やヘリカル型の受動混合構造が必要になる。現代ではCOVID-19迅速診断キットやDNAシーケンサーのマイクロチャンネルにこの技術が活用されており、CFDによるチャンネル設計最適化が製品開発の核心となっている。
マイクロ流体力学の数値計算手法
マイクロ流体のCFD手法
マイクロ流体のシミュレーションには、普通のCFDソルバーが使えるんですか？
連続体近似が成立する範囲（Kn < 0.01）では通常のNavier-StokesベースのCFDが使える。ただし、マイクロスケール特有の数値的課題がある。
二相流（液滴・気泡）の界面追跡
マイクロ流体デバイスでは液滴生成やT字合流が頻出する。界面追跡の主要手法を比較しよう。
| 手法 | 長所 | 短所 | 適用例 |
|---|---|---|---|
| VOF (Volume of Fluid) | 質量保存良好、コスト低 | 界面が拡散、接触線問題 | 液滴生成、スラグ流 |
| Level-Set | 界面形状が滑らか | 質量非保存 | 合流、分裂 |
| CLSVOF | VOF+LSの長所を統合 | 実装が複雑 | 高精度が必要な問題 |
| Phase-Field | 接触線の自然な記述 | メッシュ要求が厳しい | 濡れ現象、接触角制御 |
| Front-Tracking | 界面の明示的追跡 | トポロジー変化に弱い | 単一液滴変形 |
表面張力の数値計算ではCSF (Continuum Surface Force) モデルが標準だが、マイクロスケールでは寄生電流（spurious current）が問題になる。これは表面張力の離散化に起因する非物理的な流れで、Ca が小さいほど顕著になる。
表面張力が大きいのにメッシュが粗いと、非物理的な流れが生じるんですね。
対策としては:
- 界面あたり10セル以上の解像度を確保
- Height Function法による曲率計算（CSFより精度が大幅に向上）
- Sharp Surface Force法
- Phase-Fieldモデルの使用（界面の熱力学的な記述で寄生電流を低減）
電気浸透流（EOF）の計算
マイクロチャネルでは電場で流体を駆動する電気浸透流がよく使われる。支配方程式は:
$\rho_e$ は電荷密度、$\varepsilon$ は誘電率、$\mathbf{E} = -\nabla\phi$ は電場。Debye長 $\lambda_D \sim 1\text{--}100\,\text{nm}$ の領域を解像する必要があり、非常に細かいメッシュが要求される。
Debye長がnmオーダーだと、メッシュのスケールが流路全体と3桁以上違いますね。
そうだ。実務的にはHelmholtz-Smoluchowskiのスリップ速度条件:
を壁面境界条件として適用し、EDL（電気二重層）の内部を解像しないアプローチが多い。$\zeta$ はゼータ電位、$E_t$ は壁面接線方向の電場だ。
マイクロ流路の「毛細管力」——重力より表面張力が支配する世界
チャネル幅が100μm以下になると、重力は流れを支配する力から外れる。代わりに表面張力が主役になり、毛細管圧力 ΔP=4γcosθ/D（Dはチャネル径）が液体を引き込む。実際、疾患診断用の紙マイクロ流体（paper-based microfluidics）では、ポンプなしで血液が毛細管力で流れる。これを数値的に扱うには、通常のNavierStokesに「VOF法」や「格子ボルツマン法」を組み合わせて気液界面を追跡する必要があり、手法選択が結果の精度を大きく左右します。
マイクロ流体力学の実務適用
マイクロ流体デバイスのCFD解析フロー
マイクロ流体デバイスのシミュレーションは、具体的にどんな手順で進めるんですか？
典型的なT字型ドロップレットジェネレータを例に説明しよう。
Step 1: 形状準備
- チャネル幅 100 μm、深さ 50 μm が典型的
- CADでの設計は精密に。面取りや角のR（製造上の丸み）を反映
- PDMSデバイスの製造誤差（±5 μm程度）も考慮
Step 2: メッシュ設計
- 界面が存在する領域にAMR（適応的メッシュ細分化）を活用
- チャネル幅方向に最低20セル、界面近傍はさらに2-3倍細かく
- 接触線近傍は特に高解像度が必要
Step 3: 物性値設定
- 連続相（例: オイル）: $\mu_c = 0.01\,\text{Pa}\cdot\text{s}$, $\rho_c = 800\,\text{kg/m}^3$
- 分散相（例: 水）: $\mu_d = 0.001\,\text{Pa}\cdot\text{s}$, $\rho_d = 1000\,\text{kg/m}^3$
- 界面張力: $\sigma = 0.005\,\text{N/m}$（界面活性剤添加時）
Step 4: 境界条件
- 入口: 流量指定（シリンジポンプに対応）。典型的には $Q = 1\text{--}10\,\mu\text{L/min}$
- 出口: 圧力指定（大気圧）
- 壁面: no-slip + 接触角指定（PDMSの水接触角 ≈ 110°）
接触角の設定がマクロのCFDにはない要素ですね。
接触角はドロップレットのサイズや生成頻度に大きく影響する。静的接触角だけでなく、動的接触角モデル（Cox-Voinov則など）を使うとより正確だ:
$\theta_D$ は動的接触角、$\theta_S$ は静的接触角、$L_H$ はマクロスケール、$L_m$ は分子スケール。
解析結果の評価指標
マイクロ流体デバイスの性能を評価する主要な指標:
| 指標 | 定義 | 用途 |
|---|---|---|
| 液滴径 | $d/w$（チャネル幅で無次元化） | 液滴生成デバイスの設計 |
| 生成頻度 | $f \cdot w / U$（無次元） | スループットの評価 |
| 混合度 | $1 - \sigma_c / \sigma_{c,0}$（濃度標準偏差） | マイクロミキサーの性能 |
| 圧力損失 | $\Delta p / (\mu U L / w^2)$（無次元） | ポンプ要件 |
| 分離効率 | 対象粒子の回収率 | ソーティングデバイス |
液滴径をチャネル幅で無次元化するのがマイクロ流体の流儀なんですね。
そうだ。T字型ドロップレットジェネレータでは、$d/w$ はCa数の関数として整理される。Squeezing regime（$\text{Ca} < 0.01$）では $d/w \propto Q_d / Q_c$（流量比に依存）、Dripping regime（$\text{Ca} > 0.1$）では $d/w \propto \text{Ca}^{-1/3}$ とスケーリングが変わる。
マイクロ流体デバイスの「失敗あるある」——流路が詰まる本当の原因
マイクロ流体実験で「流れが止まった」と報告されるトラブルの多くは、実はバブル（気泡）の詰まりです。1μm幅のチャネルに1μmの気泡が入ると、毛細管圧力がポンプ圧力を上回り流れが完全に止まる。解決策は、チャネル内壁の親水化処理（接触角θを5°以下に）か、流路の出口にバブルトラップを設けること。CFD解析でこれを事前に予測するには、動的接触角モデルを含む多相流解析が必要です。「実験が上手くいかない」と悩む前に、表面ぬれ性の数値パラメータを確認するのが実践の鉄則です。
マイクロ流体力学のソフトウェア比較
マイクロ流体対応ツールの比較
マイクロ流体のシミュレーションに使えるソフトを教えてください。
マイクロ流体は多くの物理が絡むので、ツール選定が重要だ。用途別に整理しよう。
COMSOL Multiphysics（Microfluidicsモジュール）
マイクロ流体分野で最も広く使われているツールの一つだ。
- 強み: EOF（電気浸透流）、DEP（誘電泳動）、AC EOF、拡散反応の連成が容易
- 界面追跡: Level-Set, Phase-Field, Moving Meshの3手法
- 接触角: 静的・動的接触角の設定が直感的
- EDL解像: Poisson-Nernst-Planck方程式の完全解像が可能
- 弱点: 大規模3D計算にはメモリ不足になりやすい
Ansys Fluent
Fluentは汎用性が高く、マイクロ流体にも適用可能。
- VOFモデル: 液滴生成のシミュレーションに実績多数
- UDF: 電場駆動、接触角のカスタマイズが可能
- AMR: 界面近傍の自動メッシュ細分化
- 弱点: EOF専用モジュールがないため、UDFでの実装が必要
OpenFOAM
オープンソースの自由度が活きる分野だ。
- interFoam: VOFベースの二相流ソルバー（液滴生成に多数の論文実績）
- interPhaseChangeFoam: 相変化（蒸発・凝縮）を伴うマイクロ流体
- カスタムソルバー: EOF、DLP、音響場などのマルチフィジックスをC++で自由に実装
- コミュニティ: microfluidics向けのカスタムソルバーがGitHubで多数公開
専門ツール
| ツール | 開発元 | 特徴 |
|---|---|---|
| FLOW-3D | Flow Science | VOF法の精度が高い。表面張力の取り扱いに定評 |
| Gerris/Basilisk | S. Popinet | AMR+VOF。学術界で非常に高い評価。無料 |
| COMSOL | COMSOL AB | マルチフィジックス連成の容易さ |
| Palabos | FlowKit | 格子ボルツマン法。複雑形状のマイクロ流体に適する |
Basiliskは聞いたことがなかったです。
BasiliskはStephane Popinetが開発した適応直交格子（Adaptive Quadtree/Octree）ベースのソルバーで、表面張力問題の精度では世界トップクラスだ。Height Function法による曲率計算と、AMRの組み合わせで寄生電流が非常に小さい。学術論文で頻繁に引用されている。
選定ガイドライン
用途別の推奨ツール:
| 用途 | 推奨ツール | 理由 |
|---|---|---|
| EOF/DEP | COMSOL | 電気化学連成が容易 |
| 液滴生成 | Basilisk, Fluent, interFoam | 表面張力精度 |
| 拡散混合 | COMSOL, Fluent | 多成分輸送 |
| Lab-on-a-chip全般 | COMSOL | GUIの操作性、学習コスト |
| 学術研究（高精度） | Basilisk | AMR + 高精度表面張力 |
| 大規模3D計算 | OpenFOAM, Fluent | 並列計算スケーラビリティ |
COMSOLが入門向きで、精度を突き詰めるならBasiliskやOpenFOAMという使い分けですね。
マイクロ流体デバイス、商用ツールの「落とし穴」
マイクロ流体チップの設計では、汎用CFDツールをそのまま使うと痛い目に遭うことがある。チャネル幅が10μm以下になると、Knudsen数が0.001を超え始め、連続体仮定が怪しくなる。ANSYSやCOMSOLはデフォルトでNavierStokes連続体を解くので、スリップ流れの補正を自分で設定しないと細管では流量が実験の2〜3倍に出てしまう。専用の薄膜流れオプションや分子動力学ハイブリッドの有無が、マイクロ流体専用ツールと汎用ツールの大きな差です。
マイクロ流体力学の先端研究
ラティスボルツマン法（LBM）によるマイクロ流体計算
格子ボルツマン法がマイクロ流体に向いていると聞いたんですが。
LBMはマイクロ流体との相性が非常に良い。理由は3つある。
1. 多相流の自然な記述: Shan-Chenモデルやfree-energyモデルで界面を自動的に追跡。VOFのような界面再構成が不要
2. スリップ流への拡張: Knudsen数が大きい領域でもスリップ境界条件を自然に導入可能
3. 複雑形状への適合: バウンスバック境界条件で任意形状の壁を簡単に表現
基本的なBGK（Bhatnagar-Gross-Krook）モデルの格子ボルツマン方程式:
$f_i$ は粒子分布関数、$\tau$ は緩和時間（粘度に関連: $\nu = c_s^2(\tau - 0.5)\Delta t$）。
粒子ベースの方法なので、ミクロな物理との親和性が高いんですね。
Digital Microfluidics（EWOD）
Electrowetting-on-Dielectric（EWOD）を使ったDigital Microfluidicsは近年急速に発展している分野だ。
電圧をかけることで局所的に接触角を変化させ、液滴を自在に操作する。Lippmann-Youngの式:
$\theta_V$ は電圧 $V$ 下での接触角、$\theta_0$ は初期接触角、$d$ は絶縁膜厚さ。
CFDでのモデリングでは、Phase-Fieldモデルに電場を連成させるアプローチが主流だ。COMSOLのAC/DC+Phase-Fieldモジュールで計算できる。
電圧で液滴を動かせるのは、PCR検査やポイントオブケア診断に使えそうですね。
Organ-on-a-Chip（臓器チップ）
マイクロ流体力学の最も注目される応用が、Organ-on-a-Chipだ。
血管内皮細胞を培養したマイクロチャネルで臓器の微小環境を再現する。CFDの役割は:
- せん断応力の制御: 血管内皮には $\tau_w = 1\text{--}10\,\text{dyne/cm}^2$ の生理的せん断が必要
- 酸素・栄養の輸送: $\text{Pe} = UL/D_m$ の評価。$D_m \sim 10^{-9}\,\text{m}^2/\text{s}$（酸素の水中拡散）
- 薬物濃度勾配の設計: 拡散律速か対流律速かをPe数で判断
生体の物理をマイクロ流体で再現するわけですね。
研究動向
マイクロ流体力学の最前線:
- ナノ流体力学: チャネル径が10 nm以下になると、連続体近似が破綻。MDシミュレーションやDPD（散逸粒子動力学）が必要
- 音響流（Acoustofluidics）: 超音波定在波で粒子・細胞をサイズ選別。BAW/SAWデバイスのCFDモデリング
- インクジェット高精度化: ピコリットルの液滴を制御するためのCFD最適化。サテライト液滴の抑制
- 3Dプリンティング（Two-Photon Polymerization）: マイクロ流体デバイスの直接造形とシミュレーション連携
- Inertial Microfluidics: Re $\sim$ 10-100の中間レジームで慣性揚力を活用した粒子分離
マイクロ流体は医療・バイオとCFDが交差する最前線なんですね。
オルガノイドチップとマイクロ流体——「体の外で体内を再現」する挑戦
マイクロ流体の最前線のひとつが「Organ-on-a-chip」（臓器チップ）です。肝臓や腸、肺の細胞を10〜100μm幅のマイクロチャネルに培養し、生体内に近い流れ環境（Re≈0.01〜0.1）で薬物反応を試験します。通常の細胞培養と違い、せん断応力・拡散・温度分布を精密にコントロールできるため、動物実験の精度を超えるケースも出ています。流体解析的には「極低Re・複雑な生体物質」という非常に難しい問題で、格子ボルツマン法と分子動力学のハイブリッドが注目されています。
マイクロ流体力学のトラブル対応
マイクロ流体計算の典型的トラブル
マイクロ流体のCFDで特に注意すべきトラブルを教えてください。
マクロスケールのCFDとは異なるトラブルが多い。パターン別に見ていこう。
1. 寄生電流（Spurious Current）で液滴が崩壊
症状: 静止しているべき液滴が変形・分裂する。界面近傍に非物理的な渦が発生する。
原因: 表面張力のCSFモデルにおける曲率計算の離散化誤差。Ca << 1 の場合に顕著。
定量的目安: 寄生電流の速度が物理的な流速の10%を超えると結果が信頼できない。
対策:
- 界面あたりのメッシュ解像度を上げる（最低10セル/液滴径）
- Height Function法による曲率計算に切り替え
- Phase-Fieldモデルの使用（Cahn-Hilliard方程式）
- OpenFOAMの
interFoamならinterfaceCompressionの設定を調整
静止液滴のテストケースで寄生電流の大きさを確認してから本計算に進むべきですね。
2. 接触線がピン止めされて液滴が動かない
症状: 壁面上の液滴が駆動力があるのに移動しない
原因: no-slip条件と界面追跡の矛盾。接触線では速度がゼロなのに界面が移動するパラドックス（Huh-Scriven paradox）が数値的に問題を起こす。
対策:
- Navier-slip条件: 壁面に微小スリップ長 $\beta$ を導入。$u_{\text{slip}} = \beta \partial u/\partial n$
- スリップ長は $\beta = 1\text{--}10\,\text{nm}$ が物理的だが、メッシュで解像できない場合は数値的に大きめの値を使う
- Phase-Fieldモデルは接触線の特異性を自然に回避できる
3. 時間刻みが極端に小さくなる
症状: Courant数制限で $\Delta t \sim 10^{-10}\,\text{s}$ のオーダーになり、計算が進まない
原因: 表面張力による毛管波の伝搬速度制限。Brackbill限界:
$\Delta x = 1\,\mu\text{m}$、水の物性で計算すると $\Delta t < 1.2 \times 10^{-8}\,\text{s}$。
対策:
- メッシュを必要以上に細かくしない
- 陰的な表面張力の取り扱い（一部のソルバーでサポート）
- AMRで界面近傍だけ細かくし、遠方は粗くする
表面張力が時間刻みの制約になるのはマイクロスケール特有の問題ですね。
4. 二相流の質量保存が悪い
症状: 液滴の体積が計算中に減少（または増大）していく
原因: Level-Set法は界面を $\phi = 0$ の等値面として追跡するが、再初期化ステップで質量が失われる。
対策:
- VOF法を使用（質量保存が良い）
- CLSVOF（Level-SetとVOFの結合）
- 質量補正ステップの追加
- Phase-Fieldモデル（Cahn-Hilliard: 質量保存型）
デバッグ手順
マイクロ流体計算のデバッグは以下の順序で進める。
1. 単相流テスト: まず単相のPoiseuille流れで圧力損失が $\Delta p = 12\mu Q L/(w h^3)$（矩形断面）と一致するか確認
2. 静止液滴テスト: 球状液滴を静置し、Laplace圧力 $\Delta p = 2\sigma/R$ が正しく再現されるか。寄生電流が十分小さいか確認
3. 接触角テスト: 壁面上の液滴の平衡接触角が設定値と一致するか確認
4. 動的テスト: T字合流での液滴生成など、実験との定量比較
段階的にテストを重ねていくのがマイクロ流体でも基本なんですね。
特にLaplace圧力テストは最低限実施すべきだ。これが合わなければ二相流の結果は全く信用できない。
マイクロ流体解析の「数値発散」——メッシュサイズが犯人の9割
マイクロ流体解析で解が収束しない場合、まずメッシュのアスペクト比を疑ってください。チャネル断面が10μm×100μmの場合、アスペクト比1:10のセルを使うと圧力勾配の数値誤差が膨らみやすい。理想は各方向5セル以上、アスペクト比3以下です。もう一つのトラブル源は入口境界条件の「急激な速度変化」で、発達流れを仮定せずプラグ流（一様流速）をセットすると入口付近で振動が発生し収束しない。発達長さ L≈0.06·Re·D（マイクロ流体だと数十μm）を入口前に設けるだけで改善することが多いです。
関連トピック
なった
詳しく
報告