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

Transient rod-climbing in an Oldroyd-B fluid

基本信息

  • 论文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

摘要

本文研究了非牛顿流体中的Weissenberg效应（杆爬升现象），即流体界面沿旋转杆上升的现象。作者采用Oldroyd-B模型，通过渐近分析和无量纲变量，研究了时间相关的界面高度。研究发现在小变形情况下，控制方程表明存在时间边界层，界面在短时间尺度上快速上升，然后逐渐接近稳态。通过引入拉伸时间变量，作者推导出了瞬态速度场和相应的界面轮廓，并重新引入小而有限的惯性效应来研究其与粘弹性的相互作用。

研究背景与动机

问题定义

杆爬升现象（Weissenberg效应）是非牛顿流体最具标志性的行为之一，当杆在粘弹性流体中旋转时，流体表面会沿着杆向上攀爬。这一现象在理论上仍缺乏深入理解。

重要性

  1. 理论意义: 该现象是理解复杂流体行为的基础概念，在教学和研究中具有重要地位
  2. 实用价值: 在流变学测量中有重要应用，可用于确定材料参数如法向应力系数和松弛时间
  3. 工程应用: 对聚合物溶液和其他复杂流体的工业应用具有指导意义

现有方法局限性

  1. 理论局限: 以往研究主要集中在稳态行为，使用二阶流体模型
  2. 模型问题: 二阶流体模型仅在小Deborah数或Weissenberg数下有效，且高度不稳定
  3. 时间依赖性: 缺乏对瞬态动力学的理论分析

研究动机

作者选择Oldroyd-B模型来避免二阶流体模型在非稳态流动中的问题，并首次系统研究杆爬升现象的时间依赖性。

核心贡献

  1. 首次理论分析: 提供了Oldroyd-B流体中杆爬升现象的首个完整瞬态理论分析
  2. 时间边界层理论: 发现并分析了时间边界层结构，揭示了快速上升和缓慢趋向稳态的两阶段过程
  3. 无量纲分析框架: 建立了基于Weissenberg数的系统性无量纲分析方法
  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

无量纲化

引入关键无量纲参数：

  • Weissenberg数: Wi=λΩWi = \lambda\Omega
  • 重力参数: G=ρga(μs+μp)ΩG = \frac{\rho ga}{(\mu_s + \mu_p)\Omega}
  • Reynolds数: 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. Weber变换

在处理惯性效应时，采用Weber变换处理半无限域上的Bessel微分算子，将偏微分方程转化为常微分方程。

实验设置

理论分析框架

本文为纯理论研究，采用渐近分析方法：

  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. 高Weissenberg数: 研究Wi≫1Wi \gg 1情况下的行为
  2. 复杂本构模型: 扩展到FENE-P、Giesekus等模型
  3. 有限变形: 处理大界面变形情况
  4. 实验验证: 与实验数据的对比验证

深度评价

优点

  1. 理论严谨性: 数学推导严密，渐近分析方法应用得当
  2. 物理洞察: 揭示了时间边界层结构和爬升机制
  3. 方法创新: Weber变换和域扰动方法的巧妙结合
  4. 实用价值: 为流变学测量提供了理论基础

不足

  1. 参数范围: 适用参数范围相对有限
  2. 实验对比: 缺乏与实验数据的直接对比
  3. 数值验证: 部分结果仅给出数值计算，缺乏解析表达式

影响力

  1. 学术贡献: 填补了瞬态杆爬升理论的空白
  2. 方法价值: 提供的分析方法可推广到其他复杂流体问题
  3. 工程应用: 为流变仪设计和材料表征提供理论支持

适用场景

  1. 稀聚合物溶液: Boger流体等恒剪切粘度粘弹性流体
  2. 低转速: 小Weissenberg数的实验条件
  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流体中杆爬升现象的瞬态行为。论文数学推导严谨，物理洞察深刻，方法具有创新性，为复杂流体理论和流变学应用做出了重要贡献。尽管存在参数范围限制，但其理论框架和分析方法具有很好的推广价值。