2025-11-19T11:46:14.009582

Mixing in viscoelastic fluids using elastic turbulence

van Buel, Stark
We investigate the influence of elastic turbulence on mixing of a scalar concentration field within a viscoelastic fluid in a two-dimensional Taylor-Couette geometry using numerical solutions of the Oldroyd-B model. The flow state is determined through the secondary-flow order parameter indicating the regime of elastic turbulence. When starting in the turbulent state and subsequently lowering the Weissenberg number, a weakly-chaotic flow occurs below $\mathrm{Wi_c}$. Advection in both the turbulent and weakly-chaotic flow states induces mixing, which we illustrate by the time evolution of the standard deviation of the solute concentration from the uniform distribution. In particular, in the elastic turbulent state mixing is strong and we quantify it by the mixing rate, the mixing time, and the mixing efficiency. All three quantities follow scaling laws. Importantly, we show that the order parameter is strongly correlated to the mixing rate and hence is also a good indication of mixing within the fluid.
academic

粘弾性流体における弾性乱流を用いた混合

基本情報

  • 論文ID: 2409.06391
  • タイトル: Mixing in viscoelastic fluids using elastic turbulence
  • 著者: Reinier van Buel, Holger Stark (ベルリン工科大学)
  • 分類: physics.flu-dyn
  • 発表時期: 2024年9月
  • 論文リンク: https://arxiv.org/abs/2409.06391

要約

本論文は、Oldroyd-Bモデルの数値解を用いて、二次元Taylor-Couette幾何における粘弾性流体内のスカラー濃度場の混合に対する弾性乱流の影響を研究している。流動状態は二次流秩序パラメータにより決定され、弾性乱流領域を示す。乱流状態から開始してWeissenberg数を低下させると、Wic以下で弱カオス流が出現する。乱流および弱カオス流状態における移流は両者とも混合を誘起し、均一分布に対する溶質濃度の標準偏差の時間発展により示される。特に、弾性乱流状態では混合が強く、混合速度、混合時間、および混合効率により定量化される。これら3つの量はすべてスケーリング則に従う。重要なことに、秩序パラメータは混合速度と強い相関を示し、したがって流体内混合の良好な指標でもある。

研究背景と動機

問題背景

  1. 低Reynolds数における混合の課題: ニュートン流体では、低Reynolds数条件下の流動は層流であり、混合過程は遅く、主に拡散により支配される。これはマイクロフルイディクス装置における重要な制限である。
  2. 弾性乱流の発見: 流体に高分子量ポリマーを添加すると、流体は粘弾性となり、慣性が無視できる場合でも弾性乱流現象を示す。
  3. 微小スケール混合の重要性: ラボオンチップ(lab-on-a-chip)装置では、特にマイクロメートルスケールにおいて、効率的な流体混合が極めて重要である。

研究動機

  • 混合効率向上における弾性乱流の可能性を探索する
  • 粘弾性流体における混合機構の物理的本質を理解する
  • マイクロフルイディクス応用に対する理論的指針を提供する

主要な貢献

  1. 転移性質の決定: 層流から弾性乱流への転移を亜臨界転移として識別し、弱カオス流状態の存在を発見した
  2. 混合表現方法の確立: 混合速度、混合時間、および混合効率の3つの定量化指標を導入し、すべてスケーリング則に従うことを示した
  3. 秩序パラメータと混合の強い相関性の発見: 二次流秩序パラメータが混合速度と強く相関し、流体混合程度の良好な指標として機能することを証明した
  4. Batchelor混合機構の検証: 弾性乱流状態下で指数減衰を観察し、Batchelor混合領域に対応し、実験結果と一致することを確認した

方法論の詳細

幾何学的設定と支配方程式

研究は二次元Taylor-Couette幾何を採用し、2つの同心円筒を含む：

  • 外筒半径: ro=10μmr_o = 10\mu m
  • 内筒半径: ri=2.5μmr_i = 2.5\mu m
  • 外筒は角速度Ω=2πs−1\Omega = 2\pi s^{-1}で反時計方向に回転

流体運動は一般化Navier-Stokes方程式により記述される： ρ(∂u∂t+u⋅∇u)=−∇p+ηs∇2u+∇⋅τ\rho\left(\frac{\partial \mathbf{u}}{\partial t} + \mathbf{u} \cdot \nabla \mathbf{u}\right) = -\nabla p + \eta_s \nabla^2 \mathbf{u} + \nabla \cdot \boldsymbol{\tau}

Oldroyd-B構成モデル

ポリマー応力テンソルτ\boldsymbol{\tau}はOldroyd-Bモデルを採用する： τ+λτ∇=ηp[∇⊗u+(∇⊗u)T]\boldsymbol{\tau} + \lambda \stackrel{\nabla}{\boldsymbol{\tau}} = \eta_p[\nabla \otimes \mathbf{u} + (\nabla \otimes \mathbf{u})^T]

ここでτ∇\stackrel{\nabla}{\boldsymbol{\tau}}は上対流導数である。

混合のモデル化

スカラー濃度場c(r,t)c(\mathbf{r}, t)の進化は移流-拡散方程式により支配される： ∂c∂t+(u⋅∇)c=D∇2c\frac{\partial c}{\partial t} + (\mathbf{u} \cdot \nabla)c = D\nabla^2 c

本研究ではD=0D = 0と設定し、分子拡散を無視する。

主要な無次元パラメータ

  • Reynolds数: Re=ρΩro2/ηs=0.628×10−3Re = \rho\Omega r_o^2/\eta_s = 0.628 \times 10^{-3}
  • Weissenberg数: Wi=λΩWi = \lambda\Omega
  • 粘度比: β=ηp/ηs=1.5\beta = \eta_p/\eta_s = 1.5

実験設定

数値手法

  • OpenFOAM®オープンソース有限体積ソルバーを使用
  • 粘弾性流動を扱うためRheoTool専用ソルバーを採用
  • メッシュ：径方向100要素、周方向120要素のスポークホイールメッシュ
  • 時間ステップ: δt=10−5s\delta t = 10^{-5}s

秩序パラメータの定義

二次流強度を秩序パラメータとして定義する： Φ=σ=⟨[u(r,t)−u0(r)]2⟩/u0max⁡\Phi = \sigma = \sqrt{\langle[\mathbf{u}(\mathbf{r},t) - \mathbf{u}_0(\mathbf{r})]^2\rangle}/u_0^{\max}

混合表現指標

正規化標準偏差： Δc=⟨[c(r,t)−c∞]2⟩⟨[c(r,0)−c∞]2⟩\Delta c = \sqrt{\frac{\langle[c(\mathbf{r},t) - c_\infty]^2\rangle}{\langle[c(\mathbf{r},0) - c_\infty]^2\rangle}}

ここでc∞c_\inftyは完全混合状態の濃度である。

実験結果

流動状態の転移

  1. 亜臨界転移: 層流から弾性乱流への転移が亜臨界性を示し、ヒステリシスが存在することを発見
  2. 臨界Weissenberg数:
    • 上昇枝: Wic=9.95Wi_c = 9.95
    • 下降枝: WicI=9.26Wi_c^I = 9.26
  3. 弱カオス状態: WicIWi_c^I以下で弱カオス流状態を観察

混合特性の分析

弱カオス流状態 (5.03<Wi≤8.175.03 < Wi \leq 8.17)

  • 標準偏差は指数がおよそ−1/6-1/6のべき乗則で減衰
  • パワースペクトルは指数が−5-5以上のべき乗則を示す

弾性乱流状態 (Wi≥9.42Wi \geq 9.42)

  • 初期段階ではべき乗則減衰(指数約1/51/5)、その後指数減衰に転移
  • 指数減衰形式: Δc∼exp⁡(−αmixt)\Delta c \sim \exp(-\alpha_{mix}t)
  • パワースペクトルはBatchelor スペクトル特性(k−1k^{-1}スケーリング)を示す

スケーリング則関係

  1. 混合時間: tmix∼(Wi−WicI)−0.52t_{mix} \sim (Wi - Wi_c^I)^{-0.52}
  2. 領域II濃度: cII/c∞∼(Wi−WicI)0.28c_{II}/c_\infty \sim (Wi - Wi_c^I)^{0.28}
  3. スケーリング関係: cII∼tmix−1/2c_{II} \sim t_{mix}^{-1/2}

秩序パラメータと混合速度の相関性

実験結果は秩序パラメータΦd\Phi_dが混合速度αmix\alpha_{mix}と強く相関することを示し、速度変動の増加が流体内の混合効率を直接向上させることを示唆している。

関連研究

弾性乱流研究

  • Groisman & Steinbergによるポリマー溶液における弾性乱流現象の初発見
  • Taylor-Couette幾何における線形および非線形安定性解析
  • 三次元数値シミュレーションで観察されたらせん流および帯状流パターン

混合増強研究

  • 四ロール混合機における伸張流混合
  • 十字スロット幾何における弾性不安定性
  • 蛇行チャネルにおける対流熱伝達増強

結論と考察

主要な結論

  1. 転移性質: Taylor-Couette流における弾性乱流転移の亜臨界性質を確認した
  2. 混合機構: 弾性乱流が極めて低いReynolds数下で効率的な混合を実現できることを証明した
  3. Batchelor領域: 弾性乱流状態下でBatchelor混合特性を観察した
  4. 予測指標: 秩序パラメータが混合効率の信頼できる予測指標として機能することを示した

制限事項

  1. 二次元制限: 研究は二次元幾何に限定され、三次元効果は異なる結果をもたらす可能性がある
  2. メッシュ依存性: ヒステリシス挙動はより細かいメッシュでは消失し、数値精度の影響を示唆している
  3. 初期条件感度: 混合特性は初期濃度分布の選択に依存する

今後の方向性

  1. 弾性波の影響: 弾性波が混合効率、混合速度、および混合時間に及ぼす影響を研究する
  2. 三次元拡張: 三次元幾何で発見されたスケーリング則を検証する
  3. 実用応用: マイクロフルイディクス装置における実用的応用の可能性を探索する

深い評価

利点

  1. 体系的研究: 層流から弾性乱流への転移過程における混合特性を包括的に分析した
  2. 定量表現: 混合程度と効率を表現するための複数の定量指標を導入した
  3. 理論検証: 数値結果はBatchelor混合理論および実験観察と一致している
  4. 実用価値: マイクロフルイディクス応用に対する重要な理論的指針を提供した

不足点

  1. 数値精度問題: ヒステリシス現象のメッシュ依存性は数値手法に潜在的な制限があることを示唆している
  2. 機構理解: 弱カオス状態下のべき乗則減衰の物理機構に対する説明が不十分である
  3. パラメータ範囲: 研究のパラメータ範囲は相対的に限定的であり、より広いパラメータ空間の探索が不足している

影響力

  1. 理論的貢献: 粘弾性流体における混合機構の理解に重要な洞察を提供した
  2. 応用前景: マイクロフルイディクス装置設計に対する重要な指導的意義を有する
  3. 方法論的価値: 確立された表現方法は他の流動幾何に一般化可能である

適用シーン

  • マイクロフルイディクス装置における混合最適化
  • ポリマー加工プロセスにおける混合制御
  • 生物医学応用における微小スケール混合

参考文献

論文は弾性乱流、粘弾性流動安定性、混合理論など複数の分野をカバーする57篇の関連文献を引用しており、研究に対する堅実な理論的基礎を提供している。