マイクロ流体力学 — 実践ガイドとベストプラクティス
マイクロ流体デバイスの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}$ とスケーリングが変わる。
レイノルズの実験(1883年)——乱流発見の瞬間
オズボーン・レイノルズは、管内の水にインクを流す実験で「層流から乱流への遷移」を発見しました。流速を上げていくと、インクの線がある瞬間にグチャグチャに乱れる。この劇的な瞬間を、レイノルズは数学的に $Re = \rho uD/\mu$ という無次元数で表現した。100年以上経った今も、CFDエンジニアが最初に確認するのはこのレイノルズ数です。
実務者のための直感的理解
この解析分野のイメージ
CFDって、要は「デジタル風洞」です。自動車メーカーが巨大な風洞実験設備に何億円もかけるところを、PCの中で再現できる。でも1つ注意——風洞実験なら「風を当てれば結果が出る」けど、CFDでは「メッシュの品質」と「乱流モデルの選択」という見えない品質要因がある。ここを手抜きすると、きれいなコンター図が出ても中身はデタラメ…なんてことになりかねません。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め(計算領域)、水の入り口と出口を設計し(境界条件)、ポンプの強さを設定する(流量条件)。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか? 意外と直感的ではありませんか?
初心者が陥りやすい落とし穴
「y+って何ですか?」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない!」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た!」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか(低速)、全開にするか(高速)。でもCFDではもう一つ——「どのくらい暴れた水を出すか」(乱流強度)も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね? CFDでも入口条件のミスは下流全体に波及します。
CFDメッシュの品質管理や乱流モデルの選定に悩む時間を、もっと創造的な設計作業に使えたら。 — Project NovaSolverはそんな実務者の声から生まれました。
次世代CAEプロジェクト:開発者と実務者をつなぐ
Project NovaSolverは、マイクロ流体力学を含む幅広い解析分野において、実務者の知見を最大限に活かせる環境の実現を探求しています。まだ道半ばですが、共に歩んでいただける方を募集しています。
開発パートナー登録 →