Numerical analysis of capillarity-driven thinning rheometry for polydisperse polymer solutions

技術要約：多分散ポリマー溶液における毛細管駆動型細線化レオメトリの数値解析

問題提起
キャピラリー・ブレイクアップ伸長レオメトリ（CaBER）は、ポリマー溶液の伸長粘性を特性評価するための標準的な手法である。Oldroyd-Bモデルによって記述される理想的な単分散系においては、弾性毛細管（EC）領域におけるフィラメントの細線化速度は流体の緩和時間（$\tau$）に反比例し、これにより、測定された減衰速度（$\dot{\epsilon}_{EC}$）から特徴的な緩和時間 $\tau_{EC} = 2/(3\dot{\epsilon}_{EC})$ を定義できる。しかし、実在する多分散ポリマー溶液において、$\tau_{EC}$ の解釈は曖昧である。先行研究では、$\tau_{EC}$ とせん断または線形粘弾性によって測定された緩和時間との間の不一致が指摘されている。具体的には、最近の二分散ポリスチレン混合物に関する実験（Calabrese et al., 2025）では、高分子量（MW）溶液に低分子量種を高濃度で添加しても、測定された $\tau_{EC}$ には無視できる程度の影響しか与えないことが明らかになった。これは、細線化のダイナミクスがアンサンブル平均ではなく、分子量分布の特定のサブセットによって支配されていることを示唆している。分子量分布、流動履歴、および結果としての $\tau_{EC}$ を結びつける包括的な枠組みが欠如していた。

手法
著者らは、多分散性を明示的に組み込んだ、円柱状ポリマーフィラメントの細線化をシミュレートするための数値モデルを開発した。このアプローチは以下の要素を組み合わせている：

1. 力バランス： 慣性、粘性、毛細管力、およびポリマー応力を考慮した、フィラメント半径 $a(t)$ に関する一次元力バランス方程式。
2. 構成則モデリング： マルチモードFENE-PM（有限伸長非線形弾性－Peterlin平均）モデル。ポリマー溶液は、異なる分子量分布（具体的にはPS7およびPS16ポリスチレン試料）を表す複数の「種（species）」の混合物として扱われる。
   • 線形領域： モデルは、線形粘弾性を記述するために、事前平均化された流体力学的相互作用を用いたZimmモデルを利用し、各種の緩和時間スペクトル $\tau_{is}$ を生成する。
   • 非線形領域： モデルは、FENEスプリングによる有限伸長性と、コンフォメーション依存型のドラッグ項を組み込んでいる。このドラッグ項は、鎖がコイルから楕円体へと引き伸ばされる際の流体力学的摩擦の変化を考慮し、最長モードの緩和時間（$\hat{\tau}_{i1}$）を変化させる。
3. シミュレーション： フィラメント半径と各種/各モードのコンフォメーションテンソルの系を同時に解く常微分方程式系を解く。モデルは、ジオクチルフタル酸（DOP）中のPS7とPS16の二分散ブレンドを用いた実験データ（Calabrese et al., 2025）に対して検証されている。

主要な貢献と結果

• 選択的伸長メカニズム： シミュレーションにより、EC領域において伸長歪速度 $\dot{\epsilon}$ が時間とともに減少することが示された。その結果、有効ワイセンベルグ数 $Wi = \dot{\epsilon}\tau > 0.5$（コイル・ストレッチ転移の閾値）を満たすポリマー鎖のみが、引き伸ばされた状態を維持し、弾性応力に有意に寄与する。緩和時間が短い（低分子量の）鎖は、この閾値を下回り、緩和してフィラメントを支えられなくなる。
• 高分子量種の支配： 緩和時間は分子量に比例するため（$\tau \sim M^{1.5}$）、高分子量種は低分子量種よりも長く $Wi > 0.5$を維持する。したがって、測定される$\tau_{EC}$は、「応力を寄与するサブアンサンブル」によって決定され、これは高分子量側の裾の部分に大きく偏っている。これが、高分子量溶液に低分子量ポリマーを添加しても、高分子量成分が総応力を担うのに十分な濃度であれば$\tau_{EC}$ が変化しない理由である。
• コンフォメーション依存型ドラッグ： 高濃度での実験データに一致させるためには、コンフォメーション依存型ドラッグの導入が必要であった。この項がない場合、モデルは $\tau_{EC}$ を低すぎると予測し、鎖が引き伸ばされ相互作用する際に観察される緩和時間の増加を捉えることができなかった。
• 実験パラメータへの依存性： 本研究は、$\tau_{EC}$ が流体の固有の性質ではなく、実験条件に依存することを裏付けている。具体的には、初期フィラメント径（$a_0$）および前伸長履歴（プレート分離速度）が、EC領域が確立される前に短モードが緩和する時間に影響を与える。初期径が大きいほど、また前伸長速度が高いほど、測定される $\tau_{EC}$ は長くなる。
• 平均の再評価： 本論文は、「伸長性加重分子量（$M_L$）」と「応力加重分子量」を対比させている。$M_L$ は高分子量鎖が応力に寄与する可能性を捉えるものであるが、シミュレーションにおける最小歪速度での応力加重平均の方が、EC領域中に実際に活性化している鎖のサブセットをより正確に反映している。

意義と主張
本論文は、多分散流体のCaBER実験から抽出される特徴的な細線化タイムスケール $\tau_{EC}$ は、流体の単一の固有の緩和時間として解釈されるべきではないと主張している。むしろ、それは、エラストキャピラリーのバランスが達成される瞬間に、流れによって能動的に引き伸ばされている特定の鎖のサブセットのリラクゼーションタイムを反映した、実験固有の量である。

著者らは、多分散溶液において $\tau_{EC}$ は以下の関数の関数であると結論付けている：

1. 分子量分布（特に高分子量側の裾）。
2. 全ポリマー濃度。
3. 初期フィラメント径や前伸長履歴などの実験パラメータ。

この知見は、複雑な流体におけるせん断/線形粘弾性緩和時間と、伸長細線化時間の間の不一致を解消するものである。これは、ポリマーフィラメントの「紡糸性（spinnability）」や安定性が、高分子量分画によって不釣り合いに影響を受けることを示唆しており、工業的な繊維紡糸から生物学的な液滴分裂に至るまでの現象を理解するための重要な洞察を提供している。本研究は、単一モードのOldroyd-B解釈の限界を超え、多分散流体のCaBERデータを分析するための計算フレームワークを提供するものである。