Abstract:In this numerical study, an original approach to simulate non-isothermal viscoelastic fluid flows at high Weissenberg numbers is presented. Stable computations over a wide range of Weissenberg numbers are assured by using the root conformation approach in a finite volume framework on general unstructured meshes. The numerical stabilization framework is extended to consider thermo-rheological properties in Oldroyd-B type viscoelastic fluids. The temperature dependence of the viscoelastic fluid is modeled with the time-temperature superposition principle. Both Arrhenius and WLF shift factors can be chosen, depending on the flow characteristics. The internal energy balance takes into account both energy and entropy elasticity. Partitioning is achieved by a constant split factor. An analytical solution of the balance equations in planar channel flow is derived to verify the results of the main field variables and to estimate the numerical error. The more complex entry flow of a polyisobutylene-based polymer solution in an axisymmetric 4:1 contraction is studied and compared to experimental data from the literature. We demonstrate the stability of the method in the experimentally relevant range of high Weissenberg numbers. The results at different imposed wall temperatures, as well as Weissenberg numbers, are found to be in good agreement with experimental data. Furthermore, the division between energy and entropy elasticity is investigated in detail with regard to the experimental setup.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop a numerical method to stably and efficiently simulate the non - isothermal viscoelastic fluid flow at high Weissenberg numbers (Wi). Specifically, the research aims to:
1. **Propose an extended numerical stabilization method**: In order to simulate non - isothermal viscoelastic fluid flow under experimental conditions, the authors propose a new numerical stabilization method, especially at high Weissenberg numbers, where such flow usually leads to numerical instability. They use the root conformation approach to solve this problem.
2. **Study complex inlet flows**: By simulating complex inlet flows at high Weissenberg numbers, predict thermo - rheological properties, such as viscous dissipation. The specific research object is a highly elastic polymer solution based on polyisobutylene (PIB - Boger fluid), and compare it with experimental data.
3. **Explore the influence of the energy - splitting factor**: Study the influence of the energy - splitting factor on the flow behavior in the two limit cases of pure energy elasticity and pure entropy elasticity. The energy - splitting factor determines the proportion of mechanical energy converted into thermal energy, which is crucial for understanding the energy storage mechanism inside the fluid.
### Main content of the paper
#### 1. Introduction
- **Importance of non - isothermal viscoelastic fluids**: In many industrial applications, such as polymer processing, non - isothermal effects cannot be ignored. Temperature changes will affect the viscosity and elasticity of the fluid, thus changing the flow behavior.
- **Time - Temperature Superposition Principle**: Describe the temperature dependence of the fluid through the Time - Temperature Superposition Principle (TTS), using Arrhenius or WLF (William - Landel - Ferry) shift factors.
- **Internal energy balance equation**: Introduce a constant splitting factor to describe the ratio between entropy elasticity and energy elasticity, simplifying the internal energy balance equation.
#### 2. Mathematical model
- **Continuity and momentum balance equations**: Basic equations for describing viscoelastic fluid dynamics.
- **Oldroyd - B model**: Constitutive equation for describing the polymer stress tensor.
- **Temperature dependence**: Describe the temperature dependence of the fluid through the Arrhenius shift factor.
- **Internal energy balance equation**: Thermal energy equation considering viscous dissipation and elastic storage.
#### 3. Numerical model
- **Finite volume method framework**: Based on the OpenFOAM library, implement the numerical model on general unstructured grids using the Finite Volume Method (FVM).
- **Numerical stabilization**: Use the root conformation approach to stabilize the calculations at high Weissenberg numbers.
- **Discretization and velocity - stress coupling**: Adopt the High - Resolution (HR) scheme and the Total Variation Diminishing (TVD) format to discretize the convection terms, and use the Deferred Correction (DC) method to handle the conformation tensor equation.
#### 4. Verification and validation
- **Analytical solution**: Derive the analytical solution for plane channel flow to verify the numerical results.
- **Grid convergence**: Verify grid convergence and numerical errors through the calculation results at different grid resolutions.
#### 5. Comparison with experimental data
- **Test case setup**: Simulate the axisymmetric contraction flow in the experiment of Yesilata et al., setting different wall temperatures and Weissenberg numbers.
- **Result discussion**: Compare the numerical results with the experimental data to verify the validity of the model. Pay special attention to the temperature distribution near the wall and the flow field characteristics.
### Conclusion
The paper has successfully developed a numerical method that can stably simulate non - isothermal viscoelastic fluid flow at high Weissenberg numbers. By comparing with experimental data, the accuracy and stability of the model are verified. In addition, the research also explores the influence of the energy - splitting factor on the flow behavior, providing a valuable reference for future research.