Inertia-driven and elastoinertial viscoelastic turbulent channel flow simulated with a hybrid pseudo-spectral/finite-difference numerical scheme
arXiv:2008.10660 · doi:10.1016/j.jnnfm.2020.104410
Abstract
Numerical simulation of viscoelastic flows is challenging because of the hyperbolic nature of viscoelastic constitutive equations. Despite their superior accuracy and efficiency, pseudo-spectral methods require the introduction of artificial diffusion (AD) for numerical stability in hyperbolic problems, which alters the physical nature of the system. This study presents a hybrid numerical procedure that integrates an upwind total variation diminishing (TVD) finite-difference scheme, which is known for its stability in hyperbolic problems, for the polymer stress convection term into an overall pseudo-spectral numerical framework. Numerically stable solutions are obtained for Weissenberg number well beyond O(100) without the need for either global or local AD. Side-by-side comparison with an existing pseudo-spectral code reveals the impact of AD, which is shown to differ drastically between flow regimes. Elastoinertial turbulence (EIT) becomes unphysically suppressed when AD, at any level necessary for stabilizing the pseudo-spectral method, is used. This is attributed to the importance of sharp stress shocks in its self-sustaining cycles. Nevertheless , in regimes dominated by the classical inertial mechanism for turbulence generation, there is still an acceptable range of AD that can be safely used to predict the statistics, dynamics, and structures of drag-reduced turbulence. Detailed numerical resolution analysis of the new hybrid method, especially for capturing the EIT states, is also presented.
arXiv commentPostprint (accepted version after peer review)
References in corpus (14)
- Visualizing the geometry of state space in plane Couette flow
- Elasto-inertial turbulence
- Turbulent drag reduction by polymer additives: Fundamentals and recent advances
- Exceeding the asymptotic limit of polymer drag reduction
- Active and hibernating turbulence in minimal channel flow of Newtonian and polymeric fluids
- Two-dimensional dynamics of elasto-inertial turbulence and its role in polymer drag reduction
- Strong polymer-turbulence interactions in viscoelastic turbulent channel flow
- Critical-layer structures and mechanisms in elastoinertial turbulence
- Dynamics of viscoelastic pipe flow in the maximum drag reduction limit
- Effect of polymer-stress diffusion in the numerical simulation of elastic turbulence
- Nonlinear traveling waves as a framework for understanding turbulent drag reduction
- Nonasymptotic elastoinertial turbulence for asymptotic drag reduction
- Vortex axis tracking by iterative propagation (VATIP): a method for analysing three-dimensional turbulent structures
- Vortex dynamics in low- and high-extent polymer drag reduction regimes revealed by vortex tracking and conformation analysis