Transient Rod-Climbing in an Oldroyd-B Fluid

Basic Information

• Paper ID: 2501.01016
• Title: Transient Rod-Climbing in an Oldroyd-B Fluid
• Authors: Tachin Ruangkriengsin, Rodolfo Brandão, Katie Wu, Jonghyun Hwang, Evgeniy Boyko, Howard A. Stone
• Classification: physics.flu-dyn (Fluid Dynamics)
• Submission Date: January 2, 2025 (submitted to arXiv)
• Journal: Journal of Fluid Mechanics (under review)
• Paper Link: https://arxiv.org/abs/2501.01016

Abstract

This paper investigates the Weissenberg effect (rod-climbing phenomenon) in non-Newtonian fluids, wherein the fluid interface climbs along a rotating rod. The authors employ the Oldroyd-B model and conduct asymptotic analysis using dimensionless variables to study the time-dependent interface height. The research reveals that under small deformation conditions, the governing equations exhibit a temporal boundary layer, with the interface rising rapidly on short time scales before gradually approaching a steady state. By introducing a stretched time variable, the authors derive the transient velocity field and corresponding interface profile, and subsequently reintroduce small but finite inertial effects to investigate their interaction with viscoelasticity.

Research Background and Motivation

Problem Definition

The rod-climbing phenomenon (Weissenberg effect) is one of the most characteristic behaviors of non-Newtonian fluids. When a rod rotates in a viscoelastic fluid, the fluid surface climbs up along the rod. This phenomenon still lacks deep theoretical understanding.

Significance

1. Theoretical Importance: This phenomenon is a fundamental concept for understanding complex fluid behavior and holds significant importance in teaching and research
2. Practical Value: Important applications in rheological measurements for determining material parameters such as normal stress coefficients and relaxation times
3. Engineering Applications: Provides guidance for industrial applications of polymer solutions and other complex fluids

Limitations of Existing Methods

1. Theoretical Limitations: Previous research has primarily focused on steady-state behavior using second-order fluid models
2. Model Issues: Second-order fluid models are only valid at small Deborah or Weissenberg numbers and exhibit high instability
3. Time Dependence: Lack of theoretical analysis of transient dynamics

Research Motivation

The authors selected the Oldroyd-B model to avoid problems encountered with second-order fluid models in non-steady flows, and systematically investigate the time-dependent behavior of rod-climbing for the first time.

Core Contributions

1. First Theoretical Analysis: Provides the first complete transient theoretical analysis of rod-climbing in Oldroyd-B fluids
2. Temporal Boundary Layer Theory: Discovers and analyzes the temporal boundary layer structure, revealing a two-stage process of rapid rise followed by slow approach to steady state
3. Dimensionless Analysis Framework: Establishes a systematic dimensionless analysis method based on the Weissenberg number
4. Inertial Effects Study: First investigation of the interaction between small but finite inertial effects and viscoelasticity
5. Climbing Criterion: Proposes a theoretical criterion for determining the conditions under which rod-climbing occurs: $Re < 4\beta_p Wi$

Methodology Details

Problem Setup

Investigation of the time evolution of the free surface height $h(r,t)$ of a viscoelastic polymer solution when an infinitely long cylindrical rod (radius $a$) rotates at angular velocity $\Omega$.

Mathematical Model

Governing Equations

The Oldroyd-B constitutive model is employed:

• Conformation Tensor Evolution Equation: $\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})$
• Stress Tensor: $\boldsymbol{\sigma} = -p\mathbf{I} + 2\mu_s\mathbf{E} + G(\mathbf{A} - \mathbf{I})$
• Momentum Equation: $\nabla \cdot \boldsymbol{\sigma} = \rho g\mathbf{e}_z$

Dimensionless Formulation

Key dimensionless parameters are introduced:

• Weissenberg Number: $Wi = \lambda\Omega$
• Gravity Parameter: $G = \frac{\rho ga}{(\mu_s + \mu_p)\Omega}$
• Reynolds Number: $Re = \frac{\rho a^2\Omega}{\mu_s + \mu_p}$

Characteristic Scaling

The authors identify a key characteristic scaling: $G = Wi \cdot \mathcal{G}$, where $\mathcal{G} = O(1)$, which simplifies the two-parameter asymptotic expansion to a single-parameter expansion.

Technical Innovations

1. Time Scale Separation

Two distinct time scales are identified:

• Outer Time Scale: $T = O(1)$, corresponding to steady-state behavior
• Inner Time Scale: $\tau = T/Wi$, corresponding to the transient boundary layer

2. Domain Perturbation Method

For small interface deformations, boundary conditions are Taylor-expanded from $Z = H(R,T)$ to $Z = 0$, simplifying the treatment of free surface boundary conditions.

3. Weber Transform

When handling inertial effects, the Weber transform is employed to treat the Bessel differential operator on a semi-infinite domain, converting partial differential equations into ordinary differential equations.

Experimental Setup

Theoretical Analysis Framework

This paper is purely theoretical, employing asymptotic analysis methods:

1. Small Parameter: $Wi \ll 1$ (weak viscoelasticity)
2. Gravity Range: $G \gg Wi^{-1}$ (strong gravitational effects)
3. Small Deformation Assumption: Relatively small interface deformation

Analysis Strategy

1. Section 4: Outer layer analysis – steady-state interface profile
2. Section 5: Inner layer analysis – transient interface profile
3. Section 6: Analysis including small inertial effects

Experimental Results

Main Results

1. Steady-State Interface Profile

$H_{steady}(R) = \frac{2\beta_p Wi^2}{\mathcal{G}R^4}$

This is consistent with results obtained by Joseph & Fosdick (1973) using a second-order fluid model.

2. Transient Interface Profile

$H_{transient}(R,\tau) = \frac{2\beta_p(1-(\tau+1)e^{-\tau})Wi^2}{\mathcal{G}R^4}$

where $\tau = T/Wi$ is the stretched time variable.

3. Climbing Criterion

When inertial effects are included, the condition for rod-climbing to occur is: $4El\beta_p > 1$ where $El = Wi/Re$ is the elasticity number.

Temporal Evolution Characteristics

1. Rapid Rise Phase: The interface rises rapidly within approximately $\tau \approx 6$ time units
2. Steady-State Approach: The transient solution converges to the steady-state solution in the long-time limit
3. Non-Monotonic Behavior: When inertia is included, the interface height may initially rise then subsequently decrease

Key Findings

1. Azimuthal Stress: The azimuthal stress generated by polymer stretching is the primary mechanism driving climbing
2. Inertial Competition: The centrifugal force generated by inertia competes with the azimuthal stress from viscoelasticity
3. Boundary Layer Structure: The existence of the temporal boundary layer explains the two-stage behavior of rapid transient and slow steady-state approach

Related Work

Historical Development

1. Serrin (1959): First theoretical attempt using a second-order fluid model
2. Giesekus (1961): More general constitutive equations neglecting inertial effects
3. Joseph & Fosdick (1973): Systematic domain perturbation method framework
4. Oldroyd (1950, 1951): Early applications of the Oldroyd-B model

Advantages of This Work

1. Time Dependence: First complete transient analysis
2. Model Stability: Avoids instability problems of second-order fluid models
3. Dimensionless Framework: Systematic analysis of dimensionless parameters
4. Inertial Effects: First investigation of the interaction between inertia and viscoelasticity

Conclusions and Discussion

Main Conclusions

1. Boundary Layer Structure: The rod-climbing phenomenon exhibits a pronounced temporal boundary layer, with the interface rising rapidly before slowly approaching steady state
2. Time Scale Separation: The transient time scale is $O(Wi)$, while the steady-state time scale is $O(1)$
3. Climbing Mechanism: The azimuthal stress generated by polymer stretching is the core mechanism driving climbing
4. Inertial Influence: Small inertial effects can significantly affect climbing behavior, with a critical elasticity number

Limitations

1. Small Deformation Assumption: Only applicable to small interface deformations
2. Weak Viscoelasticity: Limited to the range $Wi \ll 1$
3. Strong Gravity: Requires the condition $G \gg Wi^{-1}$
4. Simplified Model: Neglects surface tension and contact angle effects

Future Directions

1. High Weissenberg Numbers: Investigation of behavior at $Wi \gg 1$
2. Complex Constitutive Models: Extension to FENE-P, Giesekus, and other models
3. Finite Deformations: Treatment of large interface deformations
4. Experimental Validation: Comparison with experimental data

In-Depth Evaluation

Strengths

1. Theoretical Rigor: Mathematically rigorous derivations with appropriate application of asymptotic analysis methods
2. Physical Insights: Reveals the temporal boundary layer structure and climbing mechanism
3. Methodological Innovation: Clever combination of Weber transform and domain perturbation methods
4. Practical Value: Provides theoretical foundation for rheological measurements

Weaknesses

1. Parameter Range: Relatively limited applicable parameter range
2. Experimental Comparison: Lacks direct comparison with experimental data
3. Numerical Verification: Some results presented only through numerical computation without analytical expressions

Impact

1. Academic Contribution: Fills the gap in transient rod-climbing theory
2. Methodological Value: The analysis methods can be generalized to other complex fluid problems
3. Engineering Applications: Provides theoretical support for rheometer design and material characterization

Applicable Scenarios

1. Dilute Polymer Solutions: Boger fluids and other viscoelastic fluids with constant shear viscosity
2. Low Rotation Speeds: Experimental conditions with small Weissenberg numbers
3. Rheological Measurements: Determination of normal stress coefficients and relaxation times

References

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.

Overall Assessment: This is a high-quality theoretical fluid mechanics paper that systematically analyzes the transient behavior of rod-climbing in Oldroyd-B fluids for the first time. The paper features rigorous mathematical derivations, profound physical insights, and innovative methodology, making important contributions to complex fluid theory and rheological applications. Despite limitations in parameter range, the theoretical framework and analytical methods have excellent potential for generalization.