Preserving large-scale features in simulations of elastic turbulence

TL;DR

This paper investigates numerical methods for simulating elastic turbulence, highlighting how certain tensor decompositions and transformations affect the accuracy of large-scale flow features, and emphasizes the importance of mathematical bounds for validation.

Contribution

It compares tensor decomposition techniques and introduces a determinant-based bound to identify accurate simulations, emphasizing the role of the logarithmic Cholesky transformation.

Findings

Cholesky-log decomposition preserves large-scale flow patterns.

Symmetric square root decomposition distorts flow structures.

Artificial diffusion significantly alters the flow dynamics.

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.