Preserving large-scale features in simulations of elastic turbulence
Fuente:
arXiv
Saved in:
| Main Authors: | , , , |
|---|---|
| Format: | Preprint |
| Published: |
2023
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866915034603454464 |
|---|---|
| author | Yerasi, Sumithra R. Picardo, Jason R. Gupta, Anupam Vincenzi, Dario |
| author_facet | Yerasi, Sumithra R. Picardo, Jason R. Gupta, Anupam Vincenzi, Dario |
| contents | 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. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2312_09165 |
| institution | arXiv |
| publishDate | 2023 |
| record_format | arxiv |
| spellingShingle | Preserving large-scale features in simulations of elastic turbulence Yerasi, Sumithra R. Picardo, Jason R. Gupta, Anupam Vincenzi, Dario Fluid Dynamics Chaotic Dynamics Computational Physics 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. |
| title | Preserving large-scale features in simulations of elastic turbulence |
| topic | Fluid Dynamics Chaotic Dynamics Computational Physics |
| url | https://arxiv.org/abs/2312.09165 |