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

Mixing in viscoelastic fluids using elastic turbulence

基本信息

  • 论文ID: 2409.06391
  • 标题: Mixing in viscoelastic fluids using elastic turbulence
  • 作者: Reinier van Buel, Holger Stark (Technische Universität Berlin)
  • 分类: physics.flu-dyn
  • 发表时间: 2024年9月
  • 论文链接: https://arxiv.org/abs/2409.06391

摘要

本文研究了弹性湍流对二维Taylor-Couette几何中粘弹性流体内标量浓度场混合的影响，使用Oldroyd-B模型的数值解。流动状态通过二次流序参量确定，指示弹性湍流区域。当从湍流状态开始并随后降低Weissenberg数时，在Wic以下出现弱混沌流。湍流和弱混沌流状态中的平流都诱导混合，通过溶质浓度相对于均匀分布的标准差的时间演化来说明。特别是，在弹性湍流状态下混合很强，通过混合速率、混合时间和混合效率来量化。这三个量都遵循标度律。重要的是，序参量与混合速率强相关，因此也是流体内混合的良好指标。

研究背景与动机

问题背景

  1. 低Reynolds数下的混合挑战: 在牛顿流体中，低Reynolds数条件下流动为层流，混合过程缓慢，主要由扩散控制。这在微流体器件中是一个重要限制。
  2. 弹性湍流的发现: 向流体中添加高分子量聚合物可使流体变为粘弹性，即使在惯性可忽略的情况下也能表现出弹性湍流现象。
  3. 微尺度混合的重要性: 在实验室芯片(lab-on-a-chip)设备中，高效的流体混合至关重要，特别是在微米尺度下。

研究动机

  • 探索弹性湍流在提高混合效率方面的潜力
  • 理解粘弹性流体中混合机制的物理本质
  • 为微流体应用提供理论指导

核心贡献

  1. 确定转变性质: 将从层流到弹性湍流的转变识别为亚临界转变，并发现了弱混沌流状态的存在
  2. 建立混合表征方法: 引入了混合速率、混合时间和混合效率三个量化指标，都遵循标度律
  3. 发现序参量与混合的强相关性: 证明了二次流序参量与混合速率强相关，可作为流体混合程度的良好指标
  4. 验证Batchelor混合机制: 在弹性湍流状态下观察到指数衰减，对应Batchelor混合区域，与实验结果一致

方法详解

几何设置与控制方程

研究采用二维Taylor-Couette几何，包含两个同心圆柱：

  • 外圆柱半径: 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篇相关文献，涵盖了弹性湍流、粘弹性流动稳定性、混合理论等多个方面，为研究提供了坚实的理论基础。