Abstract:Simulations of elastic turbulence, the chaotic flow of highly elastic and inertialess polymer solutions, are plagued by numerical difficulties: The chaotically advected polymer conformation tensor develops extremely large gradients and can loose its positive definiteness, which triggers numerical instabilities. While efforts to tackle these issues have produced a plethora of specialized techniques -- tensor decompositions, artificial diffusion, and shock-capturing advection schemes -- we still lack an unambiguous route to accurate and efficient simulations. In this work, we show that even when a simulation is numerically stable, maintaining positive-definiteness and displaying the expected chaotic fluctuations, it can still suffer from errors significant enough to distort the large-scale dynamics and flow-structures. Focusing on two-dimensional simulations of the Oldroyd-B and FENE-P equations, we first compare two decompositions of the conformation tensor: symmetric square root (SSR) and Cholesky with a logarithmic transformation (Cholesky-log). While both simulations yield chaotic flows, only the Cholesky-log preserves the pattern of the forcing, i.e., its vortical cells remain ordered in a lattice as opposed to the vortices of the SSR simulations which shrink, expand and reorient constantly. To identify the accurate simulation, we appeal to a hitherto overlooked mathematical bound on the determinant of the conformation tensor, which unequivocally rejects the SSR simulation. Importantly, the accuracy of the Cholesky-log simulation is shown to arise from the logarithmic transformation. We then consider local artificial diffusion, a potential low-cost alternative to high-order advection schemes, and find unfortunately that it significantly modifies the dynamics. We end with an example, showing how the spurious large-scale motions identified here contaminate predictions of scalar mixing.

Effect of polymer-stress diffusion in the numerical simulation of elastic turbulence

Conserving elastic turbulence numerically using artificial diffusivity

Polymer diffusive instability leading to elastic turbulence in plane Couette flow

Preserving large-scale features in simulations of elastic turbulence

Polymeric diffusive instability leading to elastic turbulence in plane Couette flow

Numerical Simulation of Interactions Between Rigid Rod-Like Polymers and Coherent Structures in a Mixing Layer

NUMERICAL RESEARCH ON THE EFFECT OF FIBERS ON THE PROPERTY OF TURBULENT FIBER SUSPENSION IN A CHANNEL

Transition route to elastic and elasto-inertial turbulence in polymer channel flows

Studying polymer diffusiophoresis with Non-Equilibrium Molecular Dynamics

The interplay of inertia and elasticity in polymeric flows

A note on stress-driven anisotropic diffusion and its role in active deformable media

Inertial enhancement of the polymer diffusive instability

Effects of polymer additives on Rayleigh-Taylor turbulence

Abstract Submitted for the DFD11 Meeting of The American Physical Society Concentration dependence of the effects of polymer additives on bulk turbulence1

High-fidelity Robust and Efficient Finite Difference Algorithm for Simulation of Polymer-Induced Turbulence in Cylindrical Coordinates

Turbulent duct flow with polymers

Less is more: modelling polymers in turbulent flows

Polymers in turbulence: stretching statistics and the role of extreme strain-rate fluctuations

Intermittency in the not-so-smooth elastic turbulence

A Novel Fractional Brownian Dynamics Method for Simulating the Dynamics of Confined Bottle-Brush Polymers in Viscoelastic Solution