2025-11-22T10:49:15.885853

Transient rod-climbing in an Oldroyd-B fluid

Ruangkriengsin, Brandão, Wu et al.
The Weissenberg effect, or rod-climbing phenomenon, occurs in non-Newtonian fluids where the fluid interface ascends along a rotating rod. Despite its prominence, theoretical insights into this phenomenon remain limited. In earlier work, Joseph \& Fosdick (\emph{Arch. Rat. Mech. Anal.}, vol. 49, 1973, pp. 321--380) employed domain perturbation methods for second-order fluids to determine the equilibrium interface height by expanding solutions based on the rotation speed. In this work, we investigate the time-dependent interface height through asymptotic analysis with dimensionless variables and equations using the Oldroyd-B model. We begin by neglecting surface tension and inertia to focus on the interaction between gravity and viscoelasticity. In the small-deformation scenario, the governing equations indicate the presence of a boundary layer in time, where the interface rises rapidly over a short time scale before gradually approaching a steady state. By employing a stretched time variable, we derive the transient velocity field and corresponding interface profile on this short time scale and recover the steady-state profile on a longer time scale. Subsequently, we reintroduce small but finite inertial effects to investigate their interplay with viscoelasticity and propose a criterion for determining the conditions under which rod-climbing occurs.
academic

Oldroyd-B流体における過渡的ロッド登攀現象

基本情報

  • 論文ID: 2501.01016
  • タイトル: Transient rod-climbing in an Oldroyd-B fluid
  • 著者: Tachin Ruangkriengsin、Rodolfo Brandão、Katie Wu、Jonghyun Hwang、Evgeniy Boyko、Howard A. Stone
  • 分類: physics.flu-dyn（流体力学）
  • 投稿日時: 2025年1月2日（arXivへの投稿）
  • ジャーナル: Journal of Fluid Mechanics（投稿中）
  • 論文リンク: https://arxiv.org/abs/2501.01016

要旨

本論文は、非ニュートン流体におけるワイセンベルク効果（ロッド登攀現象）、すなわち流体界面が回転するロッドに沿って上昇する現象を研究している。著者らはOldroyd-Bモデルを採用し、漸近解析と無次元変数を通じて、時間依存する界面高さを調査した。小変形の場合、支配方程式は時間境界層の存在を示唆しており、界面は短時間スケールで急速に上昇した後、徐々に定常状態に接近することが明らかになった。拡張時間変数を導入することで、著者らは過渡速度場と対応する界面輪郭を導出し、小さいが有限の慣性効果を再導入して粘弾性との相互作用を調査した。

研究背景と動機

問題の定義

ロッド登攀現象（ワイセンベルク効果）は、非ニュートン流体の最も特徴的な挙動の一つであり、ロッドが粘弾性流体中で回転するとき、流体表面はロッドに沿って上昇する。この現象は理論的にはまだ十分に理解されていない。

重要性

  1. 理論的意義: この現象は複雑流体の挙動を理解するための基礎概念であり、教育と研究において重要な位置を占めている
  2. 実用的価値: 流動曲線測定において重要な応用があり、法線応力係数と緩和時間などの材料パラメータの決定に使用できる
  3. 工学的応用: ポリマー溶液および他の複雑流体の産業応用に対する指針となる

既存方法の限界

  1. 理論的限界: 従来の研究は主に定常挙動に焦点を当て、二次流体モデルを使用していた
  2. モデルの問題: 二次流体モデルは小さなデボラ数またはワイセンベルク数でのみ有効であり、高度に不安定である
  3. 時間依存性: 過渡ダイナミクスの理論的分析が不足している

研究動機

著者らは、二次流体モデルの非定常流動における問題を回避するためにOldroyd-Bモデルを選択し、ロッド登攀現象の時間依存性を初めて系統的に研究した。

核心的貢献

  1. 初の理論分析: Oldroyd-B流体におけるロッド登攀現象の初の完全な過渡理論分析を提供
  2. 時間境界層理論: 時間境界層構造を発見・分析し、急速上昇と緩やかな定常状態への接近という二段階プロセスを解明
  3. 無次元分析フレームワーク: ワイセンベルク数に基づく系統的な無次元分析方法を確立
  4. 慣性効果の研究: 小さいが有限の慣性効果と粘弾性の相互作用を初めて研究
  5. 登攀判定基準: ロッド登攀が発生する条件を決定する理論的判定基準を提案：Re<4βpWiRe < 4\beta_p Wi

方法の詳細

問題設定

無限長円形ロッド（半径aa）が角速度Ω\Omegaで粘弾性ポリマー溶液中で回転するとき、流体自由表面高さh(r,t)h(r,t)の時間発展則を研究する。

数学モデル

支配方程式

Oldroyd-B構成モデルを採用：

  • 構成テンソル進化方程式: ∂A∂t+u⋅∇A−A⋅∇u−(∇u)T⋅A=−1λ(A−I)\frac{\partial\mathbf{A}}{\partial t} + \mathbf{u} \cdot \nabla\mathbf{A} - \mathbf{A} \cdot \nabla\mathbf{u} - (\nabla\mathbf{u})^T \cdot \mathbf{A} = -\frac{1}{\lambda}(\mathbf{A} - \mathbf{I})
  • 応力テンソル: σ=−pI+2μsE+G(A−I)\boldsymbol{\sigma} = -p\mathbf{I} + 2\mu_s\mathbf{E} + G(\mathbf{A} - \mathbf{I})
  • 運動量方程式: ∇⋅σ=ρgez\nabla \cdot \boldsymbol{\sigma} = \rho g\mathbf{e}_z

無次元化

主要な無次元パラメータを導入：

  • ワイセンベルク数: Wi=λΩWi = \lambda\Omega
  • 重力パラメータ: G=ρga(μs+μp)ΩG = \frac{\rho ga}{(\mu_s + \mu_p)\Omega}
  • レイノルズ数: Re=ρa2Ωμs+μpRe = \frac{\rho a^2\Omega}{\mu_s + \mu_p}

特性的制限

著者らは主要な特性的制限を特定した：G=Wi⋅GG = Wi \cdot \mathcal{G}、ここでG=O(1)\mathcal{G} = O(1)であり、これにより二パラメータ漸近展開が単一パラメータ展開に簡略化される。

技術的革新点

1. 時間スケール分離

二つの時間スケールの存在を発見：

  • 外層時間: T=O(1)T = O(1)、定常挙動に対応
  • 内層時間: τ=T/Wi\tau = T/Wi、過渡境界層に対応

2. 領域摂動法

小さな界面変形に対して、境界条件をZ=H(R,T)Z = H(R,T)からZ=0Z = 0にテイラー展開し、自由表面境界条件の処理を簡略化。

3. ウェーバー変換

慣性効果を処理する際に、半無限領域上のベッセル微分演算子を処理するためにウェーバー変換を採用し、偏微分方程式を常微分方程式に変換。

実験設定

理論分析フレームワーク

本論文は純粋な理論研究であり、漸近解析方法を採用：

  1. 小パラメータ: Wi≪1Wi \ll 1（弱粘弾性）
  2. 重力範囲: G≫Wi−1G \gg Wi^{-1}（強い重力効果）
  3. 小変形仮定: 界面変形は相対的に小さい

分析戦略

  1. 第4節: 外層分析 - 定常界面輪郭
  2. 第5節: 内層分析 - 過渡界面輪郭
  3. 第6節: 小慣性効果を含む分析

実験結果

主要な結果

1. 定常界面輪郭

Hsteady(R)=2βpWi2GR4H_{steady}(R) = \frac{2\beta_p Wi^2}{\mathcal{G}R^4}

これはJoseph & Fosdick (1973)が二次流体モデルを使用して得た結果と一致している。

2. 過渡界面輪郭

Htransient(R,τ)=2βp(1−(τ+1)e−τ)Wi2GR4H_{transient}(R,\tau) = \frac{2\beta_p(1-(\tau+1)e^{-\tau})Wi^2}{\mathcal{G}R^4}

ここでτ=T/Wi\tau = T/Wiは拡張時間変数である。

3. 登攀判定基準

慣性効果を含める場合、ロッド登攀が発生する条件は： 4Elβp>14El\beta_p > 1 ここでEl=Wi/ReEl = Wi/Reは弾性数である。

時間発展の特性

  1. 急速上昇段階: 界面はτ≈6\tau \approx 6の時間内に急速に上昇
  2. 定常状態への接近: 過渡解は長時間極限で定常解に収束
  3. 非単調挙動: 慣性を含める場合、界面高さは先に上昇してから低下する可能性がある

主要な発見

  1. 周方向応力: ポリマー伸張により生成される周方向応力が登攀を駆動する主要なメカニズム
  2. 慣性競争: 慣性により生成される遠心力が粘弾性の周方向応力と競争
  3. 境界層構造: 時間境界層の存在は、急速な過渡と緩やかな定常の二段階挙動を説明

関連研究

歴史的発展

  1. Serrin (1959): 初の理論的試み、二次流体モデルを使用
  2. Giesekus (1961): 慣性効果を無視した、より広い構成方程式
  3. Joseph & Fosdick (1973): 領域摂動法フレームワークの系統的研究
  4. Oldroyd (1950, 1951): 初期のOldroyd-Bモデル応用

本論文の利点

  1. 時間依存性: 初の完全な過渡分析
  2. モデル安定性: 二次流体モデルの不安定性問題を回避
  3. 無次元フレームワーク: 系統的な無次元パラメータ分析
  4. 慣性効果: 慣性と粘弾性の相互作用を初めて研究

結論と考察

主要な結論

  1. 境界層構造: ロッド登攀現象には明らかな時間境界層が存在し、界面は急速に上昇した後、徐々に定常状態に接近する
  2. スケール分離: 過渡時間スケールはO(Wi)O(Wi)、定常時間スケールはO(1)O(1)
  3. 登攀メカニズム: ポリマー伸張により生成される周方向応力が登攀を駆動する核心メカニズム
  4. 慣性の影響: 小さな慣性効果は登攀挙動に大きな影響を与える可能性があり、臨界弾性数が存在

限界

  1. 小変形仮定: 小さな界面変形の場合にのみ適用可能
  2. 弱粘弾性: Wi≪1Wi \ll 1の範囲に限定
  3. 強い重力: G≫Wi−1G \gg Wi^{-1}条件が必要
  4. 簡略化モデル: 表面張力と接触角効果を無視

今後の方向性

  1. 高ワイセンベルク数: Wi≫1Wi \gg 1の場合の挙動を研究
  2. 複雑な構成モデル: FENE-P、Giesekusなどのモデルへの拡張
  3. 有限変形: 大きな界面変形の場合を処理
  4. 実験的検証: 実験データとの比較検証

深い評価

利点

  1. 理論的厳密性: 数学的導出が厳密で、漸近解析方法が適切に適用されている
  2. 物理的洞察: 時間境界層構造と登攀メカニズムを解明
  3. 方法的革新: ウェーバー変換と領域摂動法の巧妙な組み合わせ
  4. 実用的価値: 流動曲線測定に理論的基礎を提供

不足点

  1. パラメータ範囲: 適用可能なパラメータ範囲が相対的に限定的
  2. 実験との比較: 実験データとの直接的な比較が不足
  3. 数値検証: 一部の結果は数値計算のみで、解析的表現が不足

影響力

  1. 学術的貢献: 過渡ロッド登攀理論の空白を埋める
  2. 方法的価値: 提供される分析方法は他の複雑流体問題に推広可能
  3. 工学的応用: 流動曲線計の設計と材料表征に理論的支援を提供

適用場面

  1. 希薄ポリマー溶液: Boger流体などの定剪断粘度粘弾性流体
  2. 低回転速度: 小ワイセンベルク数の実験条件
  3. 流動曲線測定: 法線応力係数と緩和時間の測定

参考文献

  1. Joseph, D.D. & Fosdick, R.L. (1973). The free surface on a liquid between cylinders rotating at different speeds. Arch. Rat. Mech. Anal., 49, 321-380.
  2. Oldroyd, J.G. (1950). On the formulation of rheological equations of state. Proc. R. Soc. A, 200, 523-541.
  3. Bird, R.B., Armstrong, R.C. & Hassager, O. (1987). Dynamics of Polymeric Liquids, Vol. 1, 2nd edn. John Wiley and Sons.

総合評価: これは高品質な理論流体力学論文であり、Oldroyd-B流体におけるロッド登攀現象の過渡挙動を初めて系統的に分析している。論文の数学的導出は厳密で、物理的洞察は深く、方法は革新的であり、複雑流体理論と流動曲線学の応用に重要な貢献をしている。パラメータ範囲の制限はあるが、その理論的フレームワークと分析方法は優れた推広価値を持つ。