Inertia-driven and elastoinertial viscoelastic turbulent channel flow simulated with a hybrid pseudo-spectral/finite-difference numerical scheme

08/24/2020
by   Lu Zhu, et al.
0

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.

READ FULL TEXT

page 11

page 12

page 21

research
01/07/2021

Pseudospectral methods with PML for nonlinear Klein-Gordon equations in classical and non-relativistic regimes

Two different Perfectly Matched Layer (PML) formulations with efficient ...
research
07/26/2020

Numerical scheme based on the spectral method for calculating nonlinear hyperbolic evolution equations

High-precision numerical scheme for nonlinear hyperbolic evolution equat...
research
02/02/2023

Stability of finite difference schemes for the hyperbolic initial boundary value problem by winding number computations

In this paper, we present a numerical strategy to check the strong stabi...
research
07/20/2023

Approximate Deconvolution Leray Reduced Order Model for Convection-Dominated Flows

In this paper, we propose a novel ROM stabilization strategy for under-r...
research
06/28/2019

An Energy Stable BDF2 Fourier Pseudo-Spectral Numerical Scheme for the Square Phase Field Crystal Equation

In this paper we propose and analyze an energy stable numerical scheme f...
research
03/24/2020

Novel, simple and robust contact-discontinuity capturing schemes for high speed compressible flows

The nonlinear convection terms in the governing equations of compressibl...
research
04/30/2022

A spectral element solution of the 2D linearized potential flow radiation problem

We present a scalable 2D Galerkin spectral element method solution to th...

Please sign up or login with your details

Forgot password? Click here to reset