マイクロ流体力学
理論と物理
マイクロ流体力学の基礎
先生、マイクロ流体力学ってマクロの流体力学と何が違うんですか？ スケールが小さいだけ？
スケールの変化は物理の「支配者」を変える。マクロスケールでは慣性力が支配的だが、マイクロスケール（代表長さ $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によるチャンネル設計最適化が製品開発の核心となっている。
各項の物理的意味
- 時間項 $\partial(\rho\phi)/\partial t$：蛇口をひねった瞬間を思い浮かべてください。最初は水がバタバタと不安定に出て、しばらくすると安定した流れになりますよね？ この「変化している最中」を記述するのが時間項です。心臓の拍動で血流が脈打つのも、エンジンのバルブが開閉するたびに流れが変動するのも、すべて非定常現象。では定常解析とは？ 「十分時間が経って流れが落ち着いた後」だけを見る——つまりこの項をゼロにする。計算コストが大幅に下がるため、まず定常で解いてみるのがCFDの基本戦略です。
- 対流項 $\nabla \cdot (\rho \mathbf{u} \phi)$：川に落ち葉を落としたらどうなりますか？ 流れに乗って下流に運ばれますよね。これが「対流」——流体の動きが物を運ぶ効果です。暖房の温風が部屋の端まで届くのも、空気という「運び屋」が熱を対流で輸送しているから。ここが面白いところ——この項は「速度×速度」を含むため非線形です。つまり、流れが速くなるとこの項が急激に強くなり、制御が難しくなる。これが乱流の根本原因です。よくある勘違い：「対流と伝導は同じようなもの」→ 全然違います！ 対流は流れが運ぶ、伝導は分子が伝える。桁違いの効率差があります。
- 拡散項 $\nabla \cdot (\Gamma \nabla \phi)$：コーヒーにミルクを入れて放置したことはありますか？ かき混ぜなくても、しばらく経つと自然に混ざりますよね。あれが分子拡散です。では次の質問——ハチミツとお水、どちらが流しやすいですか？ 当然お水ですよね。ハチミツは粘性（$\mu$）が高いから流れにくい。粘性が大きいと拡散項が強くなり、流体は「もったりした」動きになります。レイノルズ数が小さい流れ（ゆっくり、ドロドロ）では拡散が支配的。逆にRe数が大きい流れでは対流が圧倒し、拡散は脇役になります。
- 圧力項 $-\nabla p$：注射器のピストンを押すと、液体が針先から勢いよく出ますよね？ なぜでしょう？ ピストン側が高圧、針先が低圧——この圧力差が流体を押す力になるからです。ダムの放水も同じ原理。天気図で等圧線がギュッと密になっている場所では？ そう、強風が吹きます。「圧力差があるところに流れが生まれる」——これがナビエ-ストークス方程式の圧力項の物理的意味。ここでの勘違いポイント：CFDの「圧力」は絶対圧ではなくゲージ圧のことが多い。圧縮性解析に切り替えたとたんに結果がおかしくなる場合、絶対圧/ゲージ圧の混同が原因かもしれません。
- ソース項 $S_\phi$：暖められた空気が上に昇る——なぜでしょう？ 周囲より軽く（密度が低く）なったから、浮力で押し上げられるのです。この浮力はソース項として方程式に追加されます。他にも、ガスコンロの炎で化学反応熱が発生する、工場の電磁ポンプで金属溶湯にローレンツ力がかかる…これらはすべて「外部から流体にエネルギーや力を注入する」作用であり、ソース項で表現します。ソース項を忘れるとどうなるか？ 自然対流の解析で浮力を入れ忘れると、流体は一切動かない——冬の部屋で暖房をつけたのに暖かい空気が上に行かない、という物理的にありえない結果になります。
仮定条件と適用限界
- 連続体仮定：クヌッセン数 Kn < 0.01（分子平均自由行程 ≪ 代表長さ）で成立
- ニュートン流体仮定：せん断応力と歪み速度が線形関係（非ニュートン流体では粘度モデルが必要）
- 非圧縮性仮定（Ma < 0.3の場合）：密度を一定として扱う。マッハ数0.3以上では圧縮性効果を考慮
- ブシネスク近似（自然対流）：密度変化を浮力項のみで考慮し、他の項では一定密度を使用
- 適用外ケース：希薄気体（Kn > 0.1）、超音速・極超音速流れ（衝撃波捕捉が必要）、自由表面流れ（VOF/Level Set等が必要）
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 速度 $u$ | m/s | 入口条件で体積流量から換算する際、断面積の単位に注意 |
| 圧力 $p$ | Pa | ゲージ圧と絶対圧の区別。圧縮性解析では絶対圧を使用 |
| 密度 $\rho$ | kg/m³ | 空気: 約1.225 kg/m³@20°C、水: 約998 kg/m³@20°C |
| 粘性係数 $\mu$ | Pa·s | 動粘性係数 $\nu = \mu/\rho$ [m²/s] との混同に注意 |
| レイノルズ数 $Re$ | 無次元 | $Re = \rho u L / \mu$。層流/乱流遷移の判定指標 |
| CFL数 | 無次元 | $CFL = u \Delta t / \Delta x$。時間刻みの安定性に直結 |
数値解法と実装
マイクロ流体の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法」や「格子ボルツマン法」を組み合わせて気液界面を追跡する必要があり、手法選択が結果の精度を大きく左右します。
風上差分（Upwind）
1次風上: 数値拡散が大きいが安定。2次風上: 精度向上するが振動のリスク。高レイノルズ数流れでは必須。
中心差分（Central Differencing）
2次精度だが、Pe数 > 2で数値振動が発生。低レイノルズ数の拡散支配流れに適する。
TVDスキーム（MUSCL、QUICK等）
リミッタ関数により数値振動を抑制しつつ高精度を維持。衝撃波や急勾配の捕捉に有効。
有限体積法 vs 有限要素法
FVM: 保存則を自然に満足。CFDの主流。FEM: 複雑形状・マルチフィジックスに有利。SPH等のメッシュフリー法も発展中。
CFL条件（クーラン数）
陽解法: CFL ≤ 1が安定条件。陰解法: CFL > 1でも安定だが、精度と反復回数に影響。LES: CFL ≈ 1を推奨。物理的意味: 1タイムステップで情報が1セル以上進まないこと。
残差モニタリング
連続の式・運動量・エネルギーの各残差が3〜4桁低下で収束と判断。質量保存の残差は特に重要。
緩和係数
圧力: 0.2〜0.3、速度: 0.5〜0.7が一般的な初期値。発散する場合は緩和係数を下げる。収束後は上げて加速。
非定常計算の内部反復
各タイムステップ内で定常解に収束するまで反復。内部反復数: 5〜20回が目安。残差がタイムステップ間で変動する場合は時間刻みを見直す。
SIMPLE法のたとえ
SIMPLE法は「交互に調整する」手法。まず速度を仮に求め（予測ステップ）、その速度で質量保存が満たされるよう圧力を補正し（補正ステップ）、補正された圧力で速度を修正する——このキャッチボールを繰り返して正解に近づく。2人で棚を水平にする作業に似ている：片方が高さを合わせ、もう片方がバランスを取り、これを交互に繰り返す。
風上差分のたとえ
風上差分は「川の流れに立って上流の情報を重視する」手法。川の中にいる人が下流を見ても水の出所は分からない——上流の情報が下流を決めるという物理を反映した離散化手法。精度は1次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
マイクロ流体デバイスの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解析でこれを事前に予測するには、動的接触角モデルを含む多相流解析が必要です。「実験が上手くいかない」と悩む前に、表面ぬれ性の数値パラメータを確認するのが実践の鉄則です。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ 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倍に出てしまう。専用の薄膜流れオプションや分子動力学ハイブリッドの有無が、マイクロ流体専用ツールと汎用ツールの大きな差です。
選定で最も重要な3つの問い
- 「何を解くか」：マイクロ流体力学に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
ラティスボルツマン法（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）を入口前に設けるだけで改善することが多いです。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——マイクロ流体力学の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告