An improved algorithm for simulating three-dimensional, viscoelastic turbulence

We present a new finite-difference formulation to update the conformation tensor in dumbbell models (e.g., Oldroyd-B, FENE-P, Giesekus) that guarantees positive eigenvalues of the tensor (i.e., the tensor remains positive definite) and prevents over-extension for finite-extensible models. The formulation is a generalization of the second-order, central difference scheme developed by Kurganov and Tadmor [A. Kurganov, E. Tadmor, New high-resolution central schemes for nonlinear conservation laws and convection-diffusion equations, J. Comput. Phys. 160 (2000) 241-282] that guarantees a scalar field remains everywhere positive. We have extended the algorithm to guarantee a tensor field remains everywhere positive definite following an update. Extensive testing of the algorithm shows that the volume average of the conformation tensor is conserved. Furthermore, volume averages of the conformation tensor in homogeneous turbulent shear flow made over the Eulerian grid are in quantitative agreement with Lagrangian averages made over fluid particles moving throughout the domain, highlighting the accuracy of the treatment of the convective terms.