Compliance-Induced Transition in Soft Microtubes: Near-Wall Measurements and Empirical Scaling

技術要約：柔軟なマイクロチューブにおけるコンプライアンス誘起遷移

問題提起
剛体円管における層流から乱流への遷移は古典的な問題（通常、$Re \approx 2100$で発生）であるが、柔軟な導管における流れの安定性は、特にマイクロスケールにおいて十分に理解されていない。壁面のコンプライアンス（柔軟性）は、流体構造相互作用（FSI）を介して、剛性管の閾値を大幅に下回るレイノルズ数で流れの不安定性を引き起こす可能性がある。しかし、現在の文献には以下の3つの決定的な欠落が存在する：

1. 空間解像度の欠如： 従来の研究は、不安定性の発生を推論するために、グローバルな指標（摩擦係数や圧力降下など）に依存していた。流体・弾性不安定性が理論的に発生すると予測されている近壁領域（$r/R \gtrsim 0.7$）において、速度変動の径方向の組織化を解明した研究は存在しない。
2. 大変形におけるスケーリングの未検証： 従来のスケーリング則（$Re_t \propto \Sigma^{3/4}$）は、弱い変形（数パーセント未満）を対象として導出されたものである。流動誘起の変形が$\sim 20\%$に達する場合があるマイクロスケール領域において、これらの漸近的なスケーリングが保持されるのかは不明である。
3. 系統的な特性評価のバイアス： 従来の研究では、硬化後の剪断弾性率の値や公称のテンプレート直径が用いられることが多く、溶媒による軟化や、テンプレートを用いた作製過程における事前膨張が無視されている。これは、計算されるレイノルズ数（$Re$）およびコンプライアンスパラメータ（$\Sigma$）にコンプライアンス依存のバイアスをもたらす。

手法
著者らは、PDMSマイクロチューブ（$D \approx 670\, \mu\text{m}$）を用い、壁面のコンプライアンスを系統的に変化させた（$G' = 14\text{--}392\, \text{kPa}$）、高速粒子トラッキング流速計（HSPT）を用いてこれらの欠落に対処する。

• 作製および特性評価： マイクロチューブは、剛性の高い入口部と変形可能なテスト部からなる二部構成のアプローチを用いて作製された（異なる硬化剤濃度を使用）。極めて重要な点として、本研究では、チューブの作製に使用されたものと同じ溶媒処理（エタノール-トルエン浸漬および制御された脱膨張）を適用した、作製一貫性のある2つの補正を実施している：
  1. 剪断弾性率 ($G'$)： モノマーの溶出とネットワークの軟化を考慮するため、溶媒処理後（エタノール-トルエン浸漬および制御された脱膨張）に測定された。
  2. 静止半径 ($R_0$)： 公称のテンプレート半径を仮定するのではなく、溶媒処理後に光学的に測定された。これにより、コンプライアンスに依存する事前膨張（最も柔らかいチューブで最大$7.8\%$）を考慮している。
• 実験プロトコル： シリンジポンプを用いてマイクロチューブ内への流れを駆動した。トレーサー粒子（$\sim 5\, \mu\text{m}$ ミカ）を10,000 FPSで追跡した。
• 遷移診断： 遷移レイノルズ数（$Re_t$）は、同一のデータセットに対して適用される3つの独立した基準を用いて同時に特定された：
  1. 摩擦係数（$f$）の層流における$64/Re$スケーリングからの逸脱。
  2. 近壁の速度勾配（壁面剪断応力 $\tau_w$）のハーゲン・ポアズイユ予測からの逸脱。
  3. 軸方向速度プロファイル（$V_z$）の局所的な放物線形状からの逸脱。
• 変動解析： 壁付近の空間解像度の限界 due to 空間解像度の制限により、本研究では、位相分解されたレイノルズ剪断応力の代わりに、結合された変動強度の指標として、速度変動の積 $v'_{r,\text{rms}} v'_{z,\text{rms}}$ を利用している。

主な結果

1. コンプライアンス依存の遷移： 遷移レイノルズ数は、壁面のコンプライアンスが増加するにつれて系統的に減少した。$Re_t$は、$392\, \text{kPa}$（剛体）の$\approx 2100$から、$76\, \text{kPa}$で$\approx 1200$、$35\, \text{kPa}$で$\approx 900$、$14\, \text{kPa}$で$\approx 600$へと低下した。
2. 近壁の局在化： コンプライアンスのあるチューブにおいて、結合された速度変動の顕著な近壁領域（$r/R \approx 0.7\text{--}0.9$）への局在化が観察された。これは剛体管には見られない特徴である。この局在化の強度は壁の柔らかさに比例して増加したが、グローバルな変動の大きさは減少した（遷移がより低い慣性強制で発生するため）。この空間的シグネチャーは、壁モード不安定性のメカニズムと一致している。
3. スケーリング指数： チューブレベルのブートストラップ分析（$N=4$ のコンプライアンスレベル）により、$Re_t \propto \Sigma^n$ の関係に対するスケーリング指数 $n = 0.48$（95% BCa CI $[0.37, 0.66]$）が得られた。
   • この結果は、理論的な壁モード予測である $\Sigma^{3/4}$ とは統計的に矛盾している。
   • また、$\Sigma^{1/2}$（非粘性モード）および $\Sigma^{5/8}$ とは一致している。
   • データには、$\Sigma^{1/2}$ と $\Sigma^{5/8}$ を決定的に区別するための統計的精度が不足している。
4. 文献の不均一性： 既存の実験文献に関する研究レベルのメタ解析により、報告されているスケーリング指数（$n = 0.41\text{--}1.24$）に大幅な不均一性（$I^2 = 96\%$）があることが明らかになった。これは、単一の普遍的な冪乗則が、すべてのコンプライアンスを持つ導管構成を記述するわけではないことを示している。

意義および主張
本論文は、コンプライアンス誘起遷移の発生時における、速度変動の最初の空間分解測定を提供したと主張している。材料および幾何学的特性の作製誘起バイアスを補正することにより、本研究は、見かけのスケーリング指数がこれらの方法論的な詳細に敏感であることを示している。

著者らは以下のように結論付けている：

• コンプライアンスのあるマイクロチューブにおける遷移は、変動の近壁局在化によって裏付けられるように、弾性・流体結合によって支配されている。
• すべての実験データベースにわたって普遍的な $Re_t \propto \Sigma^n$ スケーリングを想定することは支持されない。観察された不均一性は、形状、作製プロトコル、および流体の性質の差異に起因すると考えられる。
• 漸近的な $\Sigma^{3/4}$ スケーリングは現在のデータセットでは支持されないが、$\Sigma^{1/2}$ および $\Sigma^{5/8}$ は依然として有力な候補である。
• コンプライアンスを有するマイクロ流体デバイスの正確な安定性予測には、公的な（as-cured）値や公称値ではなく、作製後の溶媒処理された状態での材料剛性と幾何学を特性評価することが不可欠である。

本研究は、摂動固有関数または位相分解されたレイノルズ応力の直接的な測定なしには、特定の不安定性モード（例：変形修飾された壁モードと非粘性モードの区別）を決定的に特定したとは主張していない。これらは今後の課題である。